Rare germline structural variants increase risk for pediatric solid tumors.

Pediatric solid tumors are rare malignancies that represent a leading cause of death by disease among children in developed countries. The early age-of-onset of these tumors suggests that germline genetic factors are involved, yet conventional germline testing for short coding variants in established predisposition genes only identifies pathogenic events in 10-15% of patients. Here, we examined the role of germline structural variants (SVs)-an underexplored form of germline variation-in pediatric extracranial solid tumors using germline genome sequencing of 1,766 affected children, their 943 unaffected relatives, and 6,665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and a four-fold increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we revealed burdens of ultra-rare SVs that cause loss-of-function of highly expressed, mutationally intolerant, neurodevelopmental genes, as well as noncoding SVs predicted to disrupt three-dimensional chromatin domains in neural crest-derived tissues. Collectively, our results implicate rare germline SVs as a predisposing factor to pediatric solid tumors that may guide future studies and clinical practice.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics.

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Development, evaluation and validation of machine learning models to predict hospitalizations of patients with coronary artery disease within the next 12 months.

BACKGROUND: Improved survival of patients after acute coronary syndromes, population growth, and overall life expectancy rise have led to a significant increase in the proportion of patients with stable coronary artery disease (CAD), creating a significant load on the entire healthcare system. The disease often progresses with the development of many complications while significantly increasing the likelihood of hospitalization. Developing and applying a machine learning model for predicting hospitalizations of patients with CAD to an inpatient medical facility will allow for close monitoring of high-risk patients, early preventive interventions, and optimized medical care. AIMS: Development and external validation of personalized models for predicting the preventable hospitalizations of patients with stable CAD and its complications using ML algorithms and data of real-world clinical practice. METHODS: 135,873 depersonalized electronic health records of 49,103 patients with stable CAD were included in the study. Anthropometric measurements, physical examination results, laboratory, instrumental, anamnestic, and socio-demographic data, widely used in routine medical practice, were considered as potential predictors, a total of 73 features. Logistic regression, decision tree-based methods including gradient boosting (AdaBoost, LightGBM, XGBoost, CatBoost) and bagging (RandomForest and ExtraTrees), discriminant analysis (LinearDiscriminant, QuadraticDiscriminant), and naive Bayes classifier were compared. External validation was performed on the data of a separate region. RESULTS: The best results and stability to external validation data were shown by the CatBoost model with an AUC of 0.875 (95% CI 0.865-0.885) for the internal testing and 0.872 (95% CI 0.856-0.886) for the external validation. The best model showed good performance evaluated through AUROC, Brier score and standardized net benefit (for the target NPV threshold) for the validation dataset that was only slightly similar to the train data. CONCLUSION: The metrics of the best model were superior to previously published studies. The results of external validation demonstrated the relative stability of the model to new data from another region that confirms the possibility of the model's application in real clinical practice.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk.

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Functional dissection of complex and molecular trait variants at single nucleotide resolution.

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

Genetic and molecular architecture of complex traits.

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

Cancer Type and Histology Influence Cutaneous Immunotherapy Toxicities: A Multi-Institutional Cohort Study.

BACKGROUND: Cutaneous immune-related adverse events (cirAEs) are the most common toxicities to occur in the setting of immune checkpoint inhibitor (ICI) therapy. Identifying patients who are at increased risk of developing cirAEs may improve quality of life and outcomes. OBJECTIVE: To investigate the influence of cancer type and histology on the development of cirAEs in the setting of ICI therapy and survival outcomes. METHODS: This retrospective cohort study included patients between December 1, 2011, and October 30, 2020. They received ICIs from 2011-2020 with follow-up of outcomes through October 2021. We identified 3,668 ICI recipients who were seen at Mass General Brigham and Dana-Farber. Of these, 669 developed cirAEs. Records that were incomplete or categories of insufficient sample size were excluded from the study cohort. Multivariate Cox proportional hazards models were utilized to investigate the impact of cancer organ system and histology on cirAE development, after adjusting for demographics, Charlson Comorbidity Index, ICI type, cancer stage at ICI initiation, and year of ICI initiation. Time-varying Cox proportional hazards modeling was utilized to examine the impact of cirAE development on mortality. RESULTS: Compared to other non-epithelial cancers (neuroendocrine, leukemia, lymphoma, myeloma, sarcoma, and central nervous system malignancies), cutaneous squamous cell carcinoma (cSCC) (HR = 3.57, p < 0.001), melanoma (HR = 2.09, p < 0.001), head and neck adenocarcinoma (HR = 2.13, p = 0.009), genitourinary transitional cell carcinoma (HR = 2.15, p < 0.001), and genitourinary adenocarcinoma (HR = 1.53, p = 0.037) were at significantly higher risk of cirAEs in multivariate analyses. The increased risk of cirAEs translated into an adjusted survival benefit for melanoma (HR = 0.37, p < 0.001) and cSCC (HR = 0.51, p = 0.011). CONCLUSIONS: The highest rate of cirAEs and subsequent survival benefits were observed in cutaneous malignancies treated with ICIs. This study improves our understanding of patients who are at highest risk of developing cirAEs and would, therefore, benefit from appropriate counseling and closer monitoring by their oncologists and dermatologists throughout their ICI therapy. Limitations include its retrospective nature and cohort from one geography.

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.

Lung cancer in patients who have never smoked - an emerging disease.

Lung cancer is the most common cause of cancer-related deaths globally. Although smoking-related lung cancers continue to account for the majority of diagnoses, smoking rates have been decreasing for several decades. Lung cancer in individuals who have never smoked (LCINS) is estimated to be the fifth most common cause of cancer-related deaths worldwide in 2023, preferentially occurring in women and Asian populations. As smoking rates continue to decline, understanding the aetiology and features of this disease, which necessitate unique diagnostic and treatment paradigms, will be imperative. New data have provided important insights into the molecular and genomic characteristics of LCINS, which are distinct from those of smoking-associated lung cancers and directly affect treatment decisions and outcomes. Herein, we review the emerging data regarding the aetiology and features of LCINS, particularly the genetic and environmental underpinnings of this disease as well as their implications for treatment. In addition, we outline the unique diagnostic and therapeutic paradigms of LCINS and discuss future directions in identifying individuals at high risk of this disease for potential screening efforts.

Obesity shapes selection for driver mutations in cancer.

Obesity is a leading risk factor for cancer, but whether obesity is linked to specific genomic subtypes of cancer is unknown. Here, we examined the relationship between obesity and tumor genotype in two large clinicogenomic corpora. Obesity was associated with specific driver mutations in lung adenocarcinoma, endometrial carcinoma, and cancers of unknown primary, independent of clinical covariates and genetic ancestry. Obesity is therefore a putative driver of etiologic heterogeneity across cancers.

Development and validation of time-to-event models to predict metastatic recurrence of localized cutaneous melanoma.

BACKGROUND: The recent expansion of immunotherapy for stage IIB/IIC melanoma highlights a growing clinical need to identify patients at high risk of metastatic recurrence and, therefore, most likely to benefit from this therapeutic modality. OBJECTIVE: To develop time-to-event risk prediction models for melanoma metastatic recurrence. METHODS: Patients diagnosed with stage I/II primary cutaneous melanoma between 2000 and 2020 at Mass General Brigham and Dana-Farber Cancer Institute were included. Melanoma recurrence date and type were determined by chart review. Thirty clinicopathologic factors were extracted from electronic health records. Three types of time-to-event machine-learning models were evaluated internally and externally in the distant versus locoregional/nonrecurrence prediction. RESULTS: This study included 954 melanomas (155 distant, 163 locoregional, and 636 1:2 matched nonrecurrences). Distant recurrences were associated with worse survival compared to locoregional/nonrecurrences (HR: 6.21, P < .001) and to locoregional recurrences only (HR: 5.79, P < .001). The Gradient Boosting Survival model achieved the best performance (concordance index: 0.816; time-dependent AUC: 0.842; Brier score: 0.103) in the external validation. LIMITATIONS: Retrospective nature and cohort from one geography. CONCLUSIONS: These results suggest that time-to-event machine-learning models can reliably predict the metastatic recurrence from localized melanoma and help identify high-risk patients who are most likely to benefit from immunotherapy.

Prediction and stratification of longitudinal risk for chronic obstructive pulmonary disease across smoking behaviors.

Smoking is the leading risk factor for chronic obstructive pulmonary disease (COPD) worldwide, yet many people who never smoke develop COPD. We perform a longitudinal analysis of COPD in the UK Biobank to derive and validate the Socioeconomic and Environmental Risk Score which captures additive and cumulative environmental, behavioral, and socioeconomic exposure risks beyond tobacco smoking. The Socioeconomic and Environmental Risk Score is more predictive of COPD than smoking status and pack-years. Individuals in the highest decile of the risk score have a greater risk for incident COPD compared to the remaining population. Never smokers in the highest decile of exposure risk are more likely to develop COPD than previous and current smokers in the lowest decile. In general, the prediction accuracy of the Social and Environmental Risk Score is lower in non-European populations. While smoking status is often considered in screening COPD, our finding highlights the importance of other non-smoking environmental and socioeconomic variables.

MSGene: Derivation and validation of a multistate model for lifetime risk of coronary artery disease using genetic risk and the electronic health record.

Currently, coronary artery disease (CAD) is the leading cause of death among adults worldwide. Accurate risk stratification can support optimal lifetime prevention. We designed a novel and general multistate model (MSGene) to estimate age-specific transitions across 10 cardiometabolic states, dependent on clinical covariates and a CAD polygenic risk score. MSGene supports decision making about CAD prevention related to any of these states. We analyzed longitudinal data from 480,638 UK Biobank participants and compared predicted lifetime risk with the 30-year Framingham risk score. MSGene improved discrimination (C-index 0.71 vs 0.66), age of high-risk detection (C-index 0.73 vs 0.52), and overall prediction (RMSE 1.1% vs 10.9%), with external validation. We also used MSGene to refine estimates of lifetime absolute risk reduction from statin initiation. Our findings underscore the potential public health value of our novel multistate model for accurate lifetime CAD risk estimation using clinical factors and increasingly available genetics.

Somatic mutations reveal hyperactive Notch signaling and racial disparities in prurigo nodularis.

Prurigo nodularis (PN) is a chronic inflammatory skin disease that disproportionately affects African Americans and is characterized by pruritic skin nodules of unknown etiology. Little is known about genetic alterations in PN pathogenesis, especially relating to somatic events which are often implicated in inflammatory conditions. We thus performed whole-exome sequencing on 54 lesional and nonlesional skin biopsies from 17 PN patients and 10 atopic dermatitis (AD) patients for comparison. Somatic mutational analysis revealed that PN lesional skin harbors pervasive somatic mutations in fibrotic, neurotropic, and cancer-associated genes. Nonsynonymous mutations were most frequent in NOTCH1 and the Notch signaling pathway, a regulator of cellular proliferation and tissue fibrosis, and NOTCH1 mutations were absent in AD. Somatic copy-number analysis, combined with expression data, showed that recurrently deleted and downregulated genes in PN lesional skin are associated with axonal guidance and extension. Follow-up immunofluorescence validation demonstrated increased NOTCH1 expression in PN lesional skin fibroblasts and increased Notch signaling in PN lesional dermis. Finally, multi-center data revealed a significantly increased risk of NOTCH1-associated diseases in PN patients. In characterizing the somatic landscape of PN, we uncover novel insights into its pathophysiology and identify a role for dysregulated Notch signaling in PN.

Scalable genetic screening for regulatory circuits using compressed Perturb-seq.

Pooled CRISPR screens with single-cell RNA sequencing readout (Perturb-seq) have emerged as a key technique in functional genomics, but they are limited in scale by cost and combinatorial complexity. In this study, we modified the design of Perturb-seq by incorporating algorithms applied to random, low-dimensional observations. Compressed Perturb-seq measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq with an order of magnitude cost reduction and greater power to learn genetic interactions. We identified known and novel regulators of immune responses and uncovered evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing genome-wide association studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

Pre-Existing Inflammatory Disease Predicts Cutaneous Immunotherapy Toxicity Development: A Multi-Institutional Cohort Study.

Background: Relationships between pre-existing inflammatory diseases (pIDs) and cutaneous immune-related adverse events (cirAEs) have not been well-studied. This study is to investigate associations between pIDs and cirAEs among immune-checkpoint inhibitor (ICI) recipients at the Mass General Brigham healthcare system. Methods: Electronic health records were reviewed to ascertain cirAE status. Patients' pID status was determined using International Classification of Diseases (ICD) codes. Cox proportional hazard, logistic regression, and linear regression models were performed. Results: Among 3607 ICI recipients, 1354 had pIDs, and 672 developed cirAEs. After covariate adjustments, patients with cutaneous pIDs (HR:1.56, p<0.001) or both cutaneous and non-cutaneous pIDs (HR:1.76, p<0.001) had increased cirAE risk in contrast to patients with non-cutaneous pIDs alone (HR:1.01, p=0.9). In adjusted ordinal logistic regression modeling, cutaneous pIDs (OR:1.55, p<0.0001) and the presence of both cutaneous pIDs and non-cutaneous pIDs (OR:1.71, p=0.002) were associated with increased cirAE severity. The time to cirAE onset was different between the cutaneous pID group and the non-cutaneous pID group (Mean: 98 vs. 146 days, p=0.021; Beta: -0.11, p=0.033). Conclusions: ICI recipients with cutaneous pIDs should have increased clinical monitoring due to their increased risk of cirAE development, severity, and earlier onset.