Integrating predicted transcriptome from multiple tissues improves association detection

Multiomic integration analysis identifies atherogenic metabolites mediating between novel immune genes and cardiovascular risk

BACKGROUND

Understanding genetic-metabolite associations has translational implications for informing cardiovascular risk assessment. Interrogating functional genetic variants enhances our understanding of disease pathogenesis and the development and optimization of targeted interventions.

METHODS

In this study, a total of 187 plasma metabolite levels were profiled in 4974 individuals of European ancestry of the GCAT| Genomes for Life cohort. Results of genetic analyses were meta-analysed with additional datasets, resulting in up to approximately 40,000 European individuals. Results of meta-analyses were integrated with reference gene expression panels from 58 tissues and cell types to identify predicted gene expression associated with metabolite levels. This approach was also performed for cardiovascular outcomes in three independent large European studies (N = 700,000) to identify predicted gene expression additionally associated with cardiovascular risk. Finally, genetically informed mediation analysis was performed to infer causal mediation in the relationship between gene expression, metabolite levels and cardiovascular risk.

RESULTS

A total of 44 genetic loci were associated with 124 metabolites. Lead genetic variants included 11 non-synonymous variants. Predicted expression of 53 fine-mapped genes was associated with 108 metabolite levels; while predicted expression of 6 of these genes was also associated with cardiovascular outcomes, highlighting a new role for regulatory gene HCG27. Additionally, we found that atherogenic metabolite levels mediate the associations between gene expression and cardiovascular risk. Some of these genes showed stronger associations in immune tissues, providing further evidence of the role of immune cells in increasing cardiovascular risk.

CONCLUSIONS

These findings propose new gene targets that could be potential candidates for drug development aimed at lowering the risk of cardiovascular events through the modulation of blood atherogenic metabolite levels.

scTWAS Atlas: an integrative knowledgebase of single-cell transcriptome-wide association studies

Single-cell transcriptome-wide association studies (scTWAS) is a new method for conducting TWAS analysis at the cellular level to identify gene-trait associations with higher precision. This approach helps overcome the challenge of interpreting cell-type heterogeneity in traditional TWAS results. As the field of scTWAS rapidly advances, there is a growing need for additional database platforms to integrate this wealth of data and knowledge effectively. To address this gap, we present scTWAS Atlas (https://ngdc.cncb.ac.cn/sctwas/), a comprehensive database of scTWAS information integrating literature curation and data analysis. The current version of scTWAS Atlas amasses 2,765,211 associations encompassing 34 traits, 30 cell types, 9 cell conditions and 16,470 genes. The database features visualization tools, including an interactive knowledge graph that integrates single-cell expression quantitative trait loci (sc-eQTL) and scTWAS associations to build a multi-omics level regulatory network at the cellular level. Additionally, scTWAS Atlas facilitates cross-cell-type analysis, highlighting cell-type-specific and shared TWAS genes. The database is designed with user-friendly interfaces and allows for easy browsing, searching, and downloading of relevant information. Overall, scTWAS Atlas is instrumental in exploring the genetic regulatory mechanisms at the cellular level and shedding light on the role of various cell types in biological processes, offering novel insights for human health research.

Novel risk loci for COVID-19 hospitalization among admixed American populations

The genetic basis of severe COVID-19 has been thoroughly studied, and many genetic risk factors shared between populations have been identified. However, reduced sample sizes from non-European groups have limited the discovery of population-specific common risk loci. In this second study nested in the SCOURGE consortium, we conducted a genome-wide association study (GWAS) for COVID-19 hospitalization in admixed Americans, comprising a total of 4702 hospitalized cases recruited by SCOURGE and seven other participating studies in the COVID-19 Host Genetic Initiative. We identified four genome-wide significant associations, two of which constitute novel loci and were first discovered in Latin American populations (BAZ2B and DDIAS). A trans-ethnic meta-analysis revealed another novel cross-population risk locus in CREBBP. Finally, we assessed the performance of a cross-ancestry polygenic risk score in the SCOURGE admixed American cohort. This study constitutes the largest GWAS for COVID-19 hospitalization in admixed Latin Americans conducted to date. This allowed to reveal novel risk loci and emphasize the need of considering the diversity of populations in genomic research.

Gene discovery and biological insights into anxiety disorders from a large-scale multi-ancestry genome-wide association study

We leveraged information from more than 1.2 million participants, including 97,383 cases, to investigate the genetics of anxiety disorders across five continental groups. Through ancestry-specific and cross-ancestry genome-wide association studies, we identified 51 anxiety-associated loci, 39 of which were novel. In addition, polygenic risk scores derived from individuals of European descent were associated with anxiety in African, admixed American and East Asian groups. The heritability of anxiety was enriched for genes expressed in the limbic system, cerebral cortex, cerebellum, metencephalon, entorhinal cortex and brain stem. Transcriptome-wide and proteome-wide analyses highlighted 115 genes associated with anxiety through brain-specific and cross-tissue regulation. Anxiety also showed global and local genetic correlations with depression, schizophrenia and bipolar disorder and widespread pleiotropy with several physical health domains. Overall, this study expands our knowledge regarding the genetic risk and pathogenesis of anxiety disorders, highlighting the importance of investigating diverse populations and integrating multi-omics information.

Genetic Insights into Externalizing and Internalizing Traits through Integration of the Research Domain Criteria and Hierarchical Taxonomy of Psychopathology

Background

There is considerable comorbidity between externalizing (EXT) and internalizing (INT) psychopathology. Understanding the shared genetic underpinnings of these spectra is crucial for advancing knowledge of their biological bases and potential health impacts, and for informing empirical models like the Research Domain Criteria (RDoC) and Hierarchical Taxonomy of Psychopathology (HiTOP).

Methods

We conducted a multivariate genome-wide association study (GWAS) of EXT and INT psychopathology by applying genomic structural equation modeling to summary statistics from 16 EXT and INT traits in European-ancestry individuals (n = 16,400 to 1,074,629). Downstream analyses explored associations across RDoC units of analysis (i.e., genes, molecules, cells, circuits, physiology, and behaviors).

Results

The GWAS identified 409 lead single nucleotide polymorphisms (SNPs) for EXT, 85 for INT, and 256 for EXT+INT (i.e., shared) traits. Bivariate causal mixture models estimated that nearly all EXT and INT causal variants overlapped, despite a genetic correlation of 0.37 (SE = 0.02). Drug repurposing analyses identified potential therapeutic targets, including perturbagens affecting dopamine and serotonin pathways. EXT genes had enriched expression in GABAergic, cortical, and hippocampal neurons, while INT genes were more narrowly linked to GABAergic neurons. EXT+INT liability was associated with reduced grey matter volumes in the amygdala and subcallosal cortex.

Conclusions

These findings reveal both genetic overlap and distinct molecular and neurobiological pathways underlying EXT and INT psychopathology. By integrating genomic insights with the RDoC and HiTOP frameworks, this study advances our understanding of the mechanisms driving these dimensions of psychopathology.

Cross-population enhancement of PrediXcan predictions with a gnomAD-based east Asian reference framework

Over the past decade, genome-wide association studies have identified thousands of variants significantly associated with complex traits. For each locus, gene expression levels are needed to further explore its biological functions. To address this, the PrediXcan algorithm leverages large-scale reference data to impute the gene expression level from single nucleotide polymorphisms, and thus the gene-trait associations can be tested to identify the candidate causal genes. However, a challenge arises due to the fact that most reference data are from subjects of European ancestry, and the accuracy and robustness of predicted gene expression in subjects of East Asian (EAS) ancestry remains unclear. Here, we first simulated a variety of scenarios to explore the impact of the level of population diversity on gene expression. Population differentiated variants were estimated by using the allele frequency information from The Genome Aggregation Database. We found that the weights of a variants was the main factor that affected the gene expression predictions, and that ~70% of variants were significantly population differentiated based on proportion tests. To provide insights into this population effect on gene expression levels, we utilized the allele frequency information to develop a gene expression reference panel, Predict Asian-Population (PredictAP), for EAS ancestry. PredictAP can be viewed as an auxiliary tool for PrediXcan when using genotype data from EAS subjects.

An efficient, not-only-linear correlation coefficient based on clustering

Identifying meaningful patterns in data is crucial for understanding complex biological processes, particularly in transcriptomics, where genes with correlated expression often share functions or contribute to disease mechanisms. Traditional correlation coefficients, which primarily capture linear relationships, may overlook important nonlinear patterns. We introduce the clustermatch correlation coefficient (CCC), a not-only-linear coefficient that utilizes clustering to efficiently detect both linear and nonlinear associations. CCC outperforms standard methods by revealing biologically meaningful patterns that linear-only coefficients miss and is faster than state-of-the-art coefficients such as the maximal information coefficient. When applied to human gene expression data from genotype-tissue expression (GTEx), CCC identified robust linear relationships and nonlinear patterns, such as sex-specific differences, that are undetectable by standard methods. Highly ranked gene pairs were enriched for interactions in integrated networks built from protein-protein interactions, transcription factor regulation, and chemical and genetic perturbations, suggesting that CCC can detect functional relationships missed by linear-only approaches. CCC is a highly efficient, next-generation, not-only-linear correlation coefficient for genome-scale data. A record of this paper's transparent peer review process is included in the supplemental information.

Transcriptome- and proteome-wide association studies identify genes associated with renal cell carcinoma

We performed a series of integrative analyses including transcriptome-wide association studies (TWASs) and proteome-wide association studies (PWASs) of renal cell carcinoma (RCC) to nominate and prioritize molecular targets for laboratory investigation. On the basis of a genome-wide association study (GWAS) of 29,020 affected individuals and 835,670 control individuals and prediction models trained in transcriptomic reference models, our TWAS across four kidney transcriptomes (GTEx kidney cortex, kidney tubules, TCGA-KIRC [The Cancer Genome Atlas kidney renal clear-cell carcinoma], and TCGA-KIRP [TCGA kidney renal papillary cell carcinoma]) identified 38 gene associations (false-discovery rate <5%) in at least two of four transcriptomic panels and identified 12 genes that were independent of GWAS susceptibility regions. Analyses combining TWAS associations across 48 tissues from GTEx identified associations that were replicable in tumor transcriptomes for 23 additional genes. Analyses by the two major histologic types (clear-cell RCC and papillary RCC) revealed subtype-specific associations, although at least three gene associations were common to both subtypes. PWAS identified 13 associated proteins, all mapping to GWAS-significant loci. TWAS-identified genes were enriched for active enhancer or promoter regions in RCC tumors and hypoxia-inducible factor binding sites in relevant cell lines. Using gene expression correlation, common cancers (breast and prostate) and RCC risk factors (e.g., hypertension and BMI) display genetic contributions shared with RCC. Our work identifies potential molecular targets for RCC susceptibility for downstream functional investigation.

Investigating the shared genetic architecture between depression and subcortical volumes

Depression, a widespread and highly heritable mental health condition, profoundly affects millions of individuals worldwide. Neuroimaging studies have consistently revealed volumetric abnormalities in subcortical structures associated with depression. However, the genetic underpinnings shared between depression and subcortical volumes remain inadequately understood. Here, we investigate the extent of polygenic overlap using the bivariate causal mixture model (MiXeR), leveraging summary statistics from the largest genome-wide association studies for depression (N = 674,452) and 14 subcortical volumetric phenotypes (N = 33,224). Additionally, we identify shared genomic loci through conditional/conjunctional FDR analyses. MiXeR shows that subcortical volumetric traits share a substantial proportion of genetic variants with depression, with 44 distinct shared loci identified by subsequent conjunctional FDR analysis. These shared loci are predominantly located in intronic regions (58.7%) and non-coding RNA intronic regions (25.4%). The 269 protein-coding genes mapped by these shared loci exhibit specific developmental trajectories, with the expression level of 55 genes linked to both depression and subcortical volumes, and 30 genes linked to cognitive abilities and behavioral symptoms. These findings highlight a shared genetic architecture between depression and subcortical volumetric phenotypes, enriching our understanding of the neurobiological underpinnings of depression.

Genome-wide association studies of coffee intake in UK/US participants of European ancestry uncover cohort-specific genetic associations

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N = 130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across hundreds of biomarkers, health, and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from the UK Biobank (UKB; N = 334,659). We observed consistent positive genetic correlations with substance use and obesity in both cohorts. Other genetic correlations were discrepant, including positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in the UKB, and genetic correlations with cognition that were negative in 23andMe but positive in the UKB. Phenome-wide association study using polygenic scores of coffee intake derived from 23andMe or UKB summary statistics also revealed consistent associations with increased odds of obesity- and red blood cell-related traits, but all other associations were cohort-specific. Our study shows that the genetics of coffee intake associate with substance use and obesity across cohorts, but also that GWAS performed in different populations could capture cultural differences in the relationship between behavior and genetics.

Association of Genetically Predicted Skipping of COL4A4 Exon 27 with Hematuria and Albuminuria

Background:

Hematuria is an established sign of glomerular disease and can be associated with kidney failure, but there has been limited scientific study of this trait.

Methods:

Here, we combined genetic data from the UK Biobank with predicted gene expression and splicing from GTEx kidney cortex samples (n = 65) in a transcriptome-wide association study (TWAS) to identify additional potential biological mechanisms influencing hematuria.

Results:

The TWAS using kidney cortex identified significant associations for 5 genes in terms of expression and 3 significant splicing events. Notably, we identified an association between the skipping of COL4A4 exon 27, which is genetically predicted by intronic rs11898094 (minor allele frequency 13%), and hematuria. Association between this variant was also found with urinary albumin excretion. We found independent evidence supporting the same variant predicting this skipping event in glomeruli-derived mRNA transcriptomics data (n = 245) from NEPTUNE. The functional significance of loss of exon 27 was demonstrated using the split NanoLuc-based α3α4α5(IV) heterotrimer assay, in which type IV collagen heterotrimer formation was quantified by luminescence. The causal splicing variant for this skipping event is yet to be identified.

Conclusions:

In summary, by integrating multiple data types, we identify a potential splicing event associated with hematuria and albuminuria.

Human genetics and epigenetics of alcohol use disorder

Alcohol use disorder (AUD) is a prominent contributor to global morbidity and mortality. Its complex etiology involves genetics, epigenetics, and environmental factors. We review progress in understanding the genetics and epigenetics of AUD, summarizing the key findings. Advancements in technology over the decades have elevated research from early candidate gene studies to present-day genome-wide scans, unveiling numerous genetic and epigenetic risk factors for AUD. The latest GWAS on more than one million participants identified more than 100 genetic variants, and the largest epigenome-wide association studies (EWAS) in blood and brain samples have revealed tissue-specific epigenetic changes. Downstream analyses revealed enriched pathways, genetic correlations with other traits, transcriptome-wide association in brain tissues, and drug-gene interactions for AUD. We also discuss limitations and future directions, including increasing the power of GWAS and EWAS studies as well as expanding the diversity of populations included in these analyses. Larger samples, novel technologies, and analytic approaches are essential; these include whole-genome sequencing, multiomics, single-cell sequencing, spatial transcriptomics, deep-learning prediction of variant function, and integrated methods for disease risk prediction.

The genetic architecture of age at menarche and its causal effects on other traits

Age at menarche (AAM) is a sign of puberty of females. It is a heritable trait associated with various adult diseases. However, the genetic mechanism that determines AAM and links it to disease risk is poorly understood. Aiming to uncover the genetic basis for AAM, we conducted a joint association study in up to 438,089 women from 3 genome-wide association studies of European and East Asian ancestries. A series of bioinformatical analyses and causal inference were then followed to explore in-depth annotations at the associated loci and infer the causal relationship between AAM and other complex traits/diseases. This largest meta-analysis identified a total of 21 novel AAM associated loci at the genome wide significance level (P < 5.0  × 10-8), 4 of which were European ancestry-specific loci. Functional annotations prioritized 33 candidate genes at newly identified loci. Significant genetic correlations were observed between AAM and 67 complex traits. Further causal inference demonstrated the effects of AAM on 13 traits, including forced vital capacity (FVC), high blood pressure, age at first live birth, etc, indicating that earlier AAM causes lower FVC, worse lung function, hypertension and earlier age at first (last) live birth. Enrichment analysis identified 5 enriched tissues, including the hypothalamus middle, hypothalamo hypophyseal system, neurosecretory systems, hypothalamus and retina. Our findings may provide useful insights that elucidate the mechanisms determining AAM and the genetic interplay between AAM and some traits of women.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Multivariate, Multi-omic Analysis in 799,429 Individuals Identifies 134 Loci Associated with Somatoform Traits

Somatoform traits, which manifest as persistent physical symptoms without a clear medical cause, are prevalent and pose challenges to clinical practice. Understanding the genetic basis of these disorders could improve diagnostic and therapeutic approaches. With publicly available summary statistics, we conducted a multivariate genome-wide association study (GWAS) and multi-omic analysis of four somatoform traits-fatigue, irritable bowel syndrome, pain intensity, and health satisfaction-in 799,429 individuals genetically similar to Europeans. Using genomic structural equation modeling, GWAS identified 134 loci significantly associated with a somatoform common factor, including 44 loci not significant in the input GWAS and 8 novel loci for somatoform traits. Gene-property analyses highlighted an enrichment of genes involved in synaptic transmission and enriched gene expression in 12 brain tissues. Six genes, including members of the CD300 family, had putatively causal effects mediated by protein abundance. There was substantial polygenic overlap (76-83%) between the somatoform and externalizing, internalizing, and general psychopathology factors. Somatoform polygenic scores were associated most strongly with obesity, Type 2 diabetes, tobacco use disorder, and mood/anxiety disorders in independent biobanks. Drug repurposing analyses suggested potential therapeutic targets, including MEK inhibitors. Mendelian randomization indicated potentially protective effects of gut microbiota, including Ruminococcus bromii. These biological insights provide promising avenues for treatment development.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

GRPa-PRS: A risk stratification method to identify genetically-regulated pathways in polygenic diseases

Background

Polygenic risk scores (PRS) are tools used to evaluate an individual's susceptibility to polygenic diseases based on their genetic profile. A considerable proportion of people carry a high genetic risk but evade the disease. On the other hand, some individuals with a low risk of eventually developing the disease. We hypothesized that unknown counterfactors might be involved in reversing the PRS prediction, which might provide new insights into the pathogenesis, prevention, and early intervention of diseases.

Methods

We built a novel computational framework to identify genetically-regulated pathways (GRPas) using PRS-based stratification for each cohort. We curated two AD cohorts with genotyping data; the discovery (disc) and the replication (rep) datasets include 2722 and 2854 individuals, respectively. First, we calculated the optimized PRS model based on the three recent AD GWAS summary statistics for each cohort. Then, we stratified the individuals by their PRS and clinical diagnosis into six biologically meaningful PRS strata, such as AD cases with low/high risk and cognitively normal (CN) with low/high risk. Lastly, we imputed individual genetically-regulated expression (GReX) and identified differential GReX and GRPas between risk strata using gene-set enrichment and variational analyses in two models, with and without APOE effects. An orthogonality test was further conducted to verify those GRPas are independent of PRS risk. To verify the generalizability of other polygenic diseases, we further applied a default model of GRPa-PRS for schizophrenia (SCZ).

Results

For each stratum, we conducted the same procedures in both the disc and rep datasets for comparison. In AD, we identified several well-known AD-related pathways, including amyloid-beta clearance, tau protein binding, and astrocyte response to oxidative stress. Additionally, we discovered resilience-related GRPs that are orthogonal to AD PRS, such as the calcium signaling pathway and divalent inorganic cation homeostasis. In SCZ, pathways related to mitochondrial function and muscle development were highlighted. Finally, our GRPa-PRS method identified more consistent differential pathways compared to another variant-based pathway PRS method.

Conclusions

We developed a framework, GRPa-PRS, to systematically explore the differential GReX and GRPas among individuals stratified by their estimated PRS. The GReX-level comparison among those strata unveiled new insights into the pathways associated with disease risk and resilience. Our framework is extendable to other polygenic complex diseases.

MICOS Complex Loss Governs Age-Associated Murine Mitochondrial Architecture and Metabolism in the Liver, While Sam50 Dictates Diet Changes

The liver, the largest internal organ and a metabolic hub, undergoes significant declines due to aging, affecting mitochondrial function and increasing the risk of systemic liver diseases. How the mitochondrial three-dimensional (3D) structure changes in the liver across aging, and the biological mechanisms regulating such changes confers remain unclear. In this study, we employed Serial Block Face-Scanning Electron Microscopy (SBF-SEM) to achieve high-resolution 3D reconstructions of murine liver mitochondria to observe diverse phenotypes and structural alterations that occur with age, marked by a reduction in size and complexity. We also show concomitant metabolomic and lipidomic changes in aged samples. Aged human samples reflected altered disease risk. To find potential regulators of this change, we examined the Mitochondrial Contact Site and Cristae Organizing System (MICOS) complex, which plays a crucial role in maintaining mitochondrial architecture. We observe that the MICOS complex is lost during aging, but not Sam50. Sam50 is a component of the sorting and assembly machinery (SAM) complex that acts in tandem with the MICOS complex to modulate cristae morphology. In murine models subjected to a high-fat diet, there is a marked depletion of the mitochondrial protein SAM50. This reduction in Sam50 expression may heighten the susceptibility to liver disease, as our human biobank studies corroborate that Sam50 plays a genetically regulated role in the predisposition to multiple liver diseases. We further show that changes in mitochondrial calcium dysregulation and oxidative stress accompany the disruption of the MICOS complex. Together, we establish that a decrease in mitochondrial complexity and dysregulated metabolism occur with murine liver aging. While these changes are partially be regulated by age-related loss of the MICOS complex, the confluence of a murine high-fat diet can also cause loss of Sam50, which contributes to liver diseases. In summary, our study reveals potential regulators that affect age-related changes in mitochondrial structure and metabolism, which can be targeted in future therapeutic techniques.

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

A transcriptome-wide association study identified susceptibility genes for hepatocellular carcinoma in East Asia

Background

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide and is prevalent in East Asia. Although genome-wide association studies (GWASs) of HCC have identified 23 risk regions, the susceptibility genes underlying these associations largely remain unclear. To identify novel candidate genes for HCC, we conducted liver single-tissue and cross-tissue transcriptome-wide association studies (TWASs) in two populations of East Asia.

Methods

GWAS summary statistics of 2,514 subjects (1,161 HCC cases and 1,353 controls) from the Chinese Qidong cohort and 161,323 subjects (2,122 HCC cases and 159,201 controls) from the BioBank Japan project were used to conduct TWAS analysis. The single-tissue and cross-tissue TWAS approaches were both used to detect the association between susceptible genes and the risk of HCC. TWAS identified genes were further annotated by Metascape, UALCAN, GEPIA2, and DepMap.

Results

We identified 22 novel genes at 16 independent loci significantly associated with HCC risk after Bonferroni correction. Of these, 13 genes were located in novel regions. Besides, we found 83 genes overlapped in the Chinese and Japanese cohorts with P < 0.05, of which, three genes (NUAK2, HLA-DQA1, and ATP6V1G2) were discerned by both single-tissue and cross-tissue TWAS approaches. Among the genes identified through TWAS, a significant proportion of them exhibit a credible role in HCC biology, such as FAM96B, HSPA5, POLRMT, MPHOSPH10, and RABL2A. HLA-DQA1, NUAK2, and HSPA5 associated with the process of carcinogenesis in HCC as previously reported.

Conclusions

Our findings highlight the value of leveraging the gene expression data to identify new candidate genes beyond the GWAS associations and could further provide a genetic insight for the biology of HCC.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

An overview of detecting gene-trait associations by integrating GWAS summary statistics and eQTLs

Detecting genes that affect specific traits (such as human diseases and crop yields) is important for treating complex diseases and improving crop quality. A genome-wide association study (GWAS) provides new insights and directions for understanding complex traits by identifying important single nucleotide polymorphisms. Many GWAS summary statistics data related to various complex traits have been gathered recently. Studies have shown that GWAS risk loci and expression quantitative trait loci (eQTLs) often have a lot of overlaps, which makes gene expression gradually become an important intermediary to reveal the regulatory role of GWAS. In this review, we review three types of gene-trait association detection methods of integrating GWAS summary statistics and eQTLs data, namely colocalization methods, transcriptome-wide association study-oriented approaches, and Mendelian randomization-related methods. At the theoretical level, we discussed the differences, relationships, advantages, and disadvantages of various algorithms in the three kinds of gene-trait association detection methods. To further discuss the performance of various methods, we summarize the significant gene sets that influence high-density lipoprotein, low-density lipoprotein, total cholesterol, and triglyceride reported in 16 studies. We discuss the performance of various algorithms using the datasets of the four lipid traits. The advantages and limitations of various algorithms are analyzed based on experimental results, and we suggest directions for follow-up studies on detecting gene-trait associations.

Transcriptomic imputation of genetic risk variants uncovers novel whole-blood biomarkers of Parkinson's disease

Blood-based gene expression signatures could potentially be used as biomarkers for PD. However, it is unclear whether genetically-regulated transcriptomic signatures can provide novel gene candidates for use as PD biomarkers. We leveraged on the Genotype-Tissue Expression (GTEx) database to impute whole-blood transcriptomic expression using summary statistics of three large-scale PD GWAS. A random forest classifier was used with the consensus whole-blood imputed gene signature (IGS) to discriminate between cases and controls. Outcome measures included Area under the Curve (AUC) of Receiver Operating Characteristic (ROC) Curve. We demonstrated that the IGS (n = 37 genes) is conserved across PD GWAS studies and brain tissues. IGS discriminated between cases and controls in an independent whole-blood RNA-sequencing study (1176 PD, 254 prodromal, and 860 healthy controls) with mean AUC and accuracy of 64.8% and 69.4% for PD cohort, and 78.8% and 74% for prodromal cohort. PATL2 was the top-performing imputed gene in both PD and prodromal PD cohorts, whose classifier performance varied with biological sex (higher performance for males and females in the PD and prodromal PD, respectively). Single-cell RNA-sequencing studies (scRNA-seq) of healthy humans and PD patients found PATL2 to be enriched in terminal effector CD8+ and cytotoxic CD4+ cells, whose proportions are both increased in PD patients. We demonstrated the utility of GWAS transcriptomic imputation in identifying novel whole-blood transcriptomic signatures which could be leveraged upon for PD biomarker derivation. We identified PATL2 as a potential biomarker in both clinical and prodromic PD.

xWAS analysis in neuropsychiatric disorders by integrating multi-molecular phenotype quantitative trait loci and GWAS summary data

BACKGROUND

Integrating quantitative trait loci (QTL) data related to molecular phenotypes with genome-wide association study (GWAS) data is an important post-GWAS strategic approach employed to identify disease-associated molecular features. Various types of molecular phenotypes have been investigated in neuropsychiatric disorders. However, these findings pertaining to distinct molecular features are often independent of each other, posing challenges for having an overview of the mapped genes.

METHODS

In this study, we comprehensively summarized published analyses focusing on four types of risk-related molecular features (gene expression, splicing transcriptome, protein abundance, and DNA methylation) across five common neuropsychiatric disorders. Subsequently, we conducted supplementary analyses with the latest GWAS dataset and corresponding deficient molecular phenotypes using Functional Summary-based Imputation (FUSION) and summary data-based Mendelian randomization (SMR). Based on the curated and supplemented results, novel reliable genes and their functions were explored.

RESULTS

Our findings revealed that eQTL exhibited superior ability in prioritizing risk genes compared to the other QTL, followed by sQTL. Approximately half of the genes associated with splicing transcriptome, protein abundance, and DNA methylation were successfully replicated by eQTL-associated genes across all five disorders. Furthermore, we identified 436 novel reliable genes, which enriched in pathways related with neurotransmitter transportation such as synaptic, dendrite, vesicles, axon along with correlations with other neuropsychiatric disorders. Finally, we identified ten multiple molecular involved regulation patterns (MMRP), which may provide valuable insights into understanding the contribution of molecular regulation network targeting these disease-associated genes.

CONCLUSIONS

The analyses prioritized novel and reliable gene sets related with five molecular features based on published and supplementary results for five common neuropsychiatric disorders, which were missed in the original GWAS analysis. Besides, the involved MMRP behind these genes could be given priority for further investigation to elucidate the pathogenic molecular mechanisms underlying neuropsychiatric disorders in future studies.

Multi-tissue transcriptome-wide association study identifies novel candidate susceptibility genes for cataract

Introduction

Cataract is the leading cause of blindness among the elderly worldwide. Twin and family studies support an important role for genetic factors in cataract susceptibility with heritability estimates up to 58%. To date, 55 loci for cataract have been identified by genome-wide association studies (GWAS), however, much work remains to identify the causal genes. Here, we conducted a transcriptome-wide association study (TWAS) of cataract to prioritize causal genes and identify novel ones, and examine the impact of their expression.

Methods

We performed tissue-specific and multi-tissue TWAS analyses to assess associations between imputed gene expression from 54 tissues (including 49 from the Genotype Tissue Expression (GTEx) Project v8) with cataract using FUSION software. Meta-analyzed GWAS summary statistics from 59,944 cataract cases and 478,571 controls, all of European ancestry and from two cohorts (GERA and UK Biobank) were used. We then examined the expression of the novel genes in the lens tissue using the iSyTE database.

Results

Across tissue-specific and multi-tissue analyses, we identified 99 genes for which genetically predicted gene expression was associated with cataract after correcting for multiple testing. Of these 99 genes, 20 (AC007773.1, ANKH, ASIP, ATP13A2, CAPZB, CEP95, COQ6, CREB1, CROCC, DDX5, EFEMP1, EIF2S2, ESRRB, GOSR2, HERC4, INSRR, NIPSNAP2, PICALM, SENP3, and SH3YL1) did not overlap with previously reported cataract-associated loci. Tissue-specific analysis identified 202 significant gene-tissue associations for cataract, of which 166 (82.2%), representing 9 unique genes, were attributed to the previously reported 11q13.3 locus. Tissue-enrichment analysis revealed that gastrointestinal tissues represented one of the highest proportions of the Bonferroni-significant gene-tissue associations (21.3%). Moreover, this gastrointestinal tissue type was the only anatomical category significantly enriched in our results, after correcting for the number of tissue donors and imputable genes for each reference panel. Finally, most of the novel cataract genes (e.g., Capzb) were robustly expressed in iSyTE lens data.

Discussion

Our results provide evidence of the utility of imputation-based TWAS approaches to characterize known GWAS risk loci and identify novel candidate genes that may increase our understanding of cataract etiology. Our findings also highlight the fact that expression of genes associated with cataract susceptibility is not necessarily restricted to lens tissue.

Subset-based method for cross-tissue transcriptome-wide association studies improves power and interpretability

Integrating results from genome-wide association studies (GWASs) and studies of molecular phenotypes such as gene expressions can improve our understanding of the biological functions of trait-associated variants and can help prioritize candidate genes for downstream analysis. Using reference expression quantitative trait locus (eQTL) studies, several methods have been proposed to identify gene-trait associations, primarily based on gene expression imputation. To increase the statistical power by leveraging substantial eQTL sharing across tissues, meta-analysis methods aggregating such gene-based test results across multiple tissues or contexts have been developed as well. However, most existing meta-analysis methods have limited power to identify associations when the gene has weaker associations in only a few tissues and cannot identify the subset of tissues in which the gene is "activated." For this, we developed a cross-tissue subset-based transcriptome-wide association study (CSTWAS) meta-analysis method that improves power under such scenarios and can extract the set of potentially associated tissues. To improve applicability, CSTWAS uses only GWAS summary statistics and pre-computed correlation matrices to identify a subset of tissues that have the maximal evidence of gene-trait association. Through numerical simulations, we found that CSTWAS can maintain a well-calibrated type-I error rate, improves power especially when there is a small number of associated tissues for a gene-trait association, and identifies an accurate associated tissue set. By analyzing GWAS summary statistics of three complex traits and diseases, we demonstrate that CSTWAS could identify biological meaningful signals while providing an interpretation of disease etiology by extracting a set of potentially associated tissues.

A Multivariate Genome-Wide Association Study Reveals Neural Correlates and Common Biological Mechanisms of Psychopathology Spectra

There is considerable comorbidity across externalizing and internalizing behavior dimensions of psychopathology. We applied genomic structural equation modeling (gSEM) to genome-wide association study (GWAS) summary statistics to evaluate the factor structure of externalizing and internalizing psychopathology across 16 traits and disorders among European-ancestry individuals (n's = 16,400 to 1,074,629). We conducted GWAS on factors derived from well-fitting models. Downstream analyses served to identify biological mechanisms, explore drug repurposing targets, estimate genetic overlap between the externalizing and internalizing spectra, and evaluate causal effects of psychopathology liability on physical health. Both a correlated factors model, comprising two factors of externalizing and internalizing risk, and a higher-order single-factor model of genetic effects contributing to both spectra demonstrated acceptable t. GWAS identified 409 lead single nucleotide polymorphisms (SNPs) associated with externalizing and 85 lead SNPs associated with internalizing, while the second-order GWAS identified 256 lead SNPs contributing to broad psychopathology risk. In bivariate causal mixture models, nearly all externalizing and internalizing causal variants overlapped, despite a genetic correlation of only 0.37 (SE = 0.02) between them. Externalizing genes showed cell-type specific expression in GABAergic, cortical, and hippocampal neurons, and internalizing genes were associated with reduced subcallosal cortical volume, providing insight into the neurobiological underpinnings of psychopathology. Genetic liability for externalizing, internalizing, and broad psychopathology exerted causal effects on pain, general health, cardiovascular diseases, and chronic illnesses. These findings underscore the complex genetic architecture of psychopathology, identify potential biological pathways for the externalizing and internalizing spectra, and highlight the physical health burden of psychiatric comorbidity.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

postGWAS: A web server for deciphering the causality post the genome-wide association studies

While genome-wide association studies (GWAS) have unequivocally identified vast disease susceptibility variants, a majority of them are situated in non-coding regions and are in high linkage disequilibrium (LD). To pave the way of translating GWAS signals to clinical drug targets, it is essential to identify the underlying causal variants and further causal genes. To this end, a myriad of post-GWAS methods have been devised, each grounded in distinct principles including fine-mapping, co-localization, and transcriptome-wide association study (TWAS) techniques. Yet, no platform currently exists that seamlessly integrates these diverse post-GWAS methodologies. In this work, we present a user-friendly web server for post-GWAS analysis, that seamlessly integrates 9 distinct methods with 12 models, categorized by fine-mapping, colocalization, and TWAS. The server mainly helps users decipher the causality hindered by complex GWAS signals, including casual variants and casual genes, without the burden of computational skills and complex environment configuration, and provides a convenient platform for post-GWAS analysis, result visualization, facilitating the understanding and interpretation of the genome-wide association studies. The postGWAS server is available at http://g2g.biographml.com/.

Integration of human organoids single-cell transcriptomic profiles and human genetics repurposes critical cell type-specific drug targets for severe COVID-19

Human organoids recapitulate the cell type diversity and function of their primary organs holding tremendous potentials for basic and translational research. Advances in single-cell RNA sequencing (scRNA-seq) technology and genome-wide association study (GWAS) have accelerated the biological and therapeutic interpretation of trait-relevant cell types or states. Here, we constructed a computational framework to integrate atlas-level organoid scRNA-seq data, GWAS summary statistics, expression quantitative trait loci, and gene-drug interaction data for distinguishing critical cell populations and drug targets relevant to coronavirus disease 2019 (COVID-19) severity. We found that 39 cell types across eight kinds of organoids were significantly associated with COVID-19 outcomes. Notably, subset of lung mesenchymal stem cells increased proximity with fibroblasts predisposed to repair COVID-19-damaged lung tissue. Brain endothelial cell subset exhibited significant associations with severe COVID-19, and this cell subset showed a notable increase in cell-to-cell interactions with other brain cell types, including microglia. We repurposed 33 druggable genes, including IFNAR2, TYK2, and VIPR2, and their interacting drugs for COVID-19 in a cell-type-specific manner. Overall, our results showcase that host genetic determinants have cellular-specific contribution to COVID-19 severity, and identification of cell type-specific drug targets may facilitate to develop effective therapeutics for treating severe COVID-19 and its complications.

Gene Discovery and Biological Insights into Anxiety Disorders from a Multi-Ancestry Genome-wide Association Study of >1.2 Million Participants

We leveraged information from more than 1.2 million participants to investigate the genetics of anxiety disorders across five continental ancestral groups. Ancestry-specific and cross-ancestry genome-wide association studies identified 51 anxiety-associated loci, 39 of which are novel. Additionally, polygenic risk scores derived from individuals of European descent were associated with anxiety in African, Admixed-American, and East Asian groups. The heritability of anxiety was enriched for genes expressed in the limbic system, the cerebral cortex, the cerebellum, the metencephalon, the entorhinal cortex, and the brain stem. Transcriptome- and proteome-wide analyses highlighted 115 genes associated with anxiety through brain-specific and cross-tissue regulation. We also observed global and local genetic correlations with depression, schizophrenia, and bipolar disorder and putative causal relationships with several physical health conditions. Overall, this study expands the knowledge regarding the genetic risk and pathogenesis of anxiety disorders, highlighting the importance of investigating diverse populations and integrating multi-omics information.

Multivariate genome-wide association study of sleep health demonstrates unity and diversity

There has been a recent push to focus sleep research less on disordered sleep and more on the dimensional sleep health. Sleep health incorporates several dimensions of sleep: chronotype, efficiency, daytime alertness, duration, regularity, and satisfaction with sleep. A previous study demonstrated sleep health domains correlate only moderately with each other at the genomic level (|rGs| = 0.11-0.51) and show unique relationships with psychiatric domains (controlling for shared variances, duration, alertness, and non-insomnia independently related to a factor for internalizing psychopathology). Of the domains assessed, circadian preference was the least genetically correlated with all other facets of sleep health. This pattern is important because it suggests sleep health should be considered a multifaceted construct rather than a unitary construct. Prior genome-wide association studies (GWASs) have vastly increased our knowledge of the biological underpinnings of specific sleep traits but have only focused on univariate analyses. We present the first multivariate GWAS of sleep and circadian health (multivariate circadian preference, efficiency, and alertness factors, and three single-indicator factors of insomnia, duration, and regularity) using genomic structural equation modeling. We replicated loci found in prior sleep GWASs, but also discovered "novel" loci for each factor and found little evidence for genomic heterogeneity. While we saw overlapping genomic enrichment in subcortical brain regions and shared associations with external traits, much of the genetic architecture (loci, mapped genes, and enriched pathways) was diverse among sleep domains. These results confirm sleep health as a family of correlated but genetically distinct domains, which has important health implications.

Adjusting for genetic confounders in transcriptome-wide association studies improves discovery of risk genes of complex traits

Many methods have been developed to leverage expression quantitative trait loci (eQTL) data to nominate candidate genes from genome-wide association studies. These methods, including colocalization, transcriptome-wide association studies (TWAS) and Mendelian randomization-based methods; however, all suffer from a key problem-when assessing the role of a gene in a trait using its eQTLs, nearby variants and genetic components of other genes' expression may be correlated with these eQTLs and have direct effects on the trait, acting as potential confounders. Our extensive simulations showed that existing methods fail to account for these 'genetic confounders', resulting in severe inflation of false positives. Our new method, causal-TWAS (cTWAS), borrows ideas from statistical fine-mapping and allows us to adjust all genetic confounders. cTWAS showed calibrated false discovery rates in simulations, and its application on several common traits discovered new candidate genes. In conclusion, cTWAS provides a robust statistical framework for gene discovery.

Conditional transcriptome-wide association study for fine-mapping candidate causal genes

Transcriptome-wide association studies (TWASs) aim to integrate genome-wide association studies with expression-mapping studies to identify genes with genetically predicted expression (GReX) associated with a complex trait. In the present report, we develop a method, GIFT (gene-based integrative fine-mapping through conditional TWAS), that performs conditional TWAS analysis by explicitly controlling for GReX of all other genes residing in a local region to fine-map putatively causal genes. GIFT is frequentist in nature, explicitly models both expression correlation and cis-single nucleotide polymorphism linkage disequilibrium across multiple genes and uses a likelihood framework to account for expression prediction uncertainty. As a result, GIFT produces calibrated P values and is effective for fine-mapping. We apply GIFT to analyze six traits in the UK Biobank, where GIFT narrows down the set size of putatively causal genes by 32.16-91.32% compared with existing TWAS fine-mapping approaches. The genes identified by GIFT highlight the importance of vessel regulation in determining blood pressures and lipid metabolism for regulating lipid levels.

Dementia with Lewy Bodies: Genomics, Transcriptomics, and Its Future with Data Science

Dementia with Lewy bodies (DLB) is a significant public health issue. It is the second most common neurodegenerative dementia and presents with severe neuropsychiatric symptoms. Genomic and transcriptomic analyses have provided some insight into disease pathology. Variants within SNCA, GBA, APOE, SNCB, and MAPT have been shown to be associated with DLB in repeated genomic studies. Transcriptomic analysis, conducted predominantly on candidate genes, has identified signatures of synuclein aggregation, protein degradation, amyloid deposition, neuroinflammation, mitochondrial dysfunction, and the upregulation of heat-shock proteins in DLB. Yet, the understanding of DLB molecular pathology is incomplete. This precipitates the current clinical position whereby there are no available disease-modifying treatments or blood-based diagnostic biomarkers. Data science methods have the potential to improve disease understanding, optimising therapeutic intervention and drug development, to reduce disease burden. Genomic prediction will facilitate the early identification of cases and the timely application of future disease-modifying treatments. Transcript-level analyses across the entire transcriptome and machine learning analysis of multi-omic data will uncover novel signatures that may provide clues to DLB pathology and improve drug development. This review will discuss the current genomic and transcriptomic understanding of DLB, highlight gaps in the literature, and describe data science methods that may advance the field.

STW-MD: a novel spatio-temporal weighting and multi-step decision tree method for considering spatial heterogeneity in brain gene expression data

Gene expression during brain development or abnormal development is a biological process that is highly dynamic in spatio and temporal. Previous studies have mainly focused on individual brain regions or a certain developmental stage. Our motivation is to address this gap by incorporating spatio-temporal information to gain a more complete understanding of brain development or abnormal brain development, such as Alzheimer's disease (AD), and to identify potential determinants of response. In this study, we propose a novel two-step framework based on spatial-temporal information weighting and multi-step decision trees. This framework can effectively exploit the spatial similarity and temporal dependence between different stages and different brain regions, and facilitate differential gene analysis in brain regions with high heterogeneity. We focus on two datasets: the AD dataset, which includes gene expression data from early, middle and late stages, and the brain development dataset, spanning fetal development to adulthood. Our findings highlight the advantages of the proposed framework in discovering gene classes and elucidating their impact on brain development and AD progression across diverse brain regions and stages. These findings align with existing studies and provide insights into the processes of normal and abnormal brain development.

PigBiobank: a valuable resource for understanding genetic and biological mechanisms of diverse complex traits in pigs

To fully unlock the potential of pigs as both agricultural species for animal-based protein food and biomedical models for human biology and disease, a comprehensive understanding of molecular and cellular mechanisms underlying various complex phenotypes in pigs and how the findings can be translated to other species, especially humans, are urgently needed. Here, within the Farm animal Genotype-Tissue Expression (FarmGTEx) project, we build the PigBiobank (http://pigbiobank.farmgtex.org) to systematically investigate the relationships among genomic variants, regulatory elements, genes, molecular networks, tissues and complex traits in pigs. This first version of the PigBiobank curates 71 885 pigs with both genotypes and phenotypes from over 100 pig breeds worldwide, covering 264 distinct complex traits. The PigBiobank has the following functions: (i) imputed sequence-based genotype-phenotype associations via a standardized and uniform pipeline, (ii) molecular and cellular mechanisms underlying trait-associations via integrating multi-omics data, (iii) cross-species gene mapping of complex traits via transcriptome-wide association studies, and (iv) high-quality results display and visualization. The PigBiobank will be updated timely with the development of the FarmGTEx-PigGTEx project, serving as an open-access and easy-to-use resource for genetically and biologically dissecting complex traits in pigs and translating the findings to other species.

Unraveling phenotypic variance in metabolic syndrome through multi-omics

Complex multi-omics effects drive the clustering of cardiometabolic risk factors, underscoring the imperative to comprehend how individual and combined omics shape phenotypic variation. Our study partitions phenotypic variance in metabolic syndrome (MetS), blood glucose (GLU), triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and blood pressure through genome, transcriptome, metabolome, and exposome (i.e., lifestyle exposome) analyses. Our analysis included a cohort of 62,822 unrelated individuals with white British ancestry, sourced from the UK biobank. We employed linear mixed models to partition phenotypic variance using the restricted maximum likelihood (REML) method, implemented in MTG2 (v2.22). We initiated the analysis by individually modeling omics, followed by subsequent integration of pairwise omics in a joint model that also accounted for the covariance and interaction between omics layers. Finally, we estimated the correlations of various omics effects between the phenotypes using bivariate REML. Significant proportions of the MetS variance were attributed to distinct data sources: genome (9.47%), transcriptome (4.24%), metabolome (14.34%), and exposome (3.77%). The phenotypic variances explained by the genome, transcriptome, metabolome, and exposome ranged from 3.28% for GLU to 25.35% for HDL-C, 0% for GLU to 19.34% for HDL-C, 4.29% for systolic blood pressure (SBP) to 35.75% for TG, and 0.89% for GLU to 10.17% for HDL-C, respectively. Significant correlations were found between genomic and transcriptomic effects for TG and HDL-C. Furthermore, significant interaction effects between omics data were detected for both MetS and its components. Interestingly, significant correlation of omics effect between the phenotypes was found. This study underscores omics' roles, interaction effects, and random-effects covariance in unveiling phenotypic variation in multi-omics domains.

A compendium of genetic regulatory effects across pig tissues

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

A Positive Causal Effect of Shrimp Allergy on Major Depressive Disorder Mediated by Allergy- and Immune-Related Pathways in the East Asian Population

BACKGROUND

Observational studies have implied a potential correlation between allergic diseases and major depressive disorder (MDD). However, the relationship is still inconclusive as it is likely to be interfered with by substantial confounding factors and potential reverse causality. The present study aimed to investigate causal correlation of the two diseases by a Mendelian randomization (MR) study and further elucidate the underlying molecular mechanisms.

METHODS

With the biggest summary datasets of a genome-wide association study (GWAS) in the East Asian population, we conducted a two-sample, bidirectional MR study to assess the causal correlation between shrimp allergy (SA) and MDD. Subsequently, we identified the pleiotropic genes' susceptibility to the two diseases at whole-genome and tissue-specific levels, respectively. Enriched GO sets and KEGG pathways were also discovered to elucidate the potential underlying mechanisms.

RESULTS

With the most suitable MR method, SA was identified as a causal risk factor for MDD based on three different groups of independent genetic instruments, respectively (p < 2.81 × 10-2). In contrast, we did not observe a significant causal effect of MDD on SA. The GWAS-pairwise program successfully identified seven pleiotropic genetic variants (PPA3 > 0.8), indicating that the two diseases indeed have a shared genetic basis. At a whole-genome level, the MAGMA program identified 44 pleiotropic genes, which were enriched in allergy-related pathways, such as antigen processing and presentation pathway (p = 1.46 × 10-2). In brain-specific tissue, the S-MultiXcan program found 17 pleiotropic genes that were significantly enriched in immune-related pathways and GO sets, including asthma-related pathway, T-cell activation-related, and major histocompatibility complex protein-related GO sets. Regarding whole-blood tissue, the program identified six pleiotropic genes that are significantly enriched in tolerance induction-related GO sets.

CONCLUSIONS

The present study for the first time indicated a significant causal effect of SA on the occurrence of MDD, but the reverse was not true. Enrichment analyses of pleiotropic genes at whole-genome and tissue-specific levels implied the involvement of allergy and immune-related pathways in the shared genetic mechanism of the two diseases. Elucidating the causal effect and the acting direction may be beneficial in reducing the incidence rate of MDD for the massive group of SA patients in the East Asian region.

Isoform-level transcriptome-wide association uncovers genetic risk mechanisms for neuropsychiatric disorders in the human brain

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Multi-omics cannot replace sample size in genome-wide association studies

The integration of multi-omics information (e.g., epigenetics and transcriptomics) can be useful for interpreting findings from genome-wide association studies (GWAS). It has been suggested that multi-omics could circumvent or greatly reduce the need to increase GWAS sample sizes for novel variant discovery. We tested whether incorporating multi-omics information in earlier and smaller-sized GWAS boosts true-positive discovery of genes that were later revealed by larger GWAS of the same/similar traits. We applied 10 different analytic approaches to integrating multi-omics data from 12 sources (e.g., Genotype-Tissue Expression project) to test whether earlier and smaller GWAS of 4 brain-related traits (alcohol use disorder/problematic alcohol use, major depression/depression, schizophrenia, and intracranial volume/brain volume) could detect genes that were revealed by a later and larger GWAS. Multi-omics data did not reliably identify novel genes in earlier less-powered GWAS (PPV <0.2; 80% false-positive associations). Machine learning predictions marginally increased the number of identified novel genes, correctly identifying 1-8 additional genes, but only for well-powered early GWAS of highly heritable traits (i.e., intracranial volume and schizophrenia). Although multi-omics, particularly positional mapping (i.e., fastBAT, MAGMA, and H-MAGMA), can help to prioritize genes within genome-wide significant loci (PPVs = 0.5-1.0) and translate them into information about disease biology, it does not reliably increase novel gene discovery in brain-related GWAS. To increase power for discovery of novel genes and loci, increasing sample size is required.

Genome-wide association analysis of plasma lipidome identifies 495 genetic associations

The human plasma lipidome captures risk for cardiometabolic diseases. To discover new lipid-associated variants and understand the link between lipid species and cardiometabolic disorders, we perform univariate and multivariate genome-wide analyses of 179 lipid species in 7174 Finnish individuals. We fine-map the associated loci, prioritize genes, and examine their disease links in 377,277 FinnGen participants. We identify 495 genome-trait associations in 56 genetic loci including 8 novel loci, with a considerable boost provided by the multivariate analysis. For 26 loci, fine-mapping identifies variants with a high causal probability, including 14 coding variants indicating likely causal genes. A phenome-wide analysis across 953 disease endpoints reveals disease associations for 40 lipid loci. For 11 coronary artery disease risk variants, we detect strong associations with lipid species. Our study demonstrates the power of multivariate genetic analysis in correlated lipidomics data and reveals genetic links between diseases and lipid species beyond the standard lipids.

The pleiotropic contribution of genes in dopaminergic and serotonergic pathways to addiction and related behavioral traits

Introduction

Co-occurrence of substance use disorders (SUD) and other behavioral conditions, such as stress-related, aggressive or risk-taking behaviors, in the same individual has been frequently described. As dopamine (DA) and serotonin (5-HT) have been previously identified as key neurotransmitters for some of these phenotypes, we explored the genetic contribution of these pathways to SUD and these comorbid phenotypes in order to better understand the genetic relationship between them.

Methods

We tested the association of 275 dopaminergic genes and 176 serotonergic genes with these phenotypes by performing gene-based, gene-set and transcriptome-wide association studies in 11 genome-wide association studies (GWAS) datasets on SUD and related behaviors.

Results

At the gene-wide level, 68 DA and 27 5-HT genes were found to be associated with at least one GWAS on SUD or related behavior. Among them, six genes had a pleiotropic effect, being associated with at least three phenotypes: ADH1C, ARNTL, CHRNA3, HPRT1, HTR1B and DRD2. Additionally, we found nominal associations between the DA gene sets and SUD, opioid use disorder, antisocial behavior, irritability and neuroticism, and between the 5-HT-core gene set and neuroticism. Predicted gene expression correlates in brain were also found for 19 DA or 5-HT genes.

Discussion

Our study shows a pleiotropic contribution of dopaminergic and serotonergic genes to addiction and related behaviors such as anxiety, irritability, neuroticism and risk-taking behavior, highlighting a role for DA genes, which could explain, in part, the co-occurrence of these phenotypes.

Building an optimal predictive model for imputing tissue-specific gene expression by combining genotype and whole-blood transcriptome data

Accurate imputation of tissue-specific gene expression can be a powerful tool for understanding the biological mechanisms underlying human complex traits. Existing imputation methods can be grouped into two categories according to the types of predictors used. The first category uses genotype data, while the second category uses whole-blood expression data. Both data types can be easily collected from blood, avoiding invasive tissue biopsies. In this study, we attempted to build an optimal predictive model for imputing tissue-specific gene expression by combining the genotype and whole-blood expression data. We first evaluated the imputation performance of each standalone model (using genotype data [GEN model] and using whole-blood expression data [WBE model]) using their respective data types across 47 human tissues. The WBE model outperformed the GEN model in most tissues by a large gain. Then, we developed several combined models that leverage both types of predictors to further improve imputation performance. We tried various strategies, including utilizing a merged dataset of the two data types (MERGED models) and integrating the imputation outcomes of the two standalone models (inverse variance-weighted [IVW] models). We found that one of the MERGED models noticeably outperformed the standalone models. This model involved a fixed ratio between the two regularization penalty factors for the two predictor types so that the contribution of the whole-blood transcriptome is upweighted compared with the genotype. Our study suggests that one can improve the imputation of tissue-specific gene expression by combining the genotype and whole-blood expression, but the improvement can be largely dependent on the combination strategy chosen.

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.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.