The mystery of missing heritability: Genetic interactions create phantom heritability

A broken network of susceptibility genes in the monocytes of Crohn's disease patients

Genome-wide association studies have identified over 200 genetic loci associated with inflammatory bowel disease; however, the mechanism of such a large amount of susceptibility genes remains uncertain. In this study, we integrated bioinformatics analysis and two independent single-cell transcriptome datasets to investigate the expression network of 232 susceptibility genes in Crohn's disease (CD) patients and healthy controls. The study revealed that most of the susceptibility genes are specifically and strictly expressed in the monocytes of the human intestinal tract. The susceptibility genes established a network within the monocytes of health control. The robustness of a gene network may prevent disease onset that is influenced by the genetic and environmental alteration in the expression of susceptibility genes. In contrast, we showed a sparse network in pediatric/adult CD patients, suggesting the broken network contributed to the CD etiology. The network status of susceptibility genes at the single-cell level of monocytes provided novel insight into the etiology.

Epistatic interaction between ERAP2 and HLA modulates HIV-1 adaptation and disease outcome in an Australian population

A strong genetic predictor of outcome following untreated HIV-1 infection is the carriage of specific alleles of human leukocyte antigens (HLAs) that present viral epitopes to T cells. Residual variation in outcome measures may be attributed, in part, to viral adaptation to HLA-restricted T cell responses. Variants of the endoplasmic reticulum aminopeptidases (ERAPs) influence the repertoire of T cell epitopes presented by HLA alleles as they trim pathogen-derived peptide precursors to optimal lengths for antigen presentation, along with other functions unrelated to antigen presentation. We investigated whether ERAP variants influence HLA-associated HIV-1 adaptation with demonstrable effects on overall HIV-1 disease outcome. Utilizing host and viral data of 249 West Australian individuals with HIV-1 subtype B infection, we identified a novel association between two linked ERAP2 single nucleotide polymorphisms (SNPs; rs2248374 and rs2549782) with plasma HIV RNA concentration (viral load) (P adjusted = 0.0024 for both SNPs). Greater HLA-associated HIV-1 adaptation in the HIV-1 Gag gene correlated significantly with higher viral load, lower CD4+ T cell count and proportion; P = 0.0103, P = 0.0061, P = 0.0061, respectively). When considered together, there was a significant interaction between the two ERAP2 SNPs and HLA-associated HIV-1 adaptation on viral load (P = 0.0111). In a comprehensive multivariate model, addition of ERAP2 haplotypes and HLA associated adaptation as an interaction term to known HLA and CCR5 determinants and demographic factors, increased the explanatory variance of population viral load from 17.67% to 45.1% in this dataset. These effects were not replicated in publicly available datasets with comparably sized cohorts, suggesting that any true global epistasis may be dependent on specific HLA-ERAP allelic combinations. Our data raises the possibility that ERAP2 variants may shape peptide repertoires presented to HLA class I-restricted T cells to modulate the degree of viral adaptation within individuals, in turn contributing to disease variability at the population level. Analyses of other populations and experimental studies, ideally with locally derived ERAP genotyping and HLA-specific viral adaptations are needed to elucidate this further.

Genetic and environmental determinants of bone quality: a cross-sectional analysis of the Hungarian Twin Registry

There is abundant evidence that bone mineral content is highly heritable, while the heritability of bone quality (i.e. trabecular bone score [TBS] and quantitative ultrasound index [QUI]) is rarely investigated. We aimed to disentangle the role of genetic, shared and unique environmental factors on TBS and QUI among Hungarian twins. Our study includes 82 twin (48 monozygotic, 33 same-sex dizygotic) pairs from the Hungarian Twin Registry. TBS was determined by DXA, QUI by calcaneal bone ultrasound. To estimate the genetic and environmental effects, we utilized ACE-variance decomposition. For the unadjusted model of TBS, an AE model provided the best fit with > 80% additive genetic heritability. Adjustment for age, sex, BMI and smoking status improved model fit with 48.0% of total variance explained by independent variables. Furthermore, there was a strong dominant genetic effect (73.7%). In contrast, unadjusted and adjusted models for QUI showed an AE structure. Adjustments improved model fit and 25.7% of the total variance was explained by independent variables. Altogether 70-90% of the variance in QUI was related to additive genetic influences. We found a strong genetic heritability of bone quality in unadjusted models. Half of the variance of TBS was explained by age, sex and BMI. Furthermore, the adjusted model suggested that the genetic component of TBS could be dominant or an epistasis could be present. In contrast, independent variables explained only a quarter of the variance of QUI and the additive heritability explained more than half of all the variance.

Single-cell atlas of ABCA7 loss-of-function reveals impaired neuronal respiration via choline-dependent lipid imbalances

Loss-of-function (LoF) variants in the lipid transporter ABCA7 significantly increase the risk of Alzheimer's disease (odds ratio ∼2), yet the pathogenic mechanisms and the neural cell types affected by these variants remain largely unknown. Here, we performed single-nuclear RNA sequencing of 36 human post-mortem samples from the prefrontal cortex of 12 ABCA7 LoF carriers and 24 matched non-carrier control individuals. ABCA7 LoF was associated with gene expression changes in all major cell types. Excitatory neurons, which expressed the highest levels of ABCA7, showed transcriptional changes related to lipid metabolism, mitochondrial function, cell cycle-related pathways, and synaptic signaling. ABCA7 LoF-associated transcriptional changes in neurons were similarly perturbed in carriers of the common AD missense variant ABCA7 p.Ala1527Gly (n = 240 controls, 135 carriers), indicating that findings from our study may extend to large portions of the at-risk population. Consistent with ABCA7's function as a lipid exporter, lipidomic analysis of isogenic iPSC-derived neurons (iNs) revealed profound intracellular triglyceride accumulation in ABCA7 LoF, which was accompanied by a relative decrease in phosphatidylcholine abundance. Metabolomic and biochemical analyses of iNs further indicated that ABCA7 LoF was associated with disrupted mitochondrial bioenergetics that suggested impaired lipid breakdown by uncoupled respiration. Treatment of ABCA7 LoF iNs with CDP-choline (a rate-limiting precursor of phosphatidylcholine synthesis) reduced triglyceride accumulation and restored mitochondrial function, indicating that ABCA7 LoF-induced phosphatidylcholine dyshomeostasis may directly disrupt mitochondrial metabolism of lipids. Treatment with CDP-choline also rescued intracellular amyloid β -42 levels in ABCA7 LoF iNs, further suggesting a link between ABCA7 LoF metabolic disruptions in neurons and AD pathology. This study provides a detailed transcriptomic atlas of ABCA7 LoF in the human brain and mechanistically links ABCA7 LoF-induced lipid perturbations to neuronal energy dyshomeostasis. In line with a growing body of evidence, our study highlights the central role of lipid metabolism in the etiology of Alzheimer's disease.

Plasminogen missense variants and their involvement in cardiovascular and inflammatory disease

Human plasminogen (PLG), the zymogen of the fibrinolytic protease, plasmin, is a polymorphic protein with two widely distributed codominant alleles, PLG/Asp453 and PLG/Asn453. About 15 other missense or non-synonymous single nucleotide polymorphisms (nsSNPs) of PLG show major, yet different, relative abundances in world populations. Although the existence of these relatively abundant allelic variants is generally acknowledged, they are often overlooked or assumed to be non-pathogenic. In fact, at least half of those major variants are classified as having conflicting pathogenicity, and it is unclear if they contribute to different molecular phenotypes. From those, PLG/K19E and PLG/A601T are examples of two relatively abundant PLG variants that have been associated with PLG deficiencies (PD), but their pathogenic mechanisms are unclear. On the other hand, approximately 50 rare and ultra-rare PLG missense variants have been reported to cause PD as homozygous or compound heterozygous variants, often leading to a debilitating disease known as ligneous conjunctivitis. The true abundance of PD-associated nsSNPs is unknown since they can remain undetected in heterozygous carriers. However, PD variants may also contribute to other diseases. Recently, the ultra-rare autosomal dominant PLG/K311E has been found to be causative of hereditary angioedema (HAE) with normal C1 inhibitor. Two other rare pathogenic PLG missense variants, PLG/R153G and PLG/V709E, appear to affect platelet function and lead to HAE, respectively. Herein, PLG missense variants that are abundant and/or clinically relevant due to association with disease are examined along with their world distribution. Proposed molecular mechanisms are discussed when known or can be reasonably assumed.

SwissGenVar: A Platform for Clinical-Grade Interpretation of Genetic Variants to Foster Personalized Healthcare in Switzerland

Large-scale next-generation sequencing (NGS) germline testing is technically feasible today, but variant interpretation represents a major bottleneck in analysis workflows. This includes extensive variant prioritization, annotation, and time-consuming evidence curation. The scale of the interpretation problem is massive, and variants of uncertain significance (VUSs) are a challenge to personalized medicine. This challenge is further compounded by the complexity and heterogeneity of the standards used to describe genetic variants and the associated phenotypes when searching for relevant information to support clinical decision making. To address this, all five Swiss academic institutions for Medical Genetics joined forces with the Swiss Institute of Bioinformatics (SIB) to create SwissGenVar as a user-friendly nationwide repository and sharing platform for genetic variant data generated during routine diagnostic procedures and research sequencing projects. Its aim is to provide a protected environment for expert evidence sharing about individual variants to harmonize and upscale their significance interpretation at the clinical grade according to international standards. To corroborate the clinical assessment, the variant-related data will be combined with consented high-quality clinical information. Broader visibility will be achieved by interfacing with international databases, thus supporting global initiatives in personalized healthcare.

Genetic Association Between Polycystic Ovary Syndrome and the APOA5 rs662799 and PLIN1 rs894160 Metabolic Variants in the Western Saudi Population: A Case-Control Study

Background

Polycystic ovary syndrome (PCOS) is a common endocrinological condition affecting women of reproductive age, associated with insulin resistance and obesity. PCOS pathogenesis is complex and multifactorial, involving genetic and environmental factors.

Objectives

This study aimed to determine and compare genotype and allele frequencies of single nucleotide polymorphisms (SNPs) in the apolipoprotein A5 (APOA5; rs662799) and perilipin 1 (PLIN1; rs894160, rs1052700 and rs6496589) genes in Western Saudi women to investigate their association with PCOS and its clinical characteristics.

Design and methods

This was a case-control study conducted on women with (n = 104) and without (n = 87) PCOS. The SNPs were genotyped using TaqMan genotyping assays.

Results

Significant and direct associations were detected between PCOS susceptibility and APOA5 SNP rs662799 and PLIN1 SNP rs894160 (P < .001). For APOA5 SNP rs662799, women with the A allele were more likely to have PCOS (relative risk [RR] = 1.348, odds ratio [OR] = 2.313, P < .001) and hypertriglyceridaemia (OR = 17.0, P = .5) than women with the G allele. For PLIN1 SNP rs894160, women with the T allele were more likely to have PCOS than women with the C allele (RR = 8.043, OR = 7.427, P < .001). For PLIN1 SNP rs1052700, women with the TT genotype were more likely to have hyperandrogenism (OR = 29.75, P = .02) and an irregular period (OR = 0.07, P = .040) than women with the AT genotype.

Conclusion

We identified novel alleles and genotypes contributing to the genetic risk of PCOS in the Western Saudi population.

Genome-wide association studies (GWAS) and post-GWAS analyses of impulsivity: A systematic review

Impulsivity is related to a host of mental and behavioral problems. It is a complex construct with many different manifestations, most of which are heritable. The genetic compositions of these impulsivity manifestations, however, remain unclear. A number of genome-wide association studies (GWAS) and post-GWAS analyses have tried to address this issue. We conducted a systematic review of all GWAS and post-GWAS analyses of impulsivity published up to December 2023. Available data suggest that single nucleotide polymorphisms (SNPs) in more than a dozen of genes (e.g., CADM2, CTNNA2, GPM6B) are associated with different measures of impulsivity at genome-wide significant levels. Post-GWAS analyses further show that different measures of impulsivity are subject to different degrees of genetic influence, share few genetic variants, and have divergent genetic overlap with basic personality traits such as extroversion and neuroticism, cognitive ability, psychiatric disorders, substance use, and obesity. These findings shed light on controversies in the conceptualization and measurement of impulsivity, while providing new insights on the underlying mechanisms that yoke impulsivity to psychopathology.

Computing linkage disequilibrium aware genome embeddings using autoencoders

MOTIVATION

The completion of the genome has paved the way for genome-wide association studies (GWAS), which explained certain proportions of heritability. GWAS are not optimally suited to detect non-linear effects in disease risk, possibly hidden in non-additive interactions (epistasis). Alternative methods for epistasis detection using, e.g. deep neural networks (DNNs) are currently under active development. However, DNNs are constrained by finite computational resources, which can be rapidly depleted due to increasing complexity with the sheer size of the genome. Besides, the curse of dimensionality complicates the task of capturing meaningful genetic patterns for DNNs; therefore necessitates dimensionality reduction.

RESULTS

We propose a method to compress single nucleotide polymorphism (SNP) data, while leveraging the linkage disequilibrium (LD) structure and preserving potential epistasis. This method involves clustering correlated SNPs into haplotype blocks and training per-block autoencoders to learn a compressed representation of the block's genetic content. We provide an adjustable autoencoder design to accommodate diverse blocks and bypass extensive hyperparameter tuning. We applied this method to genotyping data from Project MinE, and achieved 99% average test reconstruction accuracy-i.e. minimal information loss-while compressing the input to nearly 10% of the original size. We demonstrate that haplotype-block based autoencoders outperform linear Principal Component Analysis (PCA) by approximately 3% chromosome-wide accuracy of reconstructed variants. To the extent of our knowledge, our approach is the first to simultaneously leverage haplotype structure and DNNs for dimensionality reduction of genetic data.

AVAILABILITY AND IMPLEMENTATION

Data are available for academic use through Project MinE at https://www.projectmine.com/research/data-sharing/, contingent upon terms and requirements specified by the source studies. Code is available at https://github.com/gizem-tas/haploblock-autoencoders.

NRN1 epistasis with BDNF and CACNA1C: mediation effects on symptom severity through neuroanatomical changes in schizophrenia

The expression of Neuritin-1 (NRN1), a neurotrophic factor crucial for neurodevelopment and synaptic plasticity, is enhanced by the Brain Derived Neurotrophic Factor (BDNF). Although the receptor of NRN1 remains unclear, it is suggested that NRN1's activation of the insulin receptor (IR) pathway promotes the transcription of the calcium voltage-gated channel subunit alpha1 C (CACNA1C). These three genes have been independently associated with schizophrenia (SZ) risk, symptomatology, and brain differences. However, research on how they synergistically modulate these phenotypes is scarce. We aimed to study whether the genetic epistasis between these genes affects the risk and clinical presentation of the disorder via its effect on brain structure. First, we tested the epistatic effect of NRN1 and BDNF or CACNA1C on (i) the risk for SZ, (ii) clinical symptoms severity and functionality (onset, PANSS, CGI and GAF), and (iii) brain cortical structure (thickness, surface area and volume measures estimated using FreeSurfer) in a sample of 86 SZ patients and 89 healthy subjects. Second, we explored whether those brain clusters influenced by epistatic effects mediate the clinical profiles. Although we did not find a direct epistatic impact on the risk, our data unveiled significant effects on the disorder's clinical presentation. Specifically, the NRN1-rs10484320 x BDNF-rs6265 interplay influenced PANSS general psychopathology, and the NRN1-rs4960155 x CACNA1C-rs1006737 interaction affected GAF scores. Moreover, several interactions between NRN1 SNPs and BDNF-rs6265 significantly influenced the surface area and cortical volume of the frontal, parietal, and temporal brain regions within patients. The NRN1-rs10484320 x BDNF-rs6265 epistasis in the left lateral orbitofrontal cortex fully mediated the effect on PANSS general psychopathology. Our study not only adds clinical significance to the well-described molecular relationship between NRN1 and BDNF but also underscores the utility of deconstructing SZ into biologically validated brain-imaging markers to explore their mediation role in the path from genetics to complex clinical manifestation.

Environmental and Microbial Factors in Inflammatory Bowel Disease Model Establishment: A Review Partly through Mendelian Randomization

Inflammatory bowel disease (IBD) is a complex condition resulting from environmental, microbial, immunologic, and genetic factors. With the advancement of Mendelian randomization research in IBD, we have gained new insights into the relationship between these factors and IBD. Many animal models of IBD have been developed using different methods, but few studies have attempted to model IBD by combining environmental factors and microbial factors. In this review, we examine how environmental factors and microbial factors affect the development and progression of IBD, and how they interact with each other and with the intestinal microbiota. We also summarize the current methods for creating animal models of IBD and compare their advantages and disadvantages. Based on the latest findings from Mendelian randomization studies on the role of environmental factors in IBD, we discuss which environmental and microbial factors could be used to construct a more realistic and reliable IBD experimental model. We propose that animal models of IBD should consider both environmental and microbial factors to better mimic human IBD pathogenesis and to reveal the underlying mechanisms of IBD at the immune and genetic levels. We highlight the importance of environmental and microbial factors in IBD pathogenesis and offer new perspectives and suggestions for improving experimental animal modeling. Our goal is to create a model that closely resembles the clinical picture of IBD.

Rare-variant association study unveils the Achilles' heel for HCC

In this issue of Cell Genomics, Wang, Liu, Zuo, Wang, et al.1 investigate rare variants in hepatocellular carcinoma (HCC) by performing the first rare-variant association study (RVAS) in a Chinese population cohort. It uncovers BRCAness phenotypes associated with the NRDE2-p.N377I variant, suggesting PARP inhibitors as a promising therapeutic approach for certain HCC patients.

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.

BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

Evaluation of non-additive genetic effects on carcass and meat quality traits in Korean Hanwoo cattle using genomic models

The traditional genetic evaluation methods generally consider additive genetic effects only and often ignore non-additive (dominance and epistasis) effects that may have contributed to genetic variation of complex traits of livestock species. The available dense single nucleotide polymorphisms (SNPs) panels offer to investigate the potential benefits of including non-additive genetic effects in the genomic evaluation models. Data from 16 971 genotyped (Illumina Bovine 50 K SNP chip) Korean Hanwoo cattle were used to estimate genetic variance components and prediction accuracy of genomic breeding values (GEBVs) for four carcass and meat quality traits: carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT) and marbling score (MS). Five different genetic models were evaluated through including additive, dominance and epistatic interactions (additive by additive, A × A; additive by dominance, A × D and dominance by dominance, D × D) successively in the models. The estimates of additive genetic variances and narrow sense heritabilities (ha2) were found similar across the evaluated models and traits except when additive interaction (A × A) was included. The dominance variance estimates relative to phenotypic variance ranged from 1.7-3.4% for CWT and MS traits, whereas, they were close to zero for EMA and BFT traits. The magnitude of A × A epistatic heritability (haa2) ranged between 14.8 and 27.7% in all traits. However, heritability estimates for A × D and D × D epistatic interactions (had2 and hdd2) were quite low compared to haa2 and were contributed only 0.0-9.7% of the total phenotypic variation. In general, broad sense heritability (hG2) estimates were almost twice (ranging between 0.54 and 0.68) the ha2 for all of the investigated traits. The inclusion of dominance effects did not improve the prediction accuracy of GEBV but improved 2.0-3.0% when epistatic effects were included in the model. More importantly, rank correlation revealed that partitioning of variance components considering dominance and epistatic effects in the model would enable to re-rank of top animals with better prediction of GEBV. The present result suggests that dominance and epistatic effects could be included in the genomic evaluation model for better estimates of variance components and more accurate prediction of GEBV for carcass and meat quality traits in Korean Hanwoo cattle.

Identifying latent genetic interactions in genome-wide association studies using multiple traits

The "missing" heritability of complex traits may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. We propose a new kernel-based method called Latent Interaction Testing (LIT) to screen for genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Using simulated data, we demonstrate that LIT increases power to detect latent genetic interactions compared to univariate methods. We then apply LIT to obesity-related traits in the UK Biobank and detect variants with interactive effects near known obesity-related genes (URL: https://CRAN.R-project.org/package=lit ).

An approach to identify gene-environment interactions and reveal new biological insight in complex traits

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the Gene-Lifestyle Interactions Working Group within the Cohorts for Heart and Aging Research in Genetic Epidemiology Consortium has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identify and confirm 5 loci (6 independent signals) interacted with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrate that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated heritability is significant (P < 0.02) for low-density lipoprotein cholesterol and triglycerides in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Distinct explanations underlie gene-environment interactions in the UK Biobank

The role of gene-environment (GxE) interaction in disease and complex trait architectures is widely hypothesized, but currently unknown. Here, we apply three statistical approaches to quantify and distinguish three different types of GxE interaction for a given trait and E variable. First, we detect locus-specific GxE interaction by testing for genetic correlation r g < 1 across E bins. Second, we detect genome-wide effects of the E variable on genetic variance by leveraging polygenic risk scores (PRS) to test for significant PRSxE in a regression of phenotypes on PRS, E, and PRSxE, together with differences in SNP-heritability across E bins. Third, we detect genome-wide proportional amplification of genetic and environmental effects as a function of the E variable by testing for significant PRSxE with no differences in SNP-heritability across E bins. Simulations show that these approaches achieve high sensitivity and specificity in distinguishing these three GxE scenarios. We applied our framework to 33 UK Biobank traits (25 quantitative traits and 8 diseases; average N = 325 K ) and 10 E variables spanning lifestyle, diet, and other environmental exposures. First, we identified 19 trait-E pairs with r g significantly < 1 (FDR<5%) (average r g = 0.95 ); for example, white blood cell count had r g = 0.95 (s.e. 0.01) between smokers and non-smokers. Second, we identified 28 trait-E pairs with significant PRSxE and significant SNP-heritability differences across E bins; for example, BMI had a significant PRSxE for physical activity (P=4.6e-5) with 5% larger SNP-heritability in the largest versus smallest quintiles of physical activity (P=7e-4). Third, we identified 15 trait-E pairs with significant PRSxE with no SNP-heritability differences across E bins; for example, waist-hip ratio adjusted for BMI had a significant PRSxE effect for time spent watching television (P=5e-3) with no SNP-heritability differences. Across the three scenarios, 8 of the trait-E pairs involved disease traits, whose interpretation is complicated by scale effects. Analyses using biological sex as the E variable produced additional significant findings in each of the three scenarios. Overall, we infer a significant contribution of GxE and GxSex effects to complex trait and disease variance.

An expression-directed linear mixed model discovering low-effect genetic variants

Detecting genetic variants with low-effect sizes using a moderate sample size is difficult, hindering downstream efforts to learn pathology and estimating heritability. In this work, by utilizing informative weights learned from training genetically predicted gene expression models, we formed an alternative approach to estimate the polygenic term in a linear mixed model. Our linear mixed model estimates the genetic background by incorporating their relevance to gene expression. Our protocol, expression-directed linear mixed model, enables the discovery of subtle signals of low-effect variants using moderate sample size. By applying expression-directed linear mixed model to cohorts of around 5,000 individuals with either binary (WTCCC) or quantitative (NFBC1966) traits, we demonstrated its power gain at the low-effect end of the genetic etiology spectrum. In aggregate, the additional low-effect variants detected by expression-directed linear mixed model substantially improved estimation of missing heritability. Expression-directed linear mixed model moves precision medicine forward by accurately detecting the contribution of low-effect genetic variants to human diseases.

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Darier disease is associated with neurodegenerative disorders and epilepsy

Darier disease (DD) is a rare monogenetic skin disorder with limited data on its potential association with neurological disorders. This study aimed to investigate the association between DD and neurological disorders, specifically Parkinson's disease, dementias, and epilepsy. Using Swedish national registers in a period spanning between 1977 and 2013, 935 individuals with DD were compared with up to 100 comparison individuals each, randomly selected from the general population based on birth year, sex, and county of residence at the time of the first diagnosis of DD. Individuals with DD had increased risks of being diagnosed with Parkinson's disease (RR 2.1, CI 1.1; 4.4), vascular dementia (RR 2.1, CI 1.0; 4.2), and epilepsy, (RR 2.5, CI 1.8; 3.5). No association of DD with other dementias were detected. This study demonstrates a new association between DD and neurodegenerative disorders and epilepsy, underlining the need for increased awareness, interdisciplinary collaboration, and further research to understand the underlying mechanisms. Early identification and management of neurological complications in DD patients could improve treatment strategies and patient outcomes. The findings also highlight the role of SERCA2 in the pathophysiology of neurological disorders, offering new targets for future research and potentials for novel treatments.

Next-Gen GWAS: full 2D epistatic interaction maps retrieve part of missing heritability and improve phenotypic prediction

The problem of missing heritability requires the consideration of genetic interactions among different loci, called epistasis. Current GWAS statistical models require years to assess the entire combinatorial epistatic space for a single phenotype. We propose Next-Gen GWAS (NGG) that evaluates over 60 billion single nucleotide polymorphism combinatorial first-order interactions within hours. We apply NGG to Arabidopsis thaliana providing two-dimensional epistatic maps at gene resolution. We demonstrate on several phenotypes that a large proportion of the missing heritability can be retrieved, that it indeed lies in epistatic interactions, and that it can be used to improve phenotype prediction.

Integration of whole-exome sequencing and structural neuroimaging analysis in major depressive disorder: a joint study

Major depressive disorder (MDD) is a common mental illness worldwide and is triggered by an intricate interplay between environmental and genetic factors. Although there are several studies on common variants in MDD, studies on rare variants are relatively limited. In addition, few studies have examined the genetic contributions to neurostructural alterations in MDD using whole-exome sequencing (WES). We performed WES in 367 patients with MDD and 161 healthy controls (HCs) to detect germline and copy number variations in the Korean population. Gene-based rare variants were analyzed to investigate the association between the genes and individuals, followed by neuroimaging-genetic analysis to explore the neural mechanisms underlying the genetic impact in 234 patients with MDD and 135 HCs using diffusion tensor imaging data. We identified 40 MDD-related genes and observed 95 recurrent regions of copy number variations. We also discovered a novel gene, FRMPD3, carrying rare variants that influence MDD. In addition, the single nucleotide polymorphism rs771995197 in the MUC6 gene was significantly associated with the integrity of widespread white matter tracts. Moreover, we identified 918 rare exonic missense variants in genes associated with MDD susceptibility. We postulate that rare variants of FRMPD3 may contribute significantly to MDD, with a mild penetration effect.

Temporal changes in brain morphology related to inflammation and schizophrenia: an omnigenic Mendelian randomization study

BACKGROUND

Over the past several decades, more research focuses have been made on the inflammation/immune hypothesis of schizophrenia. Building upon synaptic plasticity hypothesis, inflammation may contribute the underlying pathophysiology of schizophrenia. Yet, pinpointing the specific inflammatory agents responsible for schizophrenia remains a complex challenge, mainly due to medication and metabolic status. Multiple lines of evidence point to a wide-spread genetic association across genome underlying the phenotypic variations of schizophrenia.

METHOD

We collected the latest genome-wide association analysis (GWAS) summary data of schizophrenia, cytokines, and longitudinal change of brain. We utilized the omnigenic model which takes into account all genomic SNPs included in the GWAS of trait, instead of traditional Mendelian randomization (MR) methods. We conducted two round MR to investigate the inflammatory triggers of schizophrenia and the resulting longitudinal changes in the brain.

RESULTS

We identified seven inflammation markers linked to schizophrenia onset, which all passed the Bonferroni correction for multiple comparisons (bNGF, GROA(CXCL1), IL-8, M-CSF, MCP-3 (CCL7), TNF-β, CRP). Moreover, CRP were found to significantly influence the linear rate of brain morphology changes, predominantly in the white matter of the cerebrum and cerebellum.

CONCLUSION

With an omnigenic approach, our study sheds light on the immune pathology of schizophrenia. Although these findings need confirmation from future studies employing different methodologies, our work provides substantial evidence that pervasive, low-level neuroinflammation may play a pivotal role in schizophrenia, potentially leading to notable longitudinal changes in brain morphology.

Uncovering genetic associations in the human diseasome using an endophenotype-augmented disease network

MOTIVATION

Many diseases, particularly cardiometabolic disorders, exhibit complex multimorbidities with one another. An intuitive way to model the connections between phenotypes is with a disease-disease network (DDN), where nodes represent diseases and edges represent associations, such as shared single-nucleotide polymorphisms (SNPs), between pairs of diseases. To gain further genetic understanding of molecular contributors to disease associations, we propose a novel version of the shared-SNP DDN (ssDDN), denoted as ssDDN+, which includes connections between diseases derived from genetic correlations with intermediate endophenotypes. We hypothesize that a ssDDN+ can provide complementary information to the disease connections in a ssDDN, yielding insight into the role of clinical laboratory measurements in disease interactions.

RESULTS

Using PheWAS summary statistics from the UK Biobank, we constructed a ssDDN+ revealing hundreds of genetic correlations between diseases and quantitative traits. Our augmented network uncovers genetic associations across different disease categories, connects relevant cardiometabolic diseases, and highlights specific biomarkers that are associated with cross-phenotype associations. Out of the 31 clinical measurements under consideration, HDL-C connects the greatest number of diseases and is strongly associated with both type 2 diabetes and heart failure. Triglycerides, another blood lipid with known genetic causes in non-mendelian diseases, also adds a substantial number of edges to the ssDDN. This work demonstrates how association with clinical biomarkers can better explain the shared genetics between cardiometabolic disorders. Our study can facilitate future network-based investigations of cross-phenotype associations involving pleiotropy and genetic heterogeneity, potentially uncovering sources of missing heritability in multimorbidities.

AVAILABILITY AND IMPLEMENTATION

The generated ssDDN+ can be explored at https://hdpm.biomedinfolab.com/ddn/biomarkerDDN.

The intergenerational transmission of risk and trust attitudes: Replicating and extending "Dohmen, Falk, Huffman and Sunde 2012" using genetically informed twin data

This replication revisits an influential contribution on the intergenerational transmission of risk and trust attitudes, which, based on data from the German Socioeconomic Panel (GSOEP), reveals a positive correlation between parents' and children's attitudes. The authors of the original study argue that socialization in the family is important in the transmission process. The replication is motivated by mounting evidence indicating that within-family transmission has a considerable genetic component, which calls into question socialization as the main transmission pathway. To consider genetic transmission in addition to social transmission, the replication relies on the German twin family panel TwinLife. The findings reveal that, first, most of the variation in children's risk and social trust attitudes is attributable to differences in the non-shared environment, followed by genetic differences, whereas differences in the shared family environment - the main candidate for social transmission - do not matter. Second, correlations between parents' and children's attitudes essentially involve genetic similarity. Third, family conditions do not moderate these relationships. Thus, the findings do not support the socialization assumption.

Extreme trait GWAS (Et-GWAS): Unraveling rare variants in the 3,000 rice genome

Identifying high-impact, rare genetic variants associated with specific traits is crucial for crop improvement. The 3,010 rice genome (3K RG) dataset offers a valuable resource for discovering genomic regions with potential applications in crop breeding. We used Extreme Trait GWAS (Et-GWAS), employing bulk pooling and allele frequency measurement to efficiently extract rare variants from the 3K RG. This innovative approach facilitates the detection of associations between genetic variants and target traits, concentrating and quantifying rare alleles. In our study, on grain yield under drought stress, Et-GWAS successfully identified five key genes (OsPP2C11, OsK5.2, OsIRO2, OsPEX1, and OsPWA1) known for enhancing yield under drought. In addition, we examined the overlap of our results with previously reported qDTY-QTLs and observed that OsUCH1 and OsUCH2 genes were located within qDTY2.2 We compared Et-GWAS with conventional GWAS, finding it effectively capturing most candidate genes associated with the target trait. Validation with resistant starch showed similar results. To enhance user-friendliness, we developed a GUI for Et-GWAS; https://et-gwas.shinyapps.io/Et-GWAS/.

Genome-wide epistasis analysis reveals gene-gene interaction network on an intermediate endophenotype P-tau/Aβ42 ratio in ADNI cohort

Alzheimer's disease (AD) is a progressive neurodegenerative disorder and the most common cause of dementia in the elderly worldwide. The exact etiology of AD, particularly its genetic mechanisms, remains incompletely understood. Traditional genome-wide association studies (GWAS), which primarily focus on single-nucleotide polymorphisms (SNPs) with main effects, provide limited explanations for the "missing heritability" of AD, while there is growing evidence supporting the important role of epistasis. In this study, we performed a genome-wide SNP-SNP interaction detection using a linear regression model and employed multiple GPUs for parallel computing, significantly enhancing the speed of whole-genome analysis. The cerebrospinal fluid (CSF) phosphorylated tau (P-tau)/amyloid-[Formula: see text] (A[Formula: see text]) ratio was used as a quantitative trait (QT) to enhance statistical power. Age, gender, and clinical diagnosis were included as covariates to control for potential non-genetic factors influencing AD. We identified 961 pairs of statistically significant SNP-SNP interactions, explaining a high-level variance of P-tau/A[Formula: see text] level, all of which exhibited marginal main effects. Additionally, we replicated 432 previously reported AD-related genes and found 11 gene-gene interaction pairs overlapping with the protein-protein interaction (PPI) network. Our findings may contribute to partially explain the "missing heritability" of AD. The identified subnetwork may be associated with synaptic dysfunction, Wnt signaling pathway, oligodendrocytes, inflammation, hippocampus, and neuronal cells.

Gut microbial diversity moderates polygenic risk of schizophrenia

Background

Schizophrenia (SCZ) is a heritable disorder with a polygenic architecture, and the gut microbiota seems to be involved in its development and outcome. In this study, we investigate the interplay between genetic risk and gut microbial markers.

Methods

We included 159 first-episode, drug-naïve SCZ patients and 86 healthy controls. The microbial composition of feces was characterized using the 16S rRNA sequencing platform, and five microbial α-diversity indices were estimated [Shannon, Simpson, Chao1, the Abundance-based Eoverage Estimator (ACE), and a phylogenetic diversity-based estimate (PD)]. Polygenic risk scores (PRS) for SCZ were constructed using data from large-scale genome-wide association studies. Effects of microbial α-diversity, microbial abundance, and PRS on SCZ were evaluated via generalized linear models.

Results

We confirmed that PRS was associated with SCZ (OR = 2.08, p = 1.22×10-5) and that scores on the Shannon (OR = 0.29, p = 1.15×10-8) and Simpson (OR = 0.29, p = 1.25×10-8) indices were inversely associated with SCZ risk. We found significant interactions (p < 0.05) between PRS and α-diversity indices (Shannon, Simpson, and PD), with the effects of PRS being larger in those exhibiting higher diversity compared to those with lower diversity. Moreover, the PRS effects were larger in individuals with a high abundance of the genera Romboutsia, Streptococcus, and Anaerostipes than in those with low abundance (p < 0.05). All three of these genera showed protective effects against SCZ.

Conclusion

The current findings suggest an interplay between the gut microbiota and polygenic risk of SCZ that warrants replication in independent samples. Experimental studies are needed to determine the underpinning mechanisms.

SumVg: Total Heritability Explained by All Variants in Genome-Wide Association Studies Based on Summary Statistics with Standard Error Estimates

Genome-wide association studies (GWAS) are commonly employed to study the genetic basis of complex traits/diseases, and a key question is how much heritability could be explained by all single nucleotide polymorphisms (SNPs) in GWAS. One widely used approach that relies on summary statistics only is linkage disequilibrium score regression (LDSC); however, this approach requires certain assumptions about the effects of SNPs (e.g., all SNPs contribute to heritability and each SNP contributes equal variance). More flexible modeling methods may be useful. We previously developed an approach recovering the "true" effect sizes from a set of observed z-statistics with an empirical Bayes approach, using only summary statistics. However, methods for standard error (SE) estimation are not available yet, limiting the interpretation of our results and the applicability of the approach. In this study, we developed several resampling-based approaches to estimate the SE of SNP-based heritability, including two jackknife and three parametric bootstrap methods. The resampling procedures are performed at the SNP level as it is most common to estimate heritability from GWAS summary statistics alone. Simulations showed that the delete-d-jackknife and parametric bootstrap approaches provide good estimates of the SE. In particular, the parametric bootstrap approaches yield the lowest root-mean-squared-error (RMSE) of the true SE. We also explored various methods for constructing confidence intervals (CIs). In addition, we applied our method to estimate the SNP-based heritability of 12 immune-related traits (levels of cytokines and growth factors) to shed light on their genetic architecture. We also implemented the methods to compute the sum of heritability explained and the corresponding SE in an R package SumVg. In conclusion, SumVg may provide a useful alternative tool for calculating SNP heritability and estimating SE/CI, which does not rely on distributional assumptions of SNP effects.

Trans-regulatory variant network contributes to missing heritability

In a recent Cell Genomics article, Tsouris et al.1 analyze the transcriptomes of a large diallel panel of hybrids from Saccharomyces cerevisiae natural isolates to study cis- and trans-regulatory changes underlying gene expression variation. Vanessa Pereira and Elena Kuzmin discuss the authors' findings and the wider context in missing heritability research in this preview.

Estimating the heritability of SARS-CoV-2 susceptibility and COVID-19 severity

SARS-CoV-2 has infected over 340 million people, prompting therapeutic research. While genetic studies can highlight potential drug targets, understanding the heritability of SARS-CoV-2 susceptibility and COVID-19 severity can contextualize their results. To date, loci from meta-analyses explain 1.2% and 5.8% of variation in susceptibility and severity respectively. Here we estimate the importance of shared environment and additive genetic variation to SARS-CoV-2 susceptibility and COVID-19 severity using pedigree data, PCR results, and hospitalization information. The relative importance of genetics and shared environment for susceptibility shifted during the study, with heritability ranging from 33% (95% CI: 20%-46%) to 70% (95% CI: 63%-74%). Heritability was greater for days hospitalized with COVID-19 (41%, 95% CI: 33%-57%) compared to shared environment (33%, 95% CI: 24%-38%). While our estimates suggest these genetic architectures are not fully understood, the shift in susceptibility estimates highlights the challenge of estimation during a pandemic, given environmental fluctuations and vaccine introduction.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value ≤ 1 × 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

Genetic signatures of suicide attempt behavior: insights and applications

INTRODUCTION

Every year about 800,000 complete suicide events occur. The identification of biologic markers to identify subjects at risk would be helpful in targeting specific support treatments.

AREA COVERED

A narrative review defines the meta-analytic level of current evidence about the biologic markers of suicide behavior (SB). The meta-analytic evidence gathered so far indicates that the hypothesis-driven research largely failed to identify the biologic markers of suicide. The most consistent and replicated result was reported for: 1) 5-HTR2A T102C, associated with SB in patients with schizophrenia (OR = 1.73 (1.11-2.69)) and 2) BDNF Val66Met (rs6265), with the Met-Val + Val-Val carriers found to be at risk for suicide in the Caucasian population (OR: 1.96 (1.58-2.43)), while Val-Val vs. Val-Met + Met carriers found to be at risk for suicide in the Asian populations (OR: 1.36 (1.04-1.78)). GWAS-based meta-analyses indicate some positive replicated findings regarding the DRD2, Neuroligin gene, estrogen-related genes, and genes involved in gene expression.

EXPERT OPINION

Most consistent results were obtained when analyzing sub-samples of patients. Some promising results come from the implementation of the polygenic risk score. There is no current consensus about an implementable biomarker for SB.

Principles and methods for transferring polygenic risk scores across global populations

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Network medicine-based epistasis detection in complex diseases: ready for quantum computing

Most heritable diseases are polygenic. To comprehend the underlying genetic architecture, it is crucial to discover the clinically relevant epistatic interactions (EIs) between genomic single nucleotide polymorphisms (SNPs)1-3. Existing statistical computational methods for EI detection are mostly limited to pairs of SNPs due to the combinatorial explosion of higher-order EIs. With NeEDL (network-based epistasis detection via local search), we leverage network medicine to inform the selection of EIs that are an order of magnitude more statistically significant compared to existing tools and consist, on average, of five SNPs. We further show that this computationally demanding task can be substantially accelerated once quantum computing hardware becomes available. We apply NeEDL to eight different diseases and discover genes (affected by EIs of SNPs) that are partly known to affect the disease, additionally, these results are reproducible across independent cohorts. EIs for these eight diseases can be interactively explored in the Epistasis Disease Atlas (https://epistasis-disease-atlas.com). In summary, NeEDL is the first application that demonstrates the potential of seamlessly integrated quantum computing techniques to accelerate biomedical research. Our network medicine approach detects higher-order EIs with unprecedented statistical and biological evidence, yielding unique insights into polygenic diseases and providing a basis for the development of improved risk scores and combination therapies.

Factorizing polygenic epistasis improves prediction and uncovers biological pathways in complex traits

Epistasis is central in many domains of biology, but it has not yet been proven useful for understanding the etiology of complex traits. This is partly because complex-trait epistasis involves polygenic interactions that are poorly captured in current models. To address this gap, we developed a model called Epistasis Factor Analysis (EFA). EFA assumes that polygenic epistasis can be factorized into interactions between a few epistasis factors (EFs), which represent latent polygenic components of the observed complex trait. The statistical goals of EFA are to improve polygenic prediction and to increase power to detect epistasis, while the biological goal is to unravel genetic effects into more-homogeneous units. We mathematically characterize EFA and use simulations to show that EFA outperforms current epistasis models when its assumptions approximately hold. Applied to predicting yeast growth rates, EFA outperforms the additive model for several traits with large epistasis heritability and uniformly outperforms the standard epistasis model. We replicate these prediction improvements in a second dataset. We then apply EFA to four previously characterized traits in the UK Biobank and find statistically significant epistasis in all four, including two that are robust to scale transformation. Moreover, we find that the inferred EFs partly recover pre-defined biological pathways for two of the traits. Our results demonstrate that more realistic models can identify biologically and statistically meaningful epistasis in complex traits, indicating that epistasis has potential for precision medicine and characterizing the biology underlying GWAS results.

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

Shared Genetic Risk in the Association of Screen Time With Psychiatric Problems in Children

Importance

Children's exposure to screen time has been associated with poor mental health outcomes, yet the role of genetic factors remains largely unknown.

Objective

To assess the extent of genetic confounding in the associations between screen time and attention problems or internalizing problems in preadolescent children.

Design, Setting, and Participants

This cohort study analyzed data obtained between 2016 and 2019 from the Adolescent Brain Cognitive Development Study at 21 sites in the US. The sample included children aged 9 to 11 years of genetically assigned European ancestry with self-reported screen time. Data were analyzed between November 2021 and September 2023.

Exposure

Child-reported daily screen time (in hours) was ascertained from questionnaires completed by the children at baseline.

Main Outcomes and Measures

Child psychiatric problems, specifically attention and internalizing problems, were measured with the parent-completed Achenbach Child Behavior Checklist at the 1-year follow-up. Genetic sensitivity analyses model (Gsens) was used, which incorporated polygenic risk scores (PRSs) of both exposure and outcomes as well as either single-nucleotide variant (SNV; formerly single-nucleotide polymorphism)-based heritability or twin-based heritability to estimate genetic confounding.

Results

The 4262 children in the sample included 2269 males (53.2%) with a mean (SD) age of 9.9 (0.6) years. Child screen time was associated with attention problems (β = 0.10 SD; 95% CI, 0.07-0.13 SD) and internalizing problems (β = 0.03 SD; 95% CI, 0.003-0.06 SD). The television time PRS was associated with child screen time (β = 0.18 SD; 95% CI, 0.14-0.23 SD), the attention-deficit/hyperactivity disorder PRS was associated with attention problems (β = 0.13 SD; 95% CI, 0.10-0.16 SD), and the depression PRS was associated with internalizing problems (β = 0.10 SD; 95% CI, 0.07-0.13 SD). These PRSs were associated with cross-traits, suggesting genetic confounding. Estimates using PRSs and SNV-based heritability showed that genetic confounding accounted for most of the association between child screen time and attention problems and for 42.7% of the association between child screen time and internalizing problems. When PRSs and twin-based heritability estimates were used, genetic confounding fully explained both associations.

Conclusions and Relevance

Results of this study suggest that genetic confounding may explain a substantial part of the associations between child screen time and psychiatric problems. Genetic confounding should be considered in sociobehavioral studies of modifiable factors for youth mental health.

Bayesian multivariate genetic analysis improves translational insights

While lipid traits are known essential mediators of cardiovascular disease, few approaches have taken advantage of their shared genetic effects. We apply a Bayesian multivariate size estimator, mash, to GWAS of four lipid traits in the Million Veterans Program (MVP) and provide posterior mean and local false sign rates for all effects. These estimates borrow information across traits to improve effect size accuracy. We show that controlling local false sign rates accurately and powerfully identifies replicable genetic associations and that multivariate control furthers the ability to explain complex diseases. Our application yields high concordance between independent datasets, more accurately prioritizes causal genes, and significantly improves polygenic prediction beyond state-of-the-art methods by up to 59% for lipid traits. The use of Bayesian multivariate genetic shrinkage has yet to be applied to human quantitative trait GWAS results, and we present a staged approach to prediction on a polygenic scale.

Large sample size and nonlinear sparse models outline epistatic effects in inflammatory bowel disease

BACKGROUND

Despite clear evidence of nonlinear interactions in the molecular architecture of polygenic diseases, linear models have so far appeared optimal in genotype-to-phenotype modeling. A key bottleneck for such modeling is that genetic data intrinsically suffers from underdetermination ([Formula: see text]). Millions of variants are present in each individual while the collection of large, homogeneous cohorts is hindered by phenotype incidence, sequencing cost, and batch effects.

RESULTS

We demonstrate that when we provide enough training data and control the complexity of nonlinear models, a neural network outperforms additive approaches in whole exome sequencing-based inflammatory bowel disease case-control prediction. To do so, we propose a biologically meaningful sparsified neural network architecture, providing empirical evidence for positive and negative epistatic effects present in the inflammatory bowel disease pathogenesis.

CONCLUSIONS

In this paper, we show that underdetermination is likely a major driver for the apparent optimality of additive modeling in clinical genetics today.

A new Approach to Identify Gene-Environment Interactions and Reveal New Biological Insight in Complex traits

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the CHARGE Gene-Lifestyle Interactions Working Group has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identified and confirmed 5 loci (6 independent signals) interacting with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrated that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated contribution ranges from 1.76% to 14.05% of SNP heritability of serum lipids in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Trans-ancestry meta-analysis of genome wide association studies of inhibitory control

Deficits in effective executive function, including inhibitory control are associated with risk for a number of psychiatric disorders and significantly impact everyday functioning. These complex traits have been proposed to serve as endophenotypes, however, their genetic architecture is not yet well understood. To identify the common genetic variation associated with inhibitory control in the general population we performed the first trans-ancestry genome wide association study (GWAS) combining data across 8 sites and four ancestries (N = 14,877) using cognitive traits derived from the stop-signal task, namely - go reaction time (GoRT), go reaction time variability (GoRT SD) and stop signal reaction time (SSRT). Although we did not identify genome wide significant associations for any of the three traits, GoRT SD and SSRT demonstrated significant and similar SNP heritability of 8.2%, indicative of an influence of genetic factors. Power analyses demonstrated that the number of common causal variants contributing to the heritability of these phenotypes is relatively high and larger sample sizes are necessary to robustly identify associations. In Europeans, the polygenic risk for ADHD was significantly associated with GoRT SD and the polygenic risk for schizophrenia was associated with GoRT, while in East Asians polygenic risk for schizophrenia was associated with SSRT. These results support the potential of executive function measures as endophenotypes of neuropsychiatric disorders. Together these findings provide the first evidence indicating the influence of common genetic variation in the genetic architecture of inhibitory control quantified using objective behavioural traits derived from the stop-signal task.

Dopaminergic signalling and behavioural alterations by Comt-Dtnbp1 genetic interaction and their clinical relevance

BACKGROUND AND PURPOSE

Cognitive and motor functions are modulated by dopaminergic signalling, which is shaped by several genetic factors. The biological effects of single genetic variants might differ depending on epistatic interactions that can be functionally multi-directional and non-linear.

EXPERIMENTAL APPROACH

We performed behavioural and neurochemical assessments in genetically modified mice and behavioural assessments and genetic screening in human patients with 22q11.2 deletion syndrome (22q11.2DS).

KEY RESULTS

Here, we confirm a genetic interaction between the Comt (catechol-O-methyltransferase, human orthologue: COMT) and Dtnbp1 (dystrobrevin binding protein 1, alias dysbindin, human orthologue: DTNBP1) genes that modulate cortical and striatal dopaminergic signalling in a manner not predictable by the effects of each single gene. In mice, Comt-by-Dtnbp1 concomitant reduction leads to a hypoactive mesocortical and a hyperactive mesostriatal dopamine pathway, associated with specific cognitive abnormalities. Like mice, in subjects with the 22q11.2DS (characterized by COMT hemideletion and dopamine alterations), COMT-by-DTNBP1 concomitant reduction was associated with analogous cognitive disturbances. We then developed an easy and inexpensive colourimetric kit for the genetic screening of common COMT and DTNBP1 functional genetic variants for clinical application.

CONCLUSIONS AND IMPLICATIONS

These findings illustrate an epistatic interaction of two dopamine-related genes and their functional effects, supporting the need to address genetic interaction mechanisms at the base of complex behavioural traits.

Complex effects of sequence variants on lipid levels and coronary artery disease

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

Identifying latent genetic interactions in genome-wide association studies using multiple traits

Genome-wide association studies of complex traits frequently find that SNP-based estimates of heritability are considerably smaller than estimates from classic family-based studies. This 'missing' heritability may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. To circumvent these challenges, we propose a new method to detect genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Our approach, Latent Interaction Testing (LIT), uses the observation that correlated traits with shared latent genetic interactions have trait variance and covariance patterns that differ by genotype. LIT examines the relationship between trait variance/covariance patterns and genotype using a flexible kernel-based framework that is computationally scalable for biobank-sized datasets with a large number of traits. We first use simulated data to demonstrate that LIT substantially increases power to detect latent genetic interactions compared to a trait-by-trait univariate method. We then apply LIT to four obesity-related traits in the UK Biobank and detect genetic variants with interactive effects near known obesity-related genes. Overall, we show that LIT, implemented in the R package lit, uses shared information across traits to improve detection of latent genetic interactions compared to standard approaches.

GWASs and polygenic scores inherit all the old problems of heritability estimates

Polygenic score (PGS) computations assume an additive model of gene action because associations between phenotypes and alleles at different loci are compounded, ignoring interactions between alleles or loci let alone between genotype and environment. Consequently, PGSs are subject to the same objections that invalidated traditional heritability analyses in the 1970s. Thus, PGSs should not be used in the social sciences.

Human genomic data have different statistical properties than the data of randomised controlled trials

Madole & Harden argue that the Mendelian reshuffling of genes and genomes is analogous to randomised controlled trials. We are not convinced by their arguments. First, their recipe for meeting the demands on randomised experiments is inherently inconsistent. Second, disequilibrium across chromosomes conflicts with their assumption of statistical independence. Third, the genome-wide association study (GWAS) method has many pitfalls, including low repeatability.

Evolutionary Action-Machine Learning Model Identifies Candidate Genes Associated With Early-Onset Coronary Artery Disease

Background Coronary artery disease is a primary cause of death around the world, with both genetic and environmental risk factors. Although genome-wide association studies have linked >100 unique loci to its genetic basis, these only explain a fraction of disease heritability. Methods and Results To find additional gene drivers of coronary artery disease, we applied machine learning to quantitative evolutionary information on the impact of coding variants in whole exomes from the Myocardial Infarction Genetics Consortium. Using ensemble-based supervised learning, the Evolutionary Action-Machine Learning framework ranked each gene's ability to classify case and control samples and identified 79 significant associations. These were connected to known risk loci; enriched in cardiovascular processes like lipid metabolism, blood clotting, and inflammation; and enriched for cardiovascular phenotypes in knockout mouse models. Among them, INPP5F and MST1R are examples of potentially novel coronary artery disease risk genes that modulate immune signaling in response to cardiac stress. Conclusions We concluded that machine learning on the functional impact of coding variants, based on a massive amount of evolutionary information, has the power to suggest novel coronary artery disease risk genes for mechanistic and therapeutic discoveries in cardiovascular biology, and should also apply in other complex polygenic diseases.