Yuan Qiu1,2†Liping Liu1,2†
Haihong Yang1,2Hanzhang Chen1,2Qiuhua Deng1,2Dakai Xiao1,2Yongping Lin1,2Changbin Zhu3Weiwei Li3
Di Shao3Wenxi Jiang3
Kui Wu4
Jianxing He1,2*- 1National Clinical Research Center of Respiratory Disease, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
- 2State Key Laboratory of Respiratory Diseases, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
- 3BGI Genomics, BGI-Shenzhen, Shenzhen, China
- 4BGI-Shenzhen, Shenzhen, China
Background: Differentiating multiple pulmonary lesions as multiple primary lung cancer (MLC) or intra-pulmonary metastasis (IPM) is critical. Lung cancer also has a high genetic heterogeneity, which influenced the treatment strategy. Genetic information may aid in tracing lineage information on multiple lung lesions. This study applied comprehensive genomic profiling to decipher the intrinsic genetics of multiple lung lesions.
Methods: Sixty-six lung adenocarcinomas (LUAD) tumor lesions (FFEP) archived from 30 patients were included in this study. The 508 cancer-related genes were evaluated by targeted next-generation sequencing (MGI-seq 2000).
Results: The study included a total of 30 LUADs (66 samples). The majority of tumors demonstrated intra-tumoral heterogeneity. Two hundred twenty-four mutations were detected by sequencing the 66 samples. We investigated the driver gene mutations of NSCLC patients with multiple lesions. EGFR was the most frequently (48/198) mutated driver gene. The codons in EGFR mainly affected by mutations were p.L858R (18/66 [27.3%]) and exon 19del (8/66 [12.1%]). In addition, additional driver genes were found, including TP53, BRAF, ERBB2, MET, and PIK3CA. We also found that the inter-component heterogeneity of different lesions and more than two different mutation types of EGFR were detected in seven patients with two lesions (P3, P10, P24, P25, P28, P29, and P30). The TMB values of different lesions in each patient were different in 26 patients (except P4, P5, P14, and P30).
Conclusions: Comprehensive genomic profiling should be applied to distinguishing the nature of multiple lung lesions irrespective of radiologic and histologic diagnoses.
Introduction
The focus of lung cancer treatment shifted significantly with the identification of specific targetable driver mutations. Epidermal growth factor receptor (EGFR) driver mutation cancers represent a distinct subset of non-small-cell lung cancer (NSCLC) with broad molecular and clinical heterogeneity (1). For NSCLCs harboring EGFR driver mutations, the current standard of treatment in the first-line setting is an EGFR tyrosine kinase inhibitor (TKI), either a first- (gefitinib or erlotinib) or second-generation TKI (afatinib). Besides EGFR, ALK, KRAS, ROS1, c-MET, and PIK3CA have been implicated as a driver of NSCLCs, such as crizotinib for ALK‐positive lung cancer patients (1, 2). The TP53 gene mutation has also been identified as a driver mutation in NSCLC. These driver mutations are highly heterogeneous, including inter-patient, and intra- and inter-tumor variability. In particular, differentiating multiple pulmonary lesions as multiple primary lung cancer (MLC) or intra-pulmonary metastasis (IPM) is critical (3). Furthermore, a high degree of genetic diversity between the primary lung tumor and corresponding metastatic lesions could play a pivotal role in the therapeutic context of lung cancer patients. Beyond heterogeneity of druggable driver mutations, previous studies have analyzed the presence of mutational signatures across human cancer types, proving that specific mutational signatures correlate with defined tumors (3, 4). The genomics from NSCLC to SCLC have been reported, and molecular characterization of SCLC has revealed an extremely high mutational rate in TP53 and RB1 genes (5).
In contrast, the tumor mutation burden (TMB) reveals the total number of mutations occurring in a tumor specimen and indicates the status of genomic mutations (6). Previous studies have shown that PD-L expression is highly heterogeneous at different locations in the lung. The PD-L1 positive agreement rate between lung primary and metastatic tumors is 63-100%. In 12%-35% of lung cancer patients, PD-L1 expression changes during treatment (7). From driver genes to TMB, what are the genetic characteristics of patients with multiple primary lesions of lung adenocarcinoma?
Multiple somatic alterations lie at the root of cancer development and tumor heterogeneity. This tumor heterogeneity further complicates the design of strategies for effective treatment. Thus, knowledge about the distribution of driver mutations in NSCLCs, particularly in early-stage NSCLCs with multiple lesions, is an area of interest. This study applied comprehensive genomic profiling deciphering intrinsic genetics of multiple lung lesions.
Methods
Patients and samples
Between May 2019 and September 2019, 30 patients (many patients with multiple lesions, and only samples meeting the criteria were eventually processed for sequencing) with at least 2 lesions (30) undergoing concurrent surgery were enrolled and sequenced using a targeted exome capture sequencing (568-gene panel) on an MGI-500 platform. Of note, for a patient with multiple lesions, each lesion was evaluated separately, and the outcome of each lesion was reported.
Comprehensive genomic profiling
Genomic DNA was extracted from FFPE and peripheral blood samples using the Qiagen DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s recommendations. First, raw data generated by the BGISEQ-500 sequencer were filtered by SOAPnuke to exclude reads with low quality. Clean reads were then mapped to the reference human genome (GRCh37/hg19) from a UCSC genome browser. Calling of single nucleotide variants (SNVs) and small insertions/deletions (Ins/Del) was performed with the Genome Analysis Tool kit (GATK) using parameters adapted to HaloPlex-generated sequences. The copy number variants (CNVs) were called using the CNVnator read-depth algorithm. TMB assessed via targeted sequencing of ~1.25Mb broadly recapitulated previous results of the whole exome TMB analysis. Tumor mutation burden was the number of all the non-synonymous mutations/0.7 Mb targeted coding region.
Statistical analysis
A correlation graph was made using the R Package (version 3.3.0; http://www.r-project.org). A chi‐squared test was used for comparisons of categorical variables across multiple factors. A P< 0.05 was considered significant.
Results
Patient characteristics
All 30 LUAD patients had multifocal lesions, which added up to a total of 66 lesions. The clinical features of patients are summarized in Tables 1, 2. The age distribution was between 30 and 69 years. Twenty-two (73.3%) patients were female. Approximately 73.3% (n=22) of the study group were female compared to 26.7% (n=8) male (Table 1). Twenty-seven (90.0%) patients had stage IA and three individuals had stage IB and IIB (Table 1). Twenty-nine patients (96.7%) had > 2 lung primary lesions, and 1 patient (3.3%) had a metastatic lung lesion. Except for P1, P10, P24, and P28, the pathologic subtypes of all lesions were the same in every patient (Table 2 and Figure 1).
Table 1 Clinical characteristics of the patients at baseline.
Table 2 The clinical and driver mutation genes characteristics of 30 patients.
Figure 1 Heatmap showing mutations for each patient.
Genomic alterations in early-stage NSCLC with multiple lesions
The genomic alterations of each lesion also were reported and the mutations are summarized in Figure 1. Altogether, 224 mutations were detected using a 508 gene panel by sequencing the 66 samples. The most common mutations in our study included mutations in EGFR (76.7%) and TP53 (26.7%). Interestingly, we did not find EML4-ALK fusion and MET mutations. The codons in EGFR mainly affected by mutations were p.L858R (18/66 [27.3%]) and exon 19del (8/66 [12.1%]). To explore the genomic origin, we further investigated the shared driver mutations of 30 cases (Figure 1).
Four patients (P5, P6, P16, and P26) had shared variations, which occurred in EGFR (n=2), TP53 (n=1), 1 HSP90AA1 (n=1), and 1 KRAS (n=1), and the pathologic outcomes of all lesions were the same in every patient (Figure 1). It is worth noting that one of the lesions in P26 was primary, and the other was a metastatic lesion to the lungs (Figure 1). In 12 patients (P1, P2, P3, P4, p7, P8, P9, P12, P13, P21, P22, and P27), the pathologic results of the lesion were the same in each patient, but the driver genes were different (Figure 1). We also found that the inter-component heterogeneity of different lesions and more than two different mutations type of EGFR was detected in seven patients with two lesions (P3, P10, P24, P25, P28, P29, and P30) (Figure 1). In P10, P24, and P28, the micro-invasive lesions were found synchronously and the other lesions were subdivided as follows: invasive, Atypical adenomatous hyperplasia (AAH); and in situ. Among the 4 lesions in P14 patients, driver mutations (EGFR L747_P753delins) were detected in two lesions with different pathologic subtypes, and CHD1 and PIK3R1 mutations occurred in two other lesions. In addition, additional driver genes were found, including ALK, BRAF, ERBB2, MET, and PIK3CA. We also observed that MED12, TP53, and TSC1 gene mutations were detected in patients P4, P8, and P1.
Based on the above results, we analyzed known high-frequency driver mutations and evaluated the maximum and minimum detection rates in different situations (Table 2). The results clearly showed that increasing the lesions in one patient increased the mutation rate.
The distribution of TMB in early-stage NSCLC with multiple lesions
To investigate the intratumor heterogeneity of TMB, we analyzed the TMB status of 30 patients. The clinical features of patients are summarized in Table 3.
Table 3 The clinical and driver mutation genes characteristics of the two groups.
We found that the TMB values of different lesions in each patient were different in 26 (except P4, P5, P14, and P30; Table 3). Among the 26 patients, 4 (P1, P10, P24, and P28; Table 3) with 2 lesions, and all lesions in each patient had a different pathologic subtype. It is worth noting that the sampling locations of each lesion in P1 and P28 were also different (Table 3). In P5 and P14 with 3 lesions, the distribution of TMB is also more diverse. Two lesions in N5 patients (P5-1 and P5-3) had the same TMB value, stage, and pathologic classification. In the remaining lesion (P5-2), a different TMB value and sampling location were observed (Table 3). A similar phenomenon had also been found in P14. The lesions with different pathologies (P14-1 and P14-2) had the same TMB value, and the lesions with different TMB values (P14-1 and P14-3) had different TMB values. The clinical features and TMB values of all lesions in each patient (P4 and P30) were the same. Among two metastatic lesions from one patient (P1), highly heterogeneous TMB values were found in two lesions with the same driver gene height (Figure 1 and Table 3).
Discussion
Tumor heterogeneity is frequently cited as a confounding factor and limitation in molecular studies of tumors (8). The intra-tumor heterogeneity of multiple primary lung cancers (MPLC) may lead to therapy failure and cancer progression (9, 10). Further study of the molecular characteristics and TMB distribution of MPLC can help guide clinical practice effectively.
Previous studies have demonstrated that approximately 85% of individuals with LUAD have known driver mutations, including EGFR, KRAS, ALK, ERBB2, ROS1, RET, MET, BRAF, NRAS, and TP53 (11, 12). Our study showed that the main common driver gene was EGFR. The available data also showed that early‐ and advanced‐stage LUAD exhibit the same EGFR mutation frequencies and types (13). The major driver mutations other than EGFR, KRAS, ERBB2, PIK3CA, and TP53 were not detected in our study, which more frequently occurs in Asian LUAD (12, 14). This may be because ALK fusion gene-positive LUAD typically occurs in young non-smokers or rare smokers (15). METex14 alterations are enriched in sarcomatoid histologies, with a prevalence ranging from 8%-22% (16, 17). The ROS1 gene was reported to be present in 1%–2% of NSCLCs, which is common in young non-smoking female patients with LUAD (17). These findings indicate that known driver gene mutations dominate early genetic events and conferred a selective growth advantage. These mutations occurred in multiple primary and metastatic lesions of the same patient. Interestingly, if only one lesion was detected, the lowest detection rate may occur. Therefore, the simultaneous detection of multiple lesions helps to increase the probability of targeting drugs. Our key finding is that the different EGFR mutations always occur in different lesions of one patient, suggesting that heterogeneous distribution of EGFR mutation is the major driver event in the development of LUAD. The intertumoral heterogeneity of EGFR-activating mutations has also been confirmed at the single-cell level, which was associated with the EGFR-TKI response in LUAD patients harboring the EGFR L858R mutation (9). The resistance led by the T790M mutation is more attributed to 19 exon del than L858R (18). When the patients have a double mutation with exon 19del and T790M or exon 19del and L858R, patients will not benefit from first-generation EGFR-targeted TKI drugs and may benefit from AZD9291 (19). Identical tumor subtypes respond differently to the same drug, which may be among many others an effect of intratumoral heterogeneity. We also identified HSP90AA1 mutations in one patient with two lesions. It has been reported that HSP90AA1 is associated with a shorter overall survival rate in patients with NSCLC (20). Molecular heterogeneity between multiple lesions in individuals with the same pathologic type represents different biological processes, resulting frequently in different treatment responses for each patient. So, we further speculated whether multi-primary tumors enrich tumor heterogeneity, and thus have a greater impact on the distribution of TMB.
TMB sketches out the status of genomic mutations (21), which is emerging as a practical biomarker for predicting the response of immune checkpoint inhibitors (ICIs) (22). Early studies have shown that PD-L1 expression was markedly different between primary tumors and paired metastatic lymph nodes (23). Different histologic components within a tumor and different pathologic features contribute to the heterogeneous PD-L1 expression in patients with NSCLC (24). By exploring the TMB distribution in different lesions of the intratumor, our results showed that in one patient, TMB distribution in different lesions was diverse. Although the pathologic phenotypes of the multiple lesions were consistent, the TMB distribution was still quite different. Thus, heterogeneity within a tumor and across multiple tumors within a patient was further demonstrated. Currently, most clinical trials do not account for intratumoral heterogeneity of TMB. Intratumor heterogeneity may contribute to the ambiguous clinical results on TMB status as a predictor for immunotherapy response in patients with LUAD.
Morphologic heterogeneity of the tumor may sometimes provide an important clue to genomic heterogeneity of the tumor, which is often associated with responses to particular molecular-targeted and immunotherapy therapy. It also likely reflects the fact that, because of the considerable tumor, patient, and treatment heterogeneity, no one “optimal” management strategy can be delineated. This shows that the simultaneous detection of multiple primary or metastasis lesions that meet the sequencing standards can help provide more clinical guidance information.
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://db.cngb.org/, CNP0001479.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Review Board (IRB) of the First Affiliated Hospital of Guangzhou Medical University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
JH, YQ, and KW conceived the study. LL, HY, HC, QD, DX, YL, and WJ collected study materials or patients. DS, WL, and CZ performed data cleaning and statistical analysis. WL, CZ, and DS wrote the manuscript. All authors read and approved the final manuscript.
Funding
This study is funded by the Foundation and Applied Basic Research Fund of Guangdong Province (2020A1515011293) and the National Natural Science Foundation of China (Grant No. 81772486).
Conflict of interest
CZ, WL, DS, and WJ are employees of BGI Genomics which produces the panel test used in this study.
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.
Publisher’s note
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.
References
1. Best SA, De Souza DP, Kersbergen A, Policheni AN, Dayalan S, Tull D, et al. Synergy between the KEAP1/NRF2 and PI3K pathways drives non-Small-Cell lung cancer with an altered immune microenvironment. Cell Metab (2018) 27:935–943 e934. doi: 10.1016/j.cmet.2018.02.006
2. 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:1693–704. doi: 10.1038/ng.3990
3. Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med (2017) 9:34. doi: 10.1186/s13073-017-0424-2
4. Chen G, Gao C, Gao X, Zhang DH, Kuan SF, Burns TF, et al. Wnt/beta-catenin pathway activation mediates adaptive resistance to BRAF inhibition in colorectal cancer. Mol Cancer Ther (2018) 17:806–13. doi: 10.1158/1535-7163.MCT-17-0561
5. Chu QS. Targeting non-small cell lung cancer: driver mutation beyond epidermal growth factor mutation and anaplastic lymphoma kinase fusion. Ther Adv Med Oncol (2020) 12:1758835919895756. doi: 10.1177/1758835919895756
6. De Sousa VML, Carvalho L. Heterogeneity in lung cancer. Pathobiology (2018) 85:96–107. doi: 10.1159/000487440
7. Gallegos Ruiz MI, Floor K, Roepman P, Rodriguez JA, Meijer GA, Mooi WJ, et al. Integration of gene dosage and gene expression in non-small cell lung cancer, identification of HSP90 as potential target. PloS One (2008) 3:e0001722.
8. George J, Lim JS, Jang SJ, Cun Y, Ozretic L, Kong G, et al. Comprehensive genomic profiles of small cell lung cancer. Nature (2015) 524:47–53. doi: 10.1038/nature14664
9. Guo L, Chen Z, Xu C, Zhang X, Yan H, Su J, et al. Intratumoral heterogeneity of EGFR-activating mutations in advanced NSCLC patients at the single-cell level. BMC Cancer (2019) 19:369. doi: 10.1186/s12885-019-5555-y
10. Joshi A, Pande N, Noronha V, Patil V, Kumar R, Chougule A, et al. ROS1 mutation non-small cell lung cancer-access to optimal treatment and outcomes. Ecancermedicalscience (2019) 13:900. doi: 10.3332/ecancer.2019.900
11. Li S, Choi YL, Gong Z, Liu X, Lira M, Kan Z, et al. Comprehensive characterization of oncogenic drivers in asian lung adenocarcinoma. J Thorac Oncol (2016) 11:2129–40. doi: 10.1016/j.jtho.2016.08.142
12. Li Y, Xu Y, Wu X, He C, Liu Q, Wang F. Comprehensive analysis of EGFR T790M detection by ddPCR and ARMS-PCR and the effect of mutant abundance on the efficacy of osimertinib in NSCLC patients. J Thorac Dis (2019) 11:3004–14. doi: 10.21037/jtd.2019.07.42
13. Liang H, Pan Z, Wang W, Guo C, Chen D, Zhang J, et al. The alteration of T790M between 19 del and L858R in NSCLC in the course of EGFR-TKIs therapy: A literature-based pooled analysis. J Thorac Dis (2018) 10:2311–20. doi: 10.21037/jtd.2018.03.150
14. Liu Y, Dong Z, Jiang T, Hou L, Wu F, Gao G, et al. Heterogeneity of PD-L1 expression among the different histological components and metastatic lymph nodes in patients with resected lung adenosquamous carcinoma. Clin Lung Cancer (2018) 19:e421–30. doi: 10.1016/j.cllc.2018.02.008
15. Park S, Ha S, Lee SH, Paeng JC, Keam B, Kim TM, et al. Intratumoral heterogeneity characterized by pretreatment PET in non-small cell lung cancer patients predicts progression-free survival on EGFR tyrosine kinase inhibitor. PloS One (2018) 13:e0189766. doi: 10.1371/journal.pone.0189766
16. Pi C, Xu CR, Zhang MF, Peng XX, Wei XW, Gao X, et al. EGFR mutations in early-stage and advanced-stage lung adenocarcinoma: Analysis based on large-scale data from china. Thorac Cancer (2018) 9:814–9.
17. Saito Y, Horiuchi S, Morooka H, Ibi T, Takahashi N, Ikeya T, et al. Inter-tumor heterogeneity of PD-L1 expression in non-small cell lung cancer. J Thorac Dis (2019) 11:4982–91.
18. Schrock AB, Frampton GM, Suh J, Chalmers ZR, Rosenzweig M, Erlich RL, et al. Characterization of 298 patients with lung cancer harboring MET exon 14 skipping alterations. J Thorac Oncol (2016) 11:1493–502.
19. Schrock AB, Li SD, Frampton GM, Suh J, Braun E, Mehra R, et al. Pulmonary sarcomatoid carcinomas commonly harbor either potentially targetable genomic alterations or high tumor mutational burden as observed by comprehensive genomic profiling. J Thorac Oncol (2017) 12:932–42.
20. Shaw AT, Yeap BY, Mino-Kenudson M, Digumarthy SR, Costa DB, Heist RS, et al. Clinical features and outcome of patients with non-small-cell lung cancer who harbor EML4-ALK. J Clin Oncol (2009) 27:4247–53.
21. Shi J, Hua X, Zhu B, Ravichandran S, Wang M, Nguyen C, et al. Somatic genomics and clinical features of lung adenocarcinoma: A retrospective study. PloS Med (2016) 13:e1002162.
22. Shukuya T, Carbone DP. Predictive markers for the efficacy of anti-PD-1/PD-L1 antibodies in lung cancer. J Thorac Oncol (2016) 11:976–88.
23. Chalmers ZR, Huang FW, Gay LM, Ali SM, Chmielecki J, Ross JS, et al. Analysis of tumor mutation burden (TMB) in> 51,000 clinical cancer patients to identify novel non-coding PMS2 promoter mutations associated with increased TMB. J Clin Oncol (2016) 34:9572–2.
Keywords: driver mutation, egfr, tmb, heterogeneity, luad
Citation: Qiu Y, Liu L, Yang H, Chen H, Deng Q, Xiao D, Lin Y, Zhu C, Li W, Shao D, Jiang W, Wu K and He J (2022) Intratumor heterogeneity of driver mutations and TMB distribution in 30 early-stage LUAD patients with multiple lesions. Front. Oncol. 12:952572. doi: 10.3389/fonc.2022.952572
Received: 25 May 2022; Accepted: 05 August 2022;
Published: 30 August 2022.
Edited by:
Manash K. Paul, California State University, Los Angeles, United StatesReviewed by:
Jing Bai, Geneplus Beijing Institute, ChinaVartika Sharma, Department of Integrative Biology and Physiology, University of California, Los Angeles, United States
Copyright © 2022 Qiu, Liu, Yang, Chen, Deng, Xiao, Lin, Zhu, Li, Shao, Jiang, Wu and He. 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.
*Correspondence: Jianxing He, hejx@vip.163.com
†These authors have contributed equally to this work