Representation of genomic intratumor heterogeneity in multi-region non-small cell lung cancer patient-derived xenograft models.

Patient-derived xenograft (PDX) models are widely used in cancer research. To investigate the genomic fidelity of non-small cell lung cancer PDX models, we established 48 PDX models from 22 patients enrolled in the TRACERx study. Multi-region tumor sampling increased successful PDX engraftment and most models were histologically similar to their parent tumor. Whole-exome sequencing enabled comparison of tumors and PDX models and we provide an adapted mouse reference genome for improved removal of NOD scid gamma (NSG) mouse-derived reads from sequencing data. PDX model establishment caused a genomic bottleneck, with models often representing a single tumor subclone. While distinct tumor subclones were represented in independent models from the same tumor, individual PDX models did not fully recapitulate intratumor heterogeneity. On-going genomic evolution in mice contributed modestly to the genomic distance between tumors and PDX models. Our study highlights the importance of considering primary tumor heterogeneity when using PDX models and emphasizes the benefit of comprehensive tumor sampling.

Research autopsy programmes in oncology: shared experience from 14 centres across the world.

While there is a great clinical need to understand the biology of metastatic cancer in order to treat it more effectively, research is hampered by limited sample availability. Research autopsy programmes can crucially advance the field through synchronous, extensive, and high-volume sample collection. However, it remains an underused strategy in translational research. Via an extensive questionnaire, we collected information on the study design, enrolment strategy, study conduct, sample and data management, and challenges and opportunities of research autopsy programmes in oncology worldwide. Fourteen programmes participated in this study. Eight programmes operated 24 h/7 days, resulting in a lower median postmortem interval (time between death and start of the autopsy, 4 h) compared with those operating during working hours (9 h). Most programmes (n = 10) succeeded in collecting all samples within a median of 12 h after death. A large number of tumour sites were sampled during each autopsy (median 15.5 per patient). The median number of samples collected per patient was 58, including different processing methods for tumour samples but also non-tumour tissues and liquid biopsies. Unique biological insights derived from these samples included metastatic progression, treatment resistance, disease heterogeneity, tumour dormancy, interactions with the tumour micro-environment, and tumour representation in liquid biopsies. Tumour patient-derived xenograft (PDX) or organoid (PDO) models were additionally established, allowing for drug discovery and treatment sensitivity assays. Apart from the opportunities and achievements, we also present the challenges related with postmortem sample collections and strategies to overcome them, based on the shared experience of these 14 programmes. Through this work, we hope to increase the transparency of postmortem tissue donation, to encourage and aid the creation of new programmes, and to foster collaborations on these unique sample collections. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Copy number architectures define treatment-mediated selection of lethal prostate cancer clones.

Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient evolutionary relationships of metastases that evade treatment, we perform genome-wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 167 metastatic regions from 11 organs harvested post-mortem from 10 men who died from prostate cancer. We identify diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number changes, we confirm a common clone of origin across metastases and diagnostic biopsies, and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. These autosome-defined clusters are characterized by cluster-specific AR gene architectures, and in two index cases are topologically more congruent than by chance (p-values 3.07 × 10-8 and 6.4 × 10-4). Integration with anatomical sites suggests patterns of spread and points of genomic divergence. Here, we show that copy number boundaries identify treatment-selected clones with putatively distinct lethal trajectories.

Evolutionary characterization of lung adenocarcinoma morphology in TRACERx.

Lung adenocarcinomas (LUADs) display a broad histological spectrum from low-grade lepidic tumors through to mid-grade acinar and papillary and high-grade solid, cribriform and micropapillary tumors. How morphology reflects tumor evolution and disease progression is poorly understood. Whole-exome sequencing data generated from 805 primary tumor regions and 121 paired metastatic samples across 248 LUADs from the TRACERx 421 cohort, together with RNA-sequencing data from 463 primary tumor regions, were integrated with detailed whole-tumor and regional histopathological analysis. Tumors with predominantly high-grade patterns showed increased chromosomal complexity, with higher burden of loss of heterozygosity and subclonal somatic copy number alterations. Individual regions in predominantly high-grade pattern tumors exhibited higher proliferation and lower clonal diversity, potentially reflecting large recent subclonal expansions. Co-occurrence of truncal loss of chromosomes 3p and 3q was enriched in predominantly low-/mid-grade tumors, while purely undifferentiated solid-pattern tumors had a higher frequency of truncal arm or focal 3q gains and SMARCA4 gene alterations compared with mixed-pattern tumors with a solid component, suggesting distinct evolutionary trajectories. Clonal evolution analysis revealed that tumors tend to evolve toward higher-grade patterns. The presence of micropapillary pattern and 'tumor spread through air spaces' were associated with intrathoracic recurrence, in contrast to the presence of solid/cribriform patterns, necrosis and preoperative circulating tumor DNA detection, which were associated with extra-thoracic recurrence. These data provide insights into the relationship between LUAD morphology, the underlying evolutionary genomic landscape, and clinical and anatomical relapse risk.

Body composition and lung cancer-associated cachexia in TRACERx.

Cancer-associated cachexia (CAC) is a major contributor to morbidity and mortality in individuals with non-small cell lung cancer. Key features of CAC include alterations in body composition and body weight. Here, we explore the association between body composition and body weight with survival and delineate potential biological processes and mediators that contribute to the development of CAC. Computed tomography-based body composition analysis of 651 individuals in the TRACERx (TRAcking non-small cell lung Cancer Evolution through therapy (Rx)) study suggested that individuals in the bottom 20th percentile of the distribution of skeletal muscle or adipose tissue area at the time of lung cancer diagnosis, had significantly shorter lung cancer-specific survival and overall survival. This finding was validated in 420 individuals in the independent Boston Lung Cancer Study. Individuals classified as having developed CAC according to one or more features at relapse encompassing loss of adipose or muscle tissue, or body mass index-adjusted weight loss were found to have distinct tumor genomic and transcriptomic profiles compared with individuals who did not develop such features. Primary non-small cell lung cancers from individuals who developed CAC were characterized by enrichment of inflammatory signaling and epithelial-mesenchymal transitional pathways, and differentially expressed genes upregulated in these tumors included cancer-testis antigen MAGEA6 and matrix metalloproteinases, such as ADAMTS3. In an exploratory proteomic analysis of circulating putative mediators of cachexia performed in a subset of 110 individuals from TRACERx, a significant association between circulating GDF15 and loss of body weight, skeletal muscle and adipose tissue was identified at relapse, supporting the potential therapeutic relevance of targeting GDF15 in the management of CAC.

Genomic-transcriptomic evolution in lung cancer and metastasis.

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

The evolution of non-small cell lung cancer metastases in TRACERx.

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.

The evolution of lung cancer and impact of subclonal selection in TRACERx.

Lung cancer is the leading cause of cancer-associated mortality worldwide1. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource.

Lung adenocarcinoma promotion by air pollutants.

A complete understanding of how exposure to environmental substances promotes cancer formation is lacking. More than 70 years ago, tumorigenesis was proposed to occur in a two-step process: an initiating step that induces mutations in healthy cells, followed by a promoter step that triggers cancer development1. Here we propose that environmental particulate matter measuring ≤2.5 µm (PM2.5), known to be associated with lung cancer risk, promotes lung cancer by acting on cells that harbour pre-existing oncogenic mutations in healthy lung tissue. Focusing on EGFR-driven lung cancer, which is more common in never-smokers or light smokers, we found a significant association between PM2.5 levels and the incidence of lung cancer for 32,957 EGFR-driven lung cancer cases in four within-country cohorts. Functional mouse models revealed that air pollutants cause an influx of macrophages into the lung and release of interleukin-1ß. This process results in a progenitor-like cell state within EGFR mutant lung alveolar type II epithelial cells that fuels tumorigenesis. Ultradeep mutational profiling of histologically normal lung tissue from 295 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 53% of healthy tissue samples, respectively. These findings collectively support a tumour-promoting role for  PM2.5 air pollutants  and provide impetus for public health policy initiatives to address air pollution to reduce disease burden.

Folate receptor alpha in ovarian cancer tissue and patient serum is associated with disease burden and treatment outcomes.

BACKGROUND: Survival rates for ovarian cancer remain poor, and monitoring and prediction of therapeutic response may benefit from additional markers. Ovarian cancers frequently overexpress Folate Receptor alpha (FRα) and the soluble receptor (sFRα) is measurable in blood. Here we investigated sFRα as a potential biomarker. METHODS: We evaluated sFRα longitudinally, before and during neo-adjuvant, adjuvant and palliative therapies, and tumour FRα expression status by immunohistrochemistry. The impact of free FRα on the efficacy of anti-FRα treatments was evaluated by an antibody-dependent cellular cytotoxicity assay. RESULTS: Membrane and/or cytoplasmic FRα staining were observed in 52.7% tumours from 316 ovarian cancer patients with diverse histotypes. Circulating sFRα levels were significantly higher in patients, compared to healthy volunteers, specifically in patients sampled prior to neoadjuvant and palliative treatments. sFRα was associated with FRα cell membrane expression in the tumour. sFRα levels decreased alongside concurrent tumour burden in patients receiving standard therapies. High concentrations of sFRα partly reduced anti-FRα antibody tumour cell killing, an effect overcome by increased antibody doses. CONCLUSIONS: sFRα may present a non-invasive marker for tumour FRα expression, with the potential for monitoring patient response to treatment. Larger, prospective studies should evaluate FRα for assessing disease burden and response to systemic treatments.

Real-time ex vivo perfusion of human lymph nodes invaded by cancer (REPLICANT): a feasibility study.

Understanding how breast cancer (BC) grows in axillary lymph nodes (ALNs), and refining how therapies might halt that process, is clinically important. However, modelling the complex ALN microenvironment is difficult, and no human models exist at present. We harvested ALNs from ten BC patients, and perfused them at 37 °C ex vivo for up to 24 h. Controlled autologous testing showed that ALNs remain viable after 24 h of ex vivo perfusion: haematoxylin and eosin-stained histological appearance and proliferation (by Ki67 immunohistochemistry) did not change significantly over time for any perfused ALN compared with a control from time-point zero. Furthermore, targeted gene expression analysis (NanoString PanCancer IO360 panel) showed that only 21/750 genes were differentially expressed between control and perfused ALNs (|log2 FC| > 1 and q < 0.1): none were involved in apoptosis and metabolism, but rather all 21 genes were involved in immune function and angiogenesis. During perfusion, tissue acid-base balance remained stable. Interestingly, the flow rate increased (p < 0.001) in cancer-replaced (i.e. metastasis occupied more than 90% of the surface area on multiple levels) compared to cancer-free nodes (i.e. nodes with no metastasis on multiple sections). CXCL11 transcripts were significantly more abundant in cancer-replaced nodes, while CXCL12 transcripts were significantly more abundant in cancer-free nodes. These cytokines were also detected in the circulating perfusate. Monoclonal antibodies (nivolumab and trastuzumab) were administered into a further three ALNs to confirm perfusion efficacy. These drugs saturated the nodes; nivolumab even induced cancer cell death. Normothermic ALN perfusion is not only feasible but sustains the tumour microenvironment ex vivo for scientific investigation. This model could facilitate the identification of actionable immuno-oncology targets. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

An innate-like Vδ1+ γδ T cell compartment in the human breast is associated with remission in triple-negative breast cancer.

Innate-like tissue-resident γδ T cell compartments capable of protecting against carcinogenesis are well established in mice. Conversely, the degree to which they exist in humans, their potential properties, and their contributions to host benefit are mostly unresolved. Here, we demonstrate that healthy human breast harbors a distinct γδ T cell compartment, primarily expressing T cell receptor (TCR) Vδ1 chains, by comparison to Vδ2 chains that predominate in peripheral blood. Breast-resident Vδ1+ cells were functionally skewed toward cytolysis and IFN-γ production, but not IL-17, which has been linked with inflammatory pathologies. Breast-resident Vδ1+ cells could be activated innately via the NKG2D receptor, whereas neighboring CD8+ αß T cells required TCR signaling. A comparable population of Vδ1+ cells was found in human breast tumors, and when paired tumor and nonmalignant samples from 11 patients with triple-negative breast cancer were analyzed, progression-free and overall survival correlated with Vδ1+ cell representation, but not with either total γδ T cells or Vδ2+ T cells. As expected, progression-free survival also correlated with αß TCRs. However, whereas in most cases TCRαß repertoires focused, typical of antigen-specific responses, this was not observed for Vδ1+ cells, consistent with their innate-like responsiveness. Thus, maximal patient benefit may accrue from the collaboration of innate-like responses mounted by tissue-resident Vδ1+ compartments and adaptive responses mounted by αß T cells.

Gene expression modules in primary breast cancers as risk factors for organotropic patterns of first metastatic spread: a case control study.

BACKGROUND: Metastases from primary breast cancers can involve single or multiple organs at metastatic disease diagnosis. Molecular risk factors for particular patterns of metastastic spread in a clinical population are limited. METHODS: A case-control design including 1357 primary breast cancers was used to study three distinct clinical patterns of metastasis, which occur within the first six months of metastatic disease: bone and visceral metasynchronous spread, bone-only, and visceral-only metastasis. Whole-genome expression profiles were obtained using whole genome (WG)-DASL assays from formalin-fixed paraffin-embedded (FFPE) samples. A systematic protocol was developed for handling FFPE samples together with stringent data quality controls to identify robust expression profiling data. A panel of published and novel gene sets were tested for association with these specific patterns of metastatic spread and odds ratios (ORs) were calculated. RESULTS: Metasynchronous metastasis to bone and viscera was found in all intrinsic breast cancer subtypes, while immunohistochemically (IHC)-defined receptor status and specific IntClust subgroups were risk factors for visceral-only or bone-only first metastases. Among gene modules, those related to proliferation increased the risk of metasynchronous metastasis (OR (95% CI) = 2.3 (1.1-4.8)) and visceral-only first metastasis (OR (95% CI) = 2.5 (1.2-5.1)) but not bone-only metastasis (OR (95% CI) = 0.97 (0.56-1.7)). A 21-gene module (BV) was identified in estrogen-receptor-positive breast cancers with metasynchronous metastasis to bone and viscera (area under the curve = 0.77), and its expression increased the risk of bone and visceral metasynchronous spread in this population. BV was further orthogonally validated with NanoString nCounter in primary breast cancers, and was reproducible in their matched lymph nodes metastases and an external cohort. CONCLUSION: This case-control study of WG-DASL global expression profiles from FFPE tumour samples, after careful quality control and RNA selection, revealed that gene modules in the primary tumour have differing risks for clinical patterns of metasynchronous first metastases. Moreover, a novel gene module was identified as a putative risk factor for metasynchronous bone and visceral first metastatic spread, with potential implications for disease monitoring and treatment planning.

High-throughput RNA interference screening using pooled shRNA libraries and next generation sequencing.

RNA interference (RNAi) screening is a state-of-the-art technology that enables the dissection of biological processes and disease-related phenotypes. The commercial availability of genome-wide, short hairpin RNA (shRNA) libraries has fueled interest in this area but the generation and analysis of these complex data remain a challenge. Here, we describe complete experimental protocols and novel open source computational methodologies, shALIGN and shRNAseq, that allow RNAi screens to be rapidly deconvoluted using next generation sequencing. Our computational pipeline offers efficient screen analysis and the flexibility and scalability to quickly incorporate future developments in shRNA library technology.

High-throughput RNAi screening reveals novel regulators of telomerase.

Telomerase is considered an attractive anticancer target on the basis of its common and specific activation in most human cancers. While direct telomerase inhibition is being explored as a therapeutic strategy, alternative strategies to target regulators of telomerase that could disrupt telomere maintenance and cancer cell proliferation are not yet available. Here, we report the findings of a high-throughput functional RNA interference screen to globally profile the contribution of kinases to telomerase activity (TA). This analysis identified a number of novel telomerase modulators, including ERK8 kinase, whose inhibition reduces TA and elicited characteristics of telomere dysfunction. Given that kinases represent attractive drug targets, we addressed the therapeutic implications of our findings, such as demonstrating how limiting TA via kinase blockade could sensitize cells to inhibition of the telomere-associated protein tankyrase. Taken together, our findings suggest novel combinatorial approaches to targeting telomere maintenance as a strategy for cancer therapy.

Interaction between CHEK2*1100delC and other low-penetrance breast-cancer susceptibility genes: a familial study.

BACKGROUND: The allele CHEK2*1100delC doubles the risk of breast cancer in unselected women, but could confer a greater risk in women with a family history of the disease, particularly of bilateral breast cancer. Our aim was to measure the risk of breast cancer in relatives of women with bilateral breast cancer who were carriers of this allele. METHODS: A population-based series of 469 bilateral breast cancer cases ascertained through English cancer registries were genotyped for CHEK2*1100delC. Standardised incidence ratios (SIRs) and cumulative risks were calculated for breast cancer, prostate cancer, and all other cancers in the first-degree relatives of carriers and non-carriers. FINDINGS: The relatives of bilateral cases who were wild-type for CHEK2 had three times the population risk of female breast cancer (145 cases: SIR 3.48 (95% CI 2.96-4.09), twice the risk of prostate cancer (34 cases: SIR 2.41, 1.67-3.36) and a large excess of male breast cancer (five cases: SIR 15.06, 4.92-35.36). Relatives of those who were carriers of CHEK2*1100delC had a substantially higher risk of breast cancer (eight cases: SIR 12.11, 5.23-23.88) and possibly prostate cancer (two cases: SIR 9.87, 1.20-35.67). INTERPRETATION: These data suggest a multiplicative interaction between CHEK2*1100delC and other unknown susceptibility genes. In women with a family history of bilateral disease, CHEK2*1100delC confers a high lifetime risk and might be useful for predictive testing. Bilateral breast cancer cases and their families are likely to provide an efficient basis for identification of additional low-penetrance breast-cancer genes.