scholarly journals scGRNom: a computational pipeline of integrative multi-omics analyses for predicting cell-type disease genes and regulatory networks

2021 ◽  
Vol 13 (1) ◽  
Author(s):  
Ting Jin ◽  
Peter Rehani ◽  
Mufang Ying ◽  
Jiawei Huang ◽  
Shuang Liu ◽  
...  

AbstractUnderstanding cell-type-specific gene regulatory mechanisms from genetic variants to diseases remains challenging. To address this, we developed a computational pipeline, scGRNom (single-cell Gene Regulatory Network prediction from multi-omics), to predict cell-type disease genes and regulatory networks including transcription factors and regulatory elements. With applications to schizophrenia and Alzheimer’s disease, we predicted disease genes and regulatory networks for excitatory and inhibitory neurons, microglia, and oligodendrocytes. Further enrichment analyses revealed cross-disease and disease-specific functions and pathways at the cell-type level. Our machine learning analysis also found that cell-type disease genes improved clinical phenotype predictions. scGRNom is a general-purpose tool available at https://github.com/daifengwanglab/scGRNom.

2021 ◽  
Vol 119 (1) ◽  
pp. e2115601119
Author(s):  
Shining Ma ◽  
Xi Chen ◽  
Xiang Zhu ◽  
Philip S. Tsao ◽  
Wing Hung Wong

Abdominal aortic aneurysm (AAA) is a common degenerative cardiovascular disease whose pathobiology is not clearly understood. The cellular heterogeneity and cell-type-specific gene regulation of vascular cells in human AAA have not been well-characterized. Here, we performed analysis of whole-genome sequencing data in AAA patients versus controls with the aim of detecting disease-associated variants that may affect gene regulation in human aortic smooth muscle cells (AoSMC) and human aortic endothelial cells (HAEC), two cell types of high relevance to AAA disease. To support this analysis, we generated H3K27ac HiChIP data for these cell types and inferred cell-type-specific gene regulatory networks. We observed that AAA-associated variants were most enriched in regulatory regions in AoSMC, compared with HAEC and CD4+ cells. The cell-type-specific regulation defined by this HiChIP data supported the importance of ERG and the KLF family of transcription factors in AAA disease. The analysis of regulatory elements that contain noncoding variants and also are differentially open between AAA patients and controls revealed the significance of the interleukin-6-mediated signaling pathway. This finding was further validated by including information from the deleteriousness effect of nonsynonymous single-nucleotide variants in AAA patients and additional control data from the Medical Genome Reference Bank dataset. These results shed important insights into AAA pathogenesis and provide a model for cell-type-specific analysis of disease-associated variants.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Rongxin Fang ◽  
Sebastian Preissl ◽  
Yang Li ◽  
Xiaomeng Hou ◽  
Jacinta Lucero ◽  
...  

AbstractIdentification of the cis-regulatory elements controlling cell-type specific gene expression patterns is essential for understanding the origin of cellular diversity. Conventional assays to map regulatory elements via open chromatin analysis of primary tissues is hindered by sample heterogeneity. Single cell analysis of accessible chromatin (scATAC-seq) can overcome this limitation. However, the high-level noise of each single cell profile and the large volume of data pose unique computational challenges. Here, we introduce SnapATAC, a software package for analyzing scATAC-seq datasets. SnapATAC dissects cellular heterogeneity in an unbiased manner and map the trajectories of cellular states. Using the Nyström method, SnapATAC can process data from up to a million cells. Furthermore, SnapATAC incorporates existing tools into a comprehensive package for analyzing single cell ATAC-seq dataset. As demonstration of its utility, SnapATAC is applied to 55,592 single-nucleus ATAC-seq profiles from the mouse secondary motor cortex. The analysis reveals ~370,000 candidate regulatory elements in 31 distinct cell populations in this brain region and inferred candidate cell-type specific transcriptional regulators.


2020 ◽  
Author(s):  
Mufang Ying ◽  
Peter Rehani ◽  
Panagiotis Roussos ◽  
Daifeng Wang

AbstractStrong phenotype-genotype associations have been reported across brain diseases. However, understanding underlying gene regulatory mechanisms remains challenging, especially at the cellular level. To address this, we integrated the multi-omics data at the cellular resolution of the human brain: cell-type chromatin interactions, epigenomics and single cell transcriptomics, and predicted cell-type gene regulatory networks linking transcription factors, distal regulatory elements and target genes (e.g., excitatory and inhibitory neurons, microglia, oligodendrocyte). Using these cell-type networks and disease risk variants, we further identified the cell-type disease genes and regulatory networks for schizophrenia and Alzheimer’s disease. The celltype regulatory elements (e.g., enhancers) in the networks were also found to be potential pleiotropic regulatory loci for a variety of diseases. Further enrichment analyses including gene ontology and KEGG pathways revealed potential novel cross-disease and disease-specific molecular functions, advancing knowledge on the interplays among genetic, transcriptional and epigenetic risks at the cellular resolution between neurodegenerative and neuropsychiatric diseases. Finally, we summarized our computational analyses as a general-purpose pipeline for predicting gene regulatory networks via multi-omics data.


2020 ◽  
Author(s):  
Nil Aygün ◽  
Angela L. Elwell ◽  
Dan Liang ◽  
Michael J. Lafferty ◽  
Kerry E. Cheek ◽  
...  

SummaryInterpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing is mainly performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements of cells present during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map e/sQTLs and allele specific expression in primary human neural progenitors (n=85) and their sorted neuronal progeny (n=74). Using colocalization and TWAS, we uncover cell-type specific regulatory mechanisms underlying risk for these traits.


Author(s):  
Jieru Li ◽  
Alexandros Pertsinidis

Establishing cell-type-specific gene expression programs relies on the action of distal enhancers, cis-regulatory elements that can activate target genes over large genomic distances — up to Mega-bases away. How distal enhancers physically relay regulatory information to target promoters has remained a mystery. Here, we review the latest developments and insights into promoter–enhancer communication mechanisms revealed by live-cell, real-time single-molecule imaging approaches.


2020 ◽  
Vol 48 (6) ◽  
pp. 2880-2896 ◽  
Author(s):  
Jun Li ◽  
Ting Zhang ◽  
Aarthi Ramakrishnan ◽  
Bernd Fritzsch ◽  
Jinshu Xu ◽  
...  

Abstract The transcription factor Six1 is essential for induction of sensory cell fate and formation of auditory sensory epithelium, but how it activates gene expression programs to generate distinct cell-types remains unknown. Here, we perform genome-wide characterization of Six1 binding at different stages of auditory sensory epithelium development and find that Six1-binding to cis-regulatory elements changes dramatically at cell-state transitions. Intriguingly, Six1 pre-occupies enhancers of cell-type-specific regulators and effectors before their expression. We demonstrate in-vivo cell-type-specific activity of Six1-bound novel enhancers of Pbx1, Fgf8, Dusp6, Vangl2, the hair-cell master regulator Atoh1 and a cascade of Atoh1’s downstream factors, including Pou4f3 and Gfi1. A subset of Six1-bound sites carry consensus-sequences for its downstream factors, including Atoh1, Gfi1, Pou4f3, Gata3 and Pbx1, all of which physically interact with Six1. Motif analysis identifies RFX/X-box as one of the most significantly enriched motifs in Six1-bound sites, and we demonstrate that Six1-RFX proteins cooperatively regulate gene expression through binding to SIX:RFX-motifs. Six1 targets a wide range of hair-bundle regulators and late Six1 deletion disrupts hair-bundle polarity. This study provides a mechanistic understanding of how Six1 cooperates with distinct cofactors in feedforward loops to control lineage-specific gene expression programs during progressive differentiation of the auditory sensory epithelium.


1992 ◽  
Vol 12 (8) ◽  
pp. 3653-3662
Author(s):  
P Lowings ◽  
U Yavuzer ◽  
C R Goding

Melanocytes are specialized cells residing in the hair follicles, the eye, and the basal layer of the human epidermis whose primary function is the production of the pigment melanin, giving rise to skin, hair, and eye color. Melanogenesis, a process unique to melanocytes that involves the processing of tyrosine by a number of melanocyte-specific enzymes, including tyrosinase and tyrosinase-related protein 1 (TRP-1), occurs only after differentiation from the melanocyte precursor, the melanoblast. In humans, melanogenesis is inducible by UV irradiation, with melanin being transferred from the melanocyte in the epidermis to the surrounding keratinocytes as protection from UV-induced damage. Excessive exposure to UV, however, is the primary cause of malignant melanoma, an increasingly common and highly aggressive disease. As an initial approach to understanding the regulation of melanocyte differentiation and melanocyte-specific transcription, we have isolated the gene encoding TRP-1 and examined the cis- and trans-acting factors required for cell-type-specific expression. We find that the TRP-1 promoter comprises both positive and negative regulatory elements which confer efficient expression in a TRP-1-expressing, pigmented melanoma cell line but not in NIH 3T3 or JEG3 cells and that a minimal promoter extending between -44 and +107 is sufficient for cell-type-specific expression. Assays for DNA-protein interactions coupled with extensive mutagenesis identified three factors, whose binding correlated with the function of two positive and one negative regulatory element. One of these factors, termed M-box-binding factor 1, binds to an 11-bp motif, the M box, which acts as a positive regulatory element both in TRP-1-expressing and -nonexpressing cell lines, despite being entirely conserved between the melanocyte-specific tyrosinase and TRP-1 promoters. The possible mechanisms underlying melanocyte-specific gene expression are discussed.


1992 ◽  
Vol 12 (8) ◽  
pp. 3653-3662 ◽  
Author(s):  
P Lowings ◽  
U Yavuzer ◽  
C R Goding

Melanocytes are specialized cells residing in the hair follicles, the eye, and the basal layer of the human epidermis whose primary function is the production of the pigment melanin, giving rise to skin, hair, and eye color. Melanogenesis, a process unique to melanocytes that involves the processing of tyrosine by a number of melanocyte-specific enzymes, including tyrosinase and tyrosinase-related protein 1 (TRP-1), occurs only after differentiation from the melanocyte precursor, the melanoblast. In humans, melanogenesis is inducible by UV irradiation, with melanin being transferred from the melanocyte in the epidermis to the surrounding keratinocytes as protection from UV-induced damage. Excessive exposure to UV, however, is the primary cause of malignant melanoma, an increasingly common and highly aggressive disease. As an initial approach to understanding the regulation of melanocyte differentiation and melanocyte-specific transcription, we have isolated the gene encoding TRP-1 and examined the cis- and trans-acting factors required for cell-type-specific expression. We find that the TRP-1 promoter comprises both positive and negative regulatory elements which confer efficient expression in a TRP-1-expressing, pigmented melanoma cell line but not in NIH 3T3 or JEG3 cells and that a minimal promoter extending between -44 and +107 is sufficient for cell-type-specific expression. Assays for DNA-protein interactions coupled with extensive mutagenesis identified three factors, whose binding correlated with the function of two positive and one negative regulatory element. One of these factors, termed M-box-binding factor 1, binds to an 11-bp motif, the M box, which acts as a positive regulatory element both in TRP-1-expressing and -nonexpressing cell lines, despite being entirely conserved between the melanocyte-specific tyrosinase and TRP-1 promoters. The possible mechanisms underlying melanocyte-specific gene expression are discussed.


2021 ◽  
Vol 7 (20) ◽  
pp. eabf1444
Author(s):  
James D. Hocker ◽  
Olivier B. Poirion ◽  
Fugui Zhu ◽  
Justin Buchanan ◽  
Kai Zhang ◽  
...  

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease–associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type–specific gene regulation in human hearts during health and disease.


Author(s):  
Vahid H. Gazestani ◽  
Austin WT Chiang ◽  
Eric Courchesne ◽  
Nathan E. Lewis

ABSTRACTNumerous genes are associated with autism spectrum disorder (ASD); however, it remains unclear how most ASD risk genes influence neurodevelopment and result in similar traits. Recent genetic models of complex traits suggest non-tissue-specific genes converge on core disease genes; so we analyzed ASD genetics in this context. We found ASD risk genes partition cleanly into broadly-expressed and brain-specific genes. The two groups show sequential roles during neurodevelopment with broadly-expressed genes modulating chromatin remodeling, proliferation, and cell fate, while brain-specific risk genes are involved in neural maturation and synapse functioning. Broadly-expressed risk genes converge onto brain-specific risk genes and core neurodevelopmental genes through regulatory networks including PI3K/AKT, RAS/ERK, and WNT/β-catenin signaling pathways. Broadly-expressed and brain-specific risk genes show unique properties, wherein the broadly-expressed risk gene network is expressed prenatally and conserved in non-neuronal cells like microglia. However, the brain-specific gene network expression is limited to excitatory and inhibitory neurons, spanning prenatal to adulthood. Furthermore, the two groups are linked differently to comorbidities associated with ASD. Collectively, we describe here the organization of the genetic architecture of ASD as a hierarchy of broadly-expressed and brain-specific genes that disrupt successive stages of core neurodevelopmental processes.


Sign in / Sign up

Export Citation Format

Share Document