Integrated analysis of genomic and transcriptomic data for the discovery of splice-associated variants in cancer.

Somatic mutations within non-coding regions and even exons may have unidentified regulatory consequences that are often overlooked in analysis workflows. Here we present RegTools ( www.regtools.org ), a computationally efficient, free, and open-source software package designed to integrate somatic variants from genomic data with splice junctions from bulk or single cell transcriptomic data to identify variants that may cause aberrant splicing. We apply RegTools to over 9000 tumor samples with both tumor DNA and RNA sequence data. RegTools discovers 235,778 events where a splice-associated variant significantly increases the splicing of a particular junction, across 158,200 unique variants and 131,212 unique junctions. To characterize these somatic variants and their associated splice isoforms, we annotate them with the Variant Effect Predictor, SpliceAI, and Genotype-Tissue Expression junction counts and compare our results to other tools that integrate genomic and transcriptomic data. While many events are corroborated by the aforementioned tools, the flexibility of RegTools also allows us to identify splice-associated variants in known cancer drivers, such as TP53, CDKN2A, and B2M, and other genes.

A dominant negative variant of RAB5B disrupts maturation of surfactant protein B and surfactant protein C.

Pathogenic variants in surfactant proteins SP-B and SP-C cause surfactant deficiency and interstitial lung disease. Surfactant proteins are synthesized as precursors (proSP-B, proSP-C), trafficked, and processed via a vesicular-regulated secretion pathway; however, control of vesicular trafficking events is not fully understood. Through the Undiagnosed Diseases Network, we evaluated a child with interstitial lung disease suggestive of surfactant deficiency. Variants in known surfactant dysfunction disorder genes were not found in trio exome sequencing. Instead, a de novo heterozygous variant in RAB5B was identified in the Ras/Rab GTPases family nucleotide binding domain, p.Asp136His. Functional studies were performed in Caenorhabditis elegans by knocking the proband variant into the conserved position (Asp135) of the ortholog, rab-5 Genetic analysis demonstrated that rab-5[Asp135His] is damaging, producing a strong dominant negative gene product. rab-5[Asp135His] heterozygotes were also defective in endocytosis and early endosome (EE) fusion. Immunostaining studies of the proband's lung biopsy revealed that RAB5B and EE marker EEA1 were significantly reduced in alveolar type II cells and that mature SP-B and SP-C were significantly reduced, while proSP-B and proSP-C were normal. Furthermore, staining normal lung showed colocalization of RAB5B and EEA1 with proSP-B and proSP-C. These findings indicate that dominant negative-acting RAB5B Asp136His and EE dysfunction cause a defect in processing/trafficking to produce mature SP-B and SP-C, resulting in interstitial lung disease, and that RAB5B and EEs normally function in the surfactant secretion pathway. Together, the data suggest a noncanonical function for RAB5B and identify RAB5B p.Asp136His as a genetic mechanism for a surfactant dysfunction disorder.

Recurrent WNT pathway alterations are frequent in relapsed small cell lung cancer.

Nearly all patients with small cell lung cancer (SCLC) eventually relapse with chemoresistant disease. The molecular mechanisms driving chemoresistance in SCLC remain un-characterized. Here, we describe whole-exome sequencing of paired SCLC tumor samples procured at diagnosis and relapse from 12 patients, and unpaired relapse samples from 18 additional patients. Multiple somatic copy number alterations, including gains in ABCC1 and deletions in MYCL, MSH2, and MSH6, are identifiable in relapsed samples. Relapse samples also exhibit recurrent mutations and loss of heterozygosity in regulators of WNT signaling, including CHD8 and APC. Analysis of RNA-sequencing data shows enrichment for an ASCL1-low expression subtype and WNT activation in relapse samples. Activation of WNT signaling in chemosensitive human SCLC cell lines through APC knockdown induces chemoresistance. Additionally, in vitro-derived chemoresistant cell lines demonstrate increased WNT activity. Overall, our results suggest WNT signaling activation as a mechanism of chemoresistance in relapsed SCLC.

Dynamic landscape and regulation of RNA editing in mammals.

Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules. Although many editing sites have recently been discovered, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

Comprehensive genomic analysis reveals FLT3 activation and a therapeutic strategy for a patient with relapsed adult B-lymphoblastic leukemia.

The genomic events responsible for the pathogenesis of relapsed adult B-lymphoblastic leukemia (B-ALL) are not yet clear. We performed integrative analysis of whole-genome, whole-exome, custom capture, whole-transcriptome (RNA-seq), and locus-specific genomic assays across nine time points from a patient with primary de novo B-ALL. Comprehensive genome and transcriptome characterization revealed a dramatic tumor evolution during progression, yielding a tumor with complex clonal architecture at second relapse. We observed and validated point mutations in EP300 and NF1, a highly expressed EP300-ZNF384 gene fusion, a microdeletion in IKZF1, a focal deletion affecting SETD2, and large deletions affecting RB1, PAX5, NF1, and ETV6. Although the genome analysis revealed events of potential biological relevance, no clinically actionable treatment options were evident at the time of the second relapse. However, transcriptome analysis identified aberrant overexpression of the targetable protein kinase encoded by the FLT3 gene. Although the patient had refractory disease after salvage therapy for the second relapse, treatment with the FLT3 inhibitor sunitinib rapidly induced a near complete molecular response, permitting the patient to proceed to a matched-unrelated donor stem cell transplantation. The patient remains in complete remission more than 4 years later. Analysis of this patient's relapse genome revealed an unexpected, actionable therapeutic target that led to a specific therapy associated with a rapid clinical response. For some patients with relapsed or refractory cancers, this approach may indicate a novel therapeutic intervention that could alter outcome.

A Phase I Trial of BKM120 (Buparlisib) in Combination with Fulvestrant in Postmenopausal Women with Estrogen Receptor-Positive Metastatic Breast Cancer.

PURPOSE: This trial was conducted to determine the maximum tolerated dose (MTD) and preliminary efficacy of buparlisib, an oral pan-class I PI3K inhibitor, plus fulvestrant in postmenopausal women with metastatic estrogen receptor positive (ER(+)) breast cancer. EXPERIMENTAL DESIGN: Phase IA employed a 3+3 design to determine the MTD of buparlisib daily plus fulvestrant. Subsequent cohorts (phase IB and cohort C) evaluated intermittent (5/7-day) and continuous dosing of buparlisib (100 mg daily). No more than 3 prior systemic treatments in the metastatic setting were allowed in these subsequent cohorts. RESULTS: Thirty-one patients were enrolled. MTD was defined as buparlisib 100 mg daily plus fulvestrant. Common adverse events (AE) included fatigue (38.7%), transaminases elevation (35.5%), rash (29%), and diarrhea (19.4%). C-peptide was significantly increased during treatment, consistent with on-target effect of buparlisib. Compared with intermittent dosing, daily buparlisib was associated with more frequent early onset AEs and higher buparlisib plasma concentrations. Among the 29 evaluable patients, the clinical benefit rate was 58.6% (95% CI, 40.7%-74.5%). Response was not associated with PIK3CA mutation or treatment cohort; however, loss of PTEN, progesterone receptor (PgR) expression, or mutation in TP53 was most common in resistant cases, and mutations inAKT1 and ESR1 did not exclude treatment response. CONCLUSIONS: Buparlisib plus fulvestrant is clinically active with manageable AEs in patients with metastatic ER(+)breast cancer. Weekend breaks in buparlisib dosing reduced toxicity. Patients with PgR negative and TP53 mutation did poorly, suggesting buparlisib plus fulvestrant may not be adequately effective against tumors with these poor prognostic molecular features.

Optimizing cancer genome sequencing and analysis.

Tumors are typically sequenced to depths of 75-100× (exome) or 30-50× (whole genome). We demonstrate that current sequencing paradigms are inadequate for tumors that are impure, aneuploid or clonally heterogeneous. To reassess optimal sequencing strategies, we performed ultra-deep (up to ~312×) whole genome sequencing (WGS) and exome capture (up to ~433×) of a primary acute myeloid leukemia, its subsequent relapse, and a matched normal skin sample. We tested multiple alignment and variant calling algorithms and validated ~200,000 putative SNVs by sequencing them to depths of ~1,000×. Additional targeted sequencing provided over 10,000× coverage and ddPCR assays provided up to ~250,000× sampling of selected sites. We evaluated the effects of different library generation approaches, depth of sequencing, and analysis strategies on the ability to effectively characterize a complex tumor. This dataset, representing the most comprehensively sequenced tumor described to date, will serve as an invaluable community resource (dbGaP accession id phs000159).

Genome Modeling System: A Knowledge Management Platform for Genomics.

In this work, we present the Genome Modeling System (GMS), an analysis information management system capable of executing automated genome analysis pipelines at a massive scale. The GMS framework provides detailed tracking of samples and data coupled with reliable and repeatable analysis pipelines. The GMS also serves as a platform for bioinformatics development, allowing a large team to collaborate on data analysis, or an individual researcher to leverage the work of others effectively within its data management system. Rather than separating ad-hoc analysis from rigorous, reproducible pipelines, the GMS promotes systematic integration between the two. As a demonstration of the GMS, we performed an integrated analysis of whole genome, exome and transcriptome sequencing data from a breast cancer cell line (HCC1395) and matched lymphoblastoid line (HCC1395BL). These data are available for users to test the software, complete tutorials and develop novel GMS pipeline configurations. The GMS is available at https://github.com/genome/gms.

Convergent loss of PTEN leads to clinical resistance to a PI(3)Kα inhibitor.

Broad and deep tumour genome sequencing has shed new light on tumour heterogeneity and provided important insights into the evolution of metastases arising from different clones. There is an additional layer of complexity, in that tumour evolution may be influenced by selective pressure provided by therapy, in a similar fashion to that occurring in infectious diseases. Here we studied tumour genomic evolution in a patient (index patient) with metastatic breast cancer bearing an activating PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha, PI(3)Kα) mutation. The patient was treated with the PI(3)Kα inhibitor BYL719, which achieved a lasting clinical response, but the patient eventually became resistant to this drug (emergence of lung metastases) and died shortly thereafter. A rapid autopsy was performed and material from a total of 14 metastatic sites was collected and sequenced. All metastatic lesions, when compared to the pre-treatment tumour, had a copy loss of PTEN (phosphatase and tensin homolog) and those lesions that became refractory to BYL719 had additional and different PTEN genetic alterations, resulting in the loss of PTEN expression. To put these results in context, we examined six other patients also treated with BYL719. Acquired bi-allelic loss of PTEN was found in one of these patients, whereas in two others PIK3CA mutations present in the primary tumour were no longer detected at the time of progression. To characterize our findings functionally, we examined the effects of PTEN knockdown in several preclinical models (both in cell lines intrinsically sensitive to BYL719 and in PTEN-null xenografts derived from our index patient), which we found resulted in resistance to BYL719, whereas simultaneous PI(3)K p110ß blockade reverted this resistance phenotype. We conclude that parallel genetic evolution of separate metastatic sites with different PTEN genomic alterations leads to a convergent PTEN-null phenotype resistant to PI(3)Kα inhibition.