NGS-determined molecular markers and disease burden metrics from ctDNA correlate with PFS in previously untreated DLBCL.

Personalized risk stratification and treatment may help improve outcomes among patients with diffuse large B-cell lymphoma (DLBCL). We developed a next-generation sequencing (NGS)-based method to assess a range of potential prognostic indicators, and evaluated it using pretreatment plasma samples from 310 patients with previously untreated DLBCL from the GOYA trial (NCT01287741). Variant calls and DLBCL subtyping with the plasma-based method were concordant with corresponding tissue-based methods. Patients with a tumor burden greater than the median (p = .003) and non-germinal center B-cell-like (non-GCB) DLBCL (p = .049) had worse progression-free survival than patients with a tumor burden less than the median or GCB DLBCL. Multi-factor assessment combining orthogonal features from a single pretreatment plasma sample has promise as a prognostic indicator in this setting (p = .085). This minimally invasive plasma-based NGS assay could enable comprehensive prognostic assessment of patients in a clinical setting, with greater accessibility than current methods.

Quantitative PCR-Based Method to Assess Cell-Free DNA Quality, Adjust Input Mass, and Improve Next-Generation Sequencing Assay Performance.

Cell-free (cf)DNA-based testing has undergone increasingly wide adoption, including assays for the detection of circulating tumor DNA. Due to nucleosome protection, cfDNA has a distinctive fragment size of 160 to 180 bp. However, cfDNA can be contaminated with high molecular weight genomic DNA from blood cells released in plasma during sample collection. Such contamination can lead to decreased sensitivity or inconsistent results in cfDNA next-generation sequencing assays. This article describes a technical advancement in which a quantitative PCR method is used for high molecular weight contamination assessment and input mass adjustment, and has been demonstrated to improve consistency of performance in a circulating tumor DNA next-generation sequencing workflow.

Biological and technical factors in the assessment of blood-based tumor mutational burden (bTMB) in patients with NSCLC.

BACKGROUND: Patients treated with immunotherapy are at risk of considerable adverse events, and the ongoing struggle is to accurately identify the subset of patients who will benefit. Tumor mutational burden (TMB) has emerged as a promising predictive biomarker but requires tumor tissue which is not always available. Blood-based TMB (bTMB) may provide a minimally invasive assessment of mutational load. However, because of the required sequencing depth, bTMB analysis is costly and prone to false negative results. This study attempted to design a minimally sized bTMB panel, examined a counting-based method for bTMB in patients with stage I to IV non-small cell lung cancer (NSCLC) and evaluated both technical factors such as bTMB and tissue-based TMB (tTMB) cut-off, as well as sample-related factors such as cell-free DNA input mass which influence the correlation between bTMB and tTMB. METHODS: Tissue, plasma, and whole blood samples collected as part of the LEMA trial (NCT02894853) were used in this study. Samples of 185 treatment naïve patients with stage I to IV NSCLC were sequenced at the Roche Sequencing Solutions with a custom panel designed for TMB, using reagents and workflows derived from the AVENIO Tumor Tissue and circulating tumor DNA Analysis Kits. RESULTS: A TMB panel of 1.1 Mb demonstrated highly accurate TMB high calls with a positive predictive value of 95% when using a tTMB cut-off of 16 mut/Mb, corresponding with 42 mut/Mb for bTMB. The positive per cent agreement (PPA) of bTMB was relatively low at 32%. In stage IV samples with at least 20 ng of cfDNA input, PPA of bTMB improved to 63% and minimizing the panel to a subset of 577 kb was possible while maintaining 63% PPA. CONCLUSION: Plasma samples with high bTMB values are highly correspondent with tTMB, whereas bTMB low results may also be the result of low tumor burden at earlier stages of disease as well as poorly shedding tumors. For advanced stages of disease, PPA (sensitivity) of bTMB is satisfactory in comparison to tTMB, even when using a panel of less than 600 kb, warranting consideration of bTMB as a predictive biomarker for patients with NSCLC eligible for immunotherapy in the future.

Concordance of Genomic Alterations by Next-Generation Sequencing in Tumor Tissue versus Cell-Free DNA in Stage I-IV Non-Small Cell Lung Cancer.

Molecular biomarkers hold promise for personalization of cancer treatment. However, a typical tumor biopsy may be difficult to acquire and may not capture genetic variations within or across tumors. Given these limitations, tumor genotyping using next-generation sequencing of plasma-derived circulating tumor (ct)-DNA has the potential to transform non-small cell lung cancer (NSCLC) management. Importantly, mutations detected in biopsied tissue must also be detected in plasma-derived ctDNA at different disease stages. Using the AVENIO ctDNA Surveillance kit (research use only), mutations in ctDNA from NSCLC subjects were compared with those identified in matched tumor tissue samples, retrospectively. Plasma and tissue samples were collected from 141 treatment-naïve NSCLC subjects (stage I, n = 48; stage II, n = 37; stage III, n = 33; stage IV, n = 23). In plasma samples, the median numbers of variants per subject were 4, 6, 8, and 9 in those with stage I, II, III, and IV disease, respectively. The corresponding values in tissue samples were 5, 5, 6, and 4. The overall tissue-plasma concordance of stage II through IV was 62.2% by AVENIO software call. On multivariate analysis, concordance was positively and significantly associated with tumor size and cancer stage. Next-generation sequencing-based analyses with the AVENIO ctDNA Surveillance kit could be an alternative approach to detecting genetic variations in plasma-derived ctDNA isolated from NSCLC subjects.

Early Detection of Molecular Residual Disease in Localized Lung Cancer by Circulating Tumor DNA Profiling.

Identifying molecular residual disease (MRD) after treatment of localized lung cancer could facilitate early intervention and personalization of adjuvant therapies. Here, we apply cancer personalized profiling by deep sequencing (CAPP-seq) circulating tumor DNA (ctDNA) analysis to 255 samples from 40 patients treated with curative intent for stage I-III lung cancer and 54 healthy adults. In 94% of evaluable patients experiencing recurrence, ctDNA was detectable in the first posttreatment blood sample, indicating reliable identification of MRD. Posttreatment ctDNA detection preceded radiographic progression in 72% of patients by a median of 5.2 months, and 53% of patients harbored ctDNA mutation profiles associated with favorable responses to tyrosine kinase inhibitors or immune checkpoint blockade. Collectively, these results indicate that ctDNA MRD in patients with lung cancer can be accurately detected using CAPP-seq and may allow personalized adjuvant treatment while disease burden is lowest.Significance: This study shows that ctDNA analysis can robustly identify posttreatment MRD in patients with localized lung cancer, identifying residual/recurrent disease earlier than standard-of-care radiologic imaging, and thus could facilitate personalized adjuvant treatment at early time points when disease burden is lowest. Cancer Discov; 7(12); 1394-403. ©2017 AACR.See related commentary by Comino-Mendez and Turner, p. 1368This article is highlighted in the In This Issue feature, p. 1355.

mRNA pseudouridylation affects RNA metabolism in the parasite Toxoplasma gondii.

RNA contains over 100 modified nucleotides that are created post-transcriptionally, among which pseudouridine (Ψ) is one of the most abundant. Although it was one of the first modifications discovered, the biological role of this modification is still not fully understood. Recently, we reported that a pseudouridine synthase (TgPUS1) is necessary for differentiation of the single-celled eukaryotic parasite Toxoplasma gondii from active to chronic infection. To better understand the biological role of pseudouridylation, we report here gel-based and deep-sequencing methods to identify TgPUS1-dependent Ψ's in Toxoplasma RNA, and the use of TgPUS1 mutants to examine the effect of this modification on mRNAs. In addition to identifying conserved sites of pseudouridylation in Toxoplasma rRNA, tRNA, and snRNA, we also report extensive pseudouridylation of Toxoplasma mRNAs, with the Ψ's being relatively depleted in the 3'-UTR but enriched at position 1 of codons. We show that many Ψ's in tRNA and mRNA are dependent on the action of TgPUS1 and that TgPUS1-dependent mRNA Ψ's are enriched in developmentally regulated transcripts. RNA-seq data obtained from wild-type and TgPUS1-mutant parasites shows that genes containing a TgPUS1-dependent Ψ are relatively more abundant in mutant parasites, while pulse/chase labeling of RNA with 4-thiouracil shows that mRNAs containing TgPUS1-dependent Ψ have a modest but statistically significant increase in half-life in the mutant parasites. These data are some of the first evidence suggesting that mRNA Ψ's play an important biological role.

Role of KEAP1/NRF2 and TP53 Mutations in Lung Squamous Cell Carcinoma Development and Radiation Resistance.

Lung squamous cell carcinoma (LSCC) pathogenesis remains incompletely understood, and biomarkers predicting treatment response remain lacking. Here, we describe novel murine LSCC models driven by loss of Trp53 and Keap1, both of which are frequently mutated in human LSCCs. Homozygous inactivation of Keap1 or Trp53 promoted airway basal stem cell (ABSC) self-renewal, suggesting that mutations in these genes lead to expansion of mutant stem cell clones. Deletion of Trp53 and Keap1 in ABSCs, but not more differentiated tracheal cells, produced tumors recapitulating histologic and molecular features of human LSCCs, indicating that they represent the likely cell of origin in this model. Deletion of Keap1 promoted tumor aggressiveness, metastasis, and resistance to oxidative stress and radiotherapy (RT). KEAP1/NRF2 mutation status predicted risk of local recurrence after RT in patients with non-small lung cancer (NSCLC) and could be noninvasively identified in circulating tumor DNA. Thus, KEAP1/NRF2 mutations could serve as predictive biomarkers for personalization of therapeutic strategies for NSCLCs. SIGNIFICANCE: We developed an LSCC mouse model involving Trp53 and Keap1, which are frequently mutated in human LSCCs. In this model, ABSCs are the cell of origin of these tumors. KEAP1/NRF2 mutations increase radioresistance and predict local tumor recurrence in radiotherapy patients. Our findings are of potential clinical relevance and could lead to personalized treatment strategies for tumors with KEAP1/NRF2 mutations. Cancer Discov; 7(1); 86-101. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 1.

Distinct biological subtypes and patterns of genome evolution in lymphoma revealed by circulating tumor DNA.

Patients with diffuse large B cell lymphoma (DLBCL) exhibit marked diversity in tumor behavior and outcomes, yet the identification of poor-risk groups remains challenging. In addition, the biology underlying these differences is incompletely understood. We hypothesized that characterization of mutational heterogeneity and genomic evolution using circulating tumor DNA (ctDNA) profiling could reveal molecular determinants of adverse outcomes. To address this hypothesis, we applied cancer personalized profiling by deep sequencing (CAPP-Seq) analysis to tumor biopsies and cell-free DNA samples from 92 lymphoma patients and 24 healthy subjects. At diagnosis, the amount of ctDNA was found to strongly correlate with clinical indices and was independently predictive of patient outcomes. We demonstrate that ctDNA genotyping can classify transcriptionally defined tumor subtypes, including DLBCL cell of origin, directly from plasma. By simultaneously tracking multiple somatic mutations in ctDNA, our approach outperformed immunoglobulin sequencing and radiographic imaging for the detection of minimal residual disease and facilitated noninvasive identification of emergent resistance mutations to targeted therapies. In addition, we identified distinct patterns of clonal evolution distinguishing indolent follicular lymphomas from those that transformed into DLBCL, allowing for potential noninvasive prediction of histological transformation. Collectively, our results demonstrate that ctDNA analysis reveals biological factors that underlie lymphoma clinical outcomes and could facilitate individualized therapy.

Circulating tumour DNA profiling reveals heterogeneity of EGFR inhibitor resistance mechanisms in lung cancer patients.

Circulating tumour DNA (ctDNA) analysis facilitates studies of tumour heterogeneity. Here we employ CAPP-Seq ctDNA analysis to study resistance mechanisms in 43 non-small cell lung cancer (NSCLC) patients treated with the third-generation epidermal growth factor receptor (EGFR) inhibitor rociletinib. We observe multiple resistance mechanisms in 46% of patients after treatment with first-line inhibitors, indicating frequent intra-patient heterogeneity. Rociletinib resistance recurrently involves MET, EGFR, PIK3CA, ERRB2, KRAS and RB1. We describe a novel EGFR L798I mutation and find that EGFR C797S, which arises in ∼33% of patients after osimertinib treatment, occurs in <3% after rociletinib. Increased MET copy number is the most frequent rociletinib resistance mechanism in this cohort and patients with multiple pre-existing mechanisms (T790M and MET) experience inferior responses. Similarly, rociletinib-resistant xenografts develop MET amplification that can be overcome with the MET inhibitor crizotinib. These results underscore the importance of tumour heterogeneity in NSCLC and the utility of ctDNA-based resistance mechanism assessment.

Integrated digital error suppression for improved detection of circulating tumor DNA.

High-throughput sequencing of circulating tumor DNA (ctDNA) promises to facilitate personalized cancer therapy. However, low quantities of cell-free DNA (cfDNA) in the blood and sequencing artifacts currently limit analytical sensitivity. To overcome these limitations, we introduce an approach for integrated digital error suppression (iDES). Our method combines in silico elimination of highly stereotypical background artifacts with a molecular barcoding strategy for the efficient recovery of cfDNA molecules. Individually, these two methods each improve the sensitivity of cancer personalized profiling by deep sequencing (CAPP-Seq) by about threefold, and synergize when combined to yield ∼15-fold improvements. As a result, iDES-enhanced CAPP-Seq facilitates noninvasive variant detection across hundreds of kilobases. Applied to non-small cell lung cancer (NSCLC) patients, our method enabled biopsy-free profiling of EGFR kinase domain mutations with 92% sensitivity and >99.99% specificity at the variant level, and with 90% sensitivity and 96% specificity at the patient level. In addition, our approach allowed monitoring of NSCLC ctDNA down to 4 in 10(5) cfDNA molecules. We anticipate that iDES will aid the noninvasive genotyping and detection of ctDNA in research and clinical settings.

Organocatalytic removal of formaldehyde adducts from RNA and DNA bases.

Formaldehyde is universally used to fix tissue specimens, where it forms hemiaminal and aminal adducts with biomolecules, hindering the ability to retrieve molecular information. Common methods for removing these adducts involve extended heating, which can cause extensive degradation of nucleic acids, particularly RNA. Here, we show that water-soluble bifunctional catalysts (anthranilates and phosphanilates) speed the reversal of formaldehyde adducts of mononucleotides over standard buffers. Studies with formaldehyde-treated RNA oligonucleotides show that the catalysts enhance adduct removal, restoring unmodified RNA at 37 °C even when extensively modified, while avoiding the high temperatures that promote RNA degradation. Experiments with formalin-fixed, paraffin-embedded cell samples show that the catalysis is compatible with common RNA extraction protocols, with detectable RNA yields increased by 1.5-2.4-fold using a catalyst under optimized conditions and by 7-25-fold compared with a commercial kit. Such catalytic strategies show promise for general use in reversing formaldehyde adducts in clinical specimens.

Transcriptome-wide mapping of pseudouridines: pseudouridine synthases modify specific mRNAs in S. cerevisiae.

We developed a novel technique, called pseudouridine site identification sequencing (PSI-seq), for the transcriptome-wide mapping of pseudouridylation sites with single-base resolution from cellular RNAs based on the induced termination of reverse transcription specifically at pseudouridines following CMCT treatment. PSI-seq analysis of RNA samples from S. cerevisiae correctly detected all of the 43 known pseudouridines in yeast 18S and 25S ribosomal RNA with high specificity. Moreover, application of PSI-seq to the yeast transcriptome revealed the presence of site-specific pseudouridylation within dozens of mRNAs, including RPL11a, TEF1, and other genes implicated in translation. To identify the mechanisms responsible for mRNA pseudouridylation, we genetically deleted candidate pseudouridine synthase (Pus) enzymes and reconstituted their activities in vitro. These experiments demonstrated that the Pus1 enzyme was necessary and sufficient for pseudouridylation of RPL11a mRNA, whereas Pus4 modified TEF1 mRNA, and Pus6 pseudouridylated KAR2 mRNA. Finally, we determined that modification of RPL11a at Ψ -68 was observed in RNA from the related yeast S. mikitae, and Ψ -239 in TEF1 mRNA was maintained in S. mikitae as well as S. pombe, indicating that these pseudouridylations are ancient, evolutionarily conserved RNA modifications. This work establishes that site-specific pseudouridylation of eukaryotic mRNAs is a genetically programmed RNA modification that naturally occurs in multiple yeast transcripts via distinct mechanisms, suggesting that mRNA pseudouridylation may provide an important novel regulatory function. The approach and strategies that we report here should be generally applicable to the discovery of pseudouridylation, or other RNA modifications, in diverse biological contexts.

U2 snRNP binds intronless histone pre-mRNAs to facilitate U7-snRNP-dependent 3' end formation.

In metazoa, pre-mRNA 3' end formation occurs via two pathways: cleavage/polyadenylation for the majority of RNA polymerase II transcripts and U7-snRNP-dependent cleavage for replication-dependent histone pre-mRNAs. An RNA element derived from a replication-dependent histone gene affects multiple steps of pre-mRNA processing. Here, we demonstrate that a portion of this RNA element, present in the majority of histone mRNAs, stimulates U7-snRNP-dependent cleavage. Surprisingly, this element binds U2 snRNP, although it is derived from an intronless mRNA. Specifically, SF3b, a U2 and U12-snRNP component, contacts the RNA element both in vitro and in vivo in conjunction with hPrp43, a DEAH-box helicase. Tethering either U2 or U12 snRNP to histone pre-mRNA substrates stimulates U7-snRNP-dependent cleavage in vitro and in vivo. Finally, we show that U2 snRNP associates with histone pre-mRNAs in vivo. We conclude that U2 snRNP plays a nonsplicing role in histone mRNA maturation.