Integrative genomic analyses of European intrahepatic cholangiocarcinoma: Novel ROS1 fusion gene and PBX1 as prognostic marker.

BACKGROUND: Cholangiocarcinoma (CCA) is a fatal cancer of the bile duct with a poor prognosis owing to limited therapeutic options. The incidence of intrahepatic CCA (iCCA) is increasing worldwide, and its molecular basis is emerging. Environmental factors may contribute to regional differences in the mutation spectrum of European patients with iCCA, which are underrepresented in systematic genomic and transcriptomic studies of the disease. METHODS: We describe an integrated whole-exome sequencing and transcriptomic study of 37 iCCAs patients in Germany. RESULTS: We observed as most frequently mutated genes ARID1A (14%), IDH1, BAP1, TP53, KRAS, and ATM in 8% of patients. We identified FGFR2::BICC1 fusions in two tumours, and FGFR2::KCTD1 and TMEM106B::ROS1 as novel fusions with potential therapeutic implications in iCCA and confirmed oncogenic properties of TMEM106B::ROS1 in vitro. Using a data integration framework, we identified PBX1 as a novel central regulatory gene in iCCA. We performed extended screening by targeted sequencing of an additional 40 CCAs. In the joint analysis, IDH1 (13%), BAP1 (10%), TP53 (9%), KRAS (7%), ARID1A (7%), NF1 (5%), and ATM (5%) were the most frequently mutated genes, and we found PBX1 to show copy gain in 20% of the tumours. According to other studies, amplifications of PBX1 tend to occur in European iCCAs in contrast to liver fluke-associated Asian iCCAs. CONCLUSIONS: By analyzing an additional European cohort of iCCA patients, we found that PBX1 protein expression was a marker of poor prognosis. Overall, our findings provide insight into key molecular alterations in iCCA, reveal new targetable fusion genes, and suggest that PBX1 is a novel modulator of this disease.

Germline homozygosity and allelic imbalance of HLA-I are common in esophagogastric adenocarcinoma and impair the repertoire of immunogenic peptides.

BACKGROUND: The individual HLA-I genotype is associated with cancer, autoimmune diseases and infections. This study elucidates the role of germline homozygosity or allelic imbalance of HLA-I loci in esophago-gastric adenocarcinoma (EGA) and determines the resulting repertoires of potentially immunogenic peptides. METHODS: HLA genotypes and sequences of either (1) 10 relevant tumor-associated antigens (TAAs) or (2) patient-specific mutation-associated neoantigens (MANAs) were used to predict good-affinity binders using an in silico approach for MHC-binding (www.iedb.org). Imbalanced or lost expression of HLA-I-A/B/C alleles was analyzed by transcriptome sequencing. FluoroSpot assays and TCR sequencing were used to determine peptide-specific T-cell responses. RESULTS: We show that germline homozygosity of HLA-I genes is significantly enriched in EGA patients (n=80) compared with an HLA-matched reference cohort (n=7605). Whereas the overall mutational burden is similar, the repertoire of potentially immunogenic peptides derived from TAAs and MANAs was lower in homozygous patients. Promiscuity of peptides binding to different HLA-I molecules was low for most TAAs and MANAs and in silico modeling of the homozygous to a heterozygous HLA genotype revealed normalized peptide repertoires. Transcriptome sequencing showed imbalanced expression of HLA-I alleles in 75% of heterozygous patients. Out of these, 33% showed complete loss of heterozygosity, whereas 66% had altered expression of only one or two HLA-I molecules. In a FluoroSpot assay, we determined that peptide-specific T-cell responses against NY-ESO-1 are derived from multiple peptides, which often exclusively bind only one HLA-I allele. CONCLUSION: The high frequency of germline homozygosity in EGA patients suggests reduced cancer immunosurveillance leading to an increased cancer risk. Therapeutic targeting of allelic imbalance of HLA-I molecules should be considered in EGA.

Development, testing and validation of a targeted NGS-panel for the detection of actionable mutations in lung cancer (NSCLC) using anchored multiplex PCR technology in a multicentric setting.

Lung cancer is a paradigm for a genetically driven tumor. A variety of drugs were developed targeting specific biomarkers requiring testing for tumor genetic alterations in relevant biomarkers. Different next-generation sequencing technologies are available for library generation: 1) anchored multiplex-, 2) amplicon based- and 3) hybrid capture-based-PCR. Anchored multiplex PCR-based sequencing was investigated for routine molecular testing within the national Network Genomic Medicine Lung Cancer (nNGM). Four centers applied the anchored multiplex ArcherDX-Variantplex nNGMv2 panel to re-analyze samples pre-tested during routine diagnostics. Data analyses were performed by each center and compiled centrally according to study design. Pre-defined standards were utilized, and panel sensitivity was determined by dilution experiments. nNGMv2 panel sequencing was successful in 98.9% of the samples (N = 90). With default filter settings, all but two potential MET exon 14 skipping variants were identified at similar allele frequencies. Both MET variants were found with an adapted calling filter. Three additional variants (KEAP1, STK11, TP53) were called that were not identified in pre-testing analyses. Only total DNA amount but not a qPCR-based DNA quality score correlated with average coverage. Analysis was successful with a DNA input as low as 6.25 ng. Anchored multiplex PCR-based sequencing (nNGMv2) and a sophisticated user-friendly Archer-Analysis pipeline is a robust and specific technology to detect tumor genetic mutations for precision medicine of lung cancer patients.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma.

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Frequency and functional characterization of fusion genes in squamous cell carcinoma of the lung.

INTRODUCTION: In contrast to lung adenocarcinoma (LUAD), targetable genetic alterations are less frequently detected in squamous cell carcinoma of the lung (LUSC). Over the last years, gene fusions have become promising targets in many solid cancers. Here, we analysed a cohort of LUSC, identified recurrent fusion genes and functionally characterised these tumour genomes. METHODS: A subset of 1608 squamous cell carcinomas of the lung was analysed by means of the FusionPlex® Lung Panel to identify potentially targetable gene fusions using targeted next-generation sequencing. Cases harbouring recurrent gene fusions were further analysed using FISH, Cytoscan HD arrays and cell culture experiments. RESULTS: We found both, known and novel gene fusions in about 3 % of the cases. Known fusions occurring in lung cancer included ALK::EML4, EGFRvIII, EZR::ROS1 and FGFR3::TACC. We further identified recurrent gene fusions of currently unknown biological function, involving EGFR::VSTM2A and NSD3::FGFR1 and showed that the occurrence of the EGFR::VSTM2A fusion is accompanied by high-level amplification of EGFR. Our analyses further revealed that the genomes of these LUSC patients are chromosomally unstable, which leads us to believe that such non-actionable genomic rearrangements may be a result of "chromosomal chaos" most probably not representing exclusive cancer-driving genes in this cancer entity. CONCLUSIONS: We emphasise that caution should be taken when novel fusions are found and that the appearance of new gene fusions should always be interpreted in the molecular context of the respective disease.

KEAP1/NFE2L2 Pathway Signature Outperforms KEAP1/NFE2L2 Mutation Status and Reveals Alternative Pathway-Activating Mutations in NSCLC.

INTRODUCTION: Activation of the antioxidant KEAP1/NFE2L2 (NRF2) pathway leads to increased glutamine dependence and an aggressive phenotype in NSCLC. Because this pathway has been explored as a clinical target, we developed a transcriptomic signature for identifying KEAP1/NFE2L2-activated tumors. METHODS: A total of 971 NSCLC samples were used to train an expression signature (K1N2-score) to predict KEAP1/NFE2L2 mutations. There were 348 in-house NSCLCs that were analyzed using a NanoString expression panel for validation. RESULTS: The 46-gene K1N2 score robustly predicted KEAP1/NFE2L2 mutations in the validation set irrespective of histology and mutation (area under the curve: 89.5, sensitivity: 90.2%), suggesting that approximately 90% of KEAP1/NFE2L2 mutations are pathway-activating. The K1N2-score outperformed KEAP1/NFE2L2 mutational status when predicting patient survival (score p = 0.047; mutation p = 0.215). In K1N2 score-positive but KEAP1/NFE2L2 wild-type samples, enrichment testing identified SMARCA4/BRG1 and CUL3 mutations as mimics of KEAP1/NFE2L2 mutations. CONCLUSIONS: The K1N2-score identified KEAP1/NFE2L2-activated NSCLC by robustly detecting KEAP1/NFE2L2mut cases and discovering alternative genomic activators. It is a potential means for selecting patients with a constitutively active KEAP1/NFE2L2 pathway.

AI-based multi-PRS models outperform classical single-PRS models.

Polygenic risk scores (PRS) calculate the risk for a specific disease based on the weighted sum of associated alleles from different genetic loci in the germline estimated by regression models. Recent advances in genetics made it possible to create polygenic predictors of complex human traits, including risks for many important complex diseases, such as cancer, diabetes, or cardiovascular diseases, typically influenced by many genetic variants, each of which has a negligible effect on overall risk. In the current study, we analyzed whether adding additional PRS from other diseases to the prediction models and replacing the regressions with machine learning models can improve overall predictive performance. Results showed that multi-PRS models outperform single-PRS models significantly on different diseases. Moreover, replacing regression models with machine learning models, i.e., deep learning, can also improve overall accuracy.

Dissecting the genetic heterogeneity of gastric cancer.

BACKGROUND: Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture. METHODS: We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO. FINDINGS: Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level. INTERPRETATION: Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO. FUNDING: German Research Foundation (DFG).

Why loss of Y? A pan-cancer genome analysis of tumors with loss of Y chromosome.

Loss of the Y chromosome (LoY) is frequently observed in somatic cells of elderly men. However, LoY is highly increased in tumor tissue and correlates with an overall worse prognosis. The underlying causes and downstream effects of LoY are widely unknown. Therefore, we analyzed genomic and transcriptomic data of 13 cancer types (2375 patients) and classified tumors of male patients according to loss or retain of the Y chromosome (LoY or RoY, average LoY fraction: 0.46). The frequencies of LoY ranged from almost absence (glioblastoma, glioma, thyroid carcinoma) to 77% (kidney renal papillary cell carcinoma). Genomic instability, aneuploidy, and mutation burden were enriched in LoY tumors. In addition, we found more frequently in LoY tumors the gate keeping tumor suppressor gene TP53 mutated in three cancer types (colon adenocarcinoma, head and neck squamous carcinoma, lung adenocarcinoma) and oncogenes MET, CDK6, KRAS, and EGFR amplified in multiple cancer types. On the transcriptomic level, we observed MMP13, known to be involved in invasion, to be up-regulated in LoY of three adenocarcinomas and down-regulation of the tumor suppressor gene GPC5 in LoY of three cancer types. Furthermore, we found enrichment of a smoking-related mutation signature in LoY tumors of head and neck and lung cancer. Strikingly, we observed a correlation between cancer type-specific sex bias in incidence rates and frequencies of LoY, in line with the hypothesis that LoY increases cancer risk in males. Overall, LoY is a frequent phenomenon in cancer that is enriched in genomically unstable tumors. It correlates with genomic features beyond the Y chromosome and might contribute to higher incidence rates in males.

The extracellular matrix landscape in salivary gland carcinomas is defined by cellular differentiation via expression of three distinct protein modules.

The extracellular matrix (ECM) is an integral part of the tumor microenvironment of carcinomas. Even though salivary gland carcinomas (SGCs) display a range of tumor cell differentiation and distinct extracellular matrices, their ECM landscape has not been characterized in depth. The ECM composition of 89 SGC primaries, 14 metastases, and 25 normal salivary gland tissues was assessed using deep proteomic profiling. Machine learning algorithms and network analysis were used to detect tumor groups and protein modules that explain specific ECM landscapes. Multimodal in situ studies to validate exploratory findings and to infer a putative cellular origin of ECM components were applied. We revealed two fundamental SGC ECM classes which align with the presence or absence of myoepithelial tumor differentiation. We describe the SGC ECM through three biologically distinct protein modules that are differentially expressed across ECM classes and cell types. The modules have a distinct prognostic impact on different SGC types. Since targeted therapy is rarely available for SGC, we used the proteomic expression profile to identify putative therapeutic targets. In summary, we provide the first extensive inventory of ECM components in SGC, a difficult-to-treat disease that encompasses tumors with distinct cellular differentiation. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Common Germline Risk Variants Impact Somatic Alterations and Clinical Features across Cancers.

Aggregation of genome-wide common risk variants, such as polygenic risk score (PRS), can measure genetic susceptibility to cancer. A better understanding of how common germline variants associate with somatic alterations and clinical features could facilitate personalized cancer prevention and early detection. We constructed PRSs from 14 genome-wide association studies (median n = 64,905) for 12 cancer types by multiple methods and calibrated them using the UK Biobank resources (n = 335,048). Meta-analyses across cancer types in The Cancer Genome Atlas (n = 7,965) revealed that higher PRS values were associated with earlier cancer onset and lower burden of somatic alterations, including total mutations, chromosome/arm somatic copy-number alterations (SCNA), and focal SCNAs. This contrasts with rare germline pathogenic variants (e.g., BRCA1/2 variants), showing heterogeneous associations with somatic alterations. Our results suggest that common germline cancer risk variants allow early tumor development before the accumulation of many somatic alterations characteristic of later stages of carcinogenesis. SIGNIFICANCE: Meta-analyses across cancers show that common germline risk variants affect not only cancer predisposition but the age of cancer onset and burden of somatic alterations, including total mutations and copy-number alterations.

GWAS meta-analysis of 16 790 patients with Barrett's oesophagus and oesophageal adenocarcinoma identifies 16 novel genetic risk loci and provides insights into disease aetiology beyond the single marker level.

OBJECTIVE: Oesophageal cancer (EC) is the sixth leading cause of cancer-related deaths. Oesophageal adenocarcinoma (EA), with Barrett's oesophagus (BE) as a precursor lesion, is the most prevalent EC subtype in the Western world. This study aims to contribute to better understand the genetic causes of BE/EA by leveraging genome wide association studies (GWAS), genetic correlation analyses and polygenic risk modelling. DESIGN: We combined data from previous GWAS with new cohorts, increasing the sample size to 16 790 BE/EA cases and 32 476 controls. We also carried out a transcriptome wide association study (TWAS) using expression data from disease-relevant tissues to identify BE/EA candidate genes. To investigate the relationship with reported BE/EA risk factors, a linkage disequilibrium score regression (LDSR) analysis was performed. BE/EA risk models were developed combining clinical/lifestyle risk factors with polygenic risk scores (PRS) derived from the GWAS meta-analysis. RESULTS: The GWAS meta-analysis identified 27 BE and/or EA risk loci, 11 of which were novel. The TWAS identified promising BE/EA candidate genes at seven GWAS loci and at five additional risk loci. The LDSR analysis led to the identification of novel genetic correlations and pointed to differences in BE and EA aetiology. Gastro-oesophageal reflux disease appeared to contribute stronger to the metaplastic BE transformation than to EA development. Finally, combining PRS with BE/EA risk factors improved the performance of the risk models. CONCLUSION: Our findings provide further insights into BE/EA aetiology and its relationship to risk factors. The results lay the foundation for future follow-up studies to identify underlying disease mechanisms and improving risk prediction.

Trophoblast Cell Surface Antigen 2 (TROP2) as a Predictive Bio-Marker for the Therapeutic Efficacy of Sacituzumab Govitecan in Adenocarcinoma of the Esophagus.

INTRODUCTION: The Trophoblast cell surface antigen 2 (TROP2) is expressed in many carcinomas and may represent a target for treatment. Sacituzumab govitecan (SG) is a TROP2-directed antibody-drug conjugate (ADC). Nearly nothing is known about the biological effectiveness of SG in esophageal adenocarcinoma (EAC). MATERIAL AND METHODS: We determined the TROP2 expression in nearly 600 human EAC. In addition, we used the EAC cell lines (ESO-26, OACM5.1C, and FLO-1) and a xenograft mouse model to investigate this relationship. RESULTS: Of 598 human EACs analyzed, 88% showed varying degrees of TROP2 positivity. High TROP2 positive ESO-26 and low TROP2 positive OACM5.1C showed high sensitivity to SG in contrast to negative FLO-1. In vivo, the ESO-26 tumor shows a significantly better response to SG than the TROP2-negative FLO-1 tumor. ESO-26 vital tumor cells show similar TROP2 expression on all carcinoma cells as before therapy initiation, FLO-1 is persistently negative. DISCUSSION: Our data suggest that sacituzumab govitecan is a new therapy option in esophageal adenocarcinoma and the TROP2 expression in irinotecan-naïve EAC correlates with the extent of treatment response by sacituzumab govitecan. TROP2 is emerging as a predictive biomarker in completely TROP2-negative tumors. This should be considered in future clinical trials.

Genome instability is associated with ethnic differences between Asians and Europeans in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the deadliest cancer types with diverse etiological factors across the world. Although large scale genomic studies have been conducted in different countries, integrative analysis of HCC genomes and ethnic comparison across cohorts are lacking. Methods: We first integrated genomes of 1,349 HCC patients from five large cohorts across the world and applied multiple statistical methods in identifying driver genes. Subsequently, we systematically compared HCC genomes and transcriptomes between Asians and Europeans using the TCGA cohort. Results: We identified 29 novel candidate driver genes, many of which are infrequent tumor suppressors driving late-stage tumor progression. When we systematically compared ethnic differences in the genomic landscape between Asian and European HCCs using the TCGA cohort (n = 348), we found little differences in driver frequencies. Through multi-modal integrative analysis, we found higher genomic instability in Asians together with a collection of molecular events ranging from tumor mutation burden (TMB), copy number alterations as well as transcriptomic subtypes segregating distinctively between two ethnic backgrounds. Strikingly, we identified an Asian specific transcriptomic subtype with multiple ethnically enriched genomic alterations, in particular chromosome 16 deletion, leading to a clinically aggressive RNA subgroup unique to Asians. Integrating multi-modal information, we found that survival models predict patient prognosis much better in Asians than in Europeans, demonstrating a higher potential for precision medicine applications in Asia. Conclusion: For the first time, we have uncovered an unprecedented amount of genomic differences segregating distinctively across ethnicities in HCC and highlighted the importance of differential disease biology and management in HCC across ethnic backgrounds.

Observations that suggest a contribution of altered dermal papilla mitochondrial function to androgenetic alopecia.

Androgenetic alopecia (AGA) is a prevalent hair loss condition in males that develops due to the influence of androgens and genetic predisposition. With the aim of elucidating genes involved in AGA pathogenesis, we modelled AGA with three-dimensional culture of keratinocyte-surrounded dermal papilla (DP) cells. We co-cultured immortalised balding and non-balding human DP cells (DPCs) derived from male AGA patients with epidermal keratinocyte (NHEK) using multi-interfacial polyelectrolyte complexation technique. We observed up-regulated mitochondria-related gene expression in balding compared with non-balding DP aggregates which indicated altered mitochondria metabolism. Further observation of significantly reduced electron transport chain complex activity (complexes I, IV and V), ATP levels and ability to uptake metabolites for ATP generation demonstrated compromised mitochondria function in balding DPC. Balding DP was also found to be under significantly higher oxidative stress than non-balding DP. Our experiments suggest that application of antioxidants lowers oxidative stress levels and improves metabolite uptake in balding DPC. We postulate that the observed up-regulation of mitochondria-related genes in balding DP aggregates resulted from an over-compensatory effort to rescue decreased mitochondrial function in balding DP through the attempted production of new functional mitochondria. In all, our three-dimensional co-culturing revealed mitochondrial dysfunction in balding DPC, suggesting a metabolic component in the aetiology of AGA.