scTEA-db: a comprehensive database of novel terminal exon isoforms identified from human single cell transcriptomes.

The usage of alternative terminal exons results in messenger RNA (mRNA) isoforms that differ in their 3' untranslated regions (3' UTRs) and often also in their protein-coding sequences. Alternative 3' UTRs contain different sets of cis-regulatory elements known to regulate mRNA stability, translation and localization, all of which are vital to cell identity and function. In previous work, we revealed that ∼25 percent of the experimentally observed RNA 3' ends are located within regions currently annotated as intronic, indicating that many 3' end isoforms remain to be uncovered. Also, the inclusion of not yet annotated terminal exons is more tissue specific compared to the already annotated ones. Here, we present the single cell-based Terminal Exon Annotation database (scTEA-db, www.scTEA-db.org) that provides the community with 12 063 so far not yet annotated terminal exons and associated transcript isoforms identified by analysing 53 069 publicly available single cell transcriptomes. Our scTEA-db web portal offers an array of features to find and explore novel terminal exons belonging to 5538 human genes, 110 of which are known cancer drivers. In summary, scTEA-db provides the foundation for studying the biological role of large numbers of so far not annotated terminal exon isoforms in cell identity and function.

Pretreatment Tissue TCR Repertoire Evenness Is Associated with Complete Pathologic Response in Patients with NSCLC Receiving Neoadjuvant Chemoimmunotherapy.

PURPOSE: Characterization of the T-cell receptor (TCR) repertoire may be a promising source for predictive biomarkers of pathologic response to immunotherapy in locally advanced non-small cell lung cancer (NSCLC). EXPERIMENTAL DESIGN: In this study, next-generation TCR sequencing was performed in peripheral blood and tissue samples of 40 patients with NSCLC, before and after neoadjuvant chemoimmunotherapy (NADIM clinical trial, NCT03081689), considering their complete pathologic response (CPR) or non-CPR. Beyond TCR metrics, tissue clones were ranked by their frequency and spatiotemporal evolution of top 1% clones was determined. RESULTS: We have found a positive association between an uneven TCR repertoire in tissue samples at diagnosis and CPR at surgery. Moreover, TCR most frequently ranked clones (top 1%) present in diagnostic biopsies occupied greater frequency in the total clonal space of CPR patients, achieving an AUC ROC to identify CPR patients of 0.967 (95% confidence interval, 0.897-1.000; P = 0.001), and improving the results of PD-L1 tumor proportion score (TPS; AUC = 0.767; P = 0.026) or tumor mutational burden (TMB; AUC = 0.550; P = 0.687). Furthermore, tumors with high pretreatment top 1% clonal space showed similar immune cell populations but a higher immune reactive gene expression profile. Finally, the selective expansion of pretreatment tissue top 1% clones in peripheral blood of CPR patients suggests also a peripheral immunosurveillance, which could explain the high survival rate of these patients. CONCLUSIONS: We have identified two parameters derived from TCR repertoire analysis that could outperform PD-L1 TPS and TMB as predictive biomarkers of CPR after neoadjuvant chemoimmunotherapy, and unraveled possible mechanisms of CPR involving enhanced tumor immunogenicity and peripheral immunosurveillance.

ctDNA from body fluids is an adequate source for EGFR biomarker testing in advanced lung adenocarcinoma.

OBJECTIVES: Epidermal growth factor receptor (EGFR) biomarker testing using blood-based liquid biopsies remains challenging due to the low concentration of circulating tumor DNA (ctDNA) in certain plasma samples. The aim of this study is to evaluate the usefulness for EGFR biomarker testing of ctDNA from pleural effusions, cerebrospinal fluids, ascites and pericardial effusions obtained during the clinical management of lung adenocarcinoma patients. METHODS: For comparison purposes, 23 paired plasma and body fluid samples were collected from 17 patients with EGFR-positive lung adenocarcinoma. After circulating free DNA (cfDNA) isolation, samples were evaluated for the initial EGFR-sensitizing mutation and the p.T790M resistance mutation by array-based digital PCR (dPCR). RESULTS: Body fluids had more cfDNA than plasma samples (1.90 vs. 0.36 ng/µL; p=0.0130), and more samples tested positive for EGFR mutations (21 vs. 16 samples), with a total of 28 vs. 22 variants detected. Furthermore, mutant allele frequencies (MAFs) observed in body fluids were significantly higher than those assessed in the paired plasma samples for EGFR-sensitizing mutations (median MAFs = 15.8 vs. 0.8%; p=0.0004) as well as for the p.T790M resistance mutation (median MAFs = 8.69 vs. 0.16%; p=0.0390). Importantly, two patients who had progressed on first-generation EGFR-tyrosine kinase inhibitors with a dubious result for p.T790M plasma (MAFs = 0.11%) had an indisputably positive result in their respective body fluid samples (MAFs = 10.25 and 9.66%). CONCLUSIONS: ctDNA derived from body fluids is an informative source for EGFR biomarker testing, with greater sensitivity than plasma samples.

ctDNA analysis reveals different molecular patterns upon disease progression in patients treated with osimertinib.

BACKGROUND: Several clinical trials have demonstrated the efficacy and safety of osimertinib in advanced non-small-cell lung cancer (NSCLC). However, there is significant unexplained variability in treatment outcome. METHODS: Observational prospective cohort of 22 pre-treated patients with stage IV NSCLC harboring the epidermal growth factor receptor (EGFR) p.T790M resistance mutation and who were treated with osimertinib. Three hundred and twenty-six serial plasma samples were collected and analyzed by digital PCR (dPCR) and next-generation sequencing (NGS). RESULTS: The median progression-free survival (PFS), since the start of osimertinib, was 8.9 [interquartile range (IQR): 4.6-18.0] months. The median treatment durations of sequential gefitinib + osimertinib, afatinib + osimertinib and erlotinib + osimertinib treatments were 30.1, 24.6 and 21.1 months, respectively. The p.T790M mutation was detected in 19 (86%) pre-treatment blood samples. Undetectable levels of the original EGFR-sensitizing mutation after 3 months of treatment were associated with superior PFS (HR: 0.2, 95% CI: 0.05-0.7). Likewise, re-emergence of the original EGFR mutation, alone or together with the p.T790M mutation was significantly associated with shorter PFS (HR: 8.8, 95% CI: 1.1-70.7 and HR: 5.9, 95% CI: 1.2-27.9, respectively). Blood-based monitoring revealed three molecular patterns upon progression to osimertinib: sensitizing+/T790M+/C797S+, sensitizing+/T790M+/C797S-, and sensitizing+/T790M-/C797S-. Median time to progression in patients showing the triplet pattern (sensitizing+/T790M+/C797S+) was 12.27 months compared with 4.87 months in patients in whom only the original EGFR sensitizing was detected, and 2.17 months in patients showing the duplet pattern (sensitizing+/T790M+). Finally, we found that mutations in exon 545 of the PIK3CA gene were the most frequent alteration detected upon disease progression in patients without acquired EGFR-resistance mutations. CONCLUSIONS: Different molecular patterns identified by plasma genotyping may be of prognostic significance, suggesting that the use of liquid biopsy is a valuable approach for tumor monitoring.

Sex is a strong prognostic factor in stage IV non-small-cell lung cancer patients and should be considered in survival rate estimation.

BACKGROUND: Biological differences between the sexes have a major impact on disease and treatment outcome. In this paper, we evaluate the prognostic value of sex in stage IV non-small-cell lung cancer (NSCLC) in the context of routine clinical data, and compare this information with other external datasets. METHODS: Clinical data from stage IV NSCLC patients from Hospital Puerta de Hierro (HPH) were retrieved from electronic health records using big data analytics (N = 397). In addition, data from the Spanish Lung Cancer Group (GECP) Tumor Registry (N = 1382) and from a published study available from the cBioPortal (MSK) (N = 601) were analyzed. Survival curves were estimated using the Kaplan-Meier method. A Cox proportional hazards regression model was used to assess the prognostic value of sex. A meta-analysis to compare the outcome for males and females in terms of overall survival (OS) and progression free survival (PFS) was performed. RESULTS: The median OS time was 12 months for males and 19 months for females (overall HR = 0.77; 95% CI: 0.68-0.87; P < 0.001). Similarly, females with stage IV NSCLC harboring an EGFR-sensitizing mutation lived significantly longer than males (median OS: males, 19 months; females, 32 months) with a lower risk of death compared with males (overall HR = 0.75; 95% CI: 0.67-0.84). In addition, female patients benefited more from EGFR inhibitors in terms of PFS and OS (overall HR = 0.45; 95% CI: 0.32-0.64, and HR = 0.62; 95% CI: 0.48-0.80, respectively). Median PFS was 21 months in females and 12 months in males (P < 0.001). CONCLUSIONS: Using routine clinical data we confirmed the previous finding that among stage IV NSCLC patients, females had a significantly better prognosis than males. The effect size of the sex was notable, highlighting the fact that survival rates are usually estimated and patients are generally managed without considering the sexes separately, which may lead to suboptimal results.

Next-generation sequencing for tumor mutation quantification using liquid biopsies.

Background Non-small cell lung cancer (NSCLC) patients benefit from targeted therapies both in first- and second-line treatment. Nevertheless, molecular profiling of lung cancer tumors after first disease progression is seldom performed. The analysis of circulating tumor DNA (ctDNA) enables not only non-invasive biomarker testing but also monitoring tumor response to treatment. Digital PCR (dPCR), although a robust approach, only enables the analysis of a limited number of mutations. Next-generation sequencing (NGS), on the other hand, enables the analysis of significantly greater numbers of mutations. Methods A total of 54 circulating free DNA (cfDNA) samples from 52 NSCLC patients and two healthy donors were analyzed by NGS using the Oncomine™ Lung cfDNA Assay kit and dPCR. Results Lin's concordance correlation coefficient and Pearson's correlation coefficient between mutant allele frequencies (MAFs) assessed by NGS and dPCR revealed a positive and linear relationship between the two data sets (ρc = 0.986; 95% confidence interval [CI] = 0.975-0.991; r = 0.987; p < 0.0001, respectively), indicating an excellent concordance between both measurements. Similarly, the agreement between NGS and dPCR for the detection of the resistance mutation p.T790M was almost perfect (K = 0.81; 95% CI = 0.62-0.99), with an excellent correlation in terms of MAFs (ρc = 0.991; 95% CI = 0.981-0.992 and Pearson's r = 0.998; p < 0.0001). Importantly, cfDNA sequencing was successful using as low as 10 ng cfDNA input. Conclusions MAFs assessed by NGS were highly correlated with MAFs assessed by dPCR, demonstrating that NGS is a robust technique for ctDNA quantification using clinical samples, thereby allowing for dynamic genomic surveillance in the era of precision medicine.

Peritoneal washing is an adequate source for somatic BRCA1/2 mutation testing in ovarian malignancies.

Genetic screening for BRCA mutations should be offered to all women diagnosed with epithelial ovarian, fallopian tube, and/or peritoneal cancers given the implications for treatment options and cancer risk assessments. Yet, while germline breast cancer susceptibility gene 1 (BRCA1) and breast cancer susceptibility gene 2 (BRCA2) testing is commonly performed, BRCA1/2 somatic mutations testing is rather challenging since the poor quality of DNA extracted from formalin fixed paraffin embedded (FFPE) samples can significantly impair this process. Peritoneal lavage is routinely performed in surgeries of suspected ovarian malignancies. We have analyzed fresh tumor, peritoneal lavage and blood samples from ten patients and we have found an excellent agreement (88%) between fresh tumor and peritoneal lavage for BRCA mutation testing. Importantly, 112 of the 114 genomic alterations detected in fresh tumor samples were also found in peritoneal lavage fluids. Our data suggest that peritoneal washings can indeed streamline BRCA genes mutation testing procedures.

Prognostic value of quantitative ctDNA levels in non small cell lung cancer patients.

BACKGROUND: Circulating tumor DNA (ctDNA) levels correlate well with tumor bulk. In this paper we aim to estimate the prognostic value of the dynamic quantification of ctDNA levels. MATERIALS AND METHODS: A total of 251 serial plasma samples from 41 non-small-cell lung cancer patients who carried an activating EGFR mutation were analysed by digital PCR. For survival analysis, ctDNA levels were computed as a time-dependent covariate. RESULTS: Dynamic ctDNA measurements had prognostic significance (hazard ratio for overall survival and progression free survival according to p.T790M mutant allele frequency; 2.676 and 2.71 respectively; P < 0.05). In the same way, patients with p.T790M-negative or unchanging or decreasing plasma levels of sensitizing EGFR mutation were 12 and 4.8 times more likely to maintain response or stable disease, respectively, than patients in which the opposite occurred (P < 0.05).On average, the p.T790M mutation was detected in plasma 51 days before the assessment of progression disease by CT-scan. Finally, ctDNA outperformed CTCs for assessing tumor progression (P = 0.021). CONCLUSIONS: The appearance or increase in a unit of the p.T790M allele frequency almost triples the risk of death and progression. This information can be used to design clinical trials aiming to estimate whether T790M positive patients should start second line treatment based on molecular data rather than imaging data.

Dynamic circulating tumor DNA quantificaton for the individualization of non-small-cell lung cancer patients treatment.

BACKGROUND: Liquid biopsy has evolved from being a promising line to becoming a validated approach for biomarker testing. However, its utility for individualization of therapy has been scarcely reported. In this study, we show how monitoring levels of EGFR mutation in plasma can be useful for the individualization of treatment. RESULTS: Longitudinal EGFR mutation levels in plasma always correlated with tumor response ascertained by RECIST criteria. Moreover, decreasing EGFR mutation levels were detected in all patients benefiting from locoregional radiotherapy, whereas the opposite occurred when a patient progressed soon after radiotherapy treatment. Similarly, increasing EGFR mutation levels anticipated disease progression after TKI dose reduction, discontinuation of treatment, or reduced bioavailability due to drug interactions. In addition, EGFR mutation levels were useful to monitor treatment outcome of new therapies and constituted a decisive factor when the clinical situation of the patient did not correlate with responses ascertained by radiologist. Finally, our results indicate that cancer associated body fluids (pleural, pericardial or cerebrospinal fluid) are certainly a suitable source for biomarker testing that can extend EGFR mutation detection to biofluids other than blood. MATERIALS AND METHODS: A total of 180 serial plasma samples from 18 non-small-cell lung cancer patients who carried an activating EGFR mutation were investigated by digital PCR. CONCLUSIONS: Monitoring levels of EGFR mutation in plasma allows resolving doubts that frequently arise in daily clinical practice and constitutes a major step towards achieving personalized medicine.

Tumor burden monitoring using cell-free tumor DNA could be limited by tumor heterogeneity in advanced breast cancer and should be evaluated together with radiographic imaging.

BACKGROUND: Accurate measurement of tumor burden in breast cancer disease is essential to improve the clinical management of patients. In this study, we evaluate whether the fluctuations in the fraction of PIK3CA mutant allele correlates with tumor response according to RECIST criteria and tumor markers quantification. METHODS: Eighty six plasma samples were analyzed by digital PCR using Rare Mutation Assays for E542K, E545K and H1047R. Mutant cfDNA and tumor markers CA15-3 and CEA were compared with radiographic imaging. RESULTS: The agreement between PIK3CA mutation status in FFPE samples and circulating tumor DNA (ctDNA) was moderate (K = 0.591; 95% IC = 0.371-0.811). Restricting the analysis to the metastatic patients, we found a good agreement between PIK3CA mutation status assessed in liquid and solid biopsy (K = 0.798 95%; IC = 0.586-1). ctDNA showed serial changes with fluctuations correlating with tumor markers 15.3 and CEA in 7 out of 8 cases with Pearson correlation coefficients ranging from 0.99 to 0.46 and from 0.99 to 0.38 respectively. Similarly, fluctuations in the fraction of PIK3CA mutant allele always correlated with changes in lesion size seen on images, although in two cases it did not correlate with treatment responses as defined by RECIST criteria. CONCLUSION: oncogenic mutation quantification in plasma samples can be useful to monitor treatment outcome. However, it might be limited by tumor heterogeneity in advanced disease and it should be evaluated together with radiographic imaging.

Alternative transcripts of the SERPINA1 gene in alpha-1 antitrypsin deficiency.

BACKGROUND: SERPINA1 is the gene for alpha-1 antitrypsin (AAT), an acute phase protein with anti-protease and immunoregulatory activities. Mutations in SERPINA1 gene cause AAT deficiency and predispose individuals to early-onset emphysema and liver diseases. Expression of the SERPINA1 gene is regulated by different promoters and alternative splicing events among non-coding exons 1A, 1B and 1C. METHODS: We have developed three quantitative PCR (QT-PCR) assays (1A, 1B and 1C). These assays were applied for the analysis of SERPINA1 alternative transcripts in: (1) 16 human tissues and (2) peripheral blood leukocytes from 33 subjects with AAT mutations and 7 controls. RESULTS: Tissue-specific expression was found for the SERPINA1 transcripts. The 1A transcripts were mainly expressed in leukocytes and lung tissue while those detected with the 1B assay were highly restricted to leukocytes. Only 1B transcripts significantly correlated with serum AAT levels. The 1C transcripts were specifically found in lung, liver, kidney and pancreas. Furthermore, the expression of transcripts was related to AAT genotypes. While deficient variants of AAT had no pronounced effect on the transcript expression, null alleles were associated with significant reduction of different transcripts. CONCLUSIONS: The possibility to discriminate between SERPINA1 alternative splicing products will help us to understand better the regulation of SERPINA1 gene and its association with SERPINA1 mutations-related diseases.