The functional consequences of variation in transcription factor binding

Regulation of meristem and hormone function revealed through analysis of directly-regulated SHOOT MERISTEMLESS target genes

The Arabidopsis Knotted1-like homeobox (KNOX) gene SHOOT MERISTEMLESS (STM) encodes a homeodomain transcription factor that operates as a central component of the gene regulatory network (GRN) controlling shoot apical meristem formation and maintenance. It regulates the expression of target genes that include transcriptional regulators associated with meristem function, particularly those involved in pluripotency and cellular differentiation, as well as genes involved in hormone metabolism and signaling. Previous studies have identified KNOX-regulated genes and their associated cis-regulatory elements in several plant species. However, little is known about STM-DNA interactions in the regulatory regions of target genes in Arabidopsis. Here, we identify and map STM binding sites in the Arabidopsis genome using global ChIP-seq analysis to reveal potential directly-regulated STM target genes. We show that in the majority of target loci, STM binds within 1 kb upstream of the TSS, with other loci showing STM binding at more distal enhancer sites, and we reveal enrichment of DNA motifs containing a TGAC and/or TGAT core in STM-bound target gene cis-regulatory elements. We further demonstrate that many STM-bound genes are transcriptionally responsive to altered levels of STM activity, and show that among these, transcriptional regulators with key roles in meristem and hormone function are highly represented. Finally, we use a subset of these target genes to perform Bayesian network analysis to infer gene regulatory associations and to construct a refined GRN for STM-mediated control of meristem function.

Binding profiles for 961 Drosophila and C. elegans transcription factors reveal tissue-specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed "metapeaks," that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent-daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF-target relationships.

Identifying Reproducible Transcription Regulator Coexpression Patterns with Single Cell Transcriptomics

The proliferation of single cell transcriptomics has potentiated our ability to unveil patterns that reflect dynamic cellular processes, rather than cell type compositional effects that emerge from bulk tissue samples. In this study, we leverage a broad collection of single cell RNA-seq data to identify the gene partners whose expression is most coordinated with each human and mouse transcription regulator (TR). We assembled 120 human and 103 mouse scRNA-seq datasets from the literature (>28 million cells), constructing a single cell coexpression network for each. We aimed to understand the consistency of TR coexpression profiles across a broad sampling of biological contexts, rather than examine the preservation of context-specific signals. Our workflow therefore explicitly prioritizes the patterns that are most reproducible across cell types. Towards this goal, we characterize the similarity of each TR's coexpression within and across species. We create single cell coexpression rankings for each TR, demonstrating that this aggregated information recovers literature curated targets on par with ChIP-seq data. We then combine the coexpression and ChIP-seq information to identify candidate regulatory interactions supported across methods and species. Finally, we highlight interactions for the important neural TR ASCL1 to demonstrate how our compiled information can be adopted for community use.

RummaGEO: Automatic mining of human and mouse gene sets from GEO

The Gene Expression Omnibus (GEO) has millions of samples from thousands of studies. While users of GEO can search the metadata describing studies, there is a need for methods to search GEO at the data level. RummaGEO is a gene expression signature search engine for human and mouse RNA sequencing perturbation studies extracted from GEO. To develop RummaGEO, we automatically identified groups of samples and computed differential expressions to extract gene sets from each study. The contents of RummaGEO are served for gene set, PubMed, and metadata search. Next, we analyzed the contents of RummaGEO to identify patterns and perform global analyses. Overall, RummaGEO provides a resource that is enabling users to identify relevant GEO studies based on their own gene expression results. Users of RummaGEO can incorporate RummaGEO into their analysis workflows for integrative analyses and hypothesis generation.

COFFEE: consensus single cell-type specific inference for gene regulatory networks

The inference of gene regulatory networks (GRNs) is crucial to understanding the regulatory mechanisms that govern biological processes. GRNs may be represented as edges in a graph, and hence, it have been inferred computationally for scRNA-seq data. A wisdom of crowds approach to integrate edges from several GRNs to create one composite GRN has demonstrated improved performance when compared with individual algorithm implementations on bulk RNA-seq and microarray data. In an effort to extend this approach to scRNA-seq data, we present COFFEE (COnsensus single cell-type speciFic inFerence for gEnE regulatory networks), a Borda voting-based consensus algorithm that integrates information from 10 established GRN inference methods. We conclude that COFFEE has improved performance across synthetic, curated, and experimental datasets when compared with baseline methods. Additionally, we show that a modified version of COFFEE can be leveraged to improve performance on newer cell-type specific GRN inference methods. Overall, our results demonstrate that consensus-based methods with pertinent modifications continue to be valuable for GRN inference at the single cell level. While COFFEE is benchmarked on 10 algorithms, it is a flexible strategy that can incorporate any set of GRN inference algorithms according to user preference. A Python implementation of COFFEE may be found on GitHub: https://github.com/lodimk2/coffee.

PoSnf1 affects cellulose utilization through interaction with cellobiose transporter in Pleurotus ostreatus

Pleurotus ostreatus is one of the most cultivated edible fungi worldwide, but its lignocellulose utilization efficiency is relatively low (<50 %), which eventually affects the biological efficiency of P. ostreatus. Improving cellulase production and activity will contribute to enhancing the lignocellulose-degrading capacity of P. ostreatus. AMP-activated/Snf1 protein kinase plays important roles in regulating carbon and energy metabolism. The Snf1 homolog (PoSnf1) in P. ostreatus was obtained and analyzed using bioinformatics. The cellulose response of PoSnf1, the effect of the phosphorylation level of PoSnf1 on the expression of cellulose degradation-related genes, the putative proteins that interact with the phosphorylated PoSnf1 (P-PoSnf1), the cellobiose transport function of two sugar transporters (STP1 and STP2), and the interactions between PoSnf1 and STP1/STP2 were studied in this research. We found that cellulose treatment improved the phosphorylation level of PoSnf1, which further affected cellulase activity and the expression of most cellulose degradation-related genes. A total of 1, 024 proteins putatively interacting with P-PoSnf1 were identified, and they were enriched mainly in the substances transport and metabolism. Most of the putative cellulose degradation-related protein-coding genes could respond to cellulose. Among the P-PoSnf1-interacting proteins, the functions of two sugar transporters (STP1 and STP2) were further studied, and the results showed that both could transport cellobiose and were indirectly regulated by P-PoSnf1, and that STP2 could directly interact with PoSnf1. The results of this study indicated that PoSnf1 plays an important role in regulating the expression of cellulose degradation genes possibly by affecting cellobiose transport.

Position-dependent function of human sequence-specific transcription factors

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

White-opaque switching in Candida albicans: cell biology, regulation, and function

SUMMARYCandida albicans remains a major fungal pathogen colonizing humans and opportunistically invading tissue when conditions are predisposing. Part of the success of C. albicans was attributed to its capacity to form hyphae that facilitate tissue invasion. However, in 1987, a second developmental program was discovered, the "white-opaque transition," a high-frequency reversible switching system that impacted most aspects of the physiology, cell architecture, virulence, and gene expression of C. albicans. For the 15 years following the discovery of white-opaque switching, its role in the biology of C. albicans remained elusive. Then in 2002, it was discovered that in order to mate, C. albicans had to switch from white to opaque, a unique step in a yeast mating program. In 2006, three laboratories simultaneously identified a putative master switch gene, which led to a major quest to elucidate the underlying mechanisms that regulate white-opaque switching. Here, the evolving discoveries related to this complicated phenotypic transition are reviewed in a quasi-chronological order not only to provide a historical perspective but also to highlight several unique characteristics of white-opaque switching, which are fascinating and may be important to the life history and virulence of this persistent pathogen. Many of these characteristics have not been fully investigated, in many cases, leaving intriguing questions unresolved. Some of these include the function of unique channeled pimples on the opaque cell wall, the capacity to form opaque cells in the absence of the master switch gene WOR1, the formation of separate "pathogenic" and "sexual" biofilms, and the possibility that a significant portion of natural strains colonizing the lower gastrointestinal tract may be in the opaque phase. This review addresses many of these characteristics with the intent of engendering interest in resolving questions that remain unanswered.

RummaGEO: Automatic Mining of Human and Mouse Gene Sets from GEO

The Gene Expression Omnibus (GEO) is a major open biomedical research repository for transcriptomics and other omics datasets. It currently contains millions of gene expression samples from tens of thousands of studies collected by many biomedical research laboratories from around the world. While users of the GEO repository can search the metadata describing studies for locating relevant datasets, there are currently no methods or resources that facilitate global search of GEO at the data level. To address this shortcoming, we developed RummaGEO, a webserver application that enables gene expression signature search of a large collection of human and mouse RNA-seq studies deposited into GEO. To develop the search engine, we performed offline automatic identification of sample conditions from the uniformly aligned GEO studies available from ARCHS4. We then computed differential expression signatures to extract gene sets from these studies. In total, RummaGEO currently contains 135,264 human and 158,062 mouse gene sets extracted from 23,395 GEO studies. Next, we analyzed the contents of the RummaGEO database to identify statistical patterns and perform various global analyses. The contents of the RummaGEO database are provided as a web-server search engine with signature search, PubMed search, and metadata search functionalities. Overall, RummaGEO provides an unprecedented resource for the biomedical research community enabling hypothesis generation for many future studies. The RummaGEO search engine is available from: https://rummageo.com/.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Single-cell multi-omics analysis identifies context-specific gene regulatory gates and mechanisms

There is a growing interest in inferring context specific gene regulatory networks from single-cell RNA sequencing (scRNA-seq) data. This involves identifying the regulatory relationships between transcription factors (TFs) and genes in individual cells, and then characterizing these relationships at the level of specific cell types or cell states. In this study, we introduce scGATE (single-cell gene regulatory gate) as a novel computational tool for inferring TF-gene interaction networks and reconstructing Boolean logic gates involving regulatory TFs using scRNA-seq data. In contrast to current Boolean models, scGATE eliminates the need for individual formulations and likelihood calculations for each Boolean rule (e.g. AND, OR, XOR). By employing a Bayesian framework, scGATE infers the Boolean rule after fitting the model to the data, resulting in significant reductions in time-complexities for logic-based studies. We have applied assay for transposase-accessible chromatin with sequencing (scATAC-seq) data and TF DNA binding motifs to filter out non-relevant TFs in gene regulations. By integrating single-cell clustering with these external cues, scGATE is able to infer context specific networks. The performance of scGATE is evaluated using synthetic and real single-cell multi-omics data from mouse tissues and human blood, demonstrating its superiority over existing tools for reconstructing TF-gene networks. Additionally, scGATE provides a flexible framework for understanding the complex combinatorial and cooperative relationships among TFs regulating target genes by inferring Boolean logic gates among them.

Single-molecule chromatin configurations link transcription factor binding to expression in human cells

The binding of multiple transcription factors (TFs) to genomic enhancers activates gene expression in mammalian cells. However, the molecular details that link enhancer sequence to TF binding, promoter state, and gene expression levels remain opaque. We applied single-molecule footprinting (SMF) to measure the simultaneous occupancy of TFs, nucleosomes, and components of the transcription machinery on engineered enhancer/promoter constructs with variable numbers of TF binding sites for both a synthetic and an endogenous TF. We find that activation domains enhance a TF's capacity to compete with nucleosomes for binding to DNA in a BAF-dependent manner, TF binding on nucleosome-free DNA is consistent with independent binding between TFs, and average TF occupancy linearly contributes to promoter activation rates. We also decompose TF strength into separable binding and activation terms, which can be tuned and perturbed independently. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the binding microstates observed at the enhancer and subsequent time-dependent gene expression. This work provides a template for quantitative dissection of distinct contributors to gene activation, including the activity of chromatin remodelers, TF activation domains, chromatin acetylation, TF concentration, TF binding affinity, and TF binding site configuration.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

COFFEE: Consensus Single Cell-Type Specific Inference for Gene Regulatory Networks

The inference of gene regulatory networks (GRNs) is crucial to understanding the regulatory mechanisms that govern biological processes. GRNs may be represented as edges in a graph, and hence have been inferred computationally for scRNA-seq data. A wisdom of crowds approach to integrate edges from several GRNs to create one composite GRN has demonstrated improved performance when compared to individual algorithm implementations on bulk RNA-seq and microarray data. In an effort to extend this approach to scRNA-seq data, we present COFFEE (COnsensus single cell-type speciFic inFerence for gEnE regulatory networks), a Borda voting based consensus algorithm that integrates information from 10 established GRN inference methods. We conclude that COFFEE has improved performance across synthetic, curated and experimental datasets when compared to baseline methods. Additionally, we show that a modified version of COFFEE can be leveraged to improve performance on newer cell-type specific GRN inference methods. Overall, our results demonstrate that consensus based methods with pertinent modifications continue to be valuable for GRN inference at the single cell level.

Atlas of nascent RNA transcripts reveals enhancer to gene linkages

Gene transcription is controlled and modulated by regulatory regions, including enhancers and promoters. These regions are abundant in unstable, non-coding bidirectional transcription. Using nascent RNA transcription data across hundreds of human samples, we identified over 800,000 regions containing bidirectional transcription. We then identify highly correlated transcription between bidirectional and gene regions. The identified correlated pairs, a bidirectional region and a gene, are enriched for disease associated SNPs and often supported by independent 3D data. We present these resources as an SQL database which serves as a resource for future studies into gene regulation, enhancer associated RNAs, and transcription factors.

Computational tools for inferring transcription factor activity

Transcription factors (TFs) are essential players in orchestrating the regulatory landscape in cells. Still, their exact modes of action and dependencies on other regulatory aspects remain elusive. Since TFs act cell type-specific and each TF has its own characteristics, untangling their regulatory interactions from an experimental point of view is laborious and convoluted. Thus, there is an ongoing development of computational tools that estimate transcription factor activity (TFA) from a variety of data modalities, either based on a mapping of TFs to their putative target genes or in a genome-wide, gene-unspecific fashion. These tools can help to gain insights into TF regulation and to prioritize candidates for experimental validation. We want to give an overview of available computational tools that estimate TFA, illustrate examples of their application, debate common result validation strategies, and discuss assumptions and concomitant limitations.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

Gene signature reveals decreased SOX10-dependent transcripts in malignant cells from immune checkpoint inhibitor-resistant cutaneous melanomas

Evidence is mounting for cross-resistance between immune checkpoint and targeted kinase inhibitor therapies in cutaneous melanoma patients. Since the loss of the transcription factor, SOX10, causes tolerance to MAPK pathway inhibitors, we used bioinformatic techniques to determine if reduced SOX10 expression/activity is associated with immune checkpoint inhibitor resistance. We integrated SOX10 ChIP-seq, knockout RNA-seq, and knockdown ATAC-seq data from melanoma cell models to develop a robust SOX10 gene signature. We used computational methods to validate this signature as a measure of SOX10-dependent activity in independent single-cell and bulk RNA-seq SOX10 knockdown, cell line panel, and MAPK inhibitor drug-resistant datasets. Evaluation of patient single-cell RNA-seq data revealed lower levels of SOX10-dependent transcripts in immune checkpoint inhibitor-resistant tumors. Our results suggest that SOX10-deficient melanoma cells are associated with cross-resistance between targeted and immune checkpoint inhibitors and highlight the need to identify therapeutic strategies that target this subpopulation.

Identification of mammalian transcription factors that bind to inaccessible chromatin

Transcription factors (TFs) are proteins that affect gene expression by binding to regulatory regions of DNA in a sequence specific manner. The binding of TFs to DNA is controlled by many factors, including the DNA sequence, concentration of TF, chromatin accessibility and co-factors. Here, we systematically investigated the binding mechanism of hundreds of TFs by analysing ChIP-seq data with our explainable statistical model, ChIPanalyser. This tool uses as inputs the DNA sequence binding motif; the capacity to distinguish between strong and weak binding sites; the concentration of TF; and chromatin accessibility. We found that approximately one third of TFs are predicted to bind the genome in a DNA accessibility independent fashion, which includes TFs that can open the chromatin, their co-factors and TFs with similar motifs. Our model predicted this to be the case when the TF binds to its strongest binding regions in the genome, and only a small number of TFs have the capacity to bind dense chromatin at their weakest binding regions, such as CTCF, USF2 and CEBPB. Our study demonstrated that the binding of hundreds of human and mouse TFs is predicted by ChIPanalyser with high accuracy and showed that many TFs can bind dense chromatin.

Gene regulatory network inference using mixed-norms regularized multivariate model with covariance selection

Despite extensive research efforts, reconstruction of gene regulatory networks (GRNs) from transcriptomics data remains a pressing challenge in systems biology. While non-linear approaches for reconstruction of GRNs show improved performance over simpler alternatives, we do not yet have understanding if joint modelling of multiple target genes may improve performance, even under linearity assumptions. To address this problem, we propose two novel approaches that cast the GRN reconstruction problem as a blend between regularized multivariate regression and graphical models that combine the L2,1-norm with classical regularization techniques. We used data and networks from the DREAM5 challenge to show that the proposed models provide consistently good performance in comparison to contenders whose performance varies with data sets from simulation and experiments from model unicellular organisms Escherichia coli and Saccharomyces cerevisiae. Since the models' formulation facilitates the prediction of master regulators, we also used the resulting findings to identify master regulators over all data sets as well as their plasticity across different environments. Our results demonstrate that the identified master regulators are in line with experimental evidence from the model bacterium E. coli. Together, our study demonstrates that simultaneous modelling of several target genes results in improved inference of GRNs and can be used as an alternative in different applications.

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets

Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

Human genetics seen through an evolutionary lens

Most DNA bases crucial for species perpetuation are marked by a dearth of sequence change among species related over long evolutionary time. Recently, Christmas et al.¹ and Sullivan et al.² cast light on human DNA and its variants through comparison with 239 other mammalian species' genomes.

Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets

Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Broad compatibility between yeast UAS elements and core promoters and identification of promoter elements that determine cofactor specificity

Three classes of yeast protein-coding genes are distinguished by their dependence on the transcription cofactors TFIID, SAGA, and Mediator (MED) Tail, but whether this dependence is determined by the core promoter, upstream activating sequences (UASs), or other gene features is unclear. Also unclear is whether UASs can broadly activate transcription from the different promoter classes. Here, we measure transcription and cofactor specificity for thousands of UAS-core promoter combinations and find that most UASs broadly activate promoters regardless of regulatory class, while few display strong promoter specificity. However, matching UASs and promoters from the same gene class is generally important for optimal expression. We find that sensitivity to rapid depletion of MED Tail or SAGA is dependent on the identity of both UAS and core promoter, while dependence on TFIID localizes to only the promoter. Finally, our results suggest the role of TATA and TATA-like promoter sequences in MED Tail function.

Escape from oncogene-induced senescence is controlled by POU2F2 and memorized by chromatin scars

Although oncogene-induced senescence (OIS) is a potent tumor-suppressor mechanism, recent studies revealed that cells could escape from OIS with features of transformed cells. However, the mechanisms that promote OIS escape remain unclear, and evidence of post-senescent cells in human cancers is missing. Here, we unravel the regulatory mechanisms underlying OIS escape using dynamic multidimensional profiling. We demonstrate a critical role for AP1 and POU2F2 transcription factors in escape from OIS and identify senescence-associated chromatin scars (SACSs) as an epigenetic memory of OIS detectable during colorectal cancer progression. POU2F2 levels are already elevated in precancerous lesions and as cells escape from OIS, and its expression and binding activity to cis-regulatory elements are associated with decreased patient survival. Our results support a model in which POU2F2 exploits a precoded enhancer landscape necessary for senescence escape and reveal POU2F2 and SACS gene signatures as valuable biomarkers with diagnostic and prognostic potential.

Computational drug repurposing by exploiting large-scale gene expression data: Strategy, methods and applications

De novo drug development is an extremely complex, time-consuming and costly task. Urgent needs for therapies of various diseases have greatly accelerated searches for more effective drug development methods. Luckily, drug repurposing provides a new and effective perspective on disease treatment. Rapidly increased large-scale transcriptome data paints a detailed prospect of gene expression during disease onset and thus has received wide attention in the field of computational drug repurposing. However, how to efficiently mine transcriptome data and identify new indications for old drugs remains a critical challenge. This review discussed the irreplaceable role of transcriptome data in computational drug repurposing and summarized some representative databases, tools and strategies. More importantly, it proposed a practical guideline through establishing the correspondence between three gene expression data types and five strategies, which would facilitate researchers to adopt appropriate strategies to deeply mine large-scale transcriptome data and discover more effective therapies.

NetProphet 3: a machine learning framework for transcription factor network mapping and multi-omics integration

MOTIVATION

Many methods have been proposed for mapping the targets of transcription factors (TFs) from gene expression data. It is known that combining outputs from multiple methods can improve performance. To date, outputs have been combined by using either simplistic formulae, such as geometric mean, or carefully hand-tuned formulae that may not generalize well to new inputs. Finally, the evaluation of accuracy has been challenging due to the lack of genome-scale, ground-truth networks.

RESULTS

We developed NetProphet3, which combines scores from multiple analyses automatically, using a tree boosting algorithm trained on TF binding location data. We also developed three independent, genome-scale evaluation metrics. By these metrics, NetProphet3 is more accurate than other commonly used packages, including NetProphet 2.0, when gene expression data from direct TF perturbations are available. Furthermore, its integration mode can forge a consensus network from gene expression data and TF binding location data.

AVAILABILITY AND IMPLEMENTATION

All data and code are available at https://zenodo.org/record/7504131#.Y7Wu3i-B2x8.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data

MOTIVATION

Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression is affected by the regulatory regions. A recent successful method is the Activity-By-Contact (ABC) model which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. However, it describes regulatory interactions entirely from a gene's perspective, and does not account for all the candidate target genes of an enhancer. In addition, the ABC model requires two types of assays to measure enhancer activity, which limits the applicability. Moreover, there is neither implementation available that could allow for an integration with transcription factor (TF) binding information nor an efficient analysis of single-cell data.

RESULTS

We demonstrate that the ABC score can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes and also by considering all annotated transcription start sites of a gene. Further, we show that the model is comparably accurate with only one assay to measure enhancer activity. We combined our generalized ABC model with TF binding information and illustrated an analysis of a single-cell ATAC-seq dataset of the human heart, where we were able to characterize cell type-specific regulatory interactions and predict gene expression based on TF affinities. All executed processing steps are incorporated into our new computational pipeline STARE.

AVAILABILITY AND IMPLEMENTATION

The software is available at https://github.com/schulzlab/STARE.

CONTACT

marcel.schulz@em.uni-frankfurt.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Limited conservation in cross-species comparison of GLK transcription factor binding suggested wide-spread cistrome divergence

Non-coding cis-regulatory variants in animal genomes are an important driving force in the evolution of transcription regulation and phenotype diversity. However, cistrome dynamics in plants remain largely underexplored. Here, we compare the binding of GOLDEN2-LIKE (GLK) transcription factors in tomato, tobacco, Arabidopsis, maize and rice. Although the function of GLKs is conserved, most of their binding sites are species-specific. Conserved binding sites are often found near photosynthetic genes dependent on GLK for expression, but sites near non-differentially expressed genes in the glk mutant are nevertheless under purifying selection. The binding sites' regulatory potential can be predicted by machine learning model using quantitative genome features and TF co-binding information. Our study show that genome cis-variation caused wide-spread TF binding divergence, and most of the TF binding sites are genetically redundant. This poses a major challenge for interpreting the effect of individual sites and highlights the importance of quantitatively measuring TF occupancy.

Detecting clusters of transcription factors based on a nonhomogeneous poisson process model

BACKGROUND

Rapidly growing genome-wide ChIP-seq data have provided unprecedented opportunities to explore transcription factor (TF) binding under various cellular conditions. Despite the rich resources, development of analytical methods for studying the interaction among TFs in gene regulation still lags behind.

RESULTS

In order to address cooperative TF binding and detect TF clusters with coordinative functions, we have developed novel computational methods based on clustering the sample paths of nonhomogeneous Poisson processes. Simulation studies demonstrated the capability of these methods to accurately detect TF clusters and uncover the hierarchy of TF interactions. A further application to the multiple-TF ChIP-seq data in mouse embryonic stem cells (ESCs) showed that our methods identified the cluster of core ESC regulators reported in the literature and provided new insights on functional implications of transcrisptional regulatory modules.

CONCLUSIONS

Effective analytical tools are essential for studying protein-DNA relations. Information derived from this research will help us better understand the orchestration of transcription factors in gene regulation processes.

Measuring transcription factor binding and gene expression using barcoded self-reporting transposon calling cards and transcriptomes

Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2 and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells.

Context-dependent CpG methylation directs cell-specific binding of transcription factor ZBTB38

DNA methylation on CpGs regulates transcription in mammals, both by decreasing the binding of methylation-repelled factors and by increasing the binding of methylation-attracted factors. Among the latter, zinc finger proteins have the potential to bind methylated CpGs in a sequence-specific context. The protein ZBTB38 is unique in that it has two independent sets of zinc fingers, which recognize two different methylated consensus sequences in vitro. Here, we identify the binding sites of ZBTB38 in a human cell line, and show that they contain the two methylated consensus sequences identified in vitro. In addition, we show that the distribution of ZBTB38 sites is highly unusual: while 10% of the ZBTB38 sites are also bound by CTCF, the other 90% of sites reside in closed chromatin and are not bound by any of the other factors mapped in our model cell line. Finally, a third of ZBTB38 sites are found upstream of long and active CpG islands. Our work therefore validates ZBTB38 as a methyl-DNA binder in vivo and identifies its unique distribution in the genome.

Systematic discovery and perturbation of regulatory genes in human T cells reveals the architecture of immune networks

Gene regulatory networks ensure that important genes are expressed at precise levels. When gene expression is sufficiently perturbed, it can lead to disease. To understand how gene expression disruptions percolate through a network, we must first map connections between regulatory genes and their downstream targets. However, we lack comprehensive knowledge of the upstream regulators of most genes. Here, we developed an approach for systematic discovery of upstream regulators of critical immune factors-IL2RA, IL-2 and CTLA4-in primary human T cells. Then, we mapped the network of the target genes of these regulators and putative cis-regulatory elements using CRISPR perturbations, RNA-seq and ATAC-seq. These regulators form densely interconnected networks with extensive feedback loops. Furthermore, this network is enriched for immune-associated disease variants and genes. These results provide insight into how immune-associated disease genes are regulated in T cells and broader principles about the structure of human gene regulatory networks.

Predicting which genes will respond to transcription factor perturbations

The ability to predict which genes will respond to the perturbation of a transcription factor serves as a benchmark for our systems-level understanding of transcriptional regulatory networks. In previous work, machine learning models have been trained to predict static gene expression levels in a biological sample by using data from the same or similar samples, including data on their transcription factor binding locations, histone marks, or DNA sequence. We report on a different challenge-training machine learning models to predict which genes will respond to the perturbation of a transcription factor without using any data from the perturbed cells. We find that existing transcription factor location data (ChIP-seq) from human cells have very little detectable utility for predicting which genes will respond to perturbation of a transcription factor. Features of genes, including their preperturbation expression level and expression variation, are very useful for predicting responses to perturbation of any transcription factor. This shows that some genes are poised to respond to transcription factor perturbations and others are resistant, shedding light on why it has been so difficult to predict responses from binding locations. Certain histone marks, including H3K4me1 and H3K4me3, have some predictive power when located downstream of the transcription start site. However, the predictive power of histone marks is much less than that of gene expression level and expression variation. Sequence-based or epigenetic properties of genes strongly influence their tendency to respond to direct transcription factor perturbations, partially explaining the oft-noted difficulty of predicting responsiveness from transcription factor binding location data. These molecular features are largely reflected in and summarized by the gene's expression level and expression variation. Code is available at https://github.com/BrentLab/TFPertRespExplainer.

Fast and accurate inference of gene regulatory networks through robust precision matrix estimation

MOTIVATION

Transcriptional regulation mechanisms allow cells to adapt and respond to external stimuli by altering gene expression. The possible cell transcriptional states are determined by the underlying gene regulatory network (GRN), and reliably inferring such network would be invaluable to understand biological processes and disease progression.

RESULTS

In this article, we present a novel method for the inference of GRNs, called PORTIA, which is based on robust precision matrix estimation, and we show that it positively compares with state-of-the-art methods while being orders of magnitude faster. We extensively validated PORTIA using the DREAM and MERLIN+P datasets as benchmarks. In addition, we propose a novel scoring metric that builds on graph-theoretical concepts.

AVAILABILITY AND IMPLEMENTATION

The code and instructions for data acquisition and full reproduction of our results are available at https://github.com/AntoinePassemiers/PORTIA-Manuscript. PORTIA is available on PyPI as a Python package (portia-grn).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Unveiling RCOR1 as a rheostat at transcriptionally permissive chromatin

RCOR1 is a known transcription repressor that recruits and positions LSD1 and HDAC1/2 on chromatin to erase histone methylation and acetylation. However, there is currently an incomplete understanding of RCOR1’s range of localization and function. Here, we probe RCOR1’s distribution on a genome-wide scale and unexpectedly find that RCOR1 is predominantly associated with transcriptionally active genes. Biochemical analysis reveals that RCOR1 associates with RNA Polymerase II (POL-II) during transcription and deacetylates its carboxy-terminal domain (CTD) at lysine 7. We provide evidence that this non-canonical RCOR1 activity is linked to dampening of POL-II productive elongation at actively transcribing genes. Thus, RCOR1 represses transcription in two ways—first, via a canonical mechanism by erasing transcriptionally permissive histone modifications through associating with HDACs and, second, via a non-canonical mechanism that deacetylates RNA POL-II’s CTD to inhibit productive elongation. We conclude that RCOR1 is a transcription rheostat.

The classical neuronal-gene corepressor RCOR1/CoREST is paradoxically enriched in transcriptionally active chromatin. Here the authors show RCOR1 is recruited during promoter-proximal pausing and negatively regulates the nascent-transcript synthesis. They also show that an RCOR1-LSD1- HDAC1 complex removes lysine acetylation from RNA polymerase II to repress transcription.

Highly enriched BEND3 prevents the premature activation of bivalent genes during differentiation

Bivalent genes are ready for activation upon the arrival of developmental cues. Here, we report that BEND3 is a CpG island (CGI)-binding protein that is enriched at regulatory elements. The cocrystal structure of BEND3 in complex with its target DNA reveals the structural basis for its DNA methylation-sensitive binding property. Mouse embryos ablated of Bend3 died at the pregastrulation stage. Bend3 null embryonic stem cells (ESCs) exhibited severe defects in differentiation, during which hundreds of CGI-containing bivalent genes were prematurely activated. BEND3 is required for the stable association of polycomb repressive complex 2 (PRC2) at bivalent genes that are highly occupied by BEND3, which suggests a reining function of BEND3 in maintaining high levels of H3K27me3 at these bivalent genes in ESCs to prevent their premature activation in the forthcoming developmental stage.

Baiap3 regulates depressive behaviors in mice via attenuating dense core vesicle trafficking in subsets of prefrontal cortex neurons

Selective serotonin reuptake inhibitors (SSRIs) are effective first line therapies for treating depression, but are plagued by undesirable side effects and are not effective in all patients. Because SSRIs effectively deplete the neuronal releasable serotonin (5-HT) pool, gaining a deeper understanding of intracellular mechanisms regulating 5-HT pools can help us understand the shortcomings of SSRIs and develop more effective therapies. In this study, we found that BAIAP3 (brain-specific angiogenesis inhibitor 1-associated protein 3) is significantly downregulated in two mouse models of depression (the IR- and CUMS-induced depressive mouse models). In BAIAP3 downregulated models (in vitro and in vivo), we discovered that trafficking of dense core vesicle (DCV), organelles that store, transport and release cargo via exocytosis, was reduced. Accordingly, 5-HT exocytosis and levels in the synapse were lowered, causing defective post-synaptic neurotransmission. In a screen of natural products, we identified eucalyptol, the active components of Eucalyptus, as uniquely capable of increasing neuronal Baiap3 expression and elevate synaptic 5-HT levels. Moreover, eucalyptol treatment relieved depressive behavioral symptoms and restored serotonin levels in mice. Mechanistically, eucalyptol restores Baiap3 expression by reducing inhibitory microRNAs (miR-329, miR-362). These findings illuminate how Baiap3 depletion propagates neurotransmission dysfunction and point to eucalyptol as a novel agent for restoring serotonin exocytosis, suggesting potential for developing eucalyptol as a therapy for treating depression.

Transcription factor enrichment analysis (TFEA) quantifies the activity of multiple transcription factors from a single experiment

Detecting changes in the activity of a transcription factor (TF) in response to a perturbation provides insights into the underlying cellular process. Transcription Factor Enrichment Analysis (TFEA) is a robust and reliable computational method that detects positional motif enrichment associated with changes in transcription observed in response to a perturbation. TFEA detects positional motif enrichment within a list of ranked regions of interest (ROIs), typically sites of RNA polymerase initiation inferred from regulatory data such as nascent transcription. Therefore, we also introduce muMerge, a statistically principled method of generating a consensus list of ROIs from multiple replicates and conditions. TFEA is broadly applicable to data that informs on transcriptional regulation including nascent transcription (eg. PRO-Seq), CAGE, histone ChIP-Seq, and accessibility data (e.g., ATAC-Seq). TFEA not only identifies the key regulators responding to a perturbation, but also temporally unravels regulatory networks with time series data. Consequently, TFEA serves as a hypothesis-generating tool that provides an easy, rigorous, and cost-effective means to broadly assess TF activity yielding new biological insights.

Epstein-Barr virus inactivates the transcriptome and disrupts the chromatin architecture of its host cell in the first phase of lytic reactivation

Epstein-Barr virus (EBV), a herpes virus also termed HHV 4 and the first identified human tumor virus, establishes a stable, long-term latent infection in human B cells, its preferred host. Upon induction of EBV's lytic phase, the latently infected cells turn into a virus factory, a process that is governed by EBV. In the lytic, productive phase, all herpes viruses ensure the efficient induction of all lytic viral genes to produce progeny, but certain of these genes also repress the ensuing antiviral responses of the virally infected host cells, regulate their apoptotic death or control the cellular transcriptome. We now find that EBV causes previously unknown massive and global alterations in the chromatin of its host cell upon induction of the viral lytic phase and prior to the onset of viral DNA replication. The viral initiator protein of the lytic cycle, BZLF1, binds to >105 binding sites with different sequence motifs in cellular chromatin in a concentration dependent manner implementing a binary molar switch probably to prevent noise-induced erroneous induction of EBV's lytic phase. Concomitant with DNA binding of BZLF1, silent chromatin opens locally as shown by ATAC-seq experiments, while previously wide-open cellular chromatin becomes inaccessible on a global scale within hours. While viral transcripts increase drastically, the induction of the lytic phase results in a massive reduction of cellular transcripts and a loss of chromatin-chromatin interactions of cellular promoters with their distal regulatory elements as shown in Capture-C experiments. Our data document that EBV's lytic cycle induces discrete early processes that disrupt the architecture of host cellular chromatin and repress the cellular epigenome and transcriptome likely supporting the efficient de novo synthesis of this herpes virus.

Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases

Variants of transcription factor binding sites (TFBSs) constitute an important part of the human genome. Current evidence demonstrates close links between nucleotides within TFBSs and gene expression. There are multiple pathways through which genomic sequences located in TFBSs regulate gene expression, and recent genome-wide association studies have shown the biological significance of TFBS variation in human phenotypes. However, numerous challenges remain in the study of TFBS polymorphisms. This article aims to cover the current state of understanding as regards the genomic features of TFBSs and TFBS variants; the mechanisms through which TFBS variants regulate gene expression; the approaches to studying the effects of nucleotide changes that create or disrupt TFBSs; the challenges faced in studies of TFBS sequence variations; the effects of natural selection on collections of TFBSs; in addition to the insights gained from the study of TFBS alleles related to gout, its associated comorbidities (increased body mass index, chronic kidney disease, diabetes, dyslipidemia, coronary artery disease, ischemic heart disease, hypertension, hyperuricemia, osteoporosis, and prostate cancer), and the treatment responses of patients.

YY1 regulated transcription-based stratification of gastric tumors and identification of potential therapeutic candidates

Gastric cancer is one of the leading causes of cancer-related death worldwide. The transcription factor YY1 regulates diverse biological processes, including cell proliferation, development, DNA damage responses, and carcinogenesis. This study was designed to explore the role of YY1 regulated transcription in gastric cancer. YY1 silencing in gastric cancer cells has resulted in the inhibition of Wnt/β-catenin, JNK/MAPK, ERK/MAPK, ER, and HIF-1α signaling pathways. Genome-wide mRNA profiling upon silencing the expression YY1 gene in gastric cancer cells and comparison with the previously identified YY1 regulated genes from other lineages revealed a moderate overlap among the YY1 regulated genes. Despite the differing genes, all the YY1 regulated gene sets were expressed in most of the intestinal subtype gastric tumors and a subset of diffuse subtype gastric tumors. Integrative functional genomic analysis of the YY1 gene sets revealed an association among the pathways Wnt/β-catenin, Rapamycin, Cyclin-D1, Myc, E2F, PDGF, and AKT. Further, the drugs capable of inhibiting YY1 mediated transcription were identified as suitable targeted therapeutic candidates for gastric tumors with activated YY1. The data emerging from the investigation would pave the way for the development of YY1-based targeted therapeutics for gastric cancer.

Comparison between instrumental variable and mediation-based methods for reconstructing causal gene networks in yeast

Causal gene networks model the flow of information within a cell. Reconstructing causal networks from omics data is challenging because correlation does not imply causation. When genomics and transcriptomics data from a segregating population are combined, genomic variants can be used to orient the direction of causality between gene expression traits. Instrumental variable methods use a local expression quantitative trait locus (eQTL) as a randomized instrument for a gene's expression level, and assign target genes based on distal eQTL associations. Mediation-based methods additionally require that distal eQTL associations are mediated by the source gene. A detailed comparison between these methods has not yet been conducted, due to the lack of a standardized implementation of different methods, the limited sample size of most multi-omics datasets, and the absence of ground-truth networks for most organisms. Here we used Findr, a software package providing uniform implementations of instrumental variable, mediation, and coexpression-based methods, a recent dataset of 1012 segregants from a cross between two budding yeast strains, and the Yeastract database of known transcriptional interactions to compare causal gene network inference methods. We found that causal inference methods result in a significant overlap with the ground-truth, whereas coexpression did not perform better than random. A subsampling analysis revealed that the performance of mediation saturates at large sample sizes, due to a loss of sensitivity when residual correlations become significant. Instrumental variable methods on the other hand contain false positive predictions, due to genomic linkage between eQTL instruments. Instrumental variable and mediation-based methods also have complementary roles for identifying causal genes underlying transcriptional hotspots. Instrumental variable methods correctly predicted STB5 targets for a hotspot centred on the transcription factor STB5, whereas mediation failed due to Stb5p auto-regulating its own expression. Mediation suggests a new candidate gene, DNM1, for a hotspot on Chr XII, whereas instrumental variable methods could not distinguish between multiple genes located within the hotspot. In conclusion, causal inference from genomics and transcriptomics data is a powerful approach for reconstructing causal gene networks, which could be further improved by the development of methods to control for residual correlations in mediation analyses, and for genomic linkage and pleiotropic effects from transcriptional hotspots in instrumental variable analyses.

Comprehensive network modeling from single cell RNA sequencing of human and mouse reveals well conserved transcription regulation of hematopoiesis

Background

Presently, there is no comprehensive analysis of the transcription regulation network in hematopoiesis. Comparison of networks arising from gene co-expression across species can facilitate an understanding of the conservation of functional gene modules in hematopoiesis.

Results

We used single-cell RNA sequencing to profile bone marrow from human and mouse, and inferred transcription regulatory networks in each species in order to characterize transcriptional programs governing hematopoietic stem cell differentiation. We designed an algorithm for network reconstruction to conduct comparative transcriptomic analysis of hematopoietic gene co-expression and transcription regulation in human and mouse bone marrow cells. Co-expression network connectivity of hematopoiesis-related genes was found to be well conserved between mouse and human. The co-expression network showed “small-world” and “scale-free” architecture. The gene regulatory network formed a hierarchical structure, and hematopoiesis transcription factors localized to the hierarchy’s middle level.

Conclusions

Transcriptional regulatory networks are well conserved between human and mouse. The hierarchical organization of transcription factors may provide insights into hematopoietic cell lineage commitment, and to signal processing, cell survival and disease initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07241-2.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Learning gene networks underlying clinical phenotypes using SNP perturbation

Availability of genome sequence, molecular, and clinical phenotype data for large patient cohorts generated by recent technological advances provides an opportunity to dissect the genetic architecture of complex diseases at system level. However, previous analyses of such data have largely focused on the co-localization of SNPs associated with clinical and expression traits, each identified from genome-wide association studies and expression quantitative trait locus mapping. Thus, their description of the molecular mechanisms behind the SNPs influencing clinical phenotypes was limited to the single gene linked to the co-localized SNP. Here we introduce PerturbNet, a statistical framework for learning gene networks that modulate the influence of genetic variants on phenotypes, using genetic variants as naturally occurring perturbation of a biological system. PerturbNet uses a probabilistic graphical model to directly model the cascade of perturbation from genetic variants to the gene network to the phenotype network along with the networks at each layer of the biological system. PerturbNet learns the entire model by solving a single optimization problem with an efficient algorithm that can analyze human genome-wide data within a few hours. PerturbNet inference procedures extract a detailed description of how the gene network modulates the genetic effects on phenotypes. Using simulated and asthma data, we demonstrate that PerturbNet improves statistical power for detecting disease-linked SNPs and identifies gene networks and network modules mediating the SNP effects on traits, providing deeper insights into the underlying molecular mechanisms.

Transcription Factors in Cartilage Homeostasis and Osteoarthritis

Osteoarthritis (OA) is the most common degenerative joint disease, and it is characterized by articular cartilage loss. In part, OA is caused by aberrant anabolic and catabolic activities of the chondrocyte, the only cell type present in cartilage. These chondrocyte activities depend on the intra- and extracellular signals that the cell receives and integrates into gene expression. The key proteins for this integration are transcription factors. A large number of transcription factors exist, and a better understanding of the transcription factors activated by the various signaling pathways active during OA can help us to better understand the complex etiology of OA. In addition, establishing such a profile can help to stratify patients in different subtypes, which can be a very useful approach towards personalized therapy. In this review, we discuss crucial transcription factors for extracellular matrix metabolism, chondrocyte hypertrophy, chondrocyte senescence, and autophagy in chondrocytes. In addition, we discuss how insight into these factors can be used for treatment purposes.

Genomic Resolution of DLX-Orchestrated Transcriptional Circuits Driving Development of Forebrain GABAergic Neurons

DLX transcription factors (TFs) are master regulators of the developing vertebrate brain, driving forebrain GABAergic neuronal differentiation. Ablation of Dlx1&2 alters expression of genes that are critical for forebrain GABAergic development. We integrated epigenomic and transcriptomic analyses, complemented with in situ hybridization (ISH), and in vivo and in vitro studies of regulatory element (RE) function. This revealed the DLX-organized gene regulatory network at genomic, cellular, and spatial levels in mouse embryonic basal ganglia. DLX TFs perform dual activating and repressing functions; the consequences of their binding were determined by the sequence and genomic context of target loci. Our results reveal and, in part, explain the paradox of widespread DLX binding contrasted with a limited subset of target loci that are sensitive at the epigenomic and transcriptomic level to Dlx1&2 ablation. The regulatory properties identified here for DLX TFs suggest general mechanisms by which TFs orchestrate dynamic expression programs underlying neurodevelopment.

Lindtner et al. reveal the regulatory wiring organized by DLX transcription factors in forebrain GABAergic neuronal specification, by integrating functional genomic, epigenomic, and genetic data on a transgenic mouse model. This network determines key sequence-encoded regulatory elements and implicates a combination of histone modifications and biophysical interactions.