Somatic mutational profiles and germline polygenic risk scores in human cancer

Germline polygenic risk scores are associated with immune gene expression signature and immune cell infiltration in breast cancer

The tumor immune microenvironment (TIME) plays key roles in tumor progression and response to immunotherapy. Previous studies have identified individual germline variants associated with differences in TIME. Here, we hypothesize that common variants associated with breast cancer risk or cancer-related traits, represented by polygenic risk scores (PRSs), may jointly influence immune features in TIME. We derived 154 immune traits from bulk gene expression profiles of 764 breast tumors and 598 adjacent normal tissue samples from 825 individuals with breast cancer in the Nurses' Health Study (NHS) and NHSII. Immunohistochemical staining of four immune cell markers were available for a subset of 205 individuals. Germline PRSs were calculated for 16 different traits including breast cancer, autoimmune diseases, type 2 diabetes, ages at menarche and menopause, body mass index (BMI), BMI-adjusted waist-to-hip ratio, alcohol intake, and tobacco smoking. Overall, we identified 44 associations between germline PRSs and immune traits at false discovery rate q < 0.25, including 3 associations with q < 0.05. We observed consistent inverse associations of inflammatory bowel disease (IBD) and Crohn disease (CD) PRSs with interferon signaling and STAT1 scores in breast tumor and adjacent normal tissue; these associations were replicated in a Norwegian cohort. Inverse associations were also consistently observed for IBD PRS and B cell abundance in normal tissue. We also observed positive associations between CD PRS and endothelial cell abundance in tumor. Our findings suggest that the genetic mechanisms that influence immune-related diseases are also associated with TIME in breast cancer.

Germline mutation rate predicts cancer mortality across 37 vertebrate species

Background and objectives

Cancer develops across nearly every species. However, cancer occurs at unexpected and widely different rates throughout the animal kingdom. The reason for this variation in cancer susceptibility remains an area of intense investigation. Cancer evolves in part through the accumulation of mutations, and therefore, we hypothesized that germline mutation rates would be associated with cancer prevalence and mortality across species.

Methodology

We collected previously published data on germline mutation rate and cancer mortality data for 37 vertebrate species.

Results

Germline mutation rate was positively correlated with cancer mortality (P-value = 0.0008; R2 = 0.13). Controlling for species' average parental age, maximum longevity, adult body mass or domestication did not improve the model fit (the change (Δ) in Akaike Information Criterion (AIC) was less than 2). However, this model fit was better than a model controlling for species trophic level (ΔAIC > 2).

Conclusions and implications

The increased death rate from cancer in animals with increased germline mutation rates may suggest underlying hereditary cancer predisposition syndromes similar to those diagnosed in human patients. Species with higher germline mutation rates may benefit from close monitoring for tumors due to increased genetic risk for cancer development. Early diagnoses of cancer in these species may increase their chances of overall survival, especially for threatened and endangered species.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

Molecular and immune characterization of Chinese early-stage non-squamous non-small cell lung cancer: a multi-omics cohort study

Background

Albeit considered with superior survival, around 30% of the early-stage non-squamous non-small cell lung cancer (Ns-NSCLC) patients relapse within 5 years, suggesting unique biology. However, the biological characteristics of early-stage Ns-NSCLC, especially in the Chinese population, are still unclear.

Methods

Multi-omics interrogation of early-stage Ns-NSCLC (stage I-III), paired blood samples and normal lung tissues (n=76) by whole-exome sequencing (WES), RNA sequencing, and T-cell receptor (TCR) sequencing were conducted.

Results

An average of 128 exonic mutations were identified, and the most frequently mutant gene was EGFR (55%), followed by TP53 (37%) and TTN (26%). Mutations in MUC17, ABCA2, PDE4DIP, and MYO18B predicted significantly unfavorable disease-free survival (DFS). Moreover, cytobands amplifications in 8q24.3, 14q13.1, 14q11.2, and deletion in 3p21.1 were highlighted in recurrent cases. Higher incidence of human leukocyte antigen loss of heterozygosity (HLA-LOH), higher tumor mutational burden (TMB) and tumor neoantigen burden (TNB) were identified in ever-smokers than never-smokers. HLA-LOH also correlated with higher TMB, TNB, intratumoral heterogeneity (ITH), and whole chromosomal instability (wCIN) scores. Interestingly, higher ITH was an independent predictor of better DFS in early-stage Ns-NSCLC. Up-regulation of immune-related genes, including CRABP2, ULBP2, IL31RA, and IL1A, independently portended a dismal prognosis. Enhanced TCR diversity of peripheral blood mononuclear cells (PBMCs) predicted better prognosis, indicative of a noninvasive method for relapse surveillance. Eventually, seven machine-learning (ML) algorithms were employed to evaluate the predictive accuracy of clinical, genomic, transcriptomic, and TCR repertoire data on DFS, showing that clinical and RNA features combination in the random forest (RF) algorithm, with area under the curve (AUC) of 97.5% and 83.3% in the training and testing cohort, respectively, significantly outperformed other methods.

Conclusions

This study comprehensively profiled the genomic, transcriptomic, and TCR repertoire spectrums of Chinese early-stage Ns-NSCLC, shedding light on biological underpinnings and candidate biomarkers for prognosis development.

A comprehensive analysis of clinical and polygenic germline influences on somatic mutational burden

Tumor mutational burden (TMB), the total number of somatic mutations in the tumor, and copy number burden (CNB), the corresponding measure of aneuploidy, are established fundamental somatic features and emerging biomarkers for immunotherapy. However, the genetic and non-genetic influences on TMB/CNB and, critically, the manner by which they influence patient outcomes remain poorly understood. Here, we present a large germline-somatic study of TMB/CNB with >23,000 individuals across 17 cancer types, of which 12,000 also have extensive clinical, treatment, and overall survival (OS) measurements available. We report dozens of clinical associations with TMB/CNB, observing older age and male sex to have a strong effect on TMB and weaker impact on CNB. We additionally identified significant germline influences on TMB/CNB, including fine-scale European ancestry and germline polygenic risk scores (PRSs) for smoking, tanning, white blood cell counts, and educational attainment. We quantify the causal effect of exposures on somatic mutational processes using Mendelian randomization. Many of the identified features associated with TMB/CNB were additionally associated with OS for individuals treated at a single tertiary cancer center. For individuals receiving immunotherapy, we observed a complex relationship between PRSs for educational attainment, self-reported college attainment, TMB, and survival, suggesting that the influence of this biomarker may be substantially modified by socioeconomic status. While the accumulation of somatic alterations is a stochastic process, our work demonstrates that it can be shaped by host characteristics including germline genetics.

A Multiethnic Germline-Somatic Association Database Deciphers Multilayered and Interconnected Genetic Mutations in Cancer

Inherited germline and acquired somatic alterations can both promote human tumor development. Elucidating the cooperation between somatic and germline genetic alterations that drive tumorigenesis could help inform precision cancer prevention and treatment strategies. Here, leveraging genomic genotyping and sequencing data from 9,029 patients with cancer with European, East Asian, and African ancestry, we performed a pan-cancer analysis to evaluate the associations between germline SNPs and somatic alterations, including single-nucleotide variant and small insertion/deletion mutations, copy-number variation, tumor mutational burden, and mutational signatures. Genome-wide significant germline-somatic pairs were abundant, and most of the associations were observed in one cancer type and one ancestry group. A user-friendly interactive Multiethnic Germline-Somatic Association (MGSA) database (http://wanglab-hkust.cn:3838/MGSA/) was developed, which can be used to query, browse, and download the results of the association analyses. Moreover, the MGSA database offers additional survival analysis and functional annotation. Together, this work provides a resource for uncovering the clinical and biological roles of associations between germline variants and somatic alterations in human cancer.

SIGNIFICANCE

Comprehensive analysis of connections between germline variants and somatic events in cancer offers a resource for investigating the functional significance of genetic mutations and exploring genetic factors contributing to racial disparities.

Genetic insights into carbohydrate sulfotransferase 8 and its impact on the immunotherapy efficacy of cancer

Immune checkpoint blockade (ICB) is a promising therapy for solid tumors, but its effectiveness depends on biomarkers that are not precise. Here, we utilized genome-wide association study to investigate the association between genetic variants and tumor mutation burden to interpret ICB response. We identified 16 variants (p < 5 × 10-8) probed to 17 genes on 9 chromosomes. Subsequent analysis of one of the most significant loci in 19q13.11 suggested that the rs111308825 locus at the enhancer is causal, as its A allele impairs KLF2 binding, leading to lower carbohydrate sulfotransferase 8 (CHST8) expression. Breast cancer cells expressing CHST8 suppress T cell activation, and Chst8 loss attenuates tumor growth in a syngeneic mouse model. Further investigation revealed that programmed death-ligand 1 (PD-L1) and its homologs could be sulfated by CHST8, resulting in M2-like macrophage enrichment in the tumor microenvironment. Finally, we confirmed that low-CHST8 tumors have better ICB response, supporting the genetic effect and clinical value of rs111308825 for ICB efficacy prediction.

The role of mitochondria-related lncRNAs in characterizing the immune landscape and supervising the prognosis of osteosarcoma

Mitochondrial damage is related to the functional properties of immune cells as well as to tumorigenesis and progression. Nevertheless, there is an absence concerning the systematic evaluation of mitochondria-associated lncRNAs (MALs) in the immune profile and tumor microenvironment of osteosarcoma patients. Based on transcriptomic and clinicopathological data from the TARGET database, MAL-related patterns were ascertained by consistent clustering, and gene set variation analysis of the different patterns was completed. Next, a MAL-derived scoring system was created using Cox and LASSO regression analyses and validated by Kaplan-Meier and ROC curves. The GSEA, ESTIMATE, and CIBERSORT algorithms were utilized to characterize the immune status and underlying biological functions in the different MAL score groups. MAL-derived risk scores were well stabilized and outperformed traditional clinicopathological features to reliably predict 5-year survival in osteosarcoma cohorts. Moreover, patients with increased MAL scores were observed to suffer from poorer prognosis, higher tumor purity, and an immunosuppressive microenvironment. Based on estimated half-maximal inhibitory concentrations, the low-MAL score group benefited more from gemcitabine and docetaxel, and less from thapsigargin and sunitinib compared to the high-MAL score group. Pan-cancer analysis demonstrated that six hub MALs were strongly correlated with clinical outcomes, immune subtypes, and tumor stemness indices in various common cancers. Finally, we verified the expression patterns of hub MALs in osteosarcoma with qRT-PCR. In summary, we identified the crosstalk between prognostic MALs and tumor-infiltrating immune cells in osteosarcoma, providing a potential strategy to ameliorate clinical stratification management.

Population-specific Mutation Patterns in Breast Tumors from African American, European American, and Kenyan Patients

Women of African descent have the highest breast cancer mortality in the United States and are more likely than women from other population groups to develop an aggressive disease. It remains uncertain to what extent breast cancer in Africa is reminiscent of breast cancer in African American or European American patients. Here, we performed whole-exome sequencing of genomic DNA from 191 breast tumor and non-cancerous adjacent tissue pairs obtained from 97 African American, 69 European American, 2 Asian American, and 23 Kenyan patients. Our analysis of the sequencing data revealed an elevated tumor mutational burden in both Kenyan and African American patients, when compared with European American patients. TP53 mutations were most prevalent, particularly in African American patients, followed by PIK3CA mutations, which showed similar frequencies in European American, African American, and the Kenyan patients. Mutations targeting TBX3 were confined to European Americans and those targeting the FBXW7 tumor suppressor to African American patients whereas mutations in the ARID1A gene that are known to confer resistance to endocrine therapy were distinctively enriched among Kenyan patients. A Kyoto Encyclopedia of Genes and Genomes pathway analysis could link FBXW7 mutations to an increased mitochondrial oxidative phosphorylation capacity in tumors carrying these mutations. Finally, Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures in tumors correlated with the occurrence of driver mutations, immune cell profiles, and neighborhood deprivation with associations ranging from being mostly modest to occasionally robust. To conclude, we found mutational profiles that were different between these patient groups. The differences concentrated among genes with low mutation frequencies in breast cancer.

SIGNIFICANCE

The study describes differences in tumor mutational profiles between African American, European American, and Kenyan breast cancer patients. It also investigates how these profiles may relate to the tumor immune environment and the neighborhood environment in which the patients had residence. Finally, it describes an overrepresentation of ARID1A gene mutations in breast tumors of the Kenyan patients.

Secondary Infection Surveillance with Metagenomic Next-Generation Sequencing in COVID-19 Patients: A Cross-Sectional Study

Background

Metagenomic next-generation sequencing (mNGS) is a promising tool for improving antimicrobial therapy and infection control decision-making in complex infections. Secondary infection surveillance using mNGS in COVID-19 patients has rarely been reported.

Methods

Respiratory pathogen and antibiotic resistance prediction were evaluated by BALF mNGS for 192 hospitalized COVID-19 patients between December 2022 and February 2023.

Results

Secondary infection was confirmed in 83.3% (160/192) of the COVID-19 patients, with bacterial infections (45%, 72/160) predominating, followed by mixed bacterial and fungal infections (20%, 32/160), and fungal infections (17.5%, 28/160). The incidence of bacterial or viral secondary infection was significantly higher in patients who were admitted to the ICU, received mechanical ventilation, or developed severe pneumonia (all p<0.05). Klebsiella pneumoniae (n=30, 8.4%) was the most prevalent pathogen associated with secondary infection followed by Acinetobacter baumannii (n=29, 8.1%), Candida albicans (n=29, 8.1%), Aspergillus fumigatus (n=27, 7.6%), human herpes simplex virus type 1 (n=23, 6.4%), Staphylococcus aureus (n=20, 5.6%) and Pneumocystis jiroveci (n=14, 3.9%). The overall concordance between the resistance genes detected by mNGS and the reported phenotypic resistance in 69 samples containing five clinically important pathogens (ie, K. pneumoniae, A. baumannii, S. aureus, P. aeruginosa and E. coli) that caused secondary infection was 85.5% (59/69).

Conclusion

mNGS can detect pathogens causing secondary infection and predict antimicrobial resistance for COVID19 patients. This is crucial for initiating targeted treatment and rapidly detect unsuspected spread of multidrug-resistant pathogens.

Carcino-Evo-Devo, A Theory of the Evolutionary Role of Hereditary Tumors

A theory of the evolutionary role of hereditary tumors, or the carcino-evo-devo theory, is being developed. The main hypothesis of the theory, the hypothesis of evolution by tumor neofunctionalization, posits that hereditary tumors provided additional cell masses during the evolution of multicellular organisms for the expression of evolutionarily novel genes. The carcino-evo-devo theory has formulated several nontrivial predictions that have been confirmed in the laboratory of the author. It also suggests several nontrivial explanations of biological phenomena previously unexplained by the existing theories or incompletely understood. By considering three major types of biological development-individual, evolutionary, and neoplastic development-within one theoretical framework, the carcino-evo-devo theory has the potential to become a unifying biological theory.

Genomic Landscape of Normal and Breast Cancer Tissues in a Hungarian Pilot Cohort

A limited number of studies have focused on the mutational landscape of breast cancer in different ethnic populations within Europe and compared the data with other ethnic groups and databases. We performed whole-genome sequencing of 63 samples from 29 Hungarian breast cancer patients. We validated a subset of the identified variants at the DNA level using the Illumina TruSight Oncology (TSO) 500 assay. Canonical breast-cancer-associated genes with pathogenic germline mutations were CHEK2 and ATM. Nearly all the observed germline mutations were as frequent in the Hungarian breast cancer cohort as in independent European populations. The majority of the detected somatic short variants were single-nucleotide polymorphisms (SNPs), and only 8% and 6% of them were deletions or insertions, respectively. The genes most frequently affected by somatic mutations were KMT2C (31%), MUC4 (34%), PIK3CA (18%), and TP53 (34%). Copy number alterations were most common in the NBN, RAD51C, BRIP1, and CDH1 genes. For many samples, the somatic mutational landscape was dominated by mutational processes associated with homologous recombination deficiency (HRD). Our study, as the first breast tumor/normal sequencing study in Hungary, revealed several aspects of the significantly mutated genes and mutational signatures, and some of the copy number variations and somatic fusion events. Multiple signs of HRD were detected, highlighting the value of the comprehensive genomic characterization of breast cancer patient populations.

Germline Cancer Gene Expression Quantitative Trait Loci Are Associated with Local and Global Tumor Mutations

Somatic mutations drive cancer development and are relevant to patient responses to treatment. Emerging evidence shows that variations in the somatic genome can be influenced by the germline genetic background. However, the mechanisms underlying these germline-somatic associations remain largely obscure. We hypothesized that germline variants can influence somatic mutations in a nearby cancer gene ("local impact") or a set of recurrently mutated cancer genes across the genome ("global impact") through their regulatory effect on gene expression. To test this hypothesis, tumor targeted sequencing data from 12,413 patients across 11 cancer types in the Dana-Farber Profile cohort were integrated with germline cancer gene expression quantitative trait loci (eQTL) from the Genotype-Tissue Expression Project. Variants that upregulate ATM expression were associated with a decreased risk of somatic ATM mutations across 8 cancer types. GLI2, WRN, and CBFB eQTL were associated with global tumor mutational burden of cancer genes in ovarian cancer, glioma, and esophagogastric carcinoma, respectively. An EPHA5 eQTL was associated with mutations in cancer genes specific to colorectal cancer, and eQTL related to expression of APC, WRN, GLI1, FANCA, and TP53 were associated with mutations in genes specific to endometrial cancer. These findings provide evidence that germline-somatic associations are mediated through expression of specific cancer genes, opening new avenues for research on the underlying biological processes.

SIGNIFICANCE

Analysis of associations between the germline genetic background and somatic mutations in patients with cancer suggests that germline variants can influence local and global tumor mutations by altering expression of cancer-related genes. See related commentary by Kar, p. 1165.