Predictive models of eukaryotic transcriptional regulation reveals changes in transcription factor roles and promoter usage between metabolic conditions

Tissue expression and promoter activity analysis of the porcine TNFSF11 gene

Tumour necrosis factor (TNF) superfamily member 11 (TNFSF11), also known as RANKL, plays a crucial role in regulating several physiological and pathological activities. Additionally, it is a vital factor in bone physiology, and the sex hormone progesterone regulates the expansion of stem cells and the proliferation of mammary epithelial cells. It is essential for animal growth and reproductive physiological processes. This study aimed to evaluate the tissue-specific expression characteristics and promoter activity of the TNFSF11 gene in pigs. As a result, the study examined the presence of TNFSF11 expression in the tissues of Xiangsu pigs at 0.6 and 12 months of age. Moreover, the core promoter region of TNFSF11 was also identified by utilizing a combination of bioinformatic prediction and dual-luciferase activity tests. Finally, the effect of transcription factors on the transcriptional activity of the core promoter region was determined using site-directed mutagenesis. TNFSF11 was uniformly expressed in all tissues; however, its expression in muscles was comparatively low. The core promoter region of TNFSF11 was located in the -555 to -1 region. The prediction of the transcription start site of TNFSF11 gene-2000 ∼ + 500bp showed that there was a CpG site in 17 ∼ + 487bp. Analysis of mutations in the transcription factor binding sites revealed that mutations in the Stat5b, Myog, Trl, and EN1 binding sites had significant effects on the transcriptional activity of the TNFSF11 gene, particularly following the EN1 binding site mutation (P < 0.001). This study provides insights into both the tissue-specific expression patterns of TNFSF11 in the tissues of Xiangsu pigs and the potential regulatory effects of transcription factors on its promoter activity. These results may be helpful for future research aimed at clarifying the expression and role of the porcine TNFSF11 gene.

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.

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.

Predicting Tissue-Specific mRNA and Protein Abundance in Maize: A Machine Learning Approach

Machine learning and modeling approaches have been used to classify protein sequences for a broad set of tasks including predicting protein function, structure, expression, and localization. Some recent studies have successfully predicted whether a given gene is expressed as mRNA or even translated to proteins potentially, but given that not all genes are expressed in every condition and tissue, the challenge remains to predict condition-specific expression. To address this gap, we developed a machine learning approach to predict tissue-specific gene expression across 23 different tissues in maize, solely based on DNA promoter and protein sequences. For class labels, we defined high and low expression levels for mRNA and protein abundance and optimized classifiers by systematically exploring various methods and combinations of k-mer sequences in a two-phase approach. In the first phase, we developed Markov model classifiers for each tissue and built a feature vector based on the predictions. In the second phase, the feature vector was used as an input to a Bayesian network for final classification. Our results show that these methods can achieve high classification accuracy of up to 95% for predicting gene expression for individual tissues. By relying on sequence alone, our method works in settings where costly experimental data are unavailable and reveals useful insights into the functional, evolutionary, and regulatory characteristics of genes.

Cancer progression as a learning process

Drug resistance and metastasis—the major complications in cancer—both entail adaptation of cancer cells to stress, whether a drug or a lethal new environment. Intriguingly, these adaptive processes share similar features that cannot be explained by a pure Darwinian scheme, including dormancy, increased heterogeneity, and stress-induced plasticity. Here, we propose that learning theory offers a framework to explain these features and may shed light on these two intricate processes. In this framework, learning is performed at the single-cell level, by stress-driven exploratory trial-and-error. Such a process is not contingent on pre-existing pathways but on a random search for a state that diminishes the stress. We review underlying mechanisms that may support this search, and show by using a learning model that such exploratory learning is feasible in a high-dimensional system as the cell. At the population level, we view the tissue as a network of exploring agents that communicate, restraining cancer formation in health. In this view, disease results from the breakdown of homeostasis between cellular exploratory drive and tissue homeostasis.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

Rational gRNA design based on transcription factor binding data

The clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system has become a standard tool in many genome engineering endeavors. The endonuclease-deficient version of Cas9 (dCas9) is also a powerful programmable tool for gene regulation. In this study, we made use of Saccharomyces cerevisiae transcription factor (TF) binding data to obtain a better understanding of the interplay between TF binding and binding of dCas9 fused to an activator domain, VPR. More specifically, we targeted dCas9–VPR toward binding sites of Gcr1–Gcr2 and Tye7 present in several promoters of genes encoding enzymes engaged in the central carbon metabolism. From our data, we observed an upregulation of gene expression when dCas9–VPR was targeted next to a TF binding motif, whereas a downregulation or no change was observed when dCas9 was bound on a TF motif. This suggests a steric competition between dCas9 and the specific TF. Integrating TF binding data, therefore, proved to be useful for designing guide RNAs for CRISPR interference or CRISPR activation applications.

Inferring TF activities and activity regulators from gene expression data with constraints from TF perturbation data

Motivation

The activity of a transcription factor (TF) in a sample of cells is the extent to which it is exerting its regulatory potential. Many methods of inferring TF activity from gene expression data have been described, but due to the lack of appropriate large-scale datasets, systematic and objective validation has not been possible until now.

Results

We systematically evaluate and optimize the approach to TF activity inference in which a gene expression matrix is factored into a condition-independent matrix of control strengths and a condition-dependent matrix of TF activity levels. We find that expression data in which the activities of individual TFs have been perturbed are both necessary and sufficient for obtaining good performance. To a considerable extent, control strengths inferred using expression data from one growth condition carry over to other conditions, so the control strength matrices derived here can be used by others. Finally, we apply these methods to gain insight into the upstream factors that regulate the activities of yeast TFs Gcr2, Gln3, Gcn4 and Msn2.

Availability and implementation

Evaluation code and data are available at https://doi.org/10.5281/zenodo.4050573.

Supplementary information

Supplementary data are available at Bioinformatics online.

Differential activation mechanisms of two isoforms of Gcr1 transcription factor generated from spliced and un-spliced transcripts in Saccharomyces cerevisiae

Gcr1, an important transcription factor for glycolytic genes in Saccharomyces cerevisiae, was recently revealed to have two isoforms, Gcr1^U and Gcr1^S, produced from un-spliced and spliced transcripts, respectively. In this study, by generating strains expressing only Gcr1^U or Gcr1^S using the CRISPR/Cas9 system, we elucidate differential activation mechanisms of these two isoforms. The Gcr1^U monomer forms an active complex with its coactivator Gcr2 homodimer, whereas Gcr1^S acts as a homodimer without Gcr2. The USS domain, 55 residues at the N-terminus existing only in Gcr1^U, inhibits dimerization of Gcr1^U and even acts in trans to inhibit Gcr1^S dimerization. The Gcr1^S monomer inhibits the metabolic switch from fermentation to respiration by directly binding to the ALD4 promoter, which can be restored by overexpression of the ALD4 gene, encoding a mitochondrial aldehyde dehydrogenase required for ethanol utilization. Gcr1^U and Gcr1^S regulate almost the same target genes, but show unique activities depending on growth phase, suggesting that these isoforms play differential roles through separate activation mechanisms depending on environmental conditions.

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

Genome-scale reconstruction of Gcn4/ATF4 networks driving a growth program

Growth and starvation are considered opposite ends of a spectrum. To sustain growth, cells use coordinated gene expression programs and manage biomolecule supply in order to match the demands of metabolism and translation. Global growth programs complement increased ribosomal biogenesis with sufficient carbon metabolism, amino acid and nucleotide biosynthesis. How these resources are collectively managed is a fundamental question. The role of the Gcn4/ATF4 transcription factor has been best studied in contexts where cells encounter amino acid starvation. However, high Gcn4 activity has been observed in contexts of rapid cell proliferation, and the roles of Gcn4 in such growth contexts are unclear. Here, using a methionine-induced growth program in yeast, we show that Gcn4/ATF4 is the fulcrum that maintains metabolic supply in order to sustain translation outputs. By integrating matched transcriptome and ChIP-Seq analysis, we decipher genome-wide direct and indirect roles for Gcn4 in this growth program. Genes that enable metabolic precursor biosynthesis indispensably require Gcn4; contrastingly ribosomal genes are partly repressed by Gcn4. Gcn4 directly binds promoter-regions and transcribes a subset of metabolic genes, particularly driving lysine and arginine biosynthesis. Gcn4 also globally represses lysine and arginine enriched transcripts, which include genes encoding the translation machinery. The Gcn4 dependent lysine and arginine supply thereby maintains the synthesis of the translation machinery. This is required to maintain translation capacity. Gcn4 consequently enables metabolic-precursor supply to bolster protein synthesis, and drive a growth program. Thus, we illustrate how growth and starvation outcomes are both controlled using the same Gcn4 transcriptional outputs that function in distinct contexts.

The transcription factor Leu3 shows differential binding behavior in response to changing leucine availability

The main transcriptional regulator of leucine biosynthesis in the yeast Saccharomyces cerevisiae is the transcription factor Leu3. It has previously been reported that Leu3 always binds to its target genes, but requires activation to induce their expression. In a recent large-scale study of high-resolution transcription factor binding site identification, we showed that Leu3 has divergent binding sites in different cultivation conditions, thereby questioning the results of earlier studies. Here, we present a follow-up study using chromatin immunoprecipitation followed by sequencing (ChIP-seq) to investigate the influence of leucine supplementation on Leu3 binding activity and strength. With this new data set we are able to show that Leu3 exhibits changes in binding activity in response to changing levels of leucine availability.

Obtaining a series of native gradient promoter-5'-UTR sequences in Corynebacterium glutamicum ATCC 13032

BACKGROUND

Corynebacterium glutamicum is an important industrial microorganism used for the production of many valuable compounds, especially amino acids and their derivatives. For fine-tuning of metabolic pathways, synthetic biological tools are largely based on the rational application of promoters. However, the limited number of promoters make it difficult.

RESULTS

In this study, according to the analysis of RNA-Seq data, 90 DNA fragments with lengths of 200-500 bp that may contain promoter-5'-UTR (PUTR) sequences were amplified and linked to a fluorescent protein gene. When compared with the common strong PUTR P_sodUTR, 17 strong PUTRs were obtained, which maintained stable expression strengths from the early to post stationary phase. Among them, P_NCgl1676UTR was the strongest and its fluorescent protein expression level was more than five times higher than that of P_sodUTR. Furthermore, nine typical chemicals related to the biosynthesis of sulfur-containing amino acids (such as L-methionine, L-cysteine) were selected as stress substances to preliminarily explore the stress on these PUTRs. The results showed that the expression of P_brnFUTR was activated by L-methionine, while that of P_NCgl1202UTR was severely inhibited by L-lysine.

CONCLUSIONS

These findings demonstrated that the selected PUTRs can stably express different genes, such as the red fluorescence protein gene, and can be useful for fine-tuning regulation of metabolic networks in C. glutamicum or for establishing high-throughput screening strategies through biosensor for the production of useful compounds.

Dual threshold optimization and network inference reveal convergent evidence from TF binding locations and TF perturbation responses

A high-confidence map of the direct, functional targets of each transcription factor (TF) requires convergent evidence from independent sources. Two significant sources of evidence are TF binding locations and the transcriptional responses to direct TF perturbations. Systematic data sets of both types exist for yeast and human, but they rarely converge on a common set of direct, functional targets for a TF. Even the few genes that are both bound and responsive may not be direct functional targets. Our analysis shows that when there are many nonfunctional binding sites and many indirect targets, nonfunctional sites are expected to occur in the cis-regulatory DNA of indirect targets by chance. To address this problem, we introduce dual threshold optimization (DTO), a new method for setting significance thresholds on binding and perturbation-response data, and show that it improves convergence. It also enables comparison of binding data to perturbation-response data that have been processed by network inference algorithms, which further improves convergence. The combination of dual threshold optimization and network inference greatly expands the high-confidence TF network map in both yeast and human. Next, we analyze a comprehensive new data set measuring the transcriptional response shortly after inducing overexpression of a yeast TF. We also present a new yeast binding location data set obtained by transposon calling cards and compare it to recent ChIP-exo data. These new data sets improve convergence and expand the high-confidence network synergistically.

Yeast systems biology in understanding principles of physiology underlying complex human diseases

Complex human diseases commonly arise from deregulation of cell growth, metabolism, and/or gene expression. Yeast is a eukaryal model organism that is widely used to study these processes. Yeast systems biology benefits from the ability to exert fine experimental control over the cell growth rate and nutrient composition, which allows orthogonal experimental design and generation of multi-omics data at high resolution. This has led to several insights on the principles of cellular physiology, including many cellular processes associated with complex human diseases. Here we review these biological insights together with experimental and modeling approaches developed in yeast to study systems biology. The role of yeast systems biology to further advance systems and personalized therapies for complex diseases is discussed.

ChIP-exo analysis highlights Fkh1 and Fkh2 transcription factors as hubs that integrate multi-scale networks in budding yeast

The understanding of the multi-scale nature of molecular networks represents a major challenge. For example, regulation of a timely cell cycle must be coordinated with growth, during which changes in metabolism occur, and integrate information from the extracellular environment, e.g. signal transduction. Forkhead transcription factors are evolutionarily conserved among eukaryotes, and coordinate a timely cell cycle progression in budding yeast. Specifically, Fkh1 and Fkh2 are expressed during a lengthy window of the cell cycle, thus are potentially able to function as hubs in the multi-scale cellular environment that interlocks various biochemical networks. Here we report on a novel ChIP-exo dataset for Fkh1 and Fkh2 in both logarithmic and stationary phases, which is analyzed by novel and existing software tools. Our analysis confirms known Forkhead targets from available ChIP-chip studies and highlights novel ones involved in the cell cycle, metabolism and signal transduction. Target genes are analyzed with respect to their function, temporal expression during the cell cycle, correlation with Fkh1 and Fkh2 as well as signaling and metabolic pathways they occur in. Furthermore, differences in targets between Fkh1 and Fkh2 are presented. Our work highlights Forkhead transcription factors as hubs that integrate multi-scale networks to achieve proper timing of cell division in budding yeast.

Big data in yeast systems biology

Systems biology uses computational and mathematical modeling to study complex interactions in a biological system. The yeast Saccharomyces cerevisiae, which has served as both an important model organism and cell factory, has pioneered both the early development of such models and modeling concepts, and the more recent integration of multi-omics big data in these models to elucidate fundamental principles of biology. Here, we review the advancement of big data technologies to gain biological insight in three aspects of yeast systems biology: gene expression dynamics, cellular metabolism and the regulation network between gene expression and metabolism. The role of big data and complementary modeling approaches, including the expansion of genome-scale metabolic models and machine learning methodologies, are discussed as key drivers in the rapid advancement of yeast systems biology.