Integration of rare expression outlier-associated variants improves polygenic risk prediction

Misexpression of inactive genes in whole blood is associated with nearby rare structural variants

Gene misexpression is the aberrant transcription of a gene in a context where it is usually inactive. Despite its known pathological consequences in specific rare diseases, we have a limited understanding of its wider prevalence and mechanisms in humans. To address this, we analyzed gene misexpression in 4,568 whole-blood bulk RNA sequencing samples from INTERVAL study blood donors. We found that while individual misexpression events occur rarely, in aggregate they were found in almost all samples and a third of inactive protein-coding genes. Using 2,821 paired whole-genome and RNA sequencing samples, we identified that misexpression events are enriched in cis for rare structural variants. We established putative mechanisms through which a subset of SVs lead to gene misexpression, including transcriptional readthrough, transcript fusions, and gene inversion. Overall, we develop misexpression as a type of transcriptomic outlier analysis and extend our understanding of the variety of mechanisms by which genetic variants can influence gene expression.

A polygenic score method boosted by non-additive models

Dominance heritability in complex traits has received increasing recognition. However, most polygenic score (PGS) approaches do not incorporate non-additive effects. Here, we present GenoBoost, a flexible PGS modeling framework capable of considering both additive and non-additive effects, specifically focusing on genetic dominance. Building on statistical boosting theory, we derive provably optimal GenoBoost scores and provide its efficient implementation for analyzing large-scale cohorts. We benchmark it against seven commonly used PGS methods and demonstrate its competitive predictive performance. GenoBoost is ranked the best for four traits and second-best for three traits among twelve tested disease outcomes in UK Biobank. We reveal that GenoBoost improves prediction for autoimmune diseases by incorporating non-additive effects localized in the MHC locus and, more broadly, works best in less polygenic traits. We further demonstrate that GenoBoost can infer the mode of genetic inheritance without requiring prior knowledge. For example, GenoBoost finds non-zero genetic dominance effects for 602 of 900 selected genetic variants, resulting in 2.5% improvements in predicting psoriasis cases. Lastly, we show that GenoBoost can prioritize genetic loci with genetic dominance not previously reported in the GWAS catalog. Our results highlight the increased accuracy and biological insights from incorporating non-additive effects in PGS models.

Donor and recipient genetics: implications for the development of post-transplant diabetes mellitus

Post-transplant diabetes mellitus (PTDM) is a prevalent complication of liver transplantation and is associated with cardiometabolic complications. We studied the consequences of genetic effects of liver donors and recipients on PTDM outcomes, focusing on the diverse genetic pathways related to insulin that play a role in the development of PTDM. 1115 liver transplant recipients without a pre-transplant diagnosis of type 2 diabetes mellitus (T2D) and their paired donors recruited from two transplant centers had polygenic risk scores (PRS) for T2D, insulin secretion, and insulin sensitivity calculated. Among recipients in the highest T2D-PRS quintile, donor T2D-PRS did not contribute significantly to PTDM. However, in recipients with the lowest T2D genetic risk, donor livers with the highest T2D-PRS contributed to the development of PTDM (OR (95% CI)=3.79 (1.10-13.1), p=0.035). Recipient risk was linked to factors associated with insulin secretion (OR (95% CI) = 0.85 (0.74-0.98), p=0.02), while donor livers contributed to PTDM via gene pathways involved in insulin sensitivity (OR (95% CI)=0.86 (0.75-0.99), p=0.03). Recipient and donor PRS independently and collectively serve as predictors of PTDM onset. The genetically influenced biological pathways in recipients primarily pertain to insulin secretion, whereas the genetic makeup of donors exerts an influence on insulin sensitivity.

Genetic modifiers of rare variants in monogenic developmental disorder loci

Rare damaging variants in a large number of genes are known to cause monogenic developmental disorders (DDs) and have also been shown to cause milder subclinical phenotypes in population cohorts. Here, we show that carrying multiple (2-5) rare damaging variants across 599 dominant DD genes has an additive adverse effect on numerous cognitive and socioeconomic traits in UK Biobank, which can be partially counterbalanced by a higher educational attainment polygenic score (EA-PGS). Phenotypic deviators from expected EA-PGS could be partly explained by the enrichment or depletion of rare DD variants. Among carriers of rare DD variants, those with a DD-related clinical diagnosis had a substantially lower EA-PGS and more severe phenotype than those without a clinical diagnosis. Our results suggest that the overall burden of both rare and common variants can modify the expressivity of a phenotype, which may then influence whether an individual reaches the threshold for clinical disease.

Evaluation of Genetic or Cellular Impairments in Type I IFN Immunity in a Cohort of Young Adults with Critical COVID-19

Several genetic and immunological risk factors for severe COVID-19 have been identified, with monogenic conditions relating to 13 genes of type I interferon (IFN) immunity proposed to explain 4.8% of critical cases. However, previous cohorts have been clinically heterogeneous and were not subjected to thorough genetic and immunological analyses. We therefore aimed to systematically investigate the prevalence of rare genetic variants causing inborn errors of immunity (IEI) and functionally interrogate the type I IFN pathway in young adults that suffered from critical COVID-19 yet lacked comorbidities. We selected and clinically characterized a cohort of 38 previously healthy individuals under 50 years of age who were treated in intensive care units due to critical COVID-19. Blood samples were collected after convalescence. Two patients had IFN-α autoantibodies. Genome sequencing revealed very rare variants in the type I IFN pathway in 31.6% of the patients, which was similar to controls. Analyses of cryopreserved leukocytes did not indicate any defect in plasmacytoid dendritic cell sensing of TLR7 and TLR9 agonists in patients carrying variants in these pathways. However, lymphocyte STAT phosphorylation and protein upregulation upon IFN-α stimulation revealed three possible cases of impaired type I IFN signaling in carriers of rare variants. Together, our results suggest a strategy of functional screening followed by genome analyses and biochemical validation to uncover undiagnosed causes of critical COVID-19.

CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions

Machine Learning-based scoring and classification of genetic variants aids the assessment of clinical findings and is employed to prioritize variants in diverse genetic studies and analyses. Combined Annotation-Dependent Depletion (CADD) is one of the first methods for the genome-wide prioritization of variants across different molecular functions and has been continuously developed and improved since its original publication. Here, we present our most recent release, CADD v1.7. We explored and integrated new annotation features, among them state-of-the-art protein language model scores (Meta ESM-1v), regulatory variant effect predictions (from sequence-based convolutional neural networks) and sequence conservation scores (Zoonomia). We evaluated the new version on data sets derived from ClinVar, ExAC/gnomAD and 1000 Genomes variants. For coding effects, we tested CADD on 31 Deep Mutational Scanning (DMS) data sets from ProteinGym and, for regulatory effect prediction, we used saturation mutagenesis reporter assay data of promoter and enhancer sequences. The inclusion of new features further improved the overall performance of CADD. As with previous releases, all data sets, genome-wide CADD v1.7 scores, scripts for on-site scoring and an easy-to-use webserver are readily provided via https://cadd.bihealth.org/ or https://cadd.gs.washington.edu/ to the community.

The functional impact of rare variation across the regulatory cascade

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Polygenic risk scores: An overview from bench to bedside for personalised medicine

Since the first polygenic risk score (PRS) in 2007, research in this area has progressed significantly. The increasing number of SNPs that have been identified by large scale GWAS analyses has fuelled the development of a myriad of PRSs for a wide variety of diseases and, more recently, to PRSs that potentially identify differential response to specific drugs. PRSs constitute a composite genomic biomarker and potential applications for PRSs in clinical practice encompass risk prediction and disease screening, early diagnosis, prognostication, and drug stratification to improve efficacy or reduce adverse drug reactions. Nevertheless, to our knowledge, no PRSs have yet been adopted into routine clinical practice. Beyond the technical considerations of PRS development, the major challenges that face PRSs include demonstrating clinical utility and circumnavigating the implementation of novel genomic technologies at scale into stretched healthcare systems. In this review, we discuss progress in developing disease susceptibility PRSs across multiple medical specialties, development of pharmacogenomic PRSs, and future directions for the field.