Epigenetic Element-Based Transcriptome-Wide Association Study Identifies Novel Genes for Bipolar Disorder

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified risk genes for bipolar disorder

OBJECTIVE

Bipolar disorder (BD) is a debilitating neuropsychiatric disorder, which is associated with genetic variation through "vast but mixed" Genome-Wide Association Studies (GWAS). Transcriptome-Wide Association Study (TWAS) is more effective in explaining genetic factors that influence complex diseases and can help identifying risk genes more reliably. So, this study aims to identify potential BD risk genes in pedigrees with TWAS.

METHODS

We conducted a TWAS analysis with expression quantitative trait loci (eQTL) analysis on extended BD pedigrees, and the BD genome-wide association study (GWAS) summary data acquired from the Psychiatric Genomics Consortium (PGC). Furthermore, the BD-associated genes identified by TWAS were validated by mRNA expression profiles from the Gene Expression Omnibus (GEO) Datasets (GSE23848 and GSE46416). Functional enrichment and annotation analysis were implemented by RStudio (version 4.2.0).

RESULTS

TWAS identified 362 genes with P value < 0.05, and 18 genes remain significant after Bonferroni correction, such as SEMA3G (PTWAS=1.07 × 10-11), ALOX5AP (PTWAS=3.12 × 10-8), and PLEC (PTWAS=1.27 × 10-7). Further 6 overlapped genes were detected in integrative analysis, such as UQCRB (PTWAS=0.0020, PmRNA=0.0000), TMPRSS9 (PTWAS=0.0405, PmRNA=0.0032), and SNX10 (PTWAS=0.0104, PmRNA=0.0015). Using genes identified by TWAS, Gene Ontology (GO) enrichment analysis identified 40 significant GO terms, such as mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory, aerobic electron transport chain, oxidative phosphorylation, mitochondrial membrane proteins, and ubiquinone activity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis identified significant 15 pathways for BD, such as Oxidative phosphorylation, endocannabinoids signaling, neurodegeneration, and reactive oxide species.

CONCLUSIONS

We found a set of BD-associated genes and pathways, validating the important role of neurodevelopmental abnormalities, inflammatory responses, and mitochondrial dysfunction in the pathology of BD, offering novel information for comprehending the genetic basis of BD.

Unraveling NEK4 as a Potential Drug Target in Schizophrenia and Bipolar I Disorder: A Proteomic and Genomic Approach

BACKGROUND AND HYPOTHESIS

Investigating the shared brain protein and genetic components of schizophrenia (SCZ) and bipolar I disorder (BD-I) presents a unique opportunity to understand the underlying pathophysiological processes and pinpoint potential drug targets.

STUDY DESIGN

To identify overlapping susceptibility brain proteins in SCZ and BD-I, we carried out proteome-wide association studies (PWAS) and Mendelian Randomization (MR) by integrating human brain protein quantitative trait loci with large-scale genome-wide association studies for both disorders. We utilized transcriptome-wide association studies (TWAS) to determine the consistency of mRNA-protein dysregulation in both disorders. We applied pleiotropy-informed conditional false discovery rate (pleioFDR) analysis to identify common risk genetic loci for SCZ and BD-I. Additionally, we performed a cell-type-specific analysis in the human brain to detect risk genes notably enriched in distinct brain cell types. The impact of risk gene overexpression on dendritic arborization and axon length in neurons was also examined.

STUDY RESULTS

Our PWAS identified 42 proteins associated with SCZ and 14 with BD-I, among which NEK4, HARS2, SUGP1, and DUS2 were common to both conditions. TWAS and MR analysis verified the significant risk gene NEK4 for both SCZ and BD-I. PleioFDR analysis further supported genetic risk loci associated with NEK4 for both conditions. The cell-type specificity analysis revealed that NEK4 is expressed on the surface of glutamatergic neurons, and its overexpression enhances dendritic arborization and axon length in cultured primary neurons.

CONCLUSIONS

These findings underscore a shared genetic origin for SCZ and BD-I, offering novel insights for potential therapeutic target identification.

Joint analysis of proteome, transcriptome, and multi-trait analysis to identify novel Parkinson's disease risk genes

Genome-wide association studies (GWAS) have identified multiple risk variants for Parkinson's disease (PD). Nevertheless, how the risk variants confer the risk of PD remains largely unknown. We conducted a proteome-wide association study (PWAS) and summary-data-based mendelian randomization (SMR) analysis by integrating PD GWAS with proteome and protein quantitative trait loci (pQTL) data from human brain, plasma and CSF. We also performed a large transcriptome-wide association study (TWAS) and Fine-mapping of causal gene sets (FOCUS), leveraging joint-tissue imputation (JTI) prediction models of 22 tissues to identify and prioritize putatively causal genes. We further conducted PWAS, SMR, TWAS, and FOCUS using a multi-trait analysis of GWAS (MTAG) to identify additional PD risk genes to boost statistical power. In this large-scale study, we identified 16 genes whose genetically regulated protein abundance levels were associated with Parkinson's disease risk. We undertook a large-scale analysis of PD and correlated traits, through TWAS and FOCUS studies, and discovered 26 casual genes related to PD that had not been reported in previous TWAS. 5 genes (CD38, GPNMB, RAB29, TMEM175, TTC19) showed significant associations with PD at both the proteome-wide and transcriptome-wide levels. Our study provides new insights into the etiology and underlying genetic architecture of PD.

Transcriptome-wide association studies: recent advances in methods, applications and available databases

Genome-wide association study has identified fruitful variants impacting heritable traits. Nevertheless, identifying critical genes underlying those significant variants has been a great task. Transcriptome-wide association study (TWAS) is an instrumental post-analysis to detect significant gene-trait associations focusing on modeling transcription-level regulations, which has made numerous progresses in recent years. Leveraging from expression quantitative loci (eQTL) regulation information, TWAS has advantages in detecting functioning genes regulated by disease-associated variants, thus providing insight into mechanisms of diseases and other phenotypes. Considering its vast potential, this review article comprehensively summarizes TWAS, including the methodology, applications and available resources.

Transcriptome-wide association study identifies novel candidate susceptibility genes for migraine

Genome-wide association studies (GWASs) have identified more than 130 genetic susceptibility loci for migraine; however, how most of these loci impact migraine development is unknown. To identify novel genes associated with migraine and interpret the transcriptional products of those genes, we conducted a transcriptome-wide association study (TWAS). We performed tissue-specific and multi-tissue TWAS analyses to assess associations between imputed gene expression from 53 tissues and migraine susceptibility using FUSION software. Meta-analyzed GWAS summary statistics from 26,052 migraine cases and 487,214 controls, all of European ancestry and from two cohorts (the Kaiser Permanente GERA and the UK Biobank), were used. We evaluated the associations for genes after conditioning on variant-level effects from GWAS, and we tested for colocalization of GWAS migraine-associated loci and expression quantitative trait loci (eQTLs). Across tissue-specific and multi-tissue analyses, we identified 53 genes for which genetically predicted gene expression was associated with migraine after correcting for multiple testing. Of these 53 genes, 10 (ATF5, CNTNAP1, KTN1-AS1, NEIL1, NEK4, NNT, PNKP, RUFY2, TUBG2, and VAT1) did not overlap known migraine-associated loci identified from GWAS. Tissue-specific analysis identified 45 gene-tissue pairs and cardiovascular tissues represented the highest proportion of the Bonferroni-significant gene-tissue pairs (n = 22 [49%]), followed by brain tissues (n = 6 [13%]), and gastrointestinal tissues (n = 4 [9%]). Colocalization analyses provided evidence of shared genetic variants underlying eQTL and GWAS signals in 18 of the gene-tissue pairs (40%). Our TWAS reports novel genes for migraine and highlights the important contribution of brain, cardiovascular, and gastrointestinal tissues in migraine susceptibility.

Epigenetic signatures relating to disease-associated genotypic burden in familial risk of bipolar disorder

Environmental factors contribute to risk of bipolar disorder (BD), but how environmental factors impact the development of psychopathology within the context of elevated genetic risk is unknown. We herein sought to identify epigenetic signatures operating in the context of polygenic risk for BD in young people at high familial risk (HR) of BD. Peripheral blood-derived DNA was assayed using Illumina PsychArray, and Methylation-450K or -EPIC BeadChips. Polygenic risk scores (PRS) were calculated using summary statistics from recent genome-wide association studies for BD, major depressive disorder (MDD) and cross-disorder (meta-analysis of eight psychiatric disorders). Unrelated HR participants of European ancestry (n = 103) were stratified based on their BD-PRS score within the HR-population distribution, and the top two quintiles (High-BD-PRS; n = 41) compared against the bottom two quintiles (Low-BD-PRS; n = 41). The High-BD-PRS stratum also had higher mean cross-disorder-PRS and MDD-PRS (ANCOVA p = 0.035 and p = 0.024, respectively). We evaluated DNA methylation differences between High-BD-PRS and Low-BD-PRS strata using linear models. One differentially methylated probe (DMP) (cg00933603; p = 3.54 × 10^-7) in VARS2, a mitochondrial aminoacyl-tRNA synthetase, remained significantly hypomethylated after multiple-testing correction. Overall, BD-PRS appeared to broadly impact epigenetic processes, with 1,183 genes mapped to nominal DMPs (p < 0.05); these displayed convergence with genes previously associated with BD, schizophrenia, chronotype, and risk taking. We tested poly-methylomic epigenetic profiles derived from nominal DMPs in two independent samples (n = 54 and n = 82, respectively), and conducted an exploratory evaluation of the effects of family environment, indexing cohesion and flexibility. This study highlights an important interplay between heritable risk and epigenetic factors, which warrant further exploration.

Multi-omic insights into Parkinson's Disease: From genetic associations to functional mechanisms

Genome-Wide Association Studies (GWAS) have elucidated the genetic components of Parkinson's Disease (PD). However, because the vast majority of GWAS association signals fall within non-coding regions, translating these results into an interpretable, mechanistic understanding of the disease etiology remains a major challenge in the field. In this review, we provide an overview of the approaches to prioritize putative causal variants and genes as well as summarise the primary findings of previous studies. We then discuss recent efforts to integrate multi-omics data to identify likely pathogenic cell types and biological pathways implicated in PD pathogenesis. We have compiled full summary statistics of cell-type, tissue, and phentoype enrichment analyses from multiple studies of PD GWAS and provided them in a standardized format as a resource for the research community (https://github.com/RajLabMSSM/PD_omics_review). Finally, we discuss the experimental, computational, and conceptual advances that will be necessary to fully elucidate the effects of functional variants and genes on cellular dysregulation and disease risk.

A transcriptome-wide association study identifies susceptibility genes for Parkinson's disease

Genome-wide association study (GWAS) has seen great strides in revealing initial insights into the genetic architecture of Parkinson's disease (PD). Since GWAS signals often reside in non-coding regions, relatively few of the associations have implicated specific biological mechanisms. Here, we aimed to integrate the GWAS results with large-scale expression quantitative trait loci (eQTL) in 13 brain tissues to identify candidate causal genes for PD. We conducted a transcriptome-wide association study (TWAS) for PD using the summary statistics of over 480,000 individuals from the most recent PD GWAS. We identified 18 genes significantly associated with PD after Bonferroni corrections. The most significant gene, LRRC37A2, was associated with PD in all 13 brain tissues, such as in the hypothalamus (P = 6.12 × 10-22) and nucleus accumbens basal ganglia (P = 5.62 × 10-21). We also identified eight conditionally independent genes, including four new genes at known PD loci: CD38, LRRC37A2, RNF40, and ZSWIM7. Through conditional analyses, we demonstrated that several of the GWAS significant signals on PD could be driven by genetically regulated gene expression. The most significant TWAS gene LRRC37A2 accounts for 0.855 of the GWAS signal at its loci, and ZSWIM7 accounts for all the GWAS signals at its loci. We further identified several phenotypes previously associated with PD by querying the single nucleotide polymorphisms (SNPs) in the final model of the identified genes in phenome databases. In conclusion, we prioritized genes that are likely to affect PD by using a TWAS approach and identified phenotypes associated with PD.