More Like this

1. Illuminating links between cis-regulators and trans-acting variants in the human prefrontal cortex

   https://doi.org/10.1186/s13073-022-01133-8  

   Liu, Shuang ; Won, Hyejung ; Clarke, Declan ; Matoba, Nana ; Khullar, Saniya ; Mu, Yudi ; Wang, Daifeng ; Gerstein, Mark ( November 2022 , Genome Medicine)

   Neuropsychiatric disorders afflict a large portion of the global population and constitute a significant source of disability worldwide. Although Genome-wide Association Studies (GWAS) have identified many disorder-associated variants, the underlying regulatory mechanisms linking them to disorders remain elusive, especially those involving distant genomic elements. Expression quantitative trait loci (eQTLs) constitute a powerful means of providing this missing link. However, most eQTL studies in human brains have focused exclusively on cis-eQTLs, which link variants to nearby genes (i.e., those within 1 Mb of a variant). A complete understanding of disease etiology requires a clearer understanding of trans-regulatory mechanisms, which, in turn, entails a detailed analysis of the relationships between variants and expression changes in distant genes.

   By leveraging large datasets from the PsychENCODE consortium, we conducted a genome-wide survey of trans-eQTLs in the human dorsolateral prefrontal cortex. We also performed colocalization and mediation analyses to identify mediators in trans-regulation and use trans-eQTLs to link GWAS loci to schizophrenia risk genes.

   We identified ~80,000 candidate trans-eQTLs (at FDR<0.25) that influence the expression of ~10K target genes (i.e., “trans-eGenes”). We found that many variants associated with these candidate trans-eQTLs overlap with known cis-eQTLs. Moreover, for >60% of these variants (by colocalization), the cis-eQTL’s target gene acts as a mediator for the trans-eQTL SNP's effect on the trans-eGene, highlighting examples of cis-mediation as essential for trans-regulation. Furthermore, many of these colocalized variants fall into a discernable pattern wherein cis-eQTL’s target is a transcription factor or RNA-binding protein, which, in turn, targets the gene associated with the candidate trans-eQTL. Finally, we show that trans-regulatory mechanisms provide valuable insights into psychiatric disorders: beyond what had been possible using only cis-eQTLs, we link an additional 23 GWAS loci and 90 risk genes (using colocalization between candidate trans-eQTLs and schizophrenia GWAS loci).

   We demonstrate that the transcriptional architecture of the human brain is orchestrated by both cis- and trans-regulatory variants and found that trans-eQTLs provide insights into brain-disease biology.

    

   more » « less
2. Integration of postmortem amygdala expression profiling, GWAS, and functional cell culture assays: neuroticism-associated synaptic vesicle glycoprotein 2A (SV2A) gene is regulated by miR-133a and miR-218

   https://doi.org/10.1038/s41398-020-00966-4  

   Jurkiewicz, Magdalena ; Moser, Dirk ; Koller, Antonius ; Yu, Lei ; Chen, Emily I. ; Bennett, David A. ; Canli, Turhan ( August 2020 , Translational Psychiatry)

   Recent genome-wide studies have begun to identify gene variants, expression profiles, and regulators associated with neuroticism, anxiety disorders, and depression. We conducted a set of experimental cell culture studies of gene regulation by micro RNAs (miRNAs), based on genome-wide transcriptome, proteome, and miRNA expression data from twentypostmortemsamples of lateral amygdala from donors with known neuroticism scores. Using Ingenuity Pathway Analysis and TargetScan, we identified a list of mRNA–protein–miRNA sets whose expression patterns were consistent with miRNA-based translational repression, as a function of trait anxiety. Here, we focused on one gene from that list, which is of particular translational significance in Psychiatry: synaptic vesicle glycoprotein 2A (SV2A) is the binding site of the anticonvulsant drug levetiracetam ((S)-α-Ethyl-2-oxo-1-pyrrolidineacetamide), which has shown promise in anxiety disorder treatments. We confirmed thatSV2Ais associated with neuroticism or anxiety using an original GWAS of a community cohort (N = 1,706), and cross-referencing a published GWAS of multiple cohorts (Ns ranging from 340,569 to 390,278).Postmortemamygdala expression profiling implicated three putative regulatory miRNAs to targetSV2A: miR-133a, miR-138, and miR-218. Moving from association to experimental causal testing in cell culture, we used a luciferase assay to demonstrate that miR-133a and miR-218, but not miR-138, significantly decreased relative luciferase activity from theSV2Adual-luciferase construct. In human neuroblastoma cells, transfection with miR-133a and miR-218 reduced both endogenousSV2AmRNA and protein levels, confirming miRNA targeting of theSV2Agene. This study illustrates the utility of combiningpostmortemgene expression data with GWAS to guide experimental cell culture assays examining gene regulatory mechanisms that may contribute to complex human traits. Identifying specific molecular mechanisms of gene regulation may be useful for future clinical applications in anxiety disorders or other forms of psychopathology.

    

   more » « less
3. Selective Serotonin Reuptake Inhibitor Pharmaco-Omics: Mechanisms and Prediction

   https://doi.org/10.3389/fphar.2020.614048  

   Nguyen, Thanh Thanh ; Liu, Duan ; Ho, Ming-Fen ; Athreya, Arjun P. ; Weinshilboum, Richard ( January 2021 , Frontiers in Pharmacology)

   null (Ed.)

   Selective serotonin reuptake inhibitors (SSRIs) are a standard of care for the pharmacotherapy of patients suffering from Major Depressive Disorder (MDD). However, only one-half to two-thirds of MDD patients respond to SSRI therapy. Recently, a “multiple omics” research strategy was applied to identify genetic differences between patients who did and did not respond to SSRI therapy. As a first step, plasma metabolites were assayed using samples from the 803 patients in the PGRN-AMPS SSRI MDD trial. The metabolomics data were then used to “inform” genomics by performing a genome-wide association study (GWAS) for plasma concentrations of the metabolite most highly associated with clinical response, serotonin (5-HT). Two genome-wide or near genome-wide significant single nucleotide polymorphism (SNP) signals were identified, one that mapped near the TSPAN5 gene and another across the ERICH3 gene, both genes that are highly expressed in the brain. Knocking down TSPAN5 and ERICH3 resulted in decreased 5-HT concentrations in neuroblastoma cell culture media and decreased expression of enzymes involved in 5-HT biosynthesis and metabolism. Functional genomic studies demonstrated that ERICH3 was involved in clathrin-mediated vesicle formation and TSPAN5 was an ethanol-responsive gene that may be a marker for response to acamprosate pharmacotherapy of alcohol use disorder (AUD), a neuropsychiatric disorder highly co-morbid with MDD. In parallel studies, kynurenine was the plasma metabolite most highly associated with MDD symptom severity and application of a metabolomics-informed pharmacogenomics approach identified DEFB1 and AHR as genes associated with variation in plasma kynurenine levels. Both genes also contributed to kynurenine-related inflammatory pathways. Finally, a multiply replicated predictive algorithm for SSRI clinical response with a balanced predictive accuracy of 76% (compared with 56% for clinical data alone) was developed by including the SNPs in TSPAN5 , ERICH3 , DEFB1 and AHR . In summary, application of a multiple omics research strategy that used metabolomics to inform genomics, followed by functional genomic studies, identified novel genes that influenced monoamine biology and made it possible to develop a predictive algorithm for SSRI clinical outcomes in MDD. A similar pharmaco-omic research strategy might be broadly applicable for the study of other neuropsychiatric diseases and their drug therapy. 

   more » « less
4. Gene-set Enrichment with Mathematical Biology (GEMB)

   https://doi.org/10.1093/gigascience/giaa091  

   Cochran, Amy L ; Nieser, Kenneth J ; Forger, Daniel B ; Zöllner, Sebastian ; McInnis, Melvin G ( October 2020 , GigaScience)

   null (Ed.)

   Abstract Background Gene-set analyses measure the association between a disease of interest and a “set" of genes related to a biological pathway. These analyses often incorporate gene network properties to account for differential contributions of each gene. We extend this concept further—defining gene contributions based on biophysical properties—by leveraging mathematical models of biology to predict the effects of genetic perturbations on a particular downstream function. Results We present a method that combines gene weights from model predictions and gene ranks from genome-wide association studies into a weighted gene-set test. We demonstrate in simulation how such a method can improve statistical power. To this effect, we identify a gene set, weighted by model-predicted contributions to intracellular calcium ion concentration, that is significantly related to bipolar disorder in a small dataset (P = 0.04; n = 544). We reproduce this finding using publicly available summary data from the Psychiatric Genomics Consortium (P = 1.7 × 10−4; n = 41,653). By contrast, an approach using a general calcium signaling pathway did not detect a significant association with bipolar disorder (P = 0.08). The weighted gene-set approach based on intracellular calcium ion concentration did not detect a significant relationship with schizophrenia (P = 0.09; n = 65,967) or major depression disorder (P = 0.30; n = 500,199). Conclusions Together, these findings show how incorporating math biology into gene-set analyses might help to identify biological functions that underlie certain polygenic disorders. 

   more » « less
5. Myasthenia gravis genome-wide association study implicates AGRN as a risk locus

   https://doi.org/10.1136/jmedgenet-2021-107953  

   Topaloudi, Apostolia ; Zagoriti, Zoi ; Flint, Alyssa Camille ; Martinez, Melanie Belle ; Yang, Zhiyu ; Tsetsos, Fotis ; Christou, Yiolanda-Panayiota ; Lagoumintzis, George ; Yannaki, Evangelia ; Zamba-Papanicolaou, Eleni ; et al ( January 2021 , Journal of Medical Genetics)

   Background Myasthenia gravis (MG) is a rare autoimmune disorder affecting the neuromuscular junction (NMJ). Here, we investigate the genetic architecture of MG via a genome-wide association study (GWAS) of the largest MG data set analysed to date. Methods We performed GWAS meta-analysis integrating three different data sets (total of 1401 cases and 3508 controls). We carried out human leucocyte antigen (HLA) fine-mapping, gene-based and tissue enrichment analyses and investigated genetic correlation with 13 other autoimmune disorders as well as pleiotropy across MG and correlated disorders. Results We confirmed the previously reported MG association with TNFRSF11A (rs4369774; p=1.09×10 −13 , OR=1.4). Furthermore, gene-based analysis revealed AGRN as a novel MG susceptibility gene. HLA fine-mapping pointed to two independent MG loci: HLA-DRB1 and HLA-B . MG onset-specific analysis reveals differences in the genetic architecture of early-onset MG (EOMG) versus late-onset MG (LOMG). Furthermore, we find MG to be genetically correlated with type 1 diabetes (T1D), rheumatoid arthritis (RA), late-onset vitiligo and autoimmune thyroid disease (ATD). Cross-disorder meta-analysis reveals multiple risk loci that appear pleiotropic across MG and correlated disorders. Discussion Our gene-based analysis identifies AGRN as a novel MG susceptibility gene, implicating for the first time a locus encoding a protein (agrin) that is directly relevant to NMJ activation. Mutations in AGRN have been found to underlie congenital myasthenic syndrome. Our results are also consistent with previous studies highlighting the role of HLA and TNFRSF11A in MG aetiology and the different risk genes in EOMG versus LOMG. Finally, we uncover the genetic correlation of MG with T1D, RA, ATD and late-onset vitiligo, pointing to shared underlying genetic mechanisms. 

   more » « less