Automated real-world data integration improves cancer outcome prediction.

The digitization of health records and growing availability of tumour DNA sequencing provide an opportunity to study the determinants of cancer outcomes with unprecedented richness. Patient data are often stored in unstructured text and siloed datasets. Here we combine natural language processing annotations1,2 with structured medication, patient-reported demographic, tumour registry and tumour genomic data from 24,950 patients at Memorial Sloan Kettering Cancer Center to generate a clinicogenomic, harmonized oncologic real-world dataset (MSK-CHORD). MSK-CHORD includes data for non-small-cell lung (n = 7,809), breast (n = 5,368), colorectal (n = 5,543), prostate (n = 3,211) and pancreatic (n = 3,109) cancers and enables discovery of clinicogenomic relationships not apparent in smaller datasets. Leveraging MSK-CHORD to train machine learning models to predict overall survival, we find that models including features derived from natural language processing, such as sites of disease, outperform those based on genomic data or stage alone as tested by cross-validation and an external, multi-institution dataset. By annotating 705,241 radiology reports, MSK-CHORD also uncovers predictors of metastasis to specific organ sites, including a relationship between SETD2 mutation and lower metastatic potential in immunotherapy-treated lung adenocarcinoma corroborated in independent datasets. We demonstrate the feasibility of automated annotation from unstructured notes and its utility in predicting patient outcomes. The resulting data are provided as a public resource for real-world oncologic research.

Tumor MHC Class I Expression Associates with Intralesional IL2 Response in Melanoma.

Cancer immunotherapy can result in lasting tumor regression, but predictive biomarkers of treatment response remain ill-defined. Here, we performed single-cell proteomics, transcriptomics, and genomics on matched untreated and IL2 injected metastases from patients with melanoma. Lesions that completely regressed following intralesional IL2 harbored increased fractions and densities of nonproliferating CD8+ T cells lacking expression of PD-1, LAG-3, and TIM-3 (PD-1-LAG-3-TIM-3-). Untreated lesions from patients who subsequently responded with complete eradication of all tumor cells in all injected lesions (individuals referred to herein as "extreme responders") were characterized by proliferating CD8+ T cells with an exhausted phenotype (PD-1+LAG-3+TIM-3+), stromal B-cell aggregates, and expression of IFNγ and IL2 response genes. Loss of membranous MHC class I expression in tumor cells of untreated lesions was associated with resistance to IL2 therapy. We validated this finding in an independent cohort of metastatic melanoma patients treated with intralesional or systemic IL2. Our study suggests that intact tumor-cell antigen presentation is required for melanoma response to IL2 and describes a multidimensional and spatial approach to develop immuno-oncology biomarker hypotheses using routinely collected clinical biospecimens.

Mutational landscape of EGFR-, MYC-, and Kras-driven genetically engineered mouse models of lung adenocarcinoma.

Genetically engineered mouse models (GEMMs) of cancer are increasingly being used to assess putative driver mutations identified by large-scale sequencing of human cancer genomes. To accurately interpret experiments that introduce additional mutations, an understanding of the somatic genetic profile and evolution of GEMM tumors is necessary. Here, we performed whole-exome sequencing of tumors from three GEMMs of lung adenocarcinoma driven by mutant epidermal growth factor receptor (EGFR), mutant Kirsten rat sarcoma viral oncogene homolog (Kras), or overexpression of MYC proto-oncogene. Tumors from EGFR- and Kras-driven models exhibited, respectively, 0.02 and 0.07 nonsynonymous mutations per megabase, a dramatically lower average mutational frequency than observed in human lung adenocarcinomas. Tumors from models driven by strong cancer drivers (mutant EGFR and Kras) harbored few mutations in known cancer genes, whereas tumors driven by MYC, a weaker initiating oncogene in the murine lung, acquired recurrent clonal oncogenic Kras mutations. In addition, although EGFR- and Kras-driven models both exhibited recurrent whole-chromosome DNA copy number alterations, the specific chromosomes altered by gain or loss were different in each model. These data demonstrate that GEMM tumors exhibit relatively simple somatic genotypes compared with human cancers of a similar type, making these autochthonous model systems useful for additive engineering approaches to assess the potential of novel mutations on tumorigenesis, cancer progression, and drug sensitivity.

Targeted massively parallel sequencing of angiosarcomas reveals frequent activation of the mitogen activated protein kinase pathway.

Angiosarcomas are rare malignant mesenchymal tumors of endothelial differentiation. The clinical behavior is usually aggressive and the prognosis for patients with advanced disease is poor with no effective therapies. The genetic bases of these tumors have been partially revealed in recent studies reporting genetic alterations such as amplifications of MYC (primarily in radiation-associated angiosarcomas), inactivating mutations in PTPRB and R707Q hotspot mutations of PLCG1. Here, we performed a comprehensive genomic analysis of 34 angiosarcomas using a clinically-approved, hybridization-based targeted next-generation sequencing assay for 341 well-established oncogenes and tumor suppressor genes. Over half of the angiosarcomas (n = 18, 53%) harbored genetic alterations affecting the MAPK pathway, involving mutations in KRAS, HRAS, NRAS, BRAF, MAPK1 and NF1, or amplifications in MAPK1/CRKL, CRAF or BRAF. The most frequently detected genetic aberrations were mutations in TP53 in 12 tumors(35%) and losses of CDKN2A in9 tumors (26%). MYC amplifications were generally mutually exclusive of TP53 alterations and CDKN2A loss and were identified in 8 tumors (24%), most of which (n = 7, 88%) arose post-irradiation. Previously reported mutations in PTPRB (n = 10, 29%) and one (3%) PLCG1 R707Q mutation were also identified. Our results demonstrate that angiosarcomas are a genetically heterogeneous group of tumors, harboring a wide range of genetic alterations. The high frequency of genetic events affecting the MAPK pathway suggests that targeted therapies inhibiting MAPK signaling may be promising therapeutic avenues in patients with advanced angiosarcomas.

Metastatic Competence Can Emerge with Selection of Preexisting Oncogenic Alleles without a Need of New Mutations.

Several experimental models faithfully recapitulate many important facets of human metastatic disease. Here, we have performed whole-exome sequencing in five widely used experimental metastasis models that were independently derived through in vivo selection from heterogeneous human cancer cell lines. In addition to providing an important characterization of these model systems, our study examines the genetic evolution of metastatic phenotypes. We found that in vivo selected highly metastatic cell populations showed little genetic divergence from the corresponding parental population. However, selection of genetic variations that preexisted in parental populations, including the well-established oncogenic mutations KRAS(G13D) and BRAF(G464V), was associated with increased metastatic capability. Conversely, expression of the wild-type BRAF allele in metastatic cells inhibited metastatic outgrowth as well as tumor initiation in mice. Our findings establish that metastatic competence can arise from heterogeneous cancer cell populations without the need for acquisition of additional mutations and that such competence can benefit from further selection of tumor-initiating mutations that seed primary tumorigenesis.

Synthetic lethality in ATM-deficient RAD50-mutant tumors underlies outlier response to cancer therapy.

UNLABELLED: Metastatic solid tumors are almost invariably fatal. Patients with disseminated small-cell cancers have a particularly unfavorable prognosis, with most succumbing to their disease within two years. Here, we report on the genetic and functional analysis of an outlier curative response of a patient with metastatic small-cell cancer to combined checkpoint kinase 1 (CHK1) inhibition and DNA-damaging chemotherapy. Whole-genome sequencing revealed a clonal hemizygous mutation in the Mre11 complex gene RAD50 that attenuated ATM signaling which in the context of CHK1 inhibition contributed, via synthetic lethality, to extreme sensitivity to irinotecan. As Mre11 mutations occur in a diversity of human tumors, the results suggest a tumor-specific combination therapy strategy in which checkpoint inhibition in combination with DNA-damaging chemotherapy is synthetically lethal in tumor cells but not normal cells with somatic mutations that impair Mre11 complex function. SIGNIFICANCE: Strategies to effect deep and lasting responses to cancer therapy in patients with metastatic disease have remained difficult to attain, especially in early-phase clinical trials. Here, we present an in-depth genomic and functional genetic analysis identifying RAD50 hypomorphism as a contributing factor to a curative response to systemic combination therapy in a patient with recurrent, metastatic small-cell cancer.

A recurrent neomorphic mutation in MYOD1 defines a clinically aggressive subset of embryonal rhabdomyosarcoma associated with PI3K-AKT pathway mutations.

Rhabdomyosarcoma, a cancer of skeletal muscle lineage, is the most common soft-tissue sarcoma in children. Major subtypes of rhabdomyosarcoma include alveolar (ARMS) and embryonal (ERMS) tumors. Whereas ARMS tumors typically contain translocations generating PAX3-FOXO1 or PAX7-FOXO1 fusions that block terminal myogenic differentiation, no functionally comparable genetic event has been found in ERMS tumors. Here we report the discovery, through whole-exome sequencing, of a recurrent somatic mutation encoding p.Leu122Arg in the myogenic transcription factor MYOD1 in a distinct subset of ERMS tumors with poor outcomes that also often contain mutations altering PI3K-AKT pathway components. Previous mutagenesis studies had shown that MYOD1 with a p.Leu122Arg substitution can block wild-type MYOD1 function and bind to MYC consensus sequences, suggesting a possible switch from differentiation to proliferation. Our functional data now confirm this prediction. Thus, MYOD1 p.Leu122Arg defines a subset of rhabdomyosarcomas eligible for high-risk protocols and the development of targeted therapeutics.

An improved canine genome and a comprehensive catalogue of coding genes and non-coding transcripts.

The domestic dog, Canis familiaris, is a well-established model system for mapping trait and disease loci. While the original draft sequence was of good quality, gaps were abundant particularly in promoter regions of the genome, negatively impacting the annotation and study of candidate genes. Here, we present an improved genome build, canFam3.1, which includes 85 MB of novel sequence and now covers 99.8% of the euchromatic portion of the genome. We also present multiple RNA-Sequencing data sets from 10 different canine tissues to catalog ∼175,000 expressed loci. While about 90% of the coding genes previously annotated by EnsEMBL have measurable expression in at least one sample, the number of transcript isoforms detected by our data expands the EnsEMBL annotations by a factor of four. Syntenic comparison with the human genome revealed an additional ∼3,000 loci that are characterized as protein coding in human and were also expressed in the dog, suggesting that those were previously not annotated in the EnsEMBL canine gene set. In addition to ∼20,700 high-confidence protein coding loci, we found ∼4,600 antisense transcripts overlapping exons of protein coding genes, ∼7,200 intergenic multi-exon transcripts without coding potential, likely candidates for long intergenic non-coding RNAs (lincRNAs) and ∼11,000 transcripts were reported by two different library construction methods but did not fit any of the above categories. Of the lincRNAs, about 6,000 have no annotated orthologs in human or mouse. Functional analysis of two novel transcripts with shRNA in a mouse kidney cell line altered cell morphology and motility. All in all, we provide a much-improved annotation of the canine genome and suggest regulatory functions for several of the novel non-coding transcripts.

Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data.

A large number of computational methods have been developed for analyzing differential gene expression in RNA-seq data. We describe a comprehensive evaluation of common methods using the SEQC benchmark dataset and ENCODE data. We consider a number of key features, including normalization, accuracy of differential expression detection and differential expression analysis when one condition has no detectable expression. We find significant differences among the methods, but note that array-based methods adapted to RNA-seq data perform comparably to methods designed for RNA-seq. Our results demonstrate that increasing the number of replicate samples significantly improves detection power over increased sequencing depth.

Genome sequencing identifies a basis for everolimus sensitivity.

Cancer drugs often induce dramatic responses in a small minority of patients. We used whole-genome sequencing to investigate the genetic basis of a durable remission of metastatic bladder cancer in a patient treated with everolimus, a drug that inhibits the mTOR (mammalian target of rapamycin) signaling pathway. Among the somatic mutations was a loss-of-function mutation in TSC1 (tuberous sclerosis complex 1), a regulator of mTOR pathway activation. Targeted sequencing revealed TSC1 mutations in about 8% of 109 additional bladder cancers examined, and TSC1 mutation correlated with everolimus sensitivity. These results demonstrate the feasibility of using whole-genome sequencing in the clinical setting to identify previously occult biomarkers of drug sensitivity that can aid in the identification of patients most likely to respond to targeted anticancer drugs.

The genomic basis of adaptive evolution in threespine sticklebacks.

Marine stickleback fish have colonized and adapted to thousands of streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high-quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of twenty additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine-freshwater divergence. Our results indicate that reuse of globally shared standing genetic variation, including chromosomal inversions, has an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine-freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.

Germline mutations in BAP1 predispose to melanocytic tumors.

Common acquired melanocytic nevi are benign neoplasms that are composed of small, uniform melanocytes and are typically present as flat or slightly elevated pigmented lesions on the skin. We describe two families with a new autosomal dominant syndrome characterized by multiple, skin-colored, elevated melanocytic tumors. In contrast to common acquired nevi, the melanocytic neoplasms in affected family members ranged histopathologically from epithelioid nevi to atypical melanocytic proliferations that showed overlapping features with melanoma. Some affected individuals developed uveal or cutaneous melanomas. Segregating with this phenotype, we found inactivating germline mutations of BAP1, which encodes a ubiquitin carboxy-terminal hydrolase. The majority of melanocytic neoplasms lost the remaining wild-type allele of BAP1 by various somatic alterations. In addition, we found BAP1 mutations in a subset of sporadic melanocytic neoplasms showing histological similarities to the familial tumors. These findings suggest that loss of BAP1 is associated with a clinically and morphologically distinct type of melanocytic neoplasm.

Template-based recognition of protein fold within the midnight and twilight zones of protein sequence similarity.

Most homologous pairs of proteins have no significant sequence similarity to each other and are not identified by direct sequence comparison or profile-based strategies. However, multiple sequence alignments of low similarity homologues typically reveal a limited number of positions that are well conserved despite diversity of function. It may be inferred that conservation at most of these positions is the result of the importance of the contribution of these amino acids to the folding and stability of the protein. As such, these amino acids and their relative positions may define a structural signature. We demonstrate that extraction of this fold template provides the basis for the sequence database to be searched for patterns consistent with the fold, enabling identification of homologs that are not recognized by global sequence analysis. The fold template method was developed to address the need for a tool that could comprehensively search the midnight and twilight zones of protein sequence similarity without reliance on global statistical significance. Manual implementations of the fold template method were performed on three folds--immunoglobulin, c-lectin and TIM barrel. Following proof of concept of the template method, an automated version of the approach was developed. This automated fold template method was used to develop fold templates for 10 of the more populated folds in the SCOP database. The fold template method developed three-dimensional structural motifs or signatures that were able to return a diverse collection of proteins, while maintaining a low false positive rate. Although the results of the manual fold template method were more comprehensive than the automated fold template method, the diversity of the results from the automated fold template method surpassed those of current methods that rely on statistical significance to infer evolutionary relationships among divergent proteins.

Osteogenic differentiation of human bone marrow stromal cells on partially demineralized bone scaffolds in vitro.

Tissue engineering has been used to enhance the utility of biomaterials for clinical bone repair by the incorporation of an osteogenic cell source into a scaffold followed by the in vitro promotion of osteogenic differentiation before host implantation. In this study, three-dimensional, partially demineralized bone scaffolds were investigated for their ability to support osteogenic differentiation of human bone marrow stromal cells (BMSCs) in vitro. Dynamic cell seeding resulted in homogeneous cell attachment and infiltration within the matrix and produced significantly higher seeding efficiencies when compared with a conventional static seeding method. Dynamically seeded scaffolds were cultured for 7 and 14 days in the presence of dexamethasone and evaluated on biochemical, molecular, and morphological levels for osteogenic differentiation. Significant elevation in alkaline phosphatase activity was observed versus controls over the 14-day culture, with a transient peak indicative of early mineralization on day 7. On the basis of RT-PCR, dexamethasone-treated samples showed elevations in alkaline phosphatase and osteocalcin expression levels at 7 and 14 days over nontreated controls, while bone sialoprotein was produced only in the presence of dexamethasone at 14 days. Scanning electron microscopy evaluation of dexamethasone-treated samples at 14 days revealed primarily cuboidal cells indicative of mature osteoblasts, in contrast to nontreated controls displaying a majority of cells with a fibroblastic cell morphology. These results demonstrate that partially demineralized bone can be successfully used with human BMSCs to support osteogenic differentiation in vitro. This osseous biomaterial may offer new potential benefits as a tool for clinical bone replacement.