Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Novel TCF4:TCF12 heterodimer inhibits glioblastoma growth

TWIST1 (TW) is a pro-oncogenic basic helix-loop-helix (bHLH) transcription factor and promotes the hallmark features of malignancy (e.g., cell invasion, cancer cell stemness, and treatment resistance), which contribute to poor prognoses of glioblastoma (GBM). We previously reported that specific TW dimerization motifs regulate unique cellular phenotypes in GBM. For example, the TW:E12 heterodimer increases periostin (POSTN) expression and promotes cell invasion. TW dimer-specific transcriptional regulation requires binding to the regulatory E-box consensus sequences, but alternative bHLH dimers that balance TW dimer activity in regulating pro-oncogenic TW target genes are unknown. We leveraged the ENCODE DNase I hypersensitivity data to identify E-box sites and tethered TW:E12 and TW:TW proteins to validate dimer binding to E-boxes in vitro. Subsequently, TW knockdown revealed a novel TCF4:TCF12 bHLH dimer occupying the same TW E-box site that, when expressed as a tethered TCF4:TCF12 dimer, markedly repressed POSTN expression and extended animal survival. These observations support TCF4:TCF12 as a novel dimer with tumor-suppressor activity in GBM that functions in part through displacement of and/or competitive inhibition of pro-oncogenic TW dimers at E-box sites.

Cochlear organoids reveal transcriptional programs of postnatal hair cell differentiation from supporting cells

We explore the changes in chromatin accessibility and transcriptional programs for cochlear hair cell differentiation from postmitotic supporting cells using organoids from postnatal cochlea. The organoids contain cells with transcriptional signatures of differentiating vestibular and cochlear hair cells. Construction of trajectories identifies Lgr5+ cells as progenitors for hair cells, and the genomic data reveal gene regulatory networks leading to hair cells. We validate these networks, demonstrating dynamic changes both in expression and predicted binding sites of transcription factors (TFs) during organoid differentiation. We identify known regulators of hair cell development, Atoh1, Pou4f3, and Gfi1, and the analysis predicts the regulatory factors Tcf4, an E-protein and heterodimerization partner of Atoh1, and Ddit3, a CCAAT/enhancer-binding protein (C/EBP) that represses Hes1 and activates transcription of Wnt-signaling-related genes. Deciphering the signals for hair cell regeneration from mammalian cochlear supporting cells reveals candidates for hair cell (HC) regeneration, which is limited in the adult.

Screening for functional regulatory variants in open chromatin using GenIE-ATAC

Understanding the effects of genetic variation in gene regulatory elements is crucial to interpreting genome function. This is particularly pertinent for the hundreds of thousands of disease-associated variants identified by GWAS, which frequently sit within gene regulatory elements but whose functional effects are often unknown. Current methods are limited in their scalability and ability to assay regulatory variants in their endogenous context, independently of other tightly linked variants. Here, we present a new medium-throughput screening system: genome engineering based interrogation of enhancers assay for transposase accessible chromatin (GenIE-ATAC), that measures the effect of individual variants on chromatin accessibility in their endogenous genomic and chromatin context. We employ this assay to screen for the effects of regulatory variants in human induced pluripotent stem cells, validating a subset of causal variants, and extend our software package (rgenie) to analyse these new data. We demonstrate that this methodology can be used to understand the impact of defined deletions and point mutations within transcription factor binding sites. We thus establish GenIE-ATAC as a method to screen for the effect of gene regulatory element variation, allowing identification and prioritisation of causal variants from GWAS for functional follow-up and understanding the mechanisms of regulatory element function.

Three-dimensional genome rewiring in loci with human accelerated regions

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Machine learning dissection of human accelerated regions in primate neurodevelopment

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Identifying enhancer properties associated with genetic risk for complex traits using regulome-wide association studies

Genetic risk for complex traits is strongly enriched in non-coding genomic regions involved in gene regulation, especially enhancers. However, we lack adequate tools to connect the characteristics of these disruptions to genetic risk. Here, we propose RWAS (Regulome Wide Association Study), a new application of the MAGMA software package to identify the characteristics of enhancers that contribute to genetic risk for disease. RWAS involves three steps: (i) assign genotyped SNPs to cell type- or tissue-specific regulatory features (e.g., enhancers); (ii) test associations of each regulatory feature with a trait of interest for which genome-wide association study (GWAS) summary statistics are available; (iii) perform enhancer-set enrichment analyses to identify quantitative or categorical features of regulatory elements that are associated with the trait. These steps are implemented as a novel application of MAGMA, a tool originally developed for gene-based GWAS analyses. Applying RWAS to interrogate genetic risk for schizophrenia, we discovered a class of risk-associated AT-rich enhancers that are active in the developing brain and harbor binding sites for multiple transcription factors with neurodevelopmental functions. RWAS utilizes open-source software, and we provide a comprehensive collection of annotations for tissue-specific enhancer locations and features, including their evolutionary conservation, AT content, and co-localization with binding sites for hundreds of TFs. RWAS will enable researchers to characterize properties of regulatory elements associated with any trait of interest for which GWAS summary statistics are available.

Enhancers are regulatory regions that influence gene expression via the binding of transcription factors. Risk for many heritable diseases is enriched in regulatory regions, including enhancers. In this study, we introduce a novel application of the MAGMA software tool that enables testing for associations between enhancer attributes and risk, and we use this method to determine the enhancer characteristics that are associated with risk for schizophrenia. We found that enhancers associated with schizophrenia risk are both evolutionarily conserved and in physical contact with mutation-intolerant genes, many of which have neurodevelopmental functions. Risk-associated enhancers are also AT-rich and contain binding sites for neurodevelopmental transcription factors.

Core transcription programs controlling injury-induced neurodegeneration of retinal ganglion cells

Regulatory programs governing neuronal death and axon regeneration in neurodegenerative diseases remain poorly understood. In adult mice, optic nerve crush (ONC) injury by severing retinal ganglion cell (RGC) axons results in massive RGC death and regenerative failure. We performed an in vivo CRISPR-Cas9-based genome-wide screen of 1,893 transcription factors (TFs) to seek repressors of RGC survival and axon regeneration following ONC. In parallel, we profiled the epigenetic and transcriptional landscapes of injured RGCs by ATAC-seq and RNA-seq to identify injury-responsive TFs and their targets. These analyses converged on four TFs as critical survival regulators, of which ATF3/CHOP preferentially regulate pathways activated by cytokines and innate immunity and ATF4/C/EBPγ regulate pathways engaged by intrinsic neuronal stressors. Manipulation of these TFs protects RGCs in a glaucoma model. Our results reveal core transcription programs that transform an initial axonal insult into a degenerative process and suggest novel strategies for treating neurodegenerative diseases.

Advances and Trends in Omics Technology Development

The human history has witnessed the rapid development of technologies such as high-throughput sequencing and mass spectrometry that led to the concept of “omics” and methodological advancement in systematically interrogating a cellular system. Yet, the ever-growing types of molecules and regulatory mechanisms being discovered have been persistently transforming our understandings on the cellular machinery. This renders cell omics seemingly, like the universe, expand with no limit and our goal toward the complete harness of the cellular system merely impossible. Therefore, it is imperative to review what has been done and is being done to predict what can be done toward the translation of omics information to disease control with minimal cell perturbation. With a focus on the “four big omics,” i.e., genomics, transcriptomics, proteomics, metabolomics, we delineate hierarchies of these omics together with their epiomics and interactomics, and review technologies developed for interrogation. We predict, among others, redoxomics as an emerging omics layer that views cell decision toward the physiological or pathological state as a fine-tuned redox balance.

Posterior cingulate cortex reveals an expression profile of resilience in cognitively intact elders

The posterior cingulate cortex, a key hub of the default mode network, underlies autobiographical memory retrieval and displays hypometabolic changes early in Alzheimer disease. To obtain an unbiased understanding of the molecular pathobiology of the aged posterior cingulate cortex, we performed RNA sequencing (RNA-seq) on tissue obtained from 26 participants of the Rush Religious Orders Study (11 males/15 females; aged 76-96 years) with a pre-mortem clinical diagnosis of no cognitive impairment and post-mortem neurofibrillary tangle Braak Stages I/II, III, and IV. Transcriptomic data were gathered using next-generation sequencing of RNA extracted from posterior cingulate cortex generating an average of 60 million paired reads per subject. Normalized expression of RNA-seq data was calculated using a global gene annotation and a microRNA profile. Differential expression (DESeq2, edgeR) using Braak staging as the comparison structure isolated genes for dimensional scaling, associative network building and functional clustering. Curated genes were correlated with the Mini-Mental State Examination and semantic, working and episodic memory, visuospatial ability, and a composite Global Cognitive Score. Regulatory mechanisms were determined by co-expression networks with microRNAs and an overlap of transcription factor binding sites. Analysis revealed 750 genes and 12 microRNAs significantly differentially expressed between Braak Stages I/II and III/IV and an associated six groups of transcription factor binding sites. Inputting significantly different gene/network data into a functional annotation clustering model revealed elevated presynaptic, postsynaptic and ATP-related expression in Braak Stages III and IV compared with Stages I/II, suggesting these pathways are integral for cognitive resilience seen in unimpaired elderly subjects. Principal component analysis and Kruskal-Wallis testing did not associate Braak stage with cognitive function. However, Spearman correlations between genes and cognitive test scores followed by network analysis revealed upregulation of classes of synaptic genes positively associated with performance on the visuospatial perceptual orientation domain. Upregulation of key synaptic genes suggests a role for these transcripts and associated synaptic pathways in cognitive resilience seen in elders despite Alzheimer disease pathology and dementia.

Systems modeling of metabolic dysregulation in neurodegenerative diseases

Neurodegenerative diseases (NDDs) encompass a wide range of conditions that arise owing to progressive degeneration and the ultimate loss of nerve cells in the brain and peripheral nervous system. NDDs such as Alzheimer's, Parkinson's, and Huntington's diseases negatively impact both length and quality of life, due to lack of effective disease-modifying treatments. Herein, we review the use of genome-scale metabolic models, network-based approaches, and integration with multiomics data to identify key biological processes that characterize NDDs. We describe powerful systems biology approaches for modeling NDD pathophysiology by leveraging in silico models that are informed by patient-derived multiomics data. These approaches can enable mechanistic insights into NDD-specific metabolic dysregulations that can be leveraged to identify potential metabolic markers of disease and predisease states.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Modeling gene regulatory networks using neural network architectures

Gene regulatory networks (GRNs) encode the complex molecular interactions that govern cell identity. Here we propose DeepSEM, a deep generative model that can jointly infer GRNs and biologically meaningful representation of single-cell RNA sequencing (scRNA-seq) data. In particular, we developed a neural network version of the structural equation model (SEM) to explicitly model the regulatory relationships among genes. Benchmark results show that DeepSEM achieves comparable or better performance on a variety of single-cell computational tasks, such as GRN inference, scRNA-seq data visualization, clustering and simulation, compared with the state-of-the-art methods. In addition, the gene regulations predicted by DeepSEM on cell-type marker genes in the mouse cortex can be validated by epigenetic data, which further demonstrates the accuracy and efficiency of our method. DeepSEM can provide a useful and powerful tool to analyze scRNA-seq data and infer a GRN.

ZNF718, HOXA4, and ZFP57 are differentially methylated in periodontitis in comparison with periodontal health: Epigenome-wide DNA methylation pilot study

OBJECTIVE

To investigate the differences in the epigenomic patterns of DNA methylation in peripheral leukocytes between patients with periodontitis and gingivally healthy controls evaluating its functional meaning by functional enrichment analysis.

BACKGROUND

The DNA methylation profiling of peripheral leukocytes as immune-related tissue potentially relevant as a source of biomarkers between periodontitis patients and gingivally healthy subjects has not been investigated.

METHODS

A DNA methylation epigenome-wide study of peripheral leukocytes was conducted using the Illumina MethylationEPIC platform in sixteen subjects, eight diagnosed with periodontitis patients and eight age-matched and sex-matched periodontally healthy controls. A trained periodontist performed the clinical evaluation. Global DNA methylation was estimated using methylation-sensitive high-resolution melting in LINE1. Routine cell count cytometry and metabolic laboratory tests were also performed. The analysis of differentially methylated positions (DMPs) and differentially methylated regions (DMRs) was made using R/Bioconductor environment considering leukocyte populations assessed in both routine cell counts and using the FlowSorted.Blood.EPIC package. Finally, a DMP and DMR intersection analysis was performed. Functional enrichment analysis was carried out with the differentially methylated genes found in DMP.

RESULTS

DMP analysis identified 81 differentially hypermethylated genes and 21 differentially hypomethylated genes. Importantly, the intersection analysis showed that zinc finger protein 718 (ZNF718) and homeobox A4 (HOXA4) were differentially hypermethylated and zinc finger protein 57 (ZFP57) was differentially hypomethylated in periodontitis. The functional enrichment analysis found clearly immune-related ontologies such as "detection of bacterium" and "antigen processing and presentation."

CONCLUSION

The results of this study propose three new periodontitis-related genes: ZNF718, HOXA4, and ZFP57 but also evidence the suitability and relevance of studying leukocytes' DNA methylome for biological interpretation of systemic immune-related epigenetic patterns in periodontitis.

Omics in Systems Biology: Current Progress and Future Outlook

Biological research has undergone tremendous changes over the past three decades. Research used to almost exclusively focus on a single aspect of a single molecule per experiment. Modern technologies have enabled thousands of molecules to be simultaneously analyzed and the way that these molecules influence each other to be discerned. The change is so dramatic that it has given rise to a whole new descriptive suffix (i.e., omics) to describe these fields of study. While genomics was arguably the initial driver of this new trend, it quickly spread to other biological entities resulting in the creation of transcriptomics, proteomics, metabolomics, etc. The development of these "big four omics" created a wave of other omic fields, such as epigenomics, glycomics, lipidomics, microbiomics, and even foodomics; all with the purpose of comprehensively studying all the molecular entities or processes within their respective domain. The large number of omic fields that are invented even led to the term "panomics" as a way to classify them all under one category. Ultimately, all of these omic fields are setting the foundation for developing systems biology; in which the focus will be on determining the complex interactions that occur within biological systems.

Metabolic Network Analysis Reveals Altered Bile Acid Synthesis and Metabolism in Alzheimer's Disease

Increasing evidence suggests Alzheimer's disease (AD) pathophysiology is influenced by primary and secondary bile acids, the end product of cholesterol metabolism. We analyze 2,114 post-mortem brain transcriptomes and identify genes in the alternative bile acid synthesis pathway to be expressed in the brain. A targeted metabolomic analysis of primary and secondary bile acids measured from post-mortem brain samples of 111 individuals supports these results. Our metabolic network analysis suggests that taurine transport, bile acid synthesis, and cholesterol metabolism differ in AD and cognitively normal individuals. We also identify putative transcription factors regulating metabolic genes and influencing altered metabolism in AD. Intriguingly, some bile acids measured in brain tissue cannot be explained by the presence of enzymes responsible for their synthesis, suggesting that they may originate from the gut microbiome and are transported to the brain. These findings motivate further research into bile acid metabolism in AD to elucidate their possible connection to cognitive decline.