Front. Oncol. Frontiers in Oncology Front. Oncol. 2234-943X Frontiers Media S.A. 10.3389/fonc.2020.00726 Oncology Original Research The Genomic Characteristics of ALK Fusion Positive Tumors in Chinese NSCLC Patients Liu Shaokun 1 2 3 † Huang Tanxiao 4 † Liu Ming 4 He Wenlong 1 2 3 Zhao YingShen 4 Yang Lizhen 1 2 3 Long Yingjiao 1 2 3 Zong Dandan 1 2 3 Zeng Huihui 1 2 3 Liu Yuanyuan 4 Liao Wenting 4 Duan Jingxian 4 Gong Subo 5 * Chen Shifu 4 * 1Department of Pulmonary and Critical Care Medicine, The Second Xiangya Hospital, Central South University, Changsha, China 2Research Unit of Respiratory Disease, Central South University, Changsha, China 3Diagnosis and Treatment Center of Respiratory Disease, Central South University, Changsha, China 4HaploX Biotechnology Co., Ltd., Shenzhen, China 5Department of Geriatrics, The Second Xiangya Hospital, Central South University, Changsha, China

Edited by: Ira Ida Skvortsova, Innsbruck Medical University, Austria

Reviewed by: Xiaofei Song, University of Texas MD Anderson Cancer Center, United States; Qi Zhao, Sun Yat-sen University Cancer Center (SYSUCC), China

 *Correspondence: Subo Gong gsb510@csu.edu.cn Shifu Chen chen@haplox.com

This article was submitted to Cancer Genetics, a section of the journal Frontiers in Oncology

†These authors have contributed equally to this work and share first authorship

 08 05 2020 2020 10 726 22 01 2020 16 04 2020 Copyright © 2020 Liu, Huang, Liu, He, Zhao, Yang, Long, Zong, Zeng, Liu, Liao, Duan, Gong and Chen. 2020 Liu, Huang, Liu, He, Zhao, Yang, Long, Zong, Zeng, Liu, Liao, Duan, Gong and Chen

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.

Anaplastic lymphoma kinase (ALK) fusion events account for ~3–7% genetic alterations in patients with non-small cell lung cancer (NSCLC). In this study, we identified the ALK fusion patterns and a novel ALK fusion partner in 44 ALK positive NSCLC patients using a customized HapOncoCDx panel, and identified ALK fusion partners. The most common partner is EML4, forming the variant 1 (v1, E13:A20, 18/44), variant 2 (v2, E20:A20, 5/44), and variant 3 (v3, E6:A20, 13/44). Moreover, we detected a new ALK fusion partner HMBOX1. At the mutation level, TP53 is the most frequently mutated gene (24%), followed by ALK (12%) and STED2 (12%). The median tumor mutation burden (TMB) of these samples is 2.29 mutations/Mb, ranging from 0.76 mut/Mb to 16.79 muts/Mb. We further elaborately portrayed the TP53 mutation sites on the peptide sequence of the encoded protein by lollipop. The mutational signature and copy number alterations (CNAs) of the samples were also analyzed. The CNA events were found in 13 (13/44) patients, and the most commonly amplified genes were MDM2 (n = 4/13) and TERT (n = 4/13). Together, these results may guide personalized clinical management of patients with ALK fusion in the era of precision medicine.

 NSCLC ALK NGS—next generation sequencing copy number aberrations genomic landscape Introduction

Anaplastic lymphoma kinase (ALK) fusion events, which are the result of ALK rearrangements, account for ~3–7% genetic alterations in non-small cell lung cancer (NSCLC) patients (1, 2). These oncogenic mutations could lead to the constitutive activation of the ALK tyrosine kinase domain, and further cause tumorigenesis (3). Hitherto, multiple ALK fusion partners have been identified, and the most normal one is echinoderm microtubule-associated protein-like 4 (EML4), which were observed in nearly 80% of all the ALK fusion cases (2). It is worth noting that more than a dozen of different EML4-ALK variants have been identified in NSCLC patients. The most common variants are variant 1 (v1, E13:A20), variant 2 (v2, E20:A20), and variant 3 (v3, E6:A20) (4).

Currently, ALK tyrosine kinase inhibitors (TKIs) are recommended for the treatment of NSCLC patients harboring ALK fusion (5). Prior to ALK TKIs treatment, the median overall survival (OS) of ALK fusion-positive metastatic NSCLC patients receiving chemotherapy was around 12 months (6). However, under the sequential treatment with ALK fusion TKIs, the OS of the patients were extended to 5 years (7). The first approved is Crizotinib which is the first generation TKI. Compared to traditional chemotherapy, it improves the PFS and OS of ALK fusion NSCLC patients significantly (8). Nevertheless, nearly all of the patients would develop drug resistance within 2 years after the initial treatment. The drug resistance was possibly caused by a secondary mutation in the kinase domain of ALK, the activation of alternative pathways, ALK amplification, or epithelial-mesenchymal transition (9). To overcome the resistance, the second-generation ALK TKIs were developed including Ceritinib, Alectinib, and Brigatinib. They were approved for the treatment of metastatic NSCLC patients with ALK fusion and had progressed on or intolerant to Crizotinib. Notably, Ceritinib and Alectinib were approved for the first-line treatment of the ALK fusion positive metastatic NSCLC patients (10–16). In addition, as the third generation ALK inhibitor, Lorlatinib has also been approved for the treatment of metastatic NSCLC patients with ALK fusion, on condition that the disease has progressed on Crizotinib or at least one other ALK inhibitor such as Alectinib or Ceritinib for metastatic disease (17). It is worth noting that different ALK inhibitors have different potencies and spectrums against acquired resistance mutations (18).

In the era of precision medicine, the genomic profiles of the patients may guide the planning of treatment strategy. For the ALK fusion positive NSCLC patients, detailed genomic profiles can elucidate the fusion partner and the rearranged breakpoint. Moreover, the proposed resistant mutations are critical for clinical treatment guidelines. Furthermore, several studies reported that molecular profiling is also associated with the prognosis of patients. Christopoulos et al. reported that the concomitant TP53 mutation is a strong indicator of poor prognosis in ALK fusion positive NSCLC patients (19). Their study also reported that EML4-ALK fusion variant V3 was associated with a more aggressive phenotype and worse overall survival due to earlier failure of several therapeutic modalities. In addition, they found that V3 and TP53 double positive patients had a very high risk of death with a median OS of around 2 years.

With the development of next-generation sequencing technologies, it is becoming more affordable to obtain the genomic landscape of cancer patients. In this study, we aim to demonstrate the genomic landscape of ALK fusion-positive tumors in 44 Chinese NSCLC patients sequenced with our customized HapOncoCDx panel which involves hybridization capture and deep sequencing of all protein-coding exons of 464 cancer-associated genes and other selected introns of other oncogenes and tumor suppressor genes, and illustrate their genomic mutation patterns and characteristics, which potentially helps to develop treatment strategy.

 Materials and Methods Patients and Samples

Forty-four patients were enrolled from 1349 NSCLC patients in this study. Tumor tissues were collected during surgery, and were formalin fixed, paraffin-embedded (FFPE) and archived. Peripheral blood (PBL) samples were also collected from each patient as control.

 DNA Extraction

DNA samples were extracted from Formalin-fixed paraffin-embedded (FFPE) samples with QIAamp DNA FFPE tissue kit (Qiagen). Extraction of PBL DNA was conducted using the RelaxGene Blood DNA system (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer's protocol. All the DNA samples were quantified both using the Qubit 2.0 fluorometer and the Qubit dsDNA HS Assay kit (Thermo Fisher Scientific, Inc., Waltham, MA, USA) according to the manufacturer's protocol.

 Library Construction and Sequencing

One hundred nanogram of DNA from each sample was sheared by the dsDNA Fragmentase (New England BioLabs, Inc., Ipswich, MA, USA), and then performed size selection (150–250 bp) using Ampure XP beads (Beckman Coulter, Inc., Brea, CA, USA). Library construction was performed using the KAPA Library Preparation kit (Kapa Biosystems, Inc., Wilmington, MA, USA) according to the manufacturer's protocol. The concentration of the library were assessed using the e Qubit dsDNA HS Assay kit, and fragment length was determined on a 4200 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). Target enrichment was carried out using the Agilent SureSelect XT HS kit (Agilent Technologies) according to the manufacturer's Protocol. DNA sequencing was then performed on the Illumina Novaseq 6000 system at an average depth of 2000X.

 Data Analysis and Variant Calling

Raw sequences were pre-processed by fastp version 0.18.0 (https://github.com/OpenGene/fastp) (20), and clean reads were aligned to the hg19 genome (GRch37) using Burrows-Wheeler Aligner maximal exact matches algorithm (21). The Gencore version 0.12.0 (https://github.com/OpenGene/gencore) (22) was used for removing duplicate reads. Pileup files with mapping quality ≥60 were generated using Samtools version 0.1.19 (http://www.htslib.org/) (23). Somatic variants were called using VarScan2 version 2.3.8 (http://varscan.sourceforge.net/) (24) [the minimum read depth 20; the variant allele frequency (VAF) threshold ≥0.01; somatic-P ≤ 0.01; strand-filter ≥1; others, default parameters]. CNV kit with version 0.9.3 (https://github.com/etal/cnvkit) (25) was used for copy number variation detection, and GeneFuse version v0.6.1 (https://github.com/OpenGene/GeneFuse) (26) for structural variation detection. Maftools was used for visualizing somatic variant analysis (27).

 Results Sample Collection and Patient Characteristics

Of the 1349 NSCLC cases, ALK rearrangements were detected in 44 cases (3.26 %). Those 44 Chinese patients with locally advanced or metastatic NSCLC were enrolled in this study, of which 20 (45.5%) were female. All patients carry ALK fusion events. Their mean age was 52.5 with ranging from 29 to 73. NGS was performed on 44 pairs of tumor and white blood cell samples. All the samples that passed the histology quality control (HQC) yielded sufficient amounts of DNA for NGS.

 Identification of ALK Rearrangements Using Targeted Sequencing

In order to identify ALK rearrangement from the DNA of patients' FFPE samples, we designed probes to cover the intron 18 and intron 19 of ALK, as well as introns of some well-known ALK fusion partners. We identified ALK rearrangements and corresponding breakpoints in the sequencing data of these patients. The statistical summary and breakpoints of the rearrangement events are listed in Table 1 and shown in Figures 1, 2, respectively. We found that 43 out of 44 patients had an EML4-ALK fusion, with variant 1 (v1, E13:A20), variant 2 (v2, E20:A20), and variant 3 (v3, E6:A20) detected in 18, 5, and 13 patients, respectively. We also identified one novel ALK fusion partner HMBOX1.

Fusion patterns of ALK.

Fusion type Counts Percent (%)
EML4-exon13-ALK-exon20 18 40.9
EML4-exon6-ALK-exon20 13 29.5
EML4-exon20-ALK-exon20 5 11.4
EML4-exon13-ALK-exon19 2 4.5
EML4-exon14-ALK-exon20 2 4.5
EML4-exon6-ALK-exon19 2 4.5
EML4-exon19-ALK-exon20 1 2.3
HMBOX1-exon4-ALK-exon20 1 2.3

The statistics of different ALK rearrangement forms. The number of each ALK fusion pattern identified in 44 NSCLC patients are shown in the barchart.

The breakpoint distribution in ALK and the respective fusion partners. Each arrowed line represents one fusion event. The exact breakpoints of ALK in GRch37 are shown in the middle panel, while the fused exons of ALK fusion partners are shown in the top and bottom panels, respectively. The sequence of ALK is exhibited reversely (from right to left), while the sequences of EML4 and HMBOX1 are represented in the forward direction. The genomic region of ALK between 299446000 and 29448500 on Chromosome 2 is divided into regions every 500 bp. Breakpoint position in ALK locate between 29446000 and 29446500 with an orange arrow, between 29446501 and 29447000 with a blue arrow, between 29447001 and 29447500 with a green arrow, between 29447501 and 29448000 with an yellow arrow, between 29448001 and 29448500 with a gray arrow.

 Mutational Profiles of ALK Fusion Positive NSCLC Patients

Genomic alterations were found in 34 (n = 34/44, 77.3%) samples with a total of 134 alterations identified including variants of non-synonymous mutations and splicing mutations. The detailed information is shown in Figure 3A. The mutation landscapes of ALK fusion positive NSCLC patients were highly heterogeneous. The median TMB was 2.29 mut/Mb with a range between 0.76 and 16.79 mut/Mb which is similar to the TMB of the TCGA NSCLC cohort.

Mutational profiles of ALK positive NSCLC patients. (A) This is the oncoprint of the somatic SNVs and Indels in 34 patients in our study. Somatic alterations included missense, nonsense, frameshift indel, in-frame indel, splice site, translation start site, multi_Hit mutations. The genes are ranked by the frequency of the mutations across all samples. (B) This is the oncoprint of the somatic SNVs and Indels in 34 patients from the MSK-IMPACT study. Somatic alterations included missense, nonsense, frameshift indel, in-frame indel, and splice-site mutations. The genes are ranked by the frequency of the mutations across all samples.

We constructed a heatmap to demonstrate the somatic mutations occurred in the tumor tissues of the patients (Figure 3A). TP53 was most commonly altered (n = 8/34, 24%), followed by SETD2 (n = 4, 12%), ALK (n = 4, 12%), SYNE1 (n = 3, 9%), SMAD4 (n = 3, 9%), SLX4 (n = 3, 9%), NOTCH3 (n = 3, 9%), LRP1B (n = 3, 9%), EP300 (n = 3, 9%), and CTNNB1 (n = 3, 9%).

Other genomic alterations of low frequencies are AMER1 (n = 2, 6%), ARID1B (n = 2, 6%), CSF3R (n = 2, 6%), FAT1 (n = 2, 6%), FOXP2 (n = 2, 6%), KDM6A (n = 2, 6%), KMT2A (n = 2, 6%),LATS1 (n = 2, 6%), MAP2K4 (n = 2, 6%), NFEL2L2 (n = 2, 6%), NOTCH1(n = 2, 6%), NTRK3 (n = 2, 6%), TERT (n = 2, 6%), and TGFBR2 (n = 2, 6%). Alterations in ABL1, ADH1B, ALDH2, APC, AR, ARID2, ATM, AURKA, BMPR1A, CACNA1C, CADM2, CAMTA1, CAPN2, CARD11, CDC73, CDK12, CREBBP, CSMD3, DNMT3A, EPHA3, ERBB4, ESR2, EWSR1, EXT1, EZH2, FGFR1, FLCN, FOXA1, FOXL2, GATA6, GPRIN2, HIF1A, HNF1B, JAK1, KDR, KMT2C, KMT2D, LATS2, MAP2K1, MDM4, MYCN, NF2, NSD1, NTRK1, PTEN, PTPRD, PZP, RARA, RIT1, RNF43, ROS1, SETBP1, SMARCA4, SMO, SOCS6, SOX2, SPEN, STAT3, STK11, SUZ12, TSHR, and U2AF1 were identified in one sample each (n = 1, 3%). We further compared our results with the MSK-IMPACT study (28), in which we extracted 45 ALK fusion positive cases that yielded 81 mutations. Overall, our results were highly consistent with the MSK-IMPACT findings, which showed that TP53 and ALK are the most frequently altered genes (Figure 3B).

We further studied their mutational signatures. We observed that C>T transition occurred most frequently, followed by C>A transition (Figure 4). This pattern is consistent with COSMIC signature 1 that had been found in most cancer samples.

Mutational signatures of ALK fusion positive NSCLC patients. SNPs are classified into transitions and transversions. Summarized data are visualized as a boxplot showing overall distribution of six different conversions (Top) and as a stacked barplot showing the fraction of conversions in each sample (Bottom).

Different driver gene mutations revealed inter-tumor heterogeneity. TP53 mutations in exon 5–8 were observed, and we further elaborately portrayed the TP53 mutation sites on the peptide sequence in a lollipop plot (Figure 5).

Protein variants resulted from TP53 mutations. Protein variants caused by TP53 mutations are displayed in the lollipop plot. These sites are considered to be mutational hot-spots.

 Copy Number Aberrations of ALK Fusion Positive NSCLC Patients

Somatic copy number alterations were found in 13 (n = 13/44, 29.5%) samples. A total of 22 alterations were identified, including gain and loss (Figure 6). MDM2 and TERT were most commonly amplified genes (n = 4/13, 31%), followed by CCND1, EPCAM, IKZF1, MET, MYCN, RICTOR (n = 1, 8%). Loss of copy number was most frequently observed in CD274, CDKN2A, JAK2 (n = 2/13, 15%), followed by FGFR1, FGFR3 (n = 1, 8%).

Copy number aberrations in 13 ALK fusion positive NSCLC patients. The names of the aberrant genes are shown in the y-axis, while each column represents a patient. The type of copy number aberrations, including gain, normal, and loss are indicated by red, yellow, and blue, respectively.

 Discussion

In this study, we identified ALK rearrangement events in 44 Chinese NSCLC patients using NGS technologies. Consistent with other studies, the most common ALK fusion partner is EML4, and the fusion occurs in the forms of the three most common variants. We also report a novel ALK fusion partner HMBOX1. It implied that NGS-based assessment for ALK fusions was accurate and comprehensive, having the unique advantages in detecting unknown ALK fusion partners, and identifying the exact breakpoints compared to the traditional methods, like FISH and IHC.

At the same time, we characterized the mutational profiles of the patients. The results were consistent with other studies, in terms of the relatively lower frequency of TP53 mutations, lower TMB, and fewer co-mutations compared to ALK-negative NSCLC patients (29). Besides, we identified the copy number alterations in their genome. Apart from the genes with a high frequency of copy number amplification, such as MET, MDM2, and TERT, we also identified some genes with copy number loss, such as CD274, CDKN2A, and JAK2. This information is important for guiding optimal clinical treatment. For instance, the copy number loss of CD274 probably indicates a low expression of PD-L1. MDM2 amplification had been reported to associate with a poor clinical outcome and significantly increased tumor growth rate with anti-PD-1/PD-L1 immunotherapies (30).

In conclusion, using our customized HapOncoCDx panel, we not only successfully detected the ALK fusion events in 44 Chinese NSCLC patients, but also explored their genomic mutational landscapes. To the best of our knowledge, this is the first study that exhibited the mutational landscape of Chinese NSCLC patients with ALK rearrangement. This result can provide genomic information for personalized clinical management for patients with ALK fusion in the era of precision medicine.

 Data Availability Statement

The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive in BIG Data Center, Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, under accession numbers HRA000138, HRA000138 that can be accessed at http://bigd.big.ac.cn/gsa-human.

 Ethics Statement

This study was approved by the ethics committee of the Second XIANGYA Hospital of Central South University and complied with Good Clinical Practices, the principles of the Declaration of Helsinki and all applicable regulatory requirements. All patients provided written informed consent prior to any study-specific procedures.

 Author Contributions

SL, TH, SG, and SC designed the study. WH, YZ, LY, YLo, DZ, HZ, and YLi carried out the sequencing experiment and collected data. TH, YZ, and WL performed the bioinformatics analysis. TH, SL, and ML wrote this manuscript. JD and SC revised this manuscript. SG and SC supervised the study.

 Conflict of Interest

TH, ML, YZ, YLi, WL, JD, and SC were employed by the company HaploX Biotechnology 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.

 References Cancer Genome Atlas Research Network Comprehensive molecular profiling of lung adenocarcinoma. Nature. (2014) 511:543–50. 10.1038/nature13385 Soda M Choi YL Enomoto M Takada S Yamashita Y Ishikawa S . Identification of the transforming EML4-ALK fusion gene in non-small-cell lung cancer. Nature. (2007) 448:561–6. 10.1038/nature0594517625570 Rikova K Guo A Zeng Q Possemato A Yu J Haack H . Global survey of phosphotyrosine signaling identifies oncogenic kinases in lung cancer. Cell. (2007) 131:1190–203. 10.1016/j.cell.2007.11.02518083107 Noh KW Lee MS Lee SE Song JY Shin HT Kim YJ . Molecular breakdown: a comprehensive view of anaplastic lymphoma kinase (ALK)-rearranged non-small cell lung cancer. J. Pathol. (2017) 243:307–19. 10.1002/path.495028741662 Kwak EL Bang YJ Camidge DR Shaw AT Solomon B Maki RG . Anaplastic lymphoma kinase inhibition in non-small-cell lung cancer. N Engl J Med. (2010) 363:1693–703. 10.1056/NEJMoa100644820979469 Pilkington G Boland A Brown T Oyee J Bagust A Dickson R. A systematic review of the clinical effectiveness of first-line chemotherapy for adult patients with locally advanced or metastatic non-small cell lung cancer. Thorax. (2015) 70:359–67. 10.1136/thoraxjnl-2014-205914 Duruisseaux M Besse B Cadranel J Pérol M Mennecier B Bigay-Game L . Overall survival with crizotinib and next-generation ALK inhibitors in ALK-positive non-small-cell lung cancer (IFCT-1302 CLINALK): a French nationwide cohort retrospective study. Oncotarget. (2017) 8:21903–17. 10.18632/oncotarget.1574628423535 Shaw AT Kim DW Nakagawa K Seto T Crinó L Ahn MJ . Crizotinib versus chemotherapy in advanced ALK-positive lung cancer, N. Engl. J. Med. (2013) 368:2385–94. 10.1056/NEJMoa1214886 Katayama R Shaw AT Khan TM Mino-Kenudson M Solomon BJ Halmos B . Mechanisms of acquired crizotinib resistance in ALK-rearranged lung Cancers. Sci Transl Med. (2012) 4:120ra17. 10.1126/scitranslmed.300331622277784 Shaw AT Kim DW Mehra R Tan DSW Felip E Laura QM . Ceritinib in ALK-rearranged non-small-cell lung cancer. N Engl J Med. (2014) 370:1189–97. 10.1056/NEJMoa1311107 Shaw AT Kim TM Crinò L Gridelli C Kiura K Liu G . Ceritinib versus chemotherapy in patients with ALK-rearranged non-small-cell lung cancer previously given chemotherapy and crizotinib (ASCEND-5): A randomised, controlled, open-label, phase 3 trial. Lancet Oncol. (2017) 18:874–86. 10.1016/S1470-2045(17)30339-X28602779 Ou SH Ahn JS De Petris L Govindan R Yang JC Hughes B . Alectinib in crizotinib-refractory ALK-rearranged non-small-cell lung cancer: A phase II global study. J Clin Oncol. (2016) 34:661–8. 10.1200/JCO.2015.63.944326598747 Shaw AT Gandhi L Gadgeel S Riely GJ Cetnar J West H . Alectinib in ALK-positive, crizotinib-resistant, non-small-cell lung cancer: a single-group, multicentre, phase 2 trial. Lancet Oncol. (2016) 17:234–42. 10.1016/S1470-2045(15)00488-X26708155 Peters S Camidge DR Shaw AT Gadgeel S Ahn JS Kim DW . Alectinib versus crizotinib in untreated ALK-positive non-small-cell lung cancer. N Engl J Med. (2017) 377:829–38. 10.1056/NEJMoa170479528586279 Gettinger SN Bazhenova LA Langer CJ Salgia R Gold KA Rosell R . Activity and safety of brigatinib in ALK-rearranged non-small-cell lung cancer and other malignancies: A single-arm, open-label, phase 1/2 trial. Lancet Oncol. (2016) 17:1683–96. 10.1016/S1470-2045(16)30392-827836716 Kim DW Tiseo M Ahn MJ Reckamp KL Hansen KH Kim SW . Brigatinib in patients with crizotinib-refractory anaplastic lymphoma kinase-positive non-small-cell lung cancer: A randomized, multicenter phase II trial. J Clin Oncol. (2017) 35:2490–8. 10.1200/JCO.2016.71.590428475456 Shaw AT Felip E Bauer TM Besse B Navarro A Postel-Vinay S . Lorlatinib in non-small-cell lung cancer with ALK or ROS1 rearrangement: An international, multicentre, open-label, single-arm first-in-man phase 1 trial. Lancet Oncol. (2017) 18:1590–9. 10.1016/S1470-2045(17)30680-029074098 Gainor JF Dardaei L Yoda S Friboulet L Leshchiner I Katayama R . Molecular mechanisms of resistance to first- and second- generation ALK inhibitors in ALK-rearranged lung cancer. Cancer Discov. (2016) 6:1118–33. 10.1158/2159-8290.CD-16-059627432227 Christopoulos P Kirchner M Bozorgmehr F Endris V Elsayed M Budczies J . Identification of a highly lethal V3+TP53+ subset in ALK+ lung adenocarcinoma. Int J Cancer. (2019) 144:190–9. 10.1002/ijc.3189330255938 Chen S Zhou Y Chen Y Gu J. FASTP: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. 10.1093/bioinformatics/bty56030423086 Li H Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. (2010) 26:589–95. 10.1093/bioinformatics/btp69820080505 Chen S Zhou Y Chen Y Huang T Liao W Xu Y . Gencore an efficient tool to generate consensus reads for error suppressing and duplicate removing of NGS data. BMC Bioinform. (2019) 20:606. 10.1186/s12859-019-3280-931881822 Li H Handsaker B Wysoker A Fennell T Ruan J Homer N . 1000 Genome Project Data Processing Subgroup: the sequence alignment/map (SAM) format and SAMtools. Bioinformatics. (2009) 25:2078–9. 10.1093/bioinformatics/btp352 Koboldt DC Zhang Q Larson DE Shen D McLellan MD Lin L . VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. (2012) 22:568–76. 10.1101/gr.129684.11122300766 Talevich E Shain AH Botton T Bastian BC. CNVkit: genome-wide copy number detection and visualization from targeted sequencing. PLoS Comput Biol. (2014) 12:e1004873. 10.1371/journal.pcbi.100487327100738 Chen S Liu M Huang T Liao W Xu M Gu J. GeneFuse: detection and visualization of target gene fusions from DNA sequencing data. Int J Biol Sci. (2018) 14:843–8. 10.7150/ijbs.2462629989075 Mayakonda A Lin DC Assenov Y Plass C Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. (2018) 28:1747–56. 10.1101/gr.239244.11830341162 Zehir A Benayed R Shah RH Syed A Middha S Kim HR . Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat Med. (2017) 23:703–13. 10.1038/nm.4333 Christopoulos P Budczies J Kirchner M Dietz S Sültmann H Thomas M . Defining molecular risk in ALK+ NSCLC. Oncotarget. (2019) 10:3093–103. 10.18632/oncotarget.2688631139322 Kato S Goodman A Walavalkar V Barkauskas DA Sharabi A Kurzrock R. Hyperprogressors after immunotherapy: analysis of genomic alterations associated with accelerated growth rate. Clin Cancer Res. (2017) 23:4242–50. 10.1158/1078-0432.CCR-16-313328351930

Funding. This work was supported by the project of Bureau of Industry and Information Technology of Shenzhen (Grant No. 20170922151538732), Shenzhen Science and Technology Innovation Committee Technical Research Project (Grant No. JSGG20180703164202084), the project of Development and Reform Commission of Shenzhen Municipality (Grant No. XMHT20190104006), and the National Key Clinical Specialty Construction Projects of China (Grant No. 2012-650).