skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Gene-set Enrichment with Mathematical Biology (GEMB)
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
Award ID(s):
1714094
PAR ID:
10300312
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
GigaScience
Volume:
9
Issue:
10
ISSN:
2047-217X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3 , raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS. 
    more » « less
  2. Abstract BackgroundUncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. MethodsWe performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. ResultsThe single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. ConclusionsVDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals. 
    more » « less
  3. Abstract Transitioning from pluripotency to differentiated cell fates is fundamental to both embryonic development and adult tissue homeostasis. Improving our understanding of this transition would facilitate our ability to manipulate pluripotent cells into tissues for therapeutic use. Here, we show that membrane voltage (Vm) regulates the exit from pluripotency and the onset of germ layer differentiation in the embryo, a process that affects both gastrulation and left-right patterning. By examining candidate genes of congenital heart disease and heterotaxy, we identifyKCNH6, a member of the ether-a-go-go class of potassium channels that hyperpolarizes the Vmand thus limits the activation of voltage gated calcium channels, lowering intracellular calcium. In pluripotent embryonic cells, depletion ofkcnh6leads to membrane depolarization, elevation of intracellular calcium levels, and the maintenance of a pluripotent state at the expense of differentiation into ectodermal and myogenic lineages. Using high-resolution temporal transcriptome analysis, we identify the gene regulatory networks downstream of membrane depolarization and calcium signaling and discover that inhibition of the mTOR pathway transitions the pluripotent cell to a differentiated fate. By manipulating Vmusing a suite of tools, we establish a bioelectric pathway that regulates pluripotency in vertebrates, including human embryonic stem cells. 
    more » « less
  4. ABSTRACT The generaParamoebaandNeoparamoeba, within the family Paramoebidae (order Dactylopodida), are distinguished by their dactylopodial pseudopodia and the presence of an intracellular eukaryotic symbiont, thePerkinsela‐like organism (PLO). Taxonomic classification within these genera has been challenging due to overlapping morphological traits and close phylogenetic relationships.Theyare marine, with some playing significant roles as parasites. Notably, they have been implicated in sea urchin mass mortality events and are known causative agents of Amoebic Gill Disease (AGD) in fish. Despite their ecological and economic importance, many aspects of their diversity, biology, evolution, and host interactions remain poorly understood. In this study, we describe a novel amoeba species,Paramoeba daytonin. sp., isolated from Daytona Beach, Florida. Morphological and molecular analyses confirm its placement within theParamoebaclade, closely related toP. eilhardi,P. karteshi, andP. aparasomata. Phylogenetic assessments using 18S rDNA (18S) and Cytochrome c Oxidase I (COI) markers demonstrate the limitations of the 18S gene for species delineation, highlighting COI as a more reliable genetic marker for this group. Additionally, observations on PLO morphology, movement, and microtubule association provide insights into the endosymbiotic relationship, reinforcing the need for further research into this unique eukaryote‐eukaryote symbiosis. 
    more » « less
  5. 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