Genomic evolutionary trajectory of metastatic squamous cell carcinoma of the lung.

BACKGROUND: The extent of inter- and intratumoral genomic heterogeneity and the clonal evolution of metastatic squamous cell carcinoma of the lung (LUSC) are poorly understood. Genomic studies of LUSC are challenged by their low tumor cell content. We sought to define the genomic landscape and evolutionary trajectories of metastatic LUSC combining nuclei-flow sorting and whole exome sequencing.
METHODS: Five patients with primary LUSC and six matched metastases were investigated. Tumor nuclei were sorted based on ploidy and expression of cytokeratin to enrich for tumor cells for whole exome sequencing.
RESULTS: Flow-sorting increased the mean tumor purity from 26% (range, 12-50%) to 73% (range, 42-93%). Overall, primary LUSCs and their matched metastases shared a median of 79% (range, 67-85%) of copy number aberrations (CNAs) and 74% (range, 65-94%) of non-synonymous mutations, including in tumor suppressor genes such as TP53. Furthermore, the ploidy of the tumors remained unchanged between primary and metastasis in 4/5 patients over time. We found differences in the mutational signatures of shared mutations compared to the private mutations in the primary or metastasis.
CONCLUSIONS: Our results demonstrate a close genomic relationship between primary LUSCs and their matched metastases, suggesting late dissemination of the metastases from the primary tumors during tumor evolution. 2021 Translational Lung Cancer Research. All rights reserved.

Introduction

Distant metastases are more common in lung adenocarcinoma (LUAD) than in lung squamous cell carcinoma (LUSC). This is one reason why metastatic LUSC has been understudied (1). Most published studies on metastatic LUSC are single case reports, often without genomic analysis of the metastasis due to a lack of sufficient material, or with a limited approach as panel sequencing or FISH (2-5). Advanced LUSC tends to affect older patients with more comorbidities than LUAD, which makes tissue sampling for research studies challenging (6). In daily practice, the biopsies of LUSC metastases are often not large enough for whole exome or whole genome sequencing. Therefore, the mechanisms and evolutionary trajectories of metastasis in LUSC remain poorly understood.

Multi-region sampling and comprehensive molecular profiling of resected tumors have greatly advanced our understanding of intratumoral heterogeneity (ITH) in lung and other cancer types (7). However, the high costs and the diagnostic reality of mostly small biopsy specimens from patients with advanced LUSC make such multi-region analyses difficult if not impossible in clinical practice. Molecular profiling of LUSC is further challenged by the often low tumor cell content of LUSC tissue (8). To overcome this hurdle, we previously showed that the tumor purity of LUAD can be greatly enhanced by an advanced nuclei flow-sorting technique (9). Briefly, cells were stained with 4’,6-Diamin-2-phenylindol (DAPI) and an antibody for TTF-1 to isolate LUAD specific cells. This flow-sorting approach allowed us to distinguish and isolate aneuploid and diploid populations enabling highly informative genomic profiling (9).

Here, we applied a modified nuclei flow-sorting approach to enrich for cytokeratin-labelled tumor nuclei from tumor tissue specimens. An antibody was used for labelling pan-cytokeratins that are highly expressed in LUSC but not in admixed benign mesenchymal cells. Whole exome sequencing (WES) was performed to investigate ITH and clonal evolution from five patients with LUSC and matched spatiotemporal metastases. We present the following article in accordance with the MDAR reporting checklist (available at http://dx.doi.org/10.21037/tlcr-21-48).

Methods

Patient cohort

The retrospective study comprised the following criteria for the inclusion of patients: (I) histologically confirmed LUSC with regional or distant metastases; (II) large enough specimens (i.e., surgical specimens or large biopsies) to allow for nuclei flow-sorting and WES; (III) at least two available tumor samples, fresh-frozen (FF) or formalin-fixed and paraffin-embedded (FFPE), one of which had to be the primary tumor (primary LUSC) and the other the matched metastasis; (IV) tumor samples that had been taken from different time points or body sites of the patient. Multiple regions within the primary tumor were not sampled. Five patients with primary LUSC and six matched metastases were investigated. Tumor samples, FF (n=10) or FFPE (n=6, ), were obtained from the archive of the Institute of Medical Genetics and Pathology of the University Hospital Basel, Switzerland. Histological slides were reviewed by a pulmonary pathologist (LB). The areas with highest tumor cell density were marked for subsequent dissection, avoiding areas rich in stromal cells or benign epithelial cells from entrapped bronchoalveolar tissue. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by the local Ethics Committee of Both Basels (EKBB) and Ethics Committee Northwest and Central Switzerland (EKNZ) (NO.: EKBB/EKNZ 31/12) and individual consent for this retrospective analysis was waived.

Table 1

Sequencing and sorting metrics

Patient	Tumor	Tissue	Sorted with	Ploidy	Mutations/Mb	Mean coverage	Total reads
P6	Primary	FF	DAPI	3.6	2.4	57.99	116,651,162
P6	Metastasis	FF	DAPI	3.4	6	69.02	160,142,677
P6	Non-tumor	FF	DAPI	2	–	62.17	137,946,635
P9	Primary	FF	DAPI	3.9	4.6	68.41	148,484,030
P9	Metastasis	FF	DAPI	4	6.2	73.91	150,098,988
P9	Non-tumor	FF	DAPI	2	–	62.55	132,798,619
P18	Primary	FF	DAPI	3.6	6.6	51.15	114,207,787
P18	Metastasis	FF	DAPI	3.3	6.6	73.17	127,941,373
P18	Non-tumor	FF	DAPI	2	–	77.15	152,991,220
P103	Primary	FFPE	DAPI/pCK	3.3	6.9	63.32	199,580,439
P103	Metastasis	FFPE	DAPI/pCK	3.2	6.1	75.91	177,816,752
P103	Non-tumor	FFPE	DAPI/pCK	2	–	32.34	127,962,411
P109	Primary	FF	DAPI	3.6	6.2	91.92	142,302,414
P109	Metastasis1	FFPE	DAPI/pCK	2	9.4	51.75	210,698,902
P109	Metastasis2	FFPE	DAPI/pCK	1.8	4.8	98.72	157,381,241
P109	Non-tumor	FFPE	DAPI/pCK	2	–	58.87	231,972,945

FF, fresh frozen; FFPE, formalin-fixed paraffin-embedded; DAPI, 4',6-diamino-2-phenylindole; pCK, pan-cytokeratin.

Nuclei isolation and multiparameter flow-sorting

Nuclei isolation from FF and FFPE samples was conducted as described previously (9-12). All nuclei were DAPI stained for sorting, DNA quantification and ploidy analysis. FF tumor nuclei were only sorted with DAPI (n=7, four patients). Pan-Cytokeratin (pCK, clone MNF116, Dako, Code-Nr. M0821) was used for FFPE tumors as an additional marker for flow-sorting. The applied pCK antibody is specific against CK5, CK6, CK8, CK17 and CK19 and highly expressed in LUSC (13,14). Cytokeratin positive cytoplasmic remnants seem to stay attached to the nuclei of FFPE tissue but are washed away in FF tissue during the procedure. Therefore, only FFPE samples were stained with pCK. This step was crucial to enable refined separation of diploid from aneuploid tumor populations and benign diploid epithelial and stromal cell nuclei in FFPE tissue (n=2, P109, Figure S1). Previously, we described the utility of sorted diploid population as germline (9). In DAPI-sorted aneuploid tumors, the matched diploid populations were used as germline for WES (P6, P9, P18, Figure S2). In samples sorted by DAPI and pCK, diploid pCK-negative populations were used as germline (P103, P109) (9). Initial purity was estimated after sorting by treating all sorted nuclei as 100% and calculating all sorted populations according to their proportions.

DNA extraction and whole exome sequencing

DNA from FFPE isolated nuclei was directly extracted using DNeasy blood and tissue kit (Qiagen, Germantown, MD, USA) according to the manufacturer’s instructions and quantified using the Qubit Fluorometer assay (Life Technologies, Carlsbad, CA, USA), whereas DNA from FF samples was first whole genome amplified using the illustra GenomiPhi V2 DNA Amplification Kit (GE Healthcare Life Sciences, Little Chalfont, Buckinghamshire, UK; product: 25660031) as previously described (9,10). In samples that were exclusively sorted with DAPI, the diploid population was used as germline, whereas diploid pCK-negative populations were considered as germline when applicable (9). Flat genomes were considered as a germline (Figure S2). According to manufacturer’s guidelines SureSelect Human All Exon V6 Kit (Agilent) was utilized for the whole exome capturing. Paired-end 100-bp reads were generated on the Illumina NovaSeq 6000. Sequencing was conducted by CeGaT (Tübingen, Germany).

Human genome GRCh37 was used as a reference for aligning the reads. Whole exome sequencing (WES) analysis, comprising sequence processing, calling of somatic single nucleotide variants (SNVs), small insertions and deletions (indels), and allele-specific copy number variations (CNAs), detecting mutational signatures and identifying clonality, was performed using state-of-the-art bioinformatic methods, as previously described with modifications (15,16).

Briefly, SNVs and indels were called using MuTect v1.1.7 and strelka v2.9.10, respectively (17,18). To normalize the FF and FFPE samples we used following approach: All SNVs C:G>T:A with variant allelic fractions (VAFs) less than 10% were discarded to reduce false positive results. Otherwise, SNVs or indels with VAFs <1% or that were supported by fewer than 3 reads were not considered. If SNVs and indels were found in one tumor of a patient, a cut-off of two reads was applied to the remaining tumors. We further excluded variants identified in at least two of a panel of 123 non-tumor samples, including the non-tumor samples included in our study. Variant annotation was performed using SnpEff software v4.1 (19).

The heatmap of non-synonymous mutations was generated using the R package maftools v.2.0.16 by selecting the genes in the Bailey et al. dataset (20) and the significant mutated genes in TCGA dataset (MutSig Q-value <0.05 based on 511 LUSC samples; TCGA, Firehose Legacy, http://www.cbioportal.org) that represents the significantly mutated genes in LUSC.

Analysis of copy number aberrations and clonality

FACETS v0.5.14 enabled the identification of allele-specific CNAs (21). Cancer cell fraction (CCF) of each mutation was inferred using ABSOLUTE v1.0.6 (22). A mutation was classified as clonal, if its probability of being clonal was >50% or if the lower bound of the 95% confidence interval of its CCF was >90%. Mutations were regarded as subclonal if they did not meet the mentioned criteria (23,24).

Phylogenetic analysis and mutational signatures

As described by Murugaesu et al., a maximum parsimony tree was built for each case using binary presence/absence matrices based on the repertoire of non-synonymous and synonymous somatic mutations in the biopsies of the tumors (25-27).

A mutation that was found in both tumor manifestations of one patient was considered as ‘trunk’. Mutations that were detected in only one tumor of the patient were considered as ‘branch’ or ‘private’.

Decomposition of mutational signatures was conducted using the R package deconstructSigs software by selecting the mutational signatures 1, 2, 4, 5 and 13, which are based on the set of 30 mutational signatures (COSMIC) that were observed in LUSC (28). To have a significant result, the amount of mutations must be >50 for the algorithm. For P6 (primary LUSC and metastasis), P109 primary LUSC and P18 primary LUSC it was <50. However, these results can serve as hints to further understand the mechanism of tumor evolution. Phylogenetic trees were drawn using FigTree v1.4.4.

Statistical analysis

All statistical analyses were performed using R software v3.6.3 (R Core Team, 2020). A Wilcoxon test was performed using the R function “wilcox.test” to check whether the primary LUSCs have a significantly different number of non-synonymous mutations than the metastases.

Results

Nuclei flow-sorting achieves high purity of tumor DNA

In this study, we investigated five patients with primary LUSC and their metastases from different sites to explore ITH and clonal evolution (). All patients were male, and three were former or current smokers, whereas the smoking status was unknown in two patients. Our tissue material cohort consisted of ten FF (Fresh Frozen) and six FFPE (Formalin-fixed paraffin-embedded) samples (). Although flow-sorting by DAPI allows separating aneuploid tumor nuclei from benign stromal tissue and epithelial cells, it does not detect diploid tumor cells (9,10). To enable genomic profiling of diploid tumor populations in LUSC, we implemented, in addition to DAPI, pCK staining to separately select diploid (2N±0.2) and/or aneuploid (>2.2N) tumor cell populations and remove pCK-negative diploid and near tetraploid benign cells. Across the five patients, we isolated nine aneuploid tumor cell populations using either DAPI alone or DAPI/pCK (). Diploid tumor cell populations were identified only by DAPI/pCK sorting (). Our flow-sorting approach enabled us to increase the initial mean purity from 26% (range, 12–50%) of unsorted tumor to a mean of 73% (range, 42–93%) after sorting (Figure S3).

Table 2

Clinical characteristics of the cohort

Patient	P6	P9	P18	P103	P109
Sex	Male	Male	Male	Male	Male
Age at diagnosis (years)	54	61	64	73	68
Smoker at diagnosis	Current	NA	Former	NA	Current
Smoking status (pack-years)	50	NA	NA	NA	NA
pTN*	pT4, pN1	pT2a, pN0	pT3, pN0	pT4, pN1	pT2a, pN1
Diameter of primary tumor	5.5 cm	3.5 cm	3.2 cm	7.8 cm	3.5 cm
Radiotherapy	NA	NA	No	No	Yes#
Chemotherapy	No	No	No	No	Yes#
Time to metastasis (months)	5	8	50	1	0 (Met1) 28 (Met2)
Site of metastasis/recurrence	Left kidney	Right kidney	Main bronchial carina	Brain	Interlobular lymph node (pN1) (Met1) thoracic spine (Met2)
Ploidy (N) primary tumor	3.6	3.9	3.6	3.3	3.6
Ploidy (N) metastasis	3.4	4	3.3	3.2	2 (Met1) 1.8 (Met2)

#Palliative radiotherapy and chemotherapy (Cisplatin/Vinorelbine and Cisplatin/Gemcitabine) prior to resection of Met 2; *according to UICC TNM classification (8th Edition). Met1, metastasis 1; Met2, metastasis 2.

Mutational landscape of primary-metastatic pairs of LUSC

We performed WES of the sorted tumor cell nuclei and normal samples to a median depth of 69x (range, 58×–99×) and 62× (range, 32×–77×) respectively (). Median somatic mutation rate was 6.2 mutations/Mb., primary LUSCs and metastases mutation rates were similar (primary LUSC median 6.2; range 2.4–6.9 vs. metastases median 6.1, range 4.8–9.4, ). The five primary LUSCs do not have fewer non-synonymous mutations than the metastases (average 325; range 144–414 vs. average 393.3; range 290–568, Wilcoxon test P value =0.79). A median of 74% (range, 65–94%) of non-synonymous mutations in the primary LUSC were also found in the corresponding metastases ( and https://cdn.amegroups.cn/static/public/tlcr-21-48-1.xlsx). On the other hand, a median of 69.7% (range, 35–93%) of the non-synonymous mutations found in the metastases were found in the matched primaries. The most commonly mutated cancer-related genes known to be mutated LUSC included TP53 (9 samples, 4 patients), FAT1 [5, 2], FBXW7 [5, 2], KMT2D [5, 2] and ARID1A [2, 2] () (29-32).

Figure 1

Clonal relationship of primary LUSC and matched metastases. (A) Venn diagram displays the number of non-synonymous mutations per patient that are common between the primary tumor (Pri) (blue) and the metastases (Met) (red, green). (B) Heatmap depicts a comparison between the non-synonymous genes in the dataset of Bailey et al. and the significant mutated genes in TCGA dataset (MutSig Q-value <0.05 based on 511 LUSC samples) that represent the most significantly mutated genes in LUSC compared to the primary tumor (pri) and the matched metastases (met) of the presented patients. Heatmap illustrates the cancer cell fraction of selected mutations. Clonal mutations are illustrated with a diagonal line.

PT and metastatic samples of LUSC shared common mutations known to drive tumorigenesis and tumor progression. Clonality analysis of primary LUSCs and matched metastases revealed that mutations in driver genes, such as the tumor suppressor genes TP53, FBXW7 and KMT2D, were predominantly clonal (). Clonal mutations between the primary LUSCs and metastases were found in each patient. Primary LUSC of P103 and metastases of P9 and P109 displayed private mutations in cancer-related gene.

In all five patients, the trunks were much longer than branches (). P109 harbored a long mutational trunk and similar number of mutations in Met1. Private mutations of P109 primary LUSC and Met2 had a much shorter branch compared to Met1. Notably, P109 indicated a common clonal trunk with its two metastases, a lymph node (Met1) and a thoracic spine metastasis (Met2).

Figure 2

Evolutionary change of mutational signatures in primary tumor and metastases. Evolution of the somatic genetic alterations illustrates the changes in mutational processes. The pie chart delineates the proportion of COSMIC mutational signatures. Black, red, and blue lines represent the trunk, the primary LUSC (Pri) and the metastasis (Met) branch, respectively. Branch lengths are proportional to the number of mutations. *This patient received a palliative radiotherapy and chemotherapy (Cisplatin/Vinorelbine and Cisplatin/Gemcitabine) prior to resection of Met 2.

Analysis of mutational signatures revealed that COSMIC signature 1 (CpG deamination, related to aging) and COSMIC signature 4 (related to smoking) were most common being present in all samples across all of our patients (33). COSMIC signature 5 (transcriptional strand bias, related to aging) was present in every primary-metastasis pair but not in all trunks and branches. The APOPEC signature was variably distributed. It was identified in trunk and both branches separately of primary LUSC and metastasis in P9. In contrast it was truncal in P109 but restricted to the branches in P18, and only found in the primary LUSC branch of P103. In P6, it was present in the trunk and the primary LUSC branch but not in the branch of the metastasis.

Taken together, the mutational profiles between the primary LUSCs and their metastasis are closely related as previously seen in LUAD (9). The long trunks of mutations can be primarily attributed to the effects of aging and smoking. The pattern of mutational signatures varied between patients and during the evolution of individual tumors.

Chromosomal aberrations are truncal events

CNAs of the primary LUSCs were highly concordant with those in matched metastases (, Figure S4). In primary LUSCs, a median of 79% (range, 67–85%) of CNAs were also observed in the matched metastases, while a median of 79% (range, 70–99%) of the aberrations in the metastases were also seen in the matched primaries. All tumor samples, primary LUSCs and metastases, shared amplified regions in 3q26.1-33 including the genes PIK3CA, SOX2 and TERC). In addition, there was truncal TERT amplification at 5p15.33 in three of the five patients. We found the following tumor suppressor genes to be lost in our cohort: ATM (2 samples, 2 patients), ARID1A [3, 2], APC [7, 4], CDKN2A [9, 4], CSMD1 [9, 4], FAT1 [3, 2], FOXP1 [8, 4], PTEN [5, 2], SETD2 [8, 4] and amplifications were detected in SOX2 [10, 5], MECOM [10, 5], TERC [10, 5], MYC [4, 2] and CCND1 [4, 2].

Figure 3

Similarity of primary tumors and metastases on chromosomal level. (A) Repertoire of copy number alterations as defined by WES of primary LUSC (pri), first metastasis (met) and second metastasis (met2) of each patient. Samples represented on the y-axis; chromosomes are represented along the x-axis. Light red: copy number loss; red: homozygous deletion; light blue: copy number gain; dark blue: amplification; (B) ploidy analysis measured using flow cytometry. Samples are represented on the x-axis; ploidy is represented in the y-axis. blue: primary tumor (pri); red: first metastasis (met); green: second metastasis (met2).

Both chromosomal profiles and ploidy of primary LUSCs and matched metastases were highly similar in most patients (, Figure S4). All primary LUSCs were aneuploid, including primary LUSC of P109 (3.6 N). Interestingly, both metastases of P109 were diploid (2N±0.2N).

Taken together, CNAs and ploidy status were very similar in the primary LUSCs and their metastases suggesting CNAs as early events with limited heterogeneity between primary LUSC and metastasis.

Discussion

In this study, we explored the spatiotemporal heterogeneity and evolution of LUSC by WES using highly enriched nuclei from five primary LUSCs and matched metastases. Intratumoral heterogeneity (ITH) and clonal evolution of primary NSCLC using comprehensive genomic profiling strategies have been analyzed in several previous studies (7,34-36). Although metastasis is the leading cause of death in NSCLC, data on matched pairs of primary LUSCs and metastases remain scarce (9,37,38).

Two other studies have recently been published on spatial and temporal heterogeneity of metastatic LUSC from a small number of patients (39,40). Leong et al. found that the interactions of the DNA repair genes with the tumor microenvironment may play a role in the acquisition of somatic mutations following the metastasis during lung cancer progression (39). Another recent study on metastatic LUSC, demonstrated that large structural variants may play a crucial role in intratumoral genetic heterogeneity (40).

Here, we demonstrate high concordance between the primary-metastasis pairs in terms of SNVs and CNAs in most of the patients. Most aberrations in LUSC were truncal indicating only limited heterogeneity between primary LUSC and metastasis and low subclonal diversity (9). This is in line with the data of a recent pan-cancer study on metastatic solid tumors that included treatment-naive and extensively treated tumors (41). This study revealed that 96% of driver mutations within metastases were clonal with only 4% of subclonal mutations (41). The potential clinical impact of subclonal diversity remains poorly studied. We could adequately address the question whether treatment might affect subclonal diversity in our study since four of the five patients had not received systemic treatment prior to tissue sampling (34).

The high concordance between primary-metastasis pairs in our study corroborates our previous data on a cohort of LUAD (9). In our previous study, more than 80% of CNAs and cancer-related gene mutations were shared between the primary-metastatic LUAD specimens. Ploidy remained stable and did not change over time (9). Here, we found that CNAs and ploidy were equally stable during the metastatic spread of LUSC. This stable aneuploidy in most LUSC and LUAD patients highlights the previously proposed critical importance of a whole genome doubling (WGD) event that might select for fitness to progress and metastasize (42,43). The only discrepant ploidy was found in P109 revealing an aneuploid primary LUSC (3.6N) and two diploid/near-diploid metastases (2N and 1.8N, respectively). Reversion from an aneuploid to a diploid state during progression is unlikely, since WGD is considered irreversible (44). Therefore, we hypothesize that the diploid metastases had emerged from a diploid population prior to WGD of the primary LUSC in this patient. We do not know whether the diploid population had been fully replaced by the aneuploid population after WGD, or whether we have missed a coexistent diploid population due to limited sampling of the primary LUSC. Although most LUSC and LUAD are aneuploid and metastasize after WGD, this example illustrates that stably diploid tumor cell populations can occasionally serve as the backbone of progression and metastasis in NSCLC (9). We have also previously observed this phenomenon in prostate cancer, where a diploid tumor cell population was consistently detectable in repeated biopsies during progression giving rise to temporary aneuploid populations (10).

Somatic mutations in cancers are caused by different mutational processes that generate a characteristic mutational signature (45). The pathogenesis of the vast majority of NSCLC is driven by the accumulation of smoking-induced alterations (7,46). This is in line with the large proportion of signature 4 in our tumor samples reflecting the contribution of smoking-related DNA damage for early acquisition and accumulation of truncal mutations as previously proposed. The smoking status of P9 and P103 was not known. However, the signature profiles showed a significant amount of smoking signature 4, suggesting that these patients were former or current smokers. Mutational signatures were previously inferred to track the evolutionary differences between the primary LUSCs and the metastases (32). In breast cancer, Ng et al. found significant differences in mutational signatures comparing private mutations of the primary tumor and metastasis (27). The discrepancy might contribute to spatial intratumoral heterogeneity. Our study could support these previous results. Mutational signatures changed during metastasis and over time revealing no obvious pattern. Despite the common trunk, the mutational signatures differed in the private branches of the primary LUSCs and metastases. Furthermore, the long trunk of shared genomic alterations in most samples indicates a linear progression model in the metastatic dissemination caused by one major subclone (41,43). In fact, a common ancestor acquired genomic alterations and remained genetically stable even after late dissemination. Nevertheless, there were substantial differences between the private branches of primary LUSCs and metastases. This mutational heterogeneity might be relevant in the context of specific treatments. Notably, Richard et al. have recently reported that different mutational signatures had an impact on the efficacy of immune checkpoint inhibitors (47).

Our study provides new insight into the genomic trajectories of genomic evolution of metastatic LUSC but did not reveal novel specific alterations beyond those that have already been reported by large-scale sequencing consortia (48). Among CNAs, it was most striking that all tumors shared amplifications at 3q26.1-29 [30–40% in cBioPortal (48)] indicating an important role of these CNAs for tumor initiation (49,50). SOX2 and TP63 are among the genes in the amplified 3p26.1-29 region and associated with squamous differentiation (51-53). Notably, SOX2 and p63 are known to be amplified in at least 40% and 30%, respectively, of all LUSC but almost never in LUAD (46).

Our study has several limitations. The low number of patients precludes strong general statements on the genomic trajectories of LUSC. Further confounding factors are the different states of tissue fixation (FF vs. FFPE) and sorting strategies (DAPI vs. DAPI/pCK). The variability in sample quality and the small cohort size is a major limitation but reflects the general difficulties to access appropriate samples from matched primary LUSCs and metastases. Nevertheless, our collection of LUSC is unique being the largest series undergoing WES with matched primary LUSC and metastasis pairs of LUSC so far. Due to incomplete clinical documentation of archived cases, we were unable to reconstruct the full clinicopathological history including smoking history (54). Using WES, we were bound to coding regions preventing any observations on translocations. Furthermore, limited by the sequencing depth and filtering, it is possible that subclonal variants were missed due to detection limits. In fact, one of the limitations of this study, due the possible presence of sequencing and FFPE artefacts, is to remove the mutations with an allele frequency lower then 1% and/or validated by at least 3 reads. This means, especially for the samples with a mean coverage of ~50–70x, failing to detect any subclone with an allele frequency lower than 4–6% even in the ideal circumstances.

Conclusions

Taken together, our study provides new insights into the intratumoral heterogeneity and genomic patterns of primary LUSC and matched metastases using purified tumor populations after nuclei flow-sorting. A close clonal relationship with only limited heterogeneity was revealed in the primary-metastatic pairs, similar to what was seen in LUAD, suggesting late dissemination of the metastases from the primary tumors during tumor evolution. Further studies with larger matched cohorts are required to better elucidate the critical genomic players and related molecular pathways of metastatic disease in NSCLC.

The article’s supplementary files as