Integrative approaches for large-scale transcriptome-wide association studies

Nominating novel proteins for anxiety via integrating human brain proteomes and genome-wide association study

BACKGROUND

The underlying pathogenesis of anxiety remain elusive, making the pinpointing of potential therapeutic and diagnostic biomarkers for anxiety paramount to its efficient treatment.

METHODS

We undertook a proteome-wide association study (PWAS), fusing human brain proteomes from both discovery (ROS/MAP; N = 376) and validation cohorts (Banner; N = 152) with anxiety genome-wide association study (GWAS) summary statistics. Complementing this, we executed transcriptome-wide association studies (TWAS) leveraging human brain transcriptomic data from the Common Mind Consortium (CMC) to discern the confluence of genetic influences spanning both proteomic and transcriptomic levels. We further scrutinized significant genes through a suite of methodologies.

RESULTS

We discerned 14 genes instrumental in the genesis of anxiety through their specific cis-regulated brain protein abundance. Out of these, 6 were corroborated in the confirmatory PWAS, with 4 also showing associations with anxiety via their cis-regulated brain mRNA levels. A heightened confidence level was attributed to 5 genes (RAB27B, CCDC92, BTN2A1, TMEM106B, and DOC2A), taking into account corroborative evidence from both the confirmatory PWAS and TWAS, coupled with insights from mendelian randomization analysis and colocalization evaluations. A majority of the identified genes manifest in brain regions intricately linked to anxiety and predominantly partake in lysosomal metabolic processes.

LIMITATIONS

The limited scope of the brain proteome reference datasets, stemming from a relatively modest sample size, potentially curtails our grasp on the entire gamut of genetic effects.

CONCLUSION

The genes pinpointed in our research present a promising groundwork for crafting therapeutic interventions and diagnostic tools for anxiety.

Novel genetic loci of inhibitory control in ADHD and healthy children and genetic correlations with ADHD

Cumulative evidence has showed the deficits of inhibitory control in patients with attention deficit hyperactivity disorder (ADHD), which is considered as an endophenotype of ADHD. Genetic study of inhibitory control could advance gene discovery and further facilitate the understanding of ADHD genetic basis, but the studies were limited in both the general population and ADHD patients. To reveal genetic risk variants of inhibitory control and its potential genetic relationship with ADHD, we conducted genome-wide association studies (GWAS) on inhibitory control using three datasets, which included 783 and 957 ADHD patients and 1350 healthy children. Subsequently, we employed polygenic risk scores (PRS) to explore the association of inhibitory control with ADHD and related psychiatric disorders. Firstly, we identified three significant loci for inhibitory control in the healthy dataset, two loci in the case dataset, and one locus in the meta-analysis of three datasets. Besides, we found more risk genes and variants by applying transcriptome-wide association study (TWAS) and conditional FDR method. Then, we constructed a network by connecting the genes identified in our study, leading to the identification of several vital genes. Lastly, we identified a potential relationship between inhibitory control and ADHD and autism by PRS analysis and found the direct and mediated contribution of the identified genetic loci on ADHD symptoms by mediation analysis. In conclusion, we revealed some genetic risk variants associated with inhibitory control and elucidated the benefit of inhibitory control as an endophenotype, providing valuable insights into the mechanisms underlying ADHD.

Comprehensive meta-analysis reveals distinct gene expression signatures of MASLD progression

Metabolic dysfunction-associated steatotic liver disease (MASLD) and its progressive form, metabolic dysfunction-associated steatohepatitis (MASH), pose significant risks of severe fibrosis, cirrhosis, and hepatocellular carcinoma. Despite their widespread prevalence, the molecular mechanisms underlying the development and progression of these common chronic hepatic conditions are not fully understood. Here, we conducted the most extensive meta-analysis of hepatic gene expression datasets from liver biopsy samples to date, integrating 10 RNA-sequencing and microarray datasets (1,058 samples). Using a random-effects meta-analysis model, we compared over 12,000 shared genes across datasets. We identified 685 genes differentially expressed in MASLD versus normal liver, 1,870 in MASH versus normal liver, and 3,284 in MASLD versus MASH. Integrating these results with genome-wide association studies and coexpression networks, we identified two functionally relevant, validated coexpression modules mainly driven by SMOC2, ITGBL1, LOXL1, MGP, SOD3, and TAT, HGD, SLC25A15, respectively, the latter not previously associated with MASLD and MASH. Our findings provide a comprehensive and robust analysis of hepatic gene expression alterations associated with MASLD and MASH and identify novel key drivers of MASLD progression.

From GWASs toward Mechanistic Understanding with Case Studies in Dermatogenetics

Many human skin diseases result from the complex interplay of genetic and environmental mechanisms that are largely unknown. GWASs have yielded insight into the genetic aspect of complex disease by highlighting regions of the genome or specific genetic variants associated with disease. Leveraging this information to identify causal genes and cell types will provide insight into fundamental biology, inform diagnostics, and aid drug discovery. However, the etiological mechanisms from genetic variant to disease are still unestablished in most cases. There now exists an unprecedented wealth of data and computational methods for variant interpretation in a functional context. It can be challenging to decide where to start owing to a lack of consensus on the best way to identify causal genetic mechanisms. This article highlights 3 key aspects of genetic variant interpretation: prioritizing causal genes, cell types, and pathways. We provide a practical overview of the main methods and datasets, giving examples from recent atopic dermatitis studies to provide a blueprint for variant interpretation. A collection of resources, including brief description and links to the packages and web tools, is provided for researchers looking to start in silico follow-up genetic analysis of associated genetic variants.

Integrative omics analysis elucidates the genetic basis underlying seed weight and oil content in soybean

Synergistic optimization of key agronomic traits by traditional breeding has dramatically enhanced crop productivity in the past decades. However, the genetic basis underlying coordinated regulation of yield- and quality-related traits remains poorly understood. Here, we dissected the genetic architectures of seed weight and oil content by combining genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) using 421 soybean (Glycine max) accessions. We identified 26 and 33 genetic loci significantly associated with seed weight and oil content by GWAS, respectively, and detected 5,276 expression quantitative trait loci (eQTLs) regulating expression of 3,347 genes based on population transcriptomes. Interestingly, a gene module (IC79), regulated by two eQTL hotspots, exhibited significant correlation with both seed weigh and oil content. Twenty-two candidate causal genes for seed traits were further prioritized by TWAS, including Regulator of Weight and Oil of Seed 1 (GmRWOS1), which encodes a sodium pump protein. GmRWOS1 was verified to pleiotropically regulate seed weight and oil content by gene knockout and overexpression. Notably, allelic variations of GmRWOS1 were strongly selected during domestication of soybean. This study uncovers the genetic basis and network underlying regulation of seed weight and oil content in soybean and provides a valuable resource for improving soybean yield and quality by molecular breeding.

Biobank-wide association scan identifies risk factors for late-onset Alzheimer's disease and endophenotypes

Rich data from large biobanks, coupled with increasingly accessible association statistics from genome-wide association studies (GWAS), provide great opportunities to dissect the complex relationships among human traits and diseases. We introduce BADGERS, a powerful method to perform polygenic score-based biobank-wide association scans. Compared to traditional approaches, BADGERS uses GWAS summary statistics as input and does not require multiple traits to be measured in the same cohort. We applied BADGERS to two independent datasets for late-onset Alzheimer's disease (AD; n=61,212). Among 1738 traits in the UK biobank, we identified 48 significant associations for AD. Family history, high cholesterol, and numerous traits related to intelligence and education showed strong and independent associations with AD. Furthermore, we identified 41 significant associations for a variety of AD endophenotypes. While family history and high cholesterol were strongly associated with AD subgroups and pathologies, only intelligence and education-related traits predicted pre-clinical cognitive phenotypes. These results provide novel insights into the distinct biological processes underlying various risk factors for AD.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

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

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

Integrating human endogenous retroviruses into transcriptome-wide association studies highlights novel risk factors for major psychiatric conditions

Human endogenous retroviruses (HERVs) are repetitive elements previously implicated in major psychiatric conditions, but their role in aetiology remains unclear. Here, we perform specialised transcriptome-wide association studies that consider HERV expression quantified to precise genomic locations, using RNA sequencing and genetic data from 792 post-mortem brain samples. In Europeans, we identify 1238 HERVs with expression regulated in cis, of which 26 represent expression signals associated with psychiatric disorders, with ten being conditionally independent from neighbouring expression signals. Of these, five are additionally significant in fine-mapping analyses and thus are considered high confidence risk HERVs. These include two HERV expression signatures specific to schizophrenia risk, one shared between schizophrenia and bipolar disorder, and one specific to major depressive disorder. No robust signatures are identified for autism spectrum conditions or attention deficit hyperactivity disorder in Europeans, or for any psychiatric trait in other ancestries, although this is likely a result of relatively limited statistical power. Ultimately, our study highlights extensive HERV expression and regulation in the adult cortex, including in association with psychiatric disorder risk, therefore providing a rationale for exploring neurological HERV expression in complex neuropsychiatric traits.

Identification of plasma protein markers of allergic disease risk: a mendelian randomization approach to proteomic analysis

BACKGROUND

While numerous allergy-related biomarkers and targeted treatment strategies have been developed and employed, there are still signifcant limitations and challenges in the early diagnosis and targeted treatment for allegic diseases. Our study aims to identify circulating proteins causally associated with allergic disease-related traits through Mendelian randomization (MR)-based analytical framework.

METHODS

Large-scale cis-MR was employed to estimate the effects of thousands of plasma proteins on five main allergic diseases. Additional analyses including MR Steiger analyzing and Bayesian colocalisation, were performed to test the robustness of the associations; These findings were further validated utilizing meta-analytical methods in the replication analysis. Both proteome- and transcriptome-wide association studies approach was applied, and then, a protein-protein interaction was conducted to examine the interplay between the identified proteins and the targets of existing medications.

RESULTS

Eleven plasma proteins were identified with links to atopic asthma (AA), atopic dermatitis (AD), and allergic rhinitis (AR). Subsequently, these proteins were classified into four distinct target groups, with a focus on tier 1 and 2 targets due to their higher potential to become drug targets. MR analysis and extra validation revealed STAT6 and TNFRSF6B to be Tier 1 and IL1RL2 and IL6R to be Tier 2 proteins with the potential for AA treatment. Two Tier 1 proteins, CRAT and TNFRSF6B, and five Tier 2 proteins, ERBB3, IL6R, MMP12, ICAM1, and IL1RL2, were linked to AD, and three Tier 2 proteins, MANF, STAT6, and TNFSF8, to AR.

CONCLUSION

Eleven Tier 1 and 2 protein targets that are promising drug target candidates were identified for AA, AD, and AR, which influence the development of allergic diseases and expose new diagnostic and therapeutic targets.

Genetic variants for head size share genes and pathways with cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Transcriptome-wide association study (TWAS) is a popular approach to dissect the functional consequence of disease associated non-coding variants. Most existing TWAS use bulk tissues and may not have the resolution to reveal cell-type specific target genes. Single-cell expression quantitative trait loci (sc-eQTL) datasets are emerging. The largest bulk- and sc-eQTL datasets are most conveniently available as summary statistics, but have not been broadly utilized in TWAS. Here, we present a new method EXPRESSO (EXpression PREdiction with Summary Statistics Only), to analyze sc-eQTL summary statistics, which also integrates 3D genomic data and epigenomic annotation to prioritize causal variants. EXPRESSO substantially improves existing methods. We apply EXPRESSO to analyze multi-ancestry GWAS datasets for 14 autoimmune diseases. EXPRESSO uniquely identifies 958 novel gene x trait associations, which is 26% more than the second-best method. Among them, 492 are unique to cell type level analysis and missed by TWAS using whole blood. We also develop a cell type aware drug repurposing pipeline, which leverages EXPRESSO results to identify drug compounds that can reverse disease gene expressions in relevant cell types. Our results point to multiple drugs with therapeutic potentials, including metformin for type 1 diabetes, and vitamin K for ulcerative colitis.

Total testosterone plays a crucial role in the pathway from hypothyroidism to broad depression in women

BACKGROUND

Depression tends to develop in correlation with hypothyroidism, however it's unclear how testosterone traits contribute to this association. We examined the causal association between depression, testosterone traits, and hypothyroidism using Mendelian randomization (MR).

METHOD

We conducted univariable and multivariable MR studies using summary-level statistics from genome-wide association studies (GWAS) of Hypothyroidism (n = 213,990), broad depression (n = 322,580), probable major depressive disorder (probable MDD) (n = 174,519), and International Classification of Diseases (ICD)-9 or ICD-10-coded MDD (n = 217,584) from European ancestry. The inverse variance weighted (IVW) method was used as the main MR analysis.

RESULTS

In univariate MR analysis, there is a positive causal relationship between hypothyroidism and broad depression (P = 0.0074; OR = 1.0066; 95%CI: 1.0018-1.0114) and probable MDD (P = 0.0242; OR = 1.0056; 95%CI: 1.0007-1.0105). In females, there is a causal relationship between hypothyroidism and decreased total testosterone (P < 0.001; OR = 0.9747; 95%CI: 0.9612-0.9885) and sex hormone binding globulin (SHBG) levels (P = 0.0418; OR = 0.9858; 95%CI: 0.9723-0.9995). In females, there is an inverse causal relationship between total testosterone and broad depression (P = 0.0349; OR = 0.9898; 95%CI: 0.9804-0.9993). Furthermore, in multivariate MR analysis, after adjusting for total testosterone in females, hypothyroidism only has a positive causal relationship with probable MDD, and the relationship with broad depression is no longer significant. Most notably, after adjusting for hypothyroidism, the inverse causal effect of female total testosterone levels on broad depression becomes more significant (P = 0.0154; OR = 0.9878; 95%CI: 0.9780-0.9977).

CONCLUSION

Hypothyroidism increases the risk of broad depression and probable MDD development. Total Testosterone appears to play an important role in the relationship between hypothyroidism and broad depression in female.

Genetic regulation of m6A RNA methylation and its contribution in human complex diseases

N6-methyladenosine (m6A) has been established as the most prevalent chemical modification in message RNA (mRNA), playing an essential role in determining the fate of RNA molecules. Dysregulation of m6A has been revealed to lead to abnormal physiological conditions and cause various types of human diseases. Recent studies have delineated the genetic regulatory maps for m6A methylation by mapping the quantitative trait loci of m6A (m6A-QTLs), thereby building up the regulatory circuits linking genetic variants, m6A, and human complex traits. Here, we review the recent discoveries concerning the genetic regulatory maps of m6A, describing the methodological and technical details of m6A-QTL identification, and introducing the key findings of the cis- and trans-acting drivers of m6A. We further delve into the tissue- and ethnicity-specificity of m6A-QTL, the association with other molecular phenotypes in light of genetic regulation, the regulators underlying m6A genetics, and importantly, the functional roles of m6A in mediating human complex diseases. Lastly, we discuss potential research avenues that can accelerate the translation of m6A genetics studies toward the development of therapies for human genetic diseases.

Multimodal analysis of RNA sequencing data powers discovery of complex trait genetics

Transcriptome data is commonly used to understand genome function via quantitative trait loci (QTL) mapping and to identify the molecular mechanisms driving genome wide association study (GWAS) signals through colocalization analysis and transcriptome-wide association studies (TWAS). While RNA sequencing (RNA-seq) has the potential to reveal many modalities of transcriptional regulation, such as various splicing phenotypes, such studies are often limited to gene expression due to the complexity of extracting and analyzing multiple RNA phenotypes. Here, we present Pantry (Pan-transcriptomic phenotyping), a framework to efficiently generate diverse RNA phenotypes from RNA-seq data and perform downstream integrative analyses with genetic data. Pantry currently generates phenotypes from six modalities of transcriptional regulation (gene expression, isoform ratios, splice junction usage, alternative TSS/polyA usage, and RNA stability) and integrates them with genetic data via QTL mapping, TWAS, and colocalization testing. We applied Pantry to Geuvadis and GTEx data, and found that 4,768 of the genes with no identified expression QTL in Geuvadis had QTLs in at least one other transcriptional modality, resulting in a 66% increase in genes over expression QTL mapping. We further found that QTLs exhibit modality-specific functional properties that are further reinforced by joint analysis of different RNA modalities. We also show that generalizing TWAS to multiple RNA modalities (xTWAS) approximately doubles the discovery of unique gene-trait associations, and enhances identification of regulatory mechanisms underlying GWAS signal in 42% of previously associated gene-trait pairs. We provide the Pantry code, RNA phenotypes from all Geuvadis and GTEx samples, and xQTL and xTWAS results on the web.

SR-TWAS: Leveraging Multiple Reference Panels to Improve TWAS Power by Ensemble Machine Learning

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for TWAS. To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies showed that SR-TWAS improved power, due to increased effective training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real application studies identified 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations were found for these significant risk genes.

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.

An integrative framework to prioritize genes in more than 500 loci associated with body mass index

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

Chronic Overlapping Pain Conditions and Nociplastic Pain

Chronic Overlapping Pain Conditions (COPCs) are a subset of chronic pain conditions commonly comorbid with one another and more prevalent in women and assigned female at birth (AFAB) individuals. Pain experience in these conditions may better fit with a new mechanistic pain descriptor, nociplastic pain, and nociplastic type pain may represent a shared underlying factor among COPCs. We applied GenomicSEM common-factor genome wide association study (GWAS) and multivariate transcriptome-wide association (TWAS) analyses to existing GWAS output for six COPCs in order to find genetic variation associated with nociplastic type pain, followed by genetic correlation (linkage-disequilibrium score regression), gene-set and tissue enrichment analyses. We found 24 independent single nucleotide polymorphisms (SNPs), and 127 unique genes significantly associated with nociplastic type pain, and showed nociplastic type pain to be a polygenic trait with significant SNP-heritability. We found significant genetic overlap between multisite chronic pain and nociplastic type pain, and to a smaller extent with rheumatoid arthritis and a neuropathic pain phenotype. Tissue enrichment analyses highlighted cardiac and thyroid tissue, and gene set enrichment analyses emphasized potential shared mechanisms in cognitive, personality, and metabolic traits and nociplastic type pain along with distinct pathology in migraine and headache. We use a well-powered network approach to investigate nociplastic type pain using existing COPC GWAS output, and show nociplastic type pain to be a complex, heritable trait, in addition to contributing to understanding of potential mechanisms in development of nociplastic pain.

Genetic polymorphisms associated with adverse pregnancy outcomes in nulliparas

Adverse pregnancy outcomes (APOs) affect a large proportion of pregnancies and represent an important cause of morbidity and mortality worldwide. Yet the pathophysiology of APOs is poorly understood, limiting our ability to prevent and treat these conditions. To search for genetic markers of maternal risk for four APOs, we performed multi-ancestry genome-wide association studies (GWAS) for pregnancy loss, gestational length, gestational diabetes, and preeclampsia. We clustered participants by their genetic ancestry and focused our analyses on three sub-cohorts with the largest sample sizes: European, African, and Admixed American. Association tests were carried out separately for each sub-cohort and then meta-analyzed together. Two novel loci were significantly associated with an increased risk of pregnancy loss: a cluster of SNPs located downstream of the TRMU gene (top SNP: rs142795512), and the SNP rs62021480 near RGMA. In the GWAS of gestational length we identified two new variants, rs2550487 and rs58548906 near WFDC1 and AC005052.1, respectively. Lastly, three new loci were significantly associated with gestational diabetes (top SNPs: rs72956265, rs10890563, rs79596863), located on or near ZBTB20, GUCY1A2, and RPL7P20, respectively. Fourteen loci previously correlated with preterm birth, gestational diabetes, and preeclampsia were found to be associated with these outcomes as well.

Genetic association of inflammatory marker GlycA with lung function and respiratory diseases

Association of circulating glycoprotein acetyls (GlycA), a systemic inflammation biomarker, with lung function and respiratory diseases remain to be investigated. We examined the genetic correlation, shared genetics, and potential causality of GlycA (N = 115,078) with lung function and respiratory diseases (N = 497,000). GlycA showed significant genetic correlation with FEV1 (rg = -0.14), FVC (rg = -0.18), asthma (rg = 0.21) and COPD (rg = 0.31). We consistently identified ten shared loci (including chr3p21.31 and chr8p23.1) at both SNP and gene level revealing potential shared biological mechanisms involving ubiquitination, immune response, Wnt/β-catenin signaling, cell growth and differentiation in tissues or cells including blood, epithelium, fibroblast, fetal thymus, and fetal intestine. Genetically elevated GlycA was significantly correlated with lung function and asthma susceptibility (354.13 ml decrement of FEV1, 442.28 ml decrement of FVC, and 144% increased risk of asthma per SD increment of GlycA) from MR analyses. Our findings provide insights into biological mechanisms of GlycA in relating to lung function, asthma, and COPD.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.

Using human genetics to understand the phenotypic association between chronotype and breast cancer

Little is known regarding the shared genetic influences underlying the observed phenotypic association between chronotype and breast cancer in women. Leveraging summary statistics from the hitherto largest genome-wide association study conducted in each trait, we investigated the genetic correlation, pleiotropic loci, and causal relationship of chronotype with overall breast cancer, and with its subtypes defined by the status of oestrogen receptor. We identified a negative genomic correlation between chronotype and overall breast cancer (r g  = -0.06, p = 3.00 × 10-4), consistent across oestrogen receptor-positive (r g  = -0.05, p = 3.30 × 10-3) and oestrogen receptor-negative subtypes (r g  = -0.05, p = 1.11 × 10-2). Five specific genomic regions were further identified as contributing a significant local genetic correlation. Cross-trait meta-analysis identified 78 loci shared between chronotype and breast cancer, of which 23 were novel. Transcriptome-wide association study revealed 13 shared genes, targeting tissues of the nervous, cardiovascular, digestive, and exocrine/endocrine systems. Mendelian randomisation demonstrated a significantly reduced risk of overall breast cancer (odds ratio 0.89, 95% confidence interval 0.83-0.94; p = 1.30 × 10-4) for genetically predicted morning chronotype. No reverse causality was found. Our work demonstrates an intrinsic link underlying chronotype and breast cancer, which may provide clues to inform management of sleep habits to improve female health.

Examining Differences in the Genetic and Functional Architecture of Attention-Deficit/Hyperactivity Disorder Diagnosed in Childhood and Adulthood

Background

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder with diagnostic criteria requiring symptoms to begin in childhood. We investigated whether individuals diagnosed as children differ from those diagnosed in adulthood with respect to shared and unique architecture at the genome-wide and gene expression level of analysis.

Methods

We used genomic structural equation modeling (SEM) to investigate differences in genetic correlations (rg) of childhood-diagnosed (ncases = 14,878) and adulthood-diagnosed (ncases = 6961) ADHD with 98 behavioral, psychiatric, cognitive, and health outcomes. We went on to apply transcriptome-wide SEM to identify functional annotations and patterns of gene expression associated with genetic risk sharing or divergence across the ADHD subgroups.

Results

Compared with the childhood subgroup, adulthood-diagnosed ADHD exhibited a significantly larger negative rg with educational attainment, the noncognitive skills of educational attainment, and age at first sexual intercourse. We observed a larger positive rg for adulthood-diagnosed ADHD with major depression, suicidal ideation, and a latent internalizing factor. At the gene expression level, transcriptome-wide SEM analyses revealed 22 genes that were significantly associated with shared genetic risk across the subtypes that reflected a mixture of coding and noncoding genes and included 15 novel genes relative to the ADHD subgroups.

Conclusions

This study demonstrated that ADHD diagnosed later in life shows much stronger genetic overlap with internalizing disorders and related traits. This may indicate the potential clinical relevance of distinguishing these subgroups or increased misdiagnosis for those diagnosed later in life. Top transcriptome-wide SEM results implicated genes related to neuronal function and clinical characteristics (e.g., sleep).

Integration of multi-omics summary data reveals the role of N6-methyladenosine in neuropsychiatric disorders

N6-methyladenosine (m6A) methylation regulates gene expression/protein by influencing numerous aspects of mRNA metabolism and contributes to neuropsychiatric diseases. Here, we integrated multi-omics data and genome-wide association study summary data of schizophrenia (SCZ), bipolar disorder (BP), attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD) to reveal the role of m6A in neuropsychiatric disorders by using transcriptome-wide association study (TWAS) tool and Summary-data-based Mendelian randomization (SMR). Our investigation identified 86 m6A sites associated with seven neuropsychiatric diseases and then revealed 7881 associations between m6A sites and gene expressions. Based on these results, we discovered 916 significant m6A-gene associations involving 82 disease-related m6A sites and 606 genes. Further integrating the 58 disease-related genes from TWAS and SMR analysis, we obtained 61, 8, 7, 3, and 2 associations linking m6A-disease, m6A-gene, and gene-disease for SCZ, BP, AD, MDD, and PD separately. Functional analysis showed the m6A mapped genes were enriched in "response to stimulus" pathway. In addition, we also analyzed the effect of gene expression on m6A and the post-transcription effect of m6A on protein. Our study provided new insights into the genetic component of m6A in neuropsychiatric disorders and unveiled potential pathogenic mechanisms where m6A exerts influences on disease through gene expression/protein regulation.

Large-scale genome-wide association studies reveal the genetic causal etiology between air pollutants and autoimmune diseases

BACKGROUND

Epidemiological evidence links a close correlation between long-term exposure to air pollutants and autoimmune diseases, while the causality remained unknown.

METHODS

Two-sample Mendelian randomization (TSMR) was used to investigate the role of PM10, PM2.5, NO2, and NOX (N = 423,796-456,380) in 15 autoimmune diseases (N = 14,890-314,995) using data from large European GWASs including UKB, FINNGEN, IMSGC, and IPSCSG. Multivariable Mendelian randomization (MVMR) was conducted to investigate the direct effect of each air pollutant and the mediating role of common factors, including body mass index (BMI), alcohol consumption, smoking status, and household income. Transcriptome-wide association studies (TWAS), two-step MR, and colocalization analyses were performed to explore underlying mechanisms between air pollution and autoimmune diseases.

RESULTS

In TSMR, after correction of multiple testing, hypothyroidism was causally associated with higher exposure to NO2 [odds ratio (OR): 1.37, p = 9.08 × 10-4] and NOX [OR: 1.34, p = 2.86 × 10-3], ulcerative colitis (UC) was causally associated with higher exposure to NOX [OR: 2.24, p = 1.23 × 10-2] and PM2.5 [OR: 2.60, p = 5.96 × 10-3], rheumatoid arthritis was causally associated with higher exposure to NOX [OR: 1.72, p = 1.50 × 10-2], systemic lupus erythematosus was causally associated with higher exposure to NOX [OR: 4.92, p = 6.89 × 10-3], celiac disease was causally associated with lower exposure to NOX [OR: 0.14, p = 6.74 × 10-4] and PM2.5 [OR: 0.17, p = 3.18 × 10-3]. The risky effects of PM2.5 on UC remained significant in MVMR analyses after adjusting for other air pollutants. MVMR revealed several common mediators between air pollutants and autoimmune diseases. Transcriptional analysis identified specific gene transcripts and pathways interconnecting air pollutants and autoimmune diseases. Two-step MR revealed that POR, HSPA1B, and BRD2 might mediate from air pollutants to autoimmune diseases. POR pQTL (rs59882870, PPH4=1.00) strongly colocalized with autoimmune diseases.

CONCLUSION

This research underscores the necessity of rigorous air pollutant surveillance within public health studies to curb the prevalence of autoimmune diseases.

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.

Prioritization of therapeutic targets for cancers using integrative multi-omics analysis

BACKGROUND

The integration of transcriptomic, proteomic, druggable genetic and metabolomic association studies facilitated a comprehensive investigation of molecular features and shared pathways for cancers' development and progression.

METHODS

Comprehensive approaches consisting of transcriptome-wide association studies (TWAS), proteome-wide association studies (PWAS), summary-data-based Mendelian randomization (SMR) and MR were performed to identify genes significantly associated with cancers. The results identified in above analyzes were subsequently involved in phenotype scanning and enrichment analyzes to explore the possible health effects and shared pathways. Additionally, we also conducted MR analysis   to investigate metabolic pathways related to cancers.

RESULTS

Totally 24 genes (18 transcriptomic, 1 proteomic and 5 druggable genetic) showed significant associations with cancers risk. All genes identified in multiple methods were mainly enriched in nuclear factor erythroid 2-related factor 2 (NRF2) pathway. Additionally, biosynthesis of ubiquinol and urate were found to play an important role in gastrointestinal tumors.

CONCLUSIONS

A set of putatively causal genes and pathways relevant to cancers were identified in this study, shedding light on the shared biological processes for tumorigenesis and providing compelling genetic evidence to prioritize anti-cancer drugs development.

Genome sequence analyses identify novel risk loci for multiple system atrophy

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.

Machine Learning of Three-Dimensional Protein Structures to Predict the Functional Impacts of Genome Variation

Research in the human genome sciences generates a substantial amount of genetic data for hundreds of thousands of individuals, which concomitantly increases the number of variants of unknown significance (VUS). Bioinformatic analyses can successfully reveal rare variants and variants with clear associations with disease-related phenotypes. These studies have had a significant impact on how clinical genetic screens are interpreted and how patients are stratified for treatment. There are few, if any, computational methods for variants comparable to biological activity predictions. To address this gap, we developed a machine learning method that uses protein three-dimensional structures from AlphaFold to predict how a variant will influence changes to a gene's downstream biological pathways. We trained state-of-the-art machine learning classifiers to predict which protein regions will most likely impact transcriptional activities of two proto-oncogenes, nuclear factor erythroid 2 (NFE2L2)-related factor 2 (NRF2) and c-Myc. We have identified classifiers that attain accuracies higher than 80%, which have allowed us to identify a set of key protein regions that lead to significant perturbations in c-Myc or NRF2 transcriptional pathway activities.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Tissue-specific atlas of trans-models for gene regulation elucidates complex regulation patterns

BACKGROUND

Deciphering gene regulation is essential for understanding the underlying mechanisms of healthy and disease states. While the regulatory networks formed by transcription factors (TFs) and their target genes has been mostly studied with relation to cis effects such as in TF binding sites, we focused on trans effects of TFs on the expression of their transcribed genes and their potential mechanisms.

RESULTS

We provide a comprehensive tissue-specific atlas, spanning 49 tissues of TF variations affecting gene expression through computational models considering two potential mechanisms, including combinatorial regulation by the expression of the TFs, and by genetic variants within the TF. We demonstrate that similarity between tissues based on our discovered genes corresponds to other types of tissue similarity. The genes affected by complex TF regulation, and their modelled TFs, were highly enriched for pharmacogenomic functions, while the TFs themselves were also enriched in several cancer and metabolic pathways. Additionally, genes that appear in multiple clusters are enriched for regulation of immune system while tissue clusters include cluster-specific genes that are enriched for biological functions and diseases previously associated with the tissues forming the cluster. Finally, our atlas exposes multilevel regulation across multiple tissues, where TFs regulate other TFs through the two tested mechanisms.

CONCLUSIONS

Our tissue-specific atlas provides hierarchical tissue-specific trans genetic regulations that can be further studied for association with human phenotypes.

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

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

DeLIVR: a deep learning approach to IV regression for testing nonlinear causal effects in transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) have been increasingly applied to identify (putative) causal genes for complex traits and diseases. TWAS can be regarded as a two-sample two-stage least squares method for instrumental variable (IV) regression for causal inference. The standard TWAS (called TWAS-L) only considers a linear relationship between a gene's expression and a trait in stage 2, which may lose statistical power when not true. Recently, an extension of TWAS (called TWAS-LQ) considers both the linear and quadratic effects of a gene on a trait, which however is not flexible enough due to its parametric nature and may be low powered for nonquadratic nonlinear effects. On the other hand, a deep learning (DL) approach, called DeepIV, has been proposed to nonparametrically model a nonlinear effect in IV regression. However, it is both slow and unstable due to the ill-posed inverse problem of solving an integral equation with Monte Carlo approximations. Furthermore, in the original DeepIV approach, statistical inference, that is, hypothesis testing, was not studied. Here, we propose a novel DL approach, called DeLIVR, to overcome the major drawbacks of DeepIV, by estimating a related but different target function and including a hypothesis testing framework. We show through simulations that DeLIVR was both faster and more stable than DeepIV. We applied both parametric and DL approaches to the GTEx and UK Biobank data, showcasing that DeLIVR detected additional 8 and 7 genes nonlinearly associated with high-density lipoprotein (HDL) cholesterol and low-density lipoprotein (LDL) cholesterol, respectively, all of which would be missed by TWAS-L, TWAS-LQ, and DeepIV; these genes include BUD13 associated with HDL, SLC44A2 and GMIP with LDL, all supported by previous studies.

Genetically Regulated Gene Expression in the Brain Associated With Chronic Pain: Relationships With Clinical Traits and Potential for Drug Repurposing

BACKGROUND

Chronic pain is a common, poorly understood condition. Genetic studies including genome-wide association studies have identified many relevant variants, which have yet to be translated into full understanding of chronic pain. Transcriptome-wide association studies using transcriptomic imputation methods such as S-PrediXcan can help bridge this genotype-phenotype gap.

METHODS

We carried out transcriptomic imputation using S-PrediXcan to identify genetically regulated gene expression associated with multisite chronic pain in 13 brain tissues and whole blood. Then, we imputed genetically regulated gene expression for over 31,000 Mount Sinai BioMe participants and performed a phenome-wide association study to investigate clinical relationships in chronic pain-associated gene expression changes.

RESULTS

We identified 95 experiment-wide significant gene-tissue associations (p < 7.97 × 10-7), including 36 unique genes and an additional 134 gene-tissue associations reaching within-tissue significance, including 53 additional unique genes. Of the 89 unique genes in total, 59 were novel for multisite chronic pain and 18 are established drug targets. Chronic pain genetically regulated gene expression for 10 unique genes was significantly associated with cardiac dysrhythmia, metabolic syndrome, disc disorders/dorsopathies, joint/ligament sprain, anemias, and neurologic disorder phecodes. Phenome-wide association study analyses adjusting for mean pain score showed that associations were not driven by mean pain score.

CONCLUSIONS

We carried out the largest transcriptomic imputation study of any chronic pain trait to date. Results highlight potential causal genes in chronic pain development and tissue and direction of effect. Several gene results were also drug targets. Phenome-wide association study results showed significant associations for phecodes including cardiac dysrhythmia and metabolic syndrome, thereby indicating potential shared mechanisms.

Exploration of the potential causative genes for inflammatory bowel disease: Transcriptome-wide association analysis, Mendelian randomization analysis and Bayesian colocalisation

Background

Inflammatory bowel disease (IBD) poses a complex challenge due to its intricate underlying mechanisms, and curative treatments remain elusive. Consequently, there is an urgent need to identify genes causally associated with IBD.

Methods

We extracted blood eQTL data from the GTExv8.ALL.Whole_Blood database, genome-wide association studies (GWAS) summary statistics of IBD from the IEU GWAS database, and performed a three-fold analysis protocol, including transcriptome-wide association analysis, Mendelian randomisation analysis, Bayesian colocalisation, and subsequent potential therapeutic agents identification.

Results

We identified four pathogenic genes, namely CARD9, RTEL1, STMN3 and ARFRP1, that promote the development of IBD, encompassing both ulcerative colitis (UC) and Crohn's disease (CD). Notably, ARFRP1 exhibited the ability to suppress IBD (encompassing UC and CD) development. Regarding drug prediction, cyclophosphamide emerged as a promising novel therapeutic option for IBD, encompassing UC and CD.

Conclusion

We identified several potential genes related to IBD (UC and CD), including CARD9, RTEL1, STMN3 and ARFRP1, warranting further investigation in functional studies to elucidate underlying disease mechanisms. Additionally, clinical studies exploring the potential of cyclophosphamide as a treatment avenue for IBD are warranted.

Shared Genetic Determinants of Schizophrenia and Autism Spectrum Disorder Implicate Opposite Risk Patterns: A Genome-Wide Analysis of Common Variants

BACKGROUND AND HYPOTHESIS

The synaptic pruning hypothesis posits that schizophrenia (SCZ) and autism spectrum disorder (ASD) may represent opposite ends of neurodevelopmental disorders: individuals with ASD exhibit an overabundance of synapses and connections while SCZ was characterized by excessive pruning of synapses and a reduction. Given the strong genetic predisposition of both disorders, we propose a shared genetic component, with certain loci having differential regulatory impacts.

STUDY DESIGN

Genome-Wide single nucleotide polymorphism (SNP) data of European descent from SCZ (N cases = 53 386, N controls = 77 258) and ASD (N cases = 18 381, N controls = 27 969) were analyzed. We used genetic correlation, bivariate causal mixture model, conditional false discovery rate method, colocalization, Transcriptome-Wide Association Study (TWAS), and Phenome-Wide Association Study (PheWAS) to investigate the genetic overlap and gene expression pattern.

STUDY RESULTS

We found a positive genetic correlation between SCZ and ASD (rg = .26, SE = 0.01, P = 7.87e-14), with 11 genomic loci jointly influencing both conditions (conjFDR <0.05). Functional analysis highlights a significant enrichment of shared genes during early to mid-fetal developmental stages. A notable genetic region on chromosome 17q21.31 (lead SNP rs2696609) showed strong evidence of colocalization (PP.H4.abf = 0.85). This SNP rs2696609 is linked to many imaging-derived brain phenotypes. TWAS indicated opposing gene expression patterns (primarily pseudogenes and long noncoding RNAs [lncRNAs]) for ASD and SCZ in the 17q21.31 region and some genes (LRRC37A4P, LINC02210, and DND1P1) exhibit considerable variation in the cerebellum across the lifespan.

CONCLUSIONS

Our findings support a shared genetic basis for SCZ and ASD. A common genetic variant, rs2696609, located in the Chr17q21.31 locus, may exert differential risk regulation on SCZ and ASD by altering brain structure. Future studies should focus on the role of pseudogenes, lncRNAs, and cerebellum in synaptic pruning and neurodevelopmental disorders.

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.

The broad impact of cell death genes on the human disease phenome

Cell death mediated by genetically defined signaling pathways influences the health and dynamics of all tissues, however the tissue specificity of cell death pathways and the relationships between these pathways and human disease are not well understood. We analyzed the expression profiles of an array of 44 cell death genes involved in apoptosis, necroptosis, and pyroptosis cell death pathways across 49 human tissues from GTEx, to elucidate the landscape of cell death gene expression across human tissues, and the relationship between tissue-specific genetically determined expression and the human phenome. We uncovered unique cell death gene expression profiles across tissue types, suggesting there are physiologically distinct cell death programs in different tissues. Using summary statistics-based transcriptome wide association studies (TWAS) on human traits in the UK Biobank (n ~ 500,000), we evaluated 513 traits encompassing ICD-10 defined diagnoses and laboratory-derived traits. Our analysis revealed hundreds of significant (FDR < 0.05) associations between genetically regulated cell death gene expression and an array of human phenotypes encompassing both clinical diagnoses and hematologic parameters, which were independently validated in another large-scale DNA biobank (BioVU) at Vanderbilt University Medical Center (n = 94,474) with matching phenotypes. Cell death genes were highly enriched for significant associations with blood traits versus non-cell-death genes, with apoptosis-associated genes enriched for leukocyte and platelet traits. Our findings are also concordant with independently published studies (e.g. associations between BCL2L11/BIM expression and platelet & lymphocyte counts). Overall, these results suggest that cell death genes play distinct roles in their contribution to human phenotypes, and that cell death genes influence a diverse array of human traits.

Phenotypic and genetic effect of carotid intima-media thickness on the risk of stroke

While carotid intima-media thickness (cIMT) as a noninvasive surrogate measure of atherosclerosis is widely considered a risk factor for stroke, the intrinsic link underlying cIMT and stroke has not been fully understood. We aimed to evaluate the clinical value of cIMT in stroke through the investigation of phenotypic and genetic relationships between cIMT and stroke. We evaluated phenotypic associations using observational data from UK Biobank (N = 21,526). We then investigated genetic relationships leveraging genomic data conducted in predominantly European ancestry for cIMT (N = 45,185) and any stroke (AS, Ncase/Ncontrol=40,585/406,111). Observational analyses suggested an increased hazard of stroke per one standard deviation increase in cIMT (cIMTmax-AS: hazard ratio (HR) = 1.39, 95%CI = 1.09-1.79; cIMTmean-AS: HR = 1.39, 95%CI = 1.09-1.78; cIMTmin-AS: HR = 1.32, 95%CI = 1.04-1.68). A positive global genetic correlation was observed (cIMTmax-AS: [Formula: see text]=0.23, P=9.44 × 10-5; cIMTmean-AS: [Formula: see text]=0.21, P=3.00 × 10-4; cIMTmin-AS: [Formula: see text]=0.16, P=6.30 × 10-3). This was further substantiated by five shared independent loci and 15 shared expression-trait associations. Mendelian randomization analyses suggested no causal effect of cIMT on stroke (cIMTmax-AS: odds ratio (OR)=1.12, 95%CI=0.97-1.28; cIMTmean-AS: OR=1.09, 95%CI=0.93-1.26; cIMTmin-AS: OR=1.03, 95%CI = 0.90-1.17). A putative association was observed for genetically predicted stroke on cIMT (AS-cIMTmax: beta=0.07, 95%CI = 0.01-0.13; AS-cIMTmean: beta=0.08, 95%CI = 0.01-0.15; AS-cIMTmin: beta = 0.08, 95%CI = 0.01-0.16) in the reverse direction MR, which attenuated to non-significant in sensitivity analysis. Our work does not find evidence supporting causal associations between cIMT and stroke. The pronounced cIMT-stroke association is intrinsic, and mostly attributed to shared genetic components. The clinical value of cIMT as a surrogate marker for stroke risk in the general population is likely limited.

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.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study

BACKGROUND

Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge.

RESULTS

In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes.

CONCLUSIONS

Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

A multi-ancestry genetic study of pain intensity in 598,339 veterans

Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects the quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids had a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well-characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 126 independent genetic loci, 69 of which are new. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level and cognitive traits. Integration of the genome-wide association studies findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, β-blockers and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.

Genetic contribution of reproductive traits to risk of uterine leiomyomata: a large-scale, genome-wide, cross-trait analysis

BACKGROUND

Although phenotypic associations between female reproductive characteristics and uterine leiomyomata have long been observed in epidemiologic investigations, the shared genetic architecture underlying these complex phenotypes remains unclear.

OBJECTIVE

We aimed to investigate the shared genetic basis, pleiotropic effects, and potential causal relationships underlying reproductive traits (age at menarche, age at natural menopause, and age at first birth) and uterine leiomyomata.

STUDY DESIGN

With the use of large-scale, genome-wide association studies conducted among women of European ancestry for age at menarche (n=329,345), age at natural menopause (n=201,323), age at first birth (n=418,758), and uterine leiomyomata (ncases/ncontrols=35,474/267,505), we performed a comprehensive, genome-wide, cross-trait analysis to examine systematically the common genetic influences between reproductive traits and uterine leiomyomata.

RESULTS

Significant global genetic correlations were identified between uterine leiomyomata and age at menarche (rg, -0.17; P=3.65×10-10), age at natural menopause (rg, 0.23; P=3.26×10-07), and age at first birth (rg, -0.16; P=1.96×10-06). Thirteen genomic regions were further revealed as contributing significant local correlations (P<.05/2353) to age at natural menopause and uterine leiomyomata. A cross-trait meta-analysis identified 23 shared loci, 3 of which were novel. A transcriptome-wide association study found 15 shared genes that target tissues of the digestive, exo- or endocrine, nervous, and cardiovascular systems. Mendelian randomization suggested causal relationships between a genetically predicted older age at menarche (odds ratio, 0.88; 95% confidence interval, 0.85-0.92; P=1.50×10-10) or older age at first birth (odds ratio, 0.95; 95% confidence interval, 0.90-0.99; P=.02) and a reduced risk for uterine leiomyomata and between a genetically predicted older age at natural menopause and an increased risk for uterine leiomyomata (odds ratio, 1.08; 95% confidence interval, 1.06-1.09; P=2.30×10-27). No causal association in the reverse direction was found.

CONCLUSION

Our work highlights that there are substantial shared genetic influences and putative causal links that underlie reproductive traits and uterine leiomyomata. The findings suggest that early identification of female reproductive risk factors may facilitate the initiation of strategies to modify potential uterine leiomyomata risk.

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

We present a gene-level regulatory model, single-cell ATAC + RNA linking (SCARlink), which predicts single-cell gene expression and links enhancers to target genes using multi-ome (scRNA-seq and scATAC-seq co-assay) sequencing data. The approach uses regularized Poisson regression on tile-level accessibility data to jointly model all regulatory effects at a gene locus, avoiding the limitations of pairwise gene-peak correlations and dependence on peak calling. SCARlink outperformed existing gene scoring methods for imputing gene expression from chromatin accessibility across high-coverage multi-ome datasets while giving comparable to improved performance on low-coverage datasets. Shapley value analysis on trained models identified cell-type-specific gene enhancers that are validated by promoter capture Hi-C and are 11× to 15× and 5× to 12× enriched in fine-mapped eQTLs and fine-mapped genome-wide association study (GWAS) variants, respectively. We further show that SCARlink-predicted and observed gene expression vectors provide a robust way to compute a chromatin potential vector field to enable developmental trajectory analysis.

Genomic SEM applied to explore etiological divergences in bipolar subtypes

BACKGROUND

Bipolar disorder (BD) is an overarching diagnostic class defined by the presence of at least one prior manic episode (BD I) or both a prior hypomanic episode and a prior depressive episode (BD II). Traditionally, BD II has been conceptualized as a less severe presentation of BD I, however, extant literature to investigate this claim has been mixed.

METHODS

We apply genomic structural equation modeling (Genomic SEM) to investigate divergent genetic pathways across BD's two major subtypes using the most recent GWAS summary statistics from the PGC. We begin by identifying divergences in genetic correlations across 98 external traits using a Bonferroni-corrected threshold. We also use a theoretically informed follow-up model to examine the extent to which the genetic variance in each subtype is explained by schizophrenia and major depression. Lastly, transcriptome-wide SEM (T-SEM) was used to identify neuronal gene expression patterns associated with BD subtypes.

RESULTS

BD II was characterized by significantly larger genetic overlap across non-psychiatric medical and internalizing traits (e.g. heart disease, neuroticism, insomnia), while stronger associations for BD I were absent. Consistent with these findings, follow-up modeling revealed a substantial major depression component for BD II. T-SEM results revealed 35 unique genes associated with shared risk across BD subtypes.

CONCLUSIONS

Divergent patterns of genetic relationships across external traits provide support for the distinction of the bipolar subtypes. However, our results also challenge the illness severity conceptualization of BD given stronger genetic overlap across BD II and a range of clinically relevant traits and disorders.

Cell subtype-specific effects of genetic variation in the Alzheimer's disease brain

The relationship between genetic variation and gene expression in brain cell types and subtypes remains understudied. Here, we generated single-nucleus RNA sequencing data from the neocortex of 424 individuals of advanced age; we assessed the effect of genetic variants on RNA expression in cis (cis-expression quantitative trait loci) for seven cell types and 64 cell subtypes using 1.5 million transcriptomes. This effort identified 10,004 eGenes at the cell type level and 8,099 eGenes at the cell subtype level. Many eGenes are only detected within cell subtypes. A new variant influences APOE expression only in microglia and is associated with greater cerebral amyloid angiopathy but not Alzheimer's disease pathology, after adjusting for APOEε4, providing mechanistic insights into both pathologies. Furthermore, only a TMEM106B variant affects the proportion of cell subtypes. Integration of these results with genome-wide association studies highlighted the targeted cell type and probable causal gene within Alzheimer's disease, schizophrenia, educational attainment and Parkinson's disease loci.