ORIGINAL RESEARCH article
Sec. Thoracic Oncology
The Landscape of Actionable Genomic Alterations by Next-Generation Sequencing in Tumor Tissue Versus Circulating Tumor DNA in Chinese Patients With Non-Small Cell Lung Cancer
- 1Department of Oncology, First Affiliated Hospital of Yangtze University, Jingzhou, China
- 2Zhangjiang Center for Translational Medicine, Shanghai Biotecan Pharmaceuticals Co., Ltd., Shanghai, China
- 3Department of General Surgery, Yucheng Hospital of Traditional Chinese Medicine, Dezhou City, China
- 4Department of Oncology, Affiliated Hospital of Chengde Medical University, Chengde, China
Background: Circulating tumor DNA (ctDNA) sequence analysis shows great potential in the management of non-small cell lung cancer (NSCLC) and the prediction of drug sensitivity or resistance in many cancers. Here, we drew and compared the somatic mutational profile using ctDNA and tumor tissue sequence analysis in lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC), and assess its potential clinical value.
Methods: In this study, 221 tumor tissues and 174 plasma samples from NSCLC patients were analyzed by hybridization capture-based next-generation sequencing (NGS) panel including 95 cancer-associated genes. Tumor response assessments were applied to 137 patients with advanced-stage (III and IV) NSCLC who first received targeted agents.
Results: Twenty significantly mutated genes were identified such as TP53, EGFR, RB1, KRAS, PIK3CA, CD3EAP, CTNNB1, ERBB2, APC, BRAF, TERT, FBXW7, and HRAS. Among them, TP53 was the most frequently mutated gene and had a higher mutation probability in male (p = 0.00124) and smoking (p < 0.0001) patients. A total of 48.35% (191/395) of NSCLC patients possessed at least one actionable alteration according to the OncoKB database. Although the sensitivity of genomic profiling from ctDNA was lower than that from tumor tissue DNA, the mutational landscape of target genes from ctDNA is similar to that from tumor tissue DNA, which led to 61.22% (30/49) of mutational concordance in NSCLC. Additionally, the mutational concordance between tissue DNA and ctDNA in LUAD differs from that in LUSC, which is 63.83% versus 46.67%, indicating that NSCLC subtypes influence the specificity of mutation detection in plasma-derived ctDNA. Lastly, patients with EGFR and TP53 co-alterations showed similar responses to Gefitinib and Icotinib, and the co-occurring TP53 mutation was most likely to be a poor prognostic factor for patients receiving Gefitinib, indicating that the distributions and types of TP53 mutations may contribute to the efficacy and prognosis of molecular targeted therapy.
Conclusions: As a promising alternative for tumor genomic profiling, ctDNA analysis is more credible in LUAD than in LUSC. Genomic subtyping has strong potential in prognostication and therapeutic decision-making for NSCLC patients, which indicated the necessity for the utility of target NGS in guiding clinical management.
Genomic analysis is gaining popularity for targeted therapy in a wide variety of tumors (1). The biomarkers for targeted therapy can be identified by sequencing the tumor-derived DNA (2–4). However, although tissue sequencing is still the gold standard for genomic profiling, tumor tissues are not always sufficient for sequencing after histological diagnosis in advanced or metastatic patients mainly because of the inaccessibility of some tumor sites and/or the unsafety of invasive biopsy procedures (5). Furthermore, the single-region tissue sampling may not uncover the complete molecular profiling, since heterogeneous tumor subclones exist in many cancers (6–8).
DNA from both primary and metastatic tumor cells is released into the blood during the process of necrosis, apoptosis, and lysis (9). Thus, circulating tumor DNA (ctDNA) sequence analysis, a non-invasive liquid biopsy approach (5), has shown enormous potential in genomic profiling, such as identifying targetable alternations (10–12), disease surveillance (13, 14), and the monitor of minimal residual disease (15, 16). Particularly, ctDNA has been regarded as a potential biomarker analysis in non-small cell lung cancer (NSCLC) patients, whose plasma cell-free DNA (cfDNA) concentration is higher than that in healthy individuals (17–19).
As a common and highly heterogeneous malignant tumor, NSCLC accounts for about 85% of all cases of lung cancer (20, 21). It remains the leading cause of cancer death worldwide, and its 5-year overall survival rate for advanced or metastatic patients was <5% (22). Fortunately, disease progression has been significantly blocked, and the 5-year overall survival rate has also been notably improved on account of the utilization of actionable genomic alterations in advanced-stage NSCLC (23). For example, targeted therapeutics have contributed to the higher survival rates compared with standard conventional chemotherapy in metastatic NSCLC (1, 23–25). Next-generation sequencing (NGS) is regarded as a comprehensive approach to identifying a large number of uncommon but actionable mutations (26), which is increasingly recommended by many national and international clinical guidelines (27). So far, the molecular biomarkers recommended by clinical guidelines mainly include EGFR, ERBB2 (HER2), and KRAS activating mutations; ALK and ROS1 fusions; BRAF V600E and MET copy number gain (CNG) and exon 14 skipping (E14skip) mutations; and RET and NTRK rearrangements (28–30). However, about 20%–30% of NSCLC patients do not show any objective responses to epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs). Moreover, despite initial significant responses, the genetically driven tumor will eventually be resistant to targeted therapy, which severely limits the effectiveness of TKIs (31).
Co-occurring genetic alterations might contribute to explaining the molecular mechanisms of tumor resistance (32, 33). In several studies, a potential role of TP53 mutation is associated with poor therapeutic responses (34–36). As the “guardian of the genome” (37), p53 is a DNA binding protein, which restrains the proliferation of cells with damaged DNA (38, 39). In addition, p53 is also a tumor-suppressor protein, which is regarded as the “coordinator of the underlying processes of the hallmarks of cancer”, since its inactivation paves the path for malignancy (40). However, the impact of TP53 mutations on clinical outcomes in NSCLC patients with EGFR-mutation treated with first-line TKIs requires further elucidation (36), and the molecular mechanism among them remains largely unknown. In this study, we explored whether ctDNA analysis could serve as a credible alternative for tumor genomic profiling in lung adenocarcinoma (LUAD) or lung squamous cell carcinoma (LUSC), and whether the poor therapeutic response to EGFR-TKIs is associated with types and distributions of TP53 mutations.
Materials and Methods
Patients and Samples Collection
We recruited 395 NSCLC patients from the Department of Oncology in Affiliated Hospital of Chengde Medical University (Chengde, China) and the Department of Oncology in First Affiliated Hospital of Yangtze University (Jingzhou, China) between February 2017 and November 2019. The pathological diagnosis was verified by independent pulmonary pathologists based on the 4th edition of the World Health Organization Classification of Lung Tumors (41), and tumors with histological components other than NSCLC were excluded. Eight fresh tumor tissues, 213 formalin-fixed and paraffin-embedded (FFPE) tumor specimens, and 174 plasma were collected for NGS analysis. For the group of targeted therapy in patients with advanced-stage (IIIB and IV) NSCLC, 69 tumor tissues and 68 plasma samples were collected before the initiation of any therapies. Clinical characteristics of 395 NSCLC patients are shown in Table 1.
DNA Extraction and Quality Control
Genomic DNA (gDNA) from the fresh tumor tissues, FFPE tumor specimens, and plasma was extracted by QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), GeneRead DNA FFPE Kit (Qiagen, Hilden, Germany), and HiPure Circulating DNA Midi Kit C (Magen, Guangzhou, China), respectively. Qubit® 3.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and NanoDrop ND-1000 (Thermo Scientific, Wilmington, DE, USA) were used to assess the quantity and purity of gDNA. Fragmentation status was evaluated by the Agilent 2100 Bioanalyzer instrument (Agilent Technologies) using the High-Sensitivity DNA Reagent (Agilent Technologies, Santa Clara, CA, USA) to produce a DNA integrity number (DIN). Additionally, the step of quality control (QC) was also performed to assess FFPE DNA integrity by multiplex polymerase chain reaction (PCR). In brief, gDNA (30 ng) was amplified using three different primers of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene, the size set of which was 200–400 base pairs. An Agilent 2100 bioanalyzer instrument (Agilent Technologies) was used to determine the concentration of multiplex PCR products. The fragmentation of gDNA from FFPE was estimated by the average yield ratio (AYR) value, which was calculated by dividing the yield ratio of reference DNA (Promega Madison, WI, USA) into each amplicon.
Library Preparation and Hybridization Capture
Three hundred nanograms of gDNA from each sample determined by Qubit quantification was mechanically fragmented via an E220 focused ultrasonicator Covaris (Covaris, Woburn, MA, USA). The targeted size of the DNA fragment was between 150 and 200 bp. Then, 10–100 ng DNA was used for library construction by the KAPA library preparation kit (Kapa Biosystems Inc., Wilmington, MA, USA), which was constructed with end repair, A-tailing, and adapter ligation without additional fragmentation following the manufacturer’s instructions. Finally, the NGS libraries were captured using the xGen Lockdown Probe pool (IDT Technologies), and the captured DNA fragments were amplified with 12–13 cycles of PCR using 1× KAPA HiFi Hot Start Ready Mix.
After QC and quantification by Agilent 2100 Bioanalyzer (Agilent Technologies) and Qubit® 3.0 Fluorometer (Invitrogen; Thermo Fisher Scientific, Inc.), the NGS libraries were sequenced on an Illumina NextSeq CN500 platform (Illumina Inc, San Diego, CA, USA) Medium flux chip.
Clean data were obtained following filtering adapter, low quality reads, and reads with length <36 bp, and they were aligned to reference human genome (University of California Santa Cruz ID: hg19) using the Burrows–Wheeler Aligner v. 0.7.12. Subsequently, the Picard and Genome Analysis Toolkit (GATK v.3.2) method was used for duplicate removal, local realignment, and base quality score recalibration, and it also generated the quality statistics, such as mapped reads, mean mapping quality, and mean coverage. Finally, VarDict was adopted for the identification of SNV and InDel.
The ANNOVAR software tool was used for annotating somatic variants. The candidate somatic variants were identified by the following filter conditions: (i) remove mutations with coverage depth <10×; (ii) remove variant sites with mutant allele frequency (MAF) >0.001 in the 1,000 Genomes databases (1,000 Genomes Project Consortium; https://www.internationalgenome.org/) and remove variant sites with MAF >0.001 in the ExAC (https://ncbiinsights.ncbi.nlm.nih.gov/tag/exac/); (iii) retain variant sites with MAF ≥0.001 and <0.1 in the 1,000 Genomes databases with COSMIC evidence; (iv) retain variations in the exonic or splicing region (10 bp upstream and downstream of splicing sites); (v) remove synonymous mutations; (vi) remove Unknown Variant_Classification; and (vii) the functional benign variant sites predicted by PolyPhen-2 and MutationTaster were removed.
CNVkit was employed to identify copy number variations (CNVs). Focal copy number events were determined with Genomic Identification of Significant Targets in Cancer (GISTIC) 2.0 (threshold: q-value <0.25). The correlation between the identified somatic variants and their clinical significance was established by OncoKB Precision Oncology Database (http://oncokb.org/).
The maftools package in R software (R 3.5.1, R Core Team; https://www.RProject.org) was used to create the mutational landscape, including somatic single nucleotide variations (SNVs) and InDels. The enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) were carried out to explore the biological significance of the candidate mutant genes by ClusterProfiler package (http://bioconductor.org/packages/release/bioc/html/clusterProfiler.Html) (42).
Tumor responses (CR, complete remission; PR, partial response; SD, stable disease; and PD, progressive disease) were assessed by the Response Evaluation Criteria in Solid Tumors (RECIST version 1.1) (43–46). The disease control rate (DCR) was defined as all positive responders (including complete, partial, and stable responders) divided by the total number of response-evaluable patients. The response rate (RR) was calculated by dividing complete and partial responders into the total number of response-evaluable patients. Meanwhile, the clinic pathological characteristics and DCR of targeted therapy were compared between two groups via chi-square tests by R software and GraphPad Prism (v. 7.0; GraphPad Software, La Jolla, CA) software, respectively. Fisher exact test was used to evaluate the statistical differences in categorical variables between two groups. A two-sided p-value <0.05 was used to evaluate statistical significance. The group of fewer than five samples was not adopted to perform statistical analyses. Tumor response assessments were only applied to patients with advanced-stage (IIIB and IV) NSCLC and receiving the first-line targeted therapy.
In this retrospective study, 395 NSCLC patients were enrolled including 127 patients from the Department of Oncology in Affiliated Hospital of Chengde Medical University and 268 patients from the Department of Oncology in First Affiliated Hospital of Yangtze University. Among them, 340 patients were diagnosed as LUAD, 54 patients were diagnosed as LUSC, and 1 patient was diagnosed as large cell carcinoma. In terms of tumor staging, the number of NSCLC patients in stages I–IV was 27, 14, 47, and 255, respectively. One hundred thirty-seven patients with advanced-stage (III and IV) NSCLC were treated with targeted therapy before any other therapies. For other clinical indicators, 273 of 395 cases had lymph node metastasis pathologically, 224 of 395 patients had never smoked, and 227 were male patients. Till the last time of our follow-up, 307 patients were still alive, 47 patients had died, and others were lost to follow-up (Table 1).
Landscape of Somatic Mutations
To delineate the mutation landscape of NSCLC, we first analyzed somatic mutations from 221 tumor tissue samples by an NGS panel of 95 known cancer genes (Supplementary Table S1). A total of 455 somatic variants in 46 genes were detected in 199 of 221 (90.05%) tumor tissues, including 39 nonsense mutations, 298 missense mutations, 7 frameshift insertions, 3 frameshift deletions, 17 in-frame insertions, 79 in-frame deletions, and 12 splice sites. The top 10 frequent mutated genes were TP53 (54.30%), EGFR (48.87%), RB1 (14.48%), KRAS (8.60%), PIK3CA (8.60%), CD3EAP (7.24%), CTNNB1 (5.43%), ERBB2 (5.43%), APC (4.07%), and BRAF (2.26%) (Figure 1A; Supplementary Table S2).
Figure 1 Tumor mutation landscape of NSCLC derived from tumor tissue DNA (n = 213) (A) and ctDNA (n = 139) (B). Concordance of mutated genes derived from tumor tissue DNA or ctDNA in NSCLC patients (C). Patients were arranged along the x-axis. Significantly mutated genes identified by VarDict were ranked by mutation frequencies and clinical characteristics. Tumor mutation burden (TMB, mutations per Mb) is shown in the upper panel. The bars on the left represent the mutation frequency of each gene.
To explore the feasibility of genomic profiling by peripheral blood, the NGS panel was also applied to the analysis of 174 plasma samples. A total of 236 somatic variants in 33 genes were identified in 125 of 174 (71.84%) plasma-derived ctDNA samples, including 11 nonsense mutations, 157 missense mutations, 4 frameshift insertions, 1 frameshift deletion, 3 in-frame insertions, 51 in-frame deletions, and 9 splice sites, which was lower than that from tumor tissue DNA (Figure 1B; Supplementary Table S2). Expectedly, the mutation frequencies of the most 10 common mutated genes were also lower, which were TP53 (32.76%), EGFR (25.86%), CD3EAP (12.64%), RB1 (8.05%), PIK3CA (7.47%), KRAS (5.75%), APC (2.87%), TERT (2.87%), FBXW7 (2.30%), and HRAS (1.72%) (Figure 1B; Supplementary Table S2). In short, the percentage of ctDNA-specific mutated genes and tissue DNA-specific mutated genes were 6.12% (3/49) and 32.65% (16/49), respectively, while 61.22% (30/49) of mutated genes were shared by ctDNA analysis and tumor tissue DNA analysis (Figure 1C; Supplementary Figure S1, Supplementary Table S2).
Somatic Copy Number Alterations Detection
A total of 213 tumor tissue samples and 139 plasma samples were profiled for somatic copy number alterations (SCNAs) in NSCLC. There were 20 significant peaks of copy number gain in 213 tumor tissue samples including 1q21.3, 3q26.31, 3q26.33(PIK3CA), 5p15.33(TERT), 5p13.3, 7p22.2, 7p15.3, 7p11.2(EGFR), 7q31.2(TP53TG1), 8p11.22, 8q21.13, 8q24.21, 11q13.3, 12q14.3, 13q14.2(RB1), 13q34, 14q13.3, 17q12, 19q13.12, and 20q13.2 (Figures 2A, left, 2B; Supplementary Table S3), while there were only six significant peaks of copy number gain in 139 plasma samples including 1q23.3, 5p15.33(TERT), 5p14.1, 7p11.2(EGFR), 7q31.2, and 8q22.3(TP53INP1) (Figures 2C, left, 2D; Supplementary Table S3). Additionally, we also identified 15 significant peaks of copy number loss including 1p36.33, 1q23.1, 4p16.3, 5p15.33(TERT), 6q22.1, 7p22.3, 7q34, 8q24.3, 9q34.3, 10q11.21, 10q26.13, 11p15.5, 13q34, 16p13.3, and 20q13.33 in 213 tissue samples (Figures 2A, right, 2B; Supplementary Table S3), but only three significant peaks of copy number loss including 6q22.1, 7q34 and 10q26.13 were found in 139 plasma samples (Figures 2C, right, 2D; Supplementary Table S3). Lastly, SCNAs of LUAD (Supplementary Figure S2) and LUSC (Supplementary Figure S3) were also profiled, which further indicated that the feasibility of genomic profiling using ctDNA analysis was lower than using tumor tissue DNA analysis.
Figure 2 Somatic copy number alterations (SCNAs) of NSCLC derived from tumor tissue DNA (n = 213) (A) and ctDNA (n = 139) (C). SCNAs in NSCLC are plotted by chromosomal location (vertical axis) by CNVkit. The 10 genes in the hybridization capture-based NGS panel of 95 genes (B, D).
ctDNA Analysis Has More Potential to Be an Alternative for Genomic Profiling of LUAD Compared With Its Application in LUSC
To compare the feasibility of genomic profiling using plasma samples in NSCLC subtypes, ctDNA from 143 LUAD patients and 31 LUSC patients was sequenced by the NGS panel of 95 known cancer genes.
In the group of LUAD, a total of 417 somatic variants in 44 genes were identified in 176 of 197 (89.34%) tumor tissue samples, including 35 nonsense mutations, 271 missense mutations, 7 frameshift insertions, 3 frameshift deletions, 17 in-frame insertions, 75 in-frame deletions, and 9 splice sites (Figure 3A; Supplementary Table S4). There were 197 somatic variants in 33 genes identified in 102 of 143 (71.33%) ctDNA samples including 8 nonsense mutations, 131 missense mutations, 2 frameshift insertions, 3 in-frame insertions, 46 in-frame deletions, and 7 splice sites (Figure 3C; Supplementary Table S4). Moreover, the mutation frequencies of target genes detected by ctDNA analysis were lower than that by tumor tissue DNA analysis (Figures 3A, C; Supplementary Table S4). In summary, the percentage of ctDNA-specific mutated genes and tumor tissue DNA-specific mutated genes were 6.38% (3/47) and 29.79% (14/47), respectively, while 63.83% (30/47) of mutated genes were detected by both ctDNA analysis and tumor tissue DNA analysis (Figure 3; Supplementary Table S4).
Figure 3 Tumor mutation landscape of NSCLC derived from LUAD tumor tissue DNA (n = 197) (A), LUSC tumor tissue DNA (n = 23) (B), LUAD ctDNA (n = 143) (C), and LUSC ctDNA (n = 31) (D). Patients were arranged along the x-axis. Significantly mutated genes identified by VarDict were ranked by mutation frequencies and clinical characteristics. Tumor mutation burden (TMB, mutations per Mb) is shown in the upper panel. The bars on the left represent the mutation frequency of each gene.
In the group of LUSC, there were 36 somatic variants in 13 genes were identified in 22 of 23 (95.65%) tissue samples, including 3 nonsense mutations, 26 missense mutations, 4 in-frame deletions, and 3 splice sites (Figure 3B; Supplementary Table S4). Meanwhile, a total of 39 somatic variants in 9 genes were identified in 23 of 31 (74.19%) plasma samples including 3 nonsense mutations, 26 missense mutations, 5 in-frame deletions, 1 frameshift deletion, 2 frameshift insertions, and 2 splice sites (Figure 3D; Supplementary Table S4). Additionally, the mutation frequencies of target genes detected by ctDNA analysis were similar to that by tumor tissue DNA analysis (Figures 3B, D; Supplementary Table S4). Collectively, the percentage of ctDNA-specific mutated genes, tissue DNA-specific mutated genes, and mutated genes detected by both biopsies were 13.33% (2/15), 40% (6/15), and 46.67% (7/15), respectively (Figure 3; Supplementary Table S4). All of the results were consistent with previous studies that the genetic heterogeneity of LUSC is more complex than that in LUAD (47, 48). In short, our data indicated that ctDNA analysis is more feasible as an alternative for genomic profiling of LUAD than that of LUSC by comparing the consistency between ctDNA analysis and tumor DNA analysis.
Enrichment of Somatic Mutations by KEGG and GO Analysis
To expound on the biological function of mutated genes in NSCLC, KEGG and GO enrichment analyses were performed. We identified 104 altered signaling pathways in 199 tumor tissues samples and 109 altered signaling pathways in 125 plasma samples by KEGG enrichment analyses (Figures 4A, B). All the mutated genes involved in these pathways are shown in Table 2 and Supplementary Table S5. Meanwhile, a total of 56 altered functional terms were enriched in 199 tumor tissue samples, and 41 altered functional terms were enriched in 125 plasma samples by GO enrichment analyses (Table 3; Supplementary Table S5). Figures 4D, E show the top 10 most abundant altered functional terms according to gene counts and p-value (Table 3). Additionally, the mutation types and distributions of all mutated genes involved in the top 10 altered functional terms in 199 tumor tissue samples and 125 plasma samples were shown in 28 and 19 lollipop plots (Supplementary File 2), respectively. Collectively, we observed a high level of consistency of altered signaling pathways and altered functional terms between ctDNA analysis and tumor tissue DNA analysis, which were 83.62% (97/116) and 59.02% (36/61), respectively (Figure 4C).
Figure 4 Signaling pathways by KEGG (A, B) and functional terms by GO (D, E) enrichment of somatic mutations in tumor tissue DNA (n = 213) (A, D) and ctDNA (n = 139) (B, E). Count means the number of mutated genes enriched in this pathway. Venn diagram of signaling pathways and functional terms enriched by KEGG and GO in NSCLC, respectively (C, up, down).
Table 2 Signaling pathways enriched by KEGG analysis of somatic mutations in 199 tumor tissue samples and 125 plasma samples.
Table 3 Functional terms enriched by GO analysis of somatic mutations in 199 tumor tissue samples and 125 plasma samples.
Additionally, Figure 5 shows the top 10 altered signaling pathways and altered functional terms in LUAD (Figures 5A–F; Supplementary Table S6) and LUSC (Figures 5G–L; Supplementary Table S6) enriched by KEGG enrichment analyses (Figures 5A, B, G, H) and GO enrichment analyses (Figures 5D, E, J, K) using tumor tissues (Figures 5A, D, G, J) and ctDNA (Figures 5B, E, H, K). A high level of consistency of altered signaling pathways between ctDNA analysis and tumor DNA analysis was observed in both LUAD and LUSC, which were 82.91% (97/117) versus 89.01% (81/91) (Figures 5C, I). However, we found about 50% consistency of altered functional terms between these two kinds of samples in both LUAD and LUSC, which were 55.07% (38/69) versus 46.51% (20/43) (Figures 5F, L). Among these significantly altered functional terms, we focused on the alteration in transmembrane receptor protein tyrosine kinases (Figures 5D, E, J, K), which play key roles in a wide range of cellular processes including growth, motility, differentiation, apoptosis, and metabolism (49, 50). As in previous studies, the abnormity of receptor tyrosine kinases (RTKs) signaling is closely related to diseases, especially malignancy (51–54). Constitutive activation, such as gain-of-function mutations, genomic amplification, chromosomal rearrangements, and/or autocrine activation, serves as excellent examples to show the oncogenic properties of RTKs (49, 55).
Figure 5 Signaling pathways by KEGG (A, B, G, H) and functional terms by GO (D, E, J, K) enrichment of somatic mutations in LUAD tumor tissue DNA (n = 197) (A, D), LUAD ctDNA (n = 143) (B, E), LUSC tumor tissue DNA (n = 23) (G, J) and LUSC ctDNA (n = 31) (H, K). Count means the number of mutated genes enriched in this pathway. Venn diagram of signaling pathways (C, I) and functional terms (F, L) in LUAD (C, F) and LUSC (I, L) enriched by KEGG and GO, respectively.
Clinical Actionability for the Therapy of Targeted Agents
To assess the clinical utility of anticipative molecular profiling, all mutations were divided into different levels according to the evidence of clinical actionability in OncoKB (Supplementary Figure S4A). As standard therapeutic biomarkers, a cluster of gene mutations was approved by the FDA. In our cohort, 48.35% (191/395) of patients possess at least one actionable alteration. In the group of NSCLC, level_1 accounted for 79.21% (343/433), including missense mutation of EGFR, PIK3CA, BRAF, IDH2, IDH1, ATM, and KRAS, in-frame insertion of EGFR, and in-frame deletion of EGFR; level_2 accounted for 16.17% (70/433), including in-frame insertion of ERBB2 and missense mutation of ERBB2, PIK3CA, KRAS, BRAF, and NRAS; level_3 accounted for 1.62% (7/433), including in-frame deletion of EGFR and missense mutation of AKT1; level_4 accounted for 3.00% (13/433), including missense mutation of EGFR, FGFR3, and mTOR, and nonsense mutation of CDKN2A (Table 4; Supplementary Figures S4B, E, H; Supplementary Table S7). However, patients with LUAD have a higher percentage of actionable alterations than patients with LUSC, which were 52.06% (177/340) versus 25.93% (14/54) (Supplementary Figures S4C, D, F, G, I, J; Supplementary Table S7), indicating that they may benefit more from targeted therapy.
Targeted Therapy Response Is Related to Genomic Subtyping
Next, we wanted to explore the impact of genomic characteristics on the clinical outcomes of patients with advanced-stage (III and IV) NSCLC. Based on EGFR mutations, 46 patients were classified as the group of EGFR alterations without TP53 mutations (called EGFR), and 68 patients were classified as the group of EGFR and TP53 co-alterations (called EGFR&TP53). Among them, 13 out of 46 patients in EGFR group and 19 out of 68 patients in the EGFR&TP53 group were first treated with Gefitinib. As shown in Figure 6, the disease control rate (DCR = CR + PR + SD) in patients with EGFR alterations alone and in patients with EGFR&TP53 co-alterations was 100% (9/9) and 68.75% (11/16), respectively (p = 0.0608) (Figure 6A; Table 5). Meanwhile, the response rate (RR = CR + PR) was 44.44% (4/9) in patients with EGFR alterations alone, while only 12.5% (2/16) in patients with EGFR&TP53 co-alterations (p = 0.0726) (Figure 6A; Table 5). Importantly, the distributions of drug-resistant EGFR mutations were also explored, and the percentage of T790M mutation in patients with PR, SD, and PD was 16.67% (1/6), 14.29% (2/14), and 20% (1/5), respectively (Table 5; Supplementary Table S8). Additionally, we also sought to investigate whether the TP53 concomitant mutations are associated with the responses to different first-generation EGFR-TKIs. As shown in Figure 6B and Table 5, the DCR of Gefitinib was similar to that of Icotinib (68.75% vs 69.23%, p > 0.9778).
Figure 6 The tumor response to Gefitinib was associated with genomic subtyping in NSCLC (A). The comparative analysis of p53 mutations for patients treated with targeted agents (B). DCR, disease control rate (complete response + partial response + stable disease). PD, progressive disease; EGFR, NSCLC with EGFR but not TP53 alterations (n = 15); EGFR&TP53, NSCLC with EGFR and TP53 co-alterations (n = 32). Note: Tumor response data were only collected from patients with advanced disease (stages III and IV) who were first treated with targeted agents. The group including more than five patients was used to do statistical analyses.
Co-occurring genetic alterations might play important roles in the mechanisms of tumor responses, and it may be helpful to explain the marked diversity of individual outcomes (32, 33). As previous studies have reported, a potential role of TP53 mutation is associated with poor therapeutic responses (34–36). Figure 7 shows mutations and protein structures of p53 and EGFR in NSCLC patients who were first treated with Gefitinib, Icotinib, or either of them. TP53 mutations were detected in 11 patients who responded well to Gefitinib, including nonsense mutations (numbered BT1812060034LNCTX), and missense mutations (numbered BT1806210029LNCBP, BT1806220195LNCTV, BT1806290337LNCTB, BT1808040195LNCTV, BT1808230036LNCBP, BT1810270069LNCBP, BT1812020489LNCTS, BT1902210293LNCTV, BT1905230269LNCTV, and BT1906210049LNCBP) (Figure 7A, up). However, TP53 mutation sites and types in five patients with PD were different from those in the above patients, including nonsense mutations (numbered BT1812160473LNCTV and BT1901200062LNCTB), missense mutations (numbered BT1902140231LNCTX and BT1905110010LNCBP), and splice site (numbered BT1805130336LNCTV) (Figure 7A, down). Furthermore, we also carried out the comparative analysis of p53 mutations for patients treated with Icotinib (Figure 7B; Supplementary Table S8) because it has similar chemical structures, molecular mechanisms, and clinical curative effects with Gefitinib but cost a much lower price. Our results also indicated that patients who responded differently to Icotinib carried different mutation sites. However, the roles of TP53 mutations on clinical outcomes in NSCLC patients with EGFR-mutation treated with Gefitinib and Icotinib require further elucidation, since the molecular mechanism among them remains largely unknown.
Figure 7 The protein structure, mutation types, and mutation distributions of TP53 (A, B) and EGFP (C, D) in the patients with EGFR&TP53 co-alterations who were first treated with Gefitinib (A, C) and Icotinib (B, D).
Here, we identified somatic mutations in 221 tumor tissue DNA and 174 plasma-derived ctDNA by an NGS panel of 95 known cancer genes. Among them, 340, 54, and 1 out of 395 patients are LUAD, LUSC, and large cell carcinoma, respectively. The majority of these patients were at stage IV (64.56%, 255/395), where both inter- and intratumor heterogeneity had occurred. The most two common mutated genes in both tumor tissue DNA and ctDNA were TP53 and EGFR. As a tumor-suppressor gene, TP53 is frequently mutated in almost every type of cancer (40, 56, 57). EGFR mutation is the most frequently used marker to guide targeted therapy in NSCLC clinically (3). To sum up, four points were mainly acquired in this study. First, the feasibility of genomic profiling in NSCLC from ctDNA was lower than that from tissue DNA. Second, ctDNA analysis could serve as an alternative for genomic profiling is more credible in LUAD than that in LUSC. Third, patients with EGFR & TP53 co-alterations showed similar responses to Gefitinib and Icotinib. Lastly, TP53 concomitant mutations might be a poor prognostic factor for patients treated with Gefitinib, and the complexity of identifying TP53 mutations may contribute to the efficacy and prognosis of molecular targeted therapy in NSCLC.
One primary purpose of this study was to explore the differences of somatic mutations detected by ctDNA analysis and tumor tissue DNA analysis in NSCLC subtypes. Comparatively speaking, the sensitivity of genetic mutation detected in plasma samples was lower than that in tumor tissue samples, 71.84% (125/174) versus 90.05% (199/221), which was consistent with previous studies (58, 59). The genomic concordance between tissue DNA and ctDNA in LUAD differs from that in LUSC, which is 63.83% versus 46.67%, indicating that NSCLC subtypes influence the specificity of mutation detection in plasma-derived ctDNA (60–62). In addition, the mutation frequencies of the 10 most common mutated genes (e.g., EGFR, TP53, PIK3CA, RB1, CD3EAP, CTNNB1, ERBB2, APC, ERBB2, FGFR2) in LUAD were different from that in LUSC, further indicating that the differences in genetic mutations exist in NSCLC subtypes (63–66). To be specific, the mutation frequency of TP53 in LUSC was higher than that in LUAD, which was 69.56% versus 52.28% in tumor tissue DNA and 45.16% versus 30.07% in ctDNA (64). A similar trend was also observed in PIK3CA. On the contrary, EGFR mutations were more often observed in LUAD than that in LUSC, which is 53.30% versus 13.04% in tissue DNA, and 27.97% versus 16.13% in ctDNA (63). This trend was also suitable for the mutation frequencies of KRAS and APC. Collectively, the discovery of genomic profiling among NSCLC subtypes (LUAD, LUSC, and large-cell carcinoma) could also benefit a lot from ctDNA.
Another primary purpose of this study was to explore the impact of TP53 mutations on clinical outcomes in NSCLC patients with EGFR-mutation treated with first-line TKIs, and the DCR in patients with EGFR&TP53 co-alterations was lower than that in patients with EGFR mutations (p = 0.0608). Our results were consistent with previous studies that patients with TP53 concomitant mutations had a lower ORR and shorter PFS of EGFR-TKI treatment compared to patients without TP53 alterations in advanced NSCLC, and none of these studies reached statistical significance (34, 35, 67–69). Meanwhile, Qin et al. showed that concurrent TP53 mutations could significantly reduce the responses to first-line EGFR-TKIs and were related to worse prognosis in advanced NSCLC in a meta-analysis (70). However, we were not able to confirm whether the poor prognosis is directly associated with the types or the distributions of TP53 mutations, and this association required further elucidation, since the molecular mechanism among them remains largely unknown. Overall, growing evidence has indicated that the observed clinical heterogeneity in NSCLC patients with EGFR mutation may be partly due to the co-occurring molecular events. These concomitant mutations or the measurements of p53 protein levels may be used for predicting clinical outcomes in EGFR-mutant patients treated with targeted therapy in the future.
Liquid biopsies are a potential way to identify actionable biomarkers when they are difficult to acquire from tumor tissue samples clinically. This is mostly due to the insufficiency of available tissue, the failure of quality control, and the inadequate depth or breadth of analysis by a biomarker panel. In this present study, plasma-derived ctDNA increases detection rates of driver oncogene by 14.43%, and 57 additional patients get the opportunity to receive targeted therapy and/or participate in clinical trials. Taken together, our finding is consistent with previous studies that the identification of targetable alterations is growing in patients with the aid of plasma-derived ctDNA (71, 72). However, our study has several limitations: (i) this is a retrospective study, and its database has a certain limitation in investigating other sources of potential bias; (ii) only two institutions of accrual were included; (iii) NGS data were obtained from tumor tissues and plasma from different NSCLC patient cohorts and without matched germline DNA; (iv) survival curve analysis was not performed; (v) NSCLC patients for erlotinib response were less than 5; and (vi) the generalizability of our conclusions should be considered carefully.
In conclusion, our study suggests that ctDNA analysis was regarded as an alternative for genomic profiling in advanced or metastatic patients and is more credible in LUAD than that in LUSC. NSCLC patients with EGFR&TP53 co-alterations showed similar responses to Gefitinib and Icotinib, and TP53 concomitant mutations might be a poor prognostic factor for patients with EGFR-mutation treated with Gefitinib. However, comprehensive molecular mechanisms about concurrent TP53 mutations on clinical outcomes are dependent on further studies.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA801361.
The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Yangtze University, the Ethics Committee of Affiliated Hospital of Chengde Medical University. The patients/participants provided their written informed consent to participate in this study.
JC and HJ drafted the article and reviewed the submitted version of the manuscript. JC, MW, NL, CZ, HD, and DW worked on the acquisition of data. SW, JZ, HX, SL, and HJ worked on the analysis and interpretation of data. QL, JZ, SW, HJ, YX, XY and HX critically revised the article and worked on the statistical analysis. JL, JZ, and HJ are responsible for administrative/technical/material support and study supervision. QL, JZ, JL, and SW worked on conception and design. JL and QL approved the final version of the manuscript on behalf of all authors. All authors contributed to the article and approved the submitted version.
Conflict of Interest
Authors HJ, XY, YX, HX, SW, JL and JZ were employed by Shanghai Biotecan Pharmaceuticals Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
We thank all the patients and investigators involved in this study.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.751106/full#supplementary-material
2. Eberhard DA, Johnson BE, Amler LC, Goddard AD, Heldens SL, Herbst RS, et al. Mutations in the Epidermal Growth Factor Receptor and in KRAS Are Predictive and Prognostic Indicators in Patients With Non-Small-Cell Lung Cancer Treated With Chemotherapy Alone and in Combination With Erlotinib. J Clin Oncol (2005) 23(25):5900–9. doi: 10.1200/JCO.2005.02.857
4. Cappuzzo F, Hirsch FR, Rossi E, Bartolini S, Ceresoli GL, Bemis L, et al. Epidermal Growth Factor Receptor Gene and Protein and Gefitinib Sensitivity in Non-Small-Cell Lung Cancer. J Natl Cancer Inst (2005) 97(9):643–55. doi: 10.1093/jnci/dji112
7. Weber B, Meldgaard P, Hager H, Wu L, Wei W, Tsai J, et al. Detection of EGFR Mutations in Plasma and Biopsies From Non-Small Cell Lung Cancer Patients by Allele-Specific PCR Assays. BMC Cancer (2014) 14:294. doi: 10.1186/1471-2407-14-294
9. Jahr S, Hentze H, Englisch S, Hardt D, Fackelmayer FO, Hesch RD, et al. DNA Fragments in the Blood Plasma of Cancer Patients: Quantitations and Evidence for Their Origin From Apoptotic and Necrotic Cells. Cancer Res (2001) 61(4):1659–65.
10. Sundaresan TK, Sequist LV, Heymach JV, Riely GJ, Janne PA, Koch WH, et al. Detection of T790M, the Acquired Resistance EGFR Mutation, by Tumor Biopsy Versus Noninvasive Blood-Based Analyses. Clin Cancer Res (2016) 22(5):1103–10. doi: 10.1158/1078-0432.CCR-15-1031
11. Murtaza M, Dawson SJ, Tsui DW, Gale D, Forshew T, Piskorz AM, et al. Non-Invasive Analysis of Acquired Resistance to Cancer Therapy by Sequencing of Plasma DNA. Nature (2013) 497(7447):108–12. doi: 10.1038/nature12065
12. Tabernero J, Lenz HJ, Siena S, Sobrero A, Falcone A, Ychou M, et al. Analysis of Circulating DNA and Protein Biomarkers to Predict the Clinical Activity of Regorafenib and Assess Prognosis in Patients With Metastatic Colorectal Cancer: A Retrospective, Exploratory Analysis of the CORRECT Trial. Lancet Oncol (2015) 16(8):937–48. doi: 10.1016/S1470-2045(15)00138-2
13. Murtaza M, Dawson SJ, Pogrebniak K, Rueda OM, Provenzano E, Grant J, et al. Multifocal Clonal Evolution Characterized Using Circulating Tumour DNA in a Case of Metastatic Breast Cancer. Nat Commun (2015) 6:8760. doi: 10.1038/ncomms9760
14. Bettegowda C, Sausen M, Leary RJ, Kinde I, Wang Y, Agrawal N, et al. Detection of Circulating Tumor DNA in Early- and Late-Stage Human Malignancies. Sci Transl Med (2014) 6(224):224ra24. doi: 10.1126/scitranslmed.3007094
15. Garcia-Murillas I, Schiavon G, Weigelt B, Ng C, Hrebien S, Cutts RJ, et al. Mutation Tracking in Circulating Tumor DNA Predicts Relapse in Early Breast Cancer. Sci Transl Med (2015) 7(302):302ra133. doi: 10.1126/scitranslmed.aab0021
16. Tie J, Wang Y, Tomasetti C, Li L, Springer S, Kinde I, et al. Circulating Tumor DNA Analysis Detects Minimal Residual Disease and Predicts Recurrence in Patients With Stage II Colon Cancer. Sci Transl Med (2016) 8(346):346ra92. doi: 10.1126/scitranslmed.aaf6219
17. van der Drift MA, Hol BE, Klaassen CH, Prinsen CF, van Aarssen YA, Donders R, et al. Circulating DNA is a non-Invasive Prognostic Factor for Survival in Non-Small Cell Lung Cancer. Lung Cancer (2010) 68(2):283–7. doi: 10.1016/j.lungcan.2009.06.021
18. Lee YJ, Yoon KA, Han JY, Kim HT, Yun T, Lee GK, et al. Circulating Cell-Free DNA in Plasma of Never Smokers With Advanced Lung Adenocarcinoma Receiving Gefitinib or Standard Chemotherapy as First-Line Therapy. Clin Cancer Res (2011) 17(15):5179–87. doi: 10.1158/1078-0432.CCR-11-0400
19. Couraud S, Vaca-Paniagua F, Villar S, Oliver J, Schuster T, Blanche H, et al. Noninvasive Diagnosis of Actionable Mutations by Deep Sequencing of Circulating Free DNA in Lung Cancer From Never-Smokers: A Proof-of-Concept Study From BioCAST/IFCT-1002. Clin Cancer Res (2014) 20(17):4613–24. doi: 10.1158/1078-0432.CCR-13-3063
20. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68(6):394–424. doi: 10.3322/caac.21492
24. Mitsudomi T, Morita S, Yatabe Y, Negoro S, Okamoto I, Tsurutani J, et al. Gefitinib Versus Cisplatin Plus Docetaxel in Patients With Non-Small-Cell Lung Cancer Harbouring Mutations of the Epidermal Growth Factor Receptor (WJTOG3405): An Open Label, Randomised Phase 3 Trial. Lancet Oncol (2010) 11(2):121–8. doi: 10.1016/S1470-2045(09)70364-X
25. Yang JC, Wu YL, Schuler M, Sebastian M, Popat S, Yamamoto N, et al. Afatinib Versus Cisplatin-Based Chemotherapy for EGFR Mutation-Positive Lung Adenocarcinoma (LUX-Lung 3 and LUX-Lung 6): Analysis of Overall Survival Data From Two Randomised, Phase 3 Trials. Lancet Oncol (2015) 16(2):141–51. doi: 10.1016/S1470-2045(14)71173-8
26. Malapelle U, Pisapia P, Rocco D, Smeraglio R, di Spirito M, Bellevicine C, et al. Next Generation Sequencing Techniques in Liquid Biopsy: Focus on Non-Small Cell Lung Cancer Patients. Transl Lung Cancer Res (2016) 5(5):505–10. doi: 10.21037/tlcr.2016.10.08
27. Lindeman NI, Cagle PT, Aisner DL, Arcila ME, Beasley MB, Bernicker EH, et al. Updated Molecular Testing Guideline for the Selection of Lung Cancer Patients for Treatment With Targeted Tyrosine Kinase Inhibitors: Guideline From the College of American Pathologists, the International Association for the Study of Lung Cancer, and the Association for Molecular Pathology. J Mol Diagn (2018) 20(2):129–59. doi: 10.1016/j.jmoldx.2017.11.004
28. Rolfo C, Mack PC, Scagliotti GV, Baas P, Barlesi F, Bivona TG, et al. Liquid Biopsy for Advanced Non-Small Cell Lung Cancer (NSCLC): A Statement Paper From. J Thorac Oncol (2018) 13(9):1248–68. doi: 10.1016/j.jtho.2018.05.030
30. Zill OA, Banks KC, Fairclough SR, Mortimer SA, Vowles JV, Mokhtari R, et al. The Landscape of Actionable Genomic Alterations in Cell-Free Circulating Tumor DNA From 21,807 Advanced Cancer Patients. Clin Cancer Res (2018) 24(15):3528–38. doi: 10.1158/1078-0432.CCR-17-3837
32. Hong S, Gao F, Fu S, Wang Y, Fang W, Huang Y, et al. Concomitant Genetic Alterations With Response to Treatment and Epidermal Growth Factor Receptor Tyrosine Kinase Inhibitors in Patients With EGFR-Mutant Advanced Non-Small Cell Lung Cancer. JAMA Oncol (2018) 4(5):739–42. doi: 10.1001/jamaoncol.2018.0049
33. Blakely CM, Watkins TBK, Wu W, Gini B, Chabon JJ, McCoach CE, et al. Evolution and Clinical Impact of Co-Occurring Genetic Alterations in Advanced-Stage EGFR-Mutant Lung Cancers. Nat Genet (2017) 49(12):1693–704. doi: 10.1038/ng.3990
34. VanderLaan PA, Rangachari D, Mockus SM, Spotlow V, Reddi HV, Malcolm J, et al. Mutations in TP53, PIK3CA, PTEN and Other Genes in EGFR Mutated Lung Cancers: Correlation With Clinical Outcomes. Lung Cancer (2017) 106:17–21. doi: 10.1016/j.lungcan.2017.01.011
35. Canale M, Petracci E, Delmonte A, Chiadini E, Dazzi C, Papi M, et al. Impact of TP53 Mutations on Outcome in EGFR-Mutated Patients Treated With First-Line Tyrosine Kinase Inhibitors. Clin Cancer Res (2017) 23(9):2195–202. doi: 10.1158/1078-0432.CCR-16-0966
36. Yu HA, Suzawa K, Jordan E, Zehir A, Ni A, Kim R, et al. Concurrent Alterations in EGFR-Mutant Lung Cancers Associated With Resistance to EGFR Kinase Inhibitors and Characterization of MTOR as a Mediator of Resistance. Clin Cancer Res (2018) 24(13):3108–18. doi: 10.1158/1078-0432.CCR-17-2961
39. Powell DJ, Hrstka R, Candeias M, Bourougaa K, Vojtesek B, Fahraeus R. Stress-Dependent Changes in the Properties of P53 Complexes by the Alternative Translation Product P53/47. Cell Cycle (2008) 7(7):950–9. doi: 10.4161/cc.7.7.5626
41. Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, et al. The 2015 World Health Organization Classification of Lung Tumors: Impact of Genetic, Clinical and Radiologic Advances Since the 2004 Classification. J Thorac Oncol (2015) 10(9):1243–60. doi: 10.1097/JTO.0000000000000630
43. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, et al. New Response Evaluation Criteria in Solid Tumours: Revised RECIST Guideline (Version 1.1). Eur J Cancer (2009) 45(2):228–47. doi: 10.1016/j.ejca.2008.10.026
44. Niu J, Maurice-Dror C, Lee DH, Kim DW, Nagrial A, Voskoboynik M, et al. First-In-Human Phase 1 Study of the Anti-TIGIT Antibody Vibostolimab as Monotherapy or With Pembrolizumab for Advanced Solid Tumors, Including non-Small Cell Lung Cancer. Ann Oncol (2022) 33(2):169–80. doi: 10.1016/j.annonc.2021.11.002
45. Le X, Cornelissen R, Garassino M, Clarke JM, Tchekmedyian N, Goldman JW, et al. Poziotinib in Non-Small-Cell Lung Cancer Harboring HER2 Exon 20 Insertion Mutations After Prior Therapies: ZENITH20-2 Trial. J Clin Oncol (2021), JCO2101323. doi: 10.1200/JCO.21.01323
46. Wu X, Zhu J, Wang J, Lin ZQ, Yin R, Sun W, et al. Pamiparib Monotherapy for Patients With Germline BRCA1/2-Mutated Ovarian Cancer Previously Treated With at Least Two Lines of Chemotherapy: A Multicenter, Open-Label, Phase II Study. Clin Cancer Res (2021). doi: 10.1158/1078-0432.CCR-21-1186
47. Qin N, Li Y, Wang C, Zhu M, Dai J, Hong T, et al. Comprehensive Functional Annotation of Susceptibility Variants Identifies Genetic Heterogeneity Between Lung Adenocarcinoma and Squamous Cell Carcinoma. Front Med (2021) 15(2):275–91. doi: 10.1007/s11684-020-0779-4
48. Grob TJ, Kannengiesser I, Tsourlakis MC, Atanackovic D, Koenig AM, Vashist YK, et al. Heterogeneity of ERBB2 Amplification in Adenocarcinoma, Squamous Cell Carcinoma and Large Cell Undifferentiated Carcinoma of the Lung. Mod Pathol (2012) 25(12):1566–73. doi: 10.1038/modpathol.2012.125
52. Arjona D, Bello MJ, Alonso ME, Aminoso C, Isla A, De Campos JM, et al. Molecular Analysis of the EGFR Gene in Astrocytic Gliomas: mRNA Expression, Quantitative-PCR Analysis of Non-Homogeneous Gene Amplification and DNA Sequence Alterations. Neuropathol Appl Neurobiol (2005) 31(4):384–94. doi: 10.1111/j.1365-2990.2005.00653.x
53. Lopez-Gines C, Gil-Benso R, Ferrer-Luna R, Benito R, Serna E, Gonzalez-Darder J, et al. New Pattern of EGFR Amplification in Glioblastoma and the Relationship of Gene Copy Number With Gene Expression Profile. Mod Pathol (2010) 23(6):856–65. doi: 10.1038/modpathol.2010.62
55. McDonell LM, Kernohan KD, Boycott KM, Sawyer SL. Receptor Tyrosine Kinase Mutations in Developmental Syndromes and Cancer: Two Sides of the Same Coin. Hum Mol Genet (2015) 24(R1):R60–6. doi: 10.1093/hmg/ddv254
56. Hong B, van den Heuvel AP, Prabhu VV, Zhang S, El-Deiry WS. Targeting Tumor Suppressor P53 for Cancer Therapy: Strategies, Challenges and Opportunities. Curr Drug Targets (2014) 15(1):80–9. doi: 10.2174/1389450114666140106101412
57. Kaur RP, Vasudeva K, Kumar R, Munshi A. Role of P53 Gene in Breast Cancer: Focus on Mutation Spectrum and Therapeutic Strategies. Curr Pharm Des (2018) 24(30):3566–75. doi: 10.2174/1381612824666180926095709
58. Zhao X, Han RB, Zhao J, Wang J, Yang F, Zhong W, et al. Comparison of Epidermal Growth Factor Receptor Mutation Statuses in Tissue and Plasma in Stage I-IV non-Small Cell Lung Cancer Patients. Respiration (2013) 85(2):119–25. doi: 10.1159/000338790
59. Ohira T, Sakai K, Matsubayashi J, Kajiwara N, Kakihana M, Hagiwara M, et al. Tumor Volume Determines the Feasibility of Cell-Free DNA Sequencing for Mutation Detection in Non-Small Cell Lung Cancer. Cancer Sci (2016) 107(11):1660–6. doi: 10.1111/cas.13068
61. Odegaard JI, Vincent JJ, Mortimer S, Vowles JV, Ulrich BC, Banks KC, et al. Validation of a Plasma-Based Comprehensive Cancer Genotyping Assay Utilizing Orthogonal Tissue- and Plasma-Based Methodologies. Clin Cancer Res (2018) 24(15):3539–49. doi: 10.1158/1078-0432.CCR-17-3831
62. Leighl NB, Page RD, Raymond VM, Daniel DB, Divers SG, Reckamp KL, et al. Clinical Utility of Comprehensive Cell-Free DNA Analysis to Identify Genomic Biomarkers in Patients With Newly Diagnosed Metastatic Non-Small Cell Lung Cancer. Clin Cancer Res (2019) 25(15):4691–700. doi: 10.1158/1078-0432.CCR-19-0624
65. Plagnol V, Woodhouse S, Howarth K, Lensing S, Smith M, Epstein M, et al. Analytical Validation of a Next Generation Sequencing Liquid Biopsy Assay for High Sensitivity Broad Molecular Profiling. PloS One (2018) 13(3):e0193802. doi: 10.1371/journal.pone.0193802
66. Veldore VH, Choughule A, Routhu T, Mandloi N, Noronha V, Joshi A, et al. Validation of Liquid Biopsy: Plasma Cell-Free DNA Testing in Clinical Management of Advanced Non-Small Cell Lung Cancer. Lung Cancer (Auckl) (2018) 9:1–11. doi: 10.2147/LCTT.S147841
67. Kim Y, Lee B, Shim JH, Lee SH, Park WY, Choi YL, et al. Concurrent Genetic Alterations Predict the Progression to Target Therapy in EGFR-Mutated Advanced NSCLC. J Thorac Oncol (2019) 14(2):193–202. doi: 10.1016/j.jtho.2018.10.150
68. Deben C, Deschoolmeester V, Lardon F, Rolfo C, Pauwels P. TP53 and MDM2 Genetic Alterations in Non-Small Cell Lung Cancer: Evaluating Their Prognostic and Predictive Value. Crit Rev Oncol Hematol (2016) 99:63–73. doi: 10.1016/j.critrevonc.2015.11.019
69. Labbe C, Cabanero M, Korpanty GJ, Tomasini P, Doherty MK, Mascaux C, et al. Prognostic and Predictive Effects of TP53 Co-Mutation in Patients With EGFR-Mutated Non-Small Cell Lung Cancer (NSCLC). Lung Cancer (2017) 111:23–9. doi: 10.1016/j.lungcan.2017.06.014
70. Qin K, Hou H, Liang Y, Zhang X. Prognostic Value of TP53 Concurrent Mutations for EGFR- TKIs and ALK-TKIs Based Targeted Therapy in Advanced Non-Small Cell Lung Cancer: A Meta-Analysis. BMC Cancer (2020) 20(1):328. doi: 10.1186/s12885-020-06805-5
71. Aggarwal C, Thompson JC, Black TA, Katz SI, Fan R, Yee SS, et al. Clinical Implications of Plasma-Based Genotyping With the Delivery of Personalized Therapy in Metastatic Non-Small Cell Lung Cancer. JAMA Oncol (2019) 5(2):173–80. doi: 10.1001/jamaoncol.2018.4305
72. Mack PC, Banks KC, Espenschied CR, Burich RA, Zill OA, Lee CE, et al. Spectrum of Driver Mutations and Clinical Impact of Circulating Tumor DNA Analysis in Non-Small Cell Lung Cancer: Analysis of Over 8000 Cases. Cancer (2020) 126(14):3219–28. doi: 10.1002/cncr.32876
Keywords: NSCLC, tissue, ctDNA, genomic subtyping, targeted therapy response
Citation: Cai J, Jiang H, Li S, Yan X, Wang M, Li N, Zhu C, Dong H, Wang D, Xu Y, Xie H, Wu S, Lou J, Zhao J and Li Q (2022) The Landscape of Actionable Genomic Alterations by Next-Generation Sequencing in Tumor Tissue Versus Circulating Tumor DNA in Chinese Patients With Non-Small Cell Lung Cancer. Front. Oncol. 11:751106. doi: 10.3389/fonc.2021.751106
Received: 31 July 2021; Accepted: 30 December 2021;
Published: 22 February 2022.
Edited by:Mohamed Rahouma, Weill Cornell Medical Center, United States
Reviewed by:Girindra Raval, Augusta University, United States
Mohamed Emam Sobeih, Cairo University, Egypt
Copyright © 2022 Cai, Jiang, Li, Yan, Wang, Li, Zhu, Dong, Wang, Xu, Xie, Wu, Lou, Zhao and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work