Genomic analysis of the hierarchical structure of regulatory networks

An R2R3-MYB transcription factor PdbMYB6 enhances drought tolerance by mediating reactive oxygen species scavenging, osmotic balance, and stomatal opening

Drought is a major environmental challenge that hinders the growth and development of plants. R2R3-MYB transcription factors (TFs) play a vital role in mediating responses to abiotic stress; however, their specific functions in Populus davidiana × Populus bolleana hybrid poplar plants remain underexplored. This study focused on PdbMYB6, a novel R2R3-MYB TF identified in P. davidiana × P. bolleana plants. We found that PdbMYB6 acts as a transcriptional activator. By conducting functional analyses of both overexpression and knockout models, we demonstrated that PdbMYB6 enhances drought tolerance in plants by improving reactive oxygen species scavenging and modulating osmotic balance. Additionally, PdbMYB6 plays a role in regulating stomatal openings to minimize water loss. The qRT-PCR and RNA sequencing results revealed that PdbMYB6 influences the expression of genes related to stress tolerance. TF-centered Yeast One-Hybrid (Y1H) and chromatin immunoprecipitation (ChIP) assays indicated that PdbMYB6 binds to two novel core sequences (C [A/G/C]TG and [T/A/G]GTA) as well as GT-1 (GGAAA) and MYBCORE (AACGG) elements, which are associated with light responses and stress resistance, thereby promoting the expression of stress-resistant genes. Furthermore, Y1H and ChIP assays identified four upstream factors that regulate PdbMYB6 expression by interacting with specific elements in its promoter. Notably, the overexpression of these four factors enhances plant drought resistance and affects the expression of stress-response genes. Our findings highlight the role of the PdbMYB6 TF in the drought regulatory mechanism and provide potential gene sources for the molecular breeding of drought-resistant plants through genetic engineering.

Chromatin-based memory as a self-stabilizing influence on cell identity

Cell types are traditionally thought to be specified and stabilized by gene regulatory networks. Here, we explore how chromatin memory contributes to the specification and stabilization of cell states. Through pervasive, local, feedback loops, chromatin memory enables cell states that were initially unstable to become stable. Deeper appreciation of this self-stabilizing role for chromatin broadens our perspective of Waddington's epigenetic landscape from a static surface with islands of stability shaped by evolution, to a plasticine surface molded by experience. With implications for the evolution of cell types, stabilization of resistant states in cancer, and the widespread plasticity of complex life.

Global regulator IrrE on stress tolerance: a review

Stress tolerance is a vital attribute for all living beings to cope with environmental adversities. IrrE (also named PprI) from Deinococcus radiodurans enhances resistance to extreme radiation stress by functioning as a global regulator, mediating the transcription of genes involved in deoxyribonucleic acid (DNA) damage response (DDR). The expression of IrrE augmented the resilience of various species to heat, radiation, oxidation, osmotic stresses and inhibitors, encompassing bacterial, fungal, plant, and mammalian cells. Moreover, IrrE was employed in a global regulator engineering strategy to broaden its applications in stress tolerance. The regulatory impacts of heterologously expressed IrrE have been investigated at the molecular and systems level, including the regulation of genes, proteins, modules, or pathways involved in DNA repair, detoxification proteins, protective molecules, native regulators and other aspects. In this review, we discuss the regulatory role and mechanism of IrrE in the antiradiation response of D. radiodurans. Furthermore, the applications and regulatory effects of heterologous expression of IrrE to enhance abiotic stress tolerance are summarized in particular.

Strategies to increase the robustness of microbial cell factories

Engineering microbial cell factories have achieved much progress in producing fuels, natural products and bulk chemicals. However, in industrial fermentation, microbial cells often face various predictable and stochastic disturbances resulting from intermediate metabolites or end product toxicity, metabolic burden and harsh environment. These perturbances can potentially decrease productivity and titer. Therefore, strain robustness is essential to ensure reliable and sustainable production efficiency. In this review, the current strategies to improve host robustness were summarized, including knowledge-based engineering approaches, such as transcription factors, membrane/transporters and stress proteins, and the traditional adaptive laboratory evolution based on natural selection. Computation-assisted (e.g. GEMs, deep learning and machine learning) design of robust industrial hosts was also introduced. Furthermore, the challenges and future perspectives on engineering microbial host robustness are proposed to promote the development of green, efficient and sustainable biomanufacturers.

Insights into the identification and evolutionary conservation of key genes in the transcriptional circuits of meiosis initiation and commitment in budding yeast

Initiation of meiosis in budding yeast does not commit the cells for meiosis. Thus, two distinct signaling cascades may differentially regulate meiosis initiation and commitment in budding yeast. To distinguish between the role of these signaling cascades, we reconstructed protein-protein interaction networks and gene regulatory networks with upregulated genes in meiosis initiation and commitment. Analyzing the integrated networks, we identified four master regulators (MRs) [Ume6p, Msn2p, Met31p, Ino2p], three transcription factors (TFs), and 279 target genes (TGs) unique for meiosis initiation, and three MRs [Ndt80p, Aro80p, Rds2p], 11 TFs, and 948 TGs unique for meiosis commitment. Functional enrichment analysis of these distinct members from the transcriptional cascades for meiosis initiation and commitment revealed that nutritional cues rewire gene expression for initiating meiosis and chromosomal recombination commits cells to meiosis. As meiotic chromosomal recombination is highly conserved in eukaryotes, we compared the evolutionary rate of unique members in the transcriptional cascade of two meiotic phases of Saccharomyces cerevisiae with members of the phylum Ascomycota, revealing that the transcriptional cascade governing chromosomal recombination during meiosis commitment has experienced greater purifying selection pressure (P value = 0.0013, 0.0382, 0.0448, 0.0369, 0.02967, 0.04937, 0.03046, 0.03357 and < 0.00001 for Ashbya gossypii, Yarrowia lipolytica, Debaryomyces hansenii, Aspergillus fumigatus, Neurospora crassa, Kluyveromyces lactis, Schizosaccharomyces pombe, Schizosaccharomyces cryophilus, and Schizosaccharomyces octosporus, respectively). This study demarcates crucial players driving meiosis initiation and commitment and demonstrates their differential rate of evolution in budding yeast.

Evolutionary Plasticity in Insular Lizard, Adapting over Reproduction, Metabolism, and Color Variation

The Italian wall lizard (Podarcis siculus) living on islets exhibits a melanic skin coloration and a suite of adaptive traits lacking in nearby mainland populations. On islets, the unpredictable environmental conditions and highly fluctuating population densities are believed to have produced reversed island syndrome (RIS). Several physiological, behavioral, and life-history changes based on the RIS could result from positive selection on increased activity of melanocortins. We hypothesize that phenotypes on islets are the product of a plastic variation depending on the regulation of specific genes. Focusing on control systems that determine the insular-adapted phenotype, we demonstrated that reproductive markers, involved in the hypothalamus-hypophysis-gonadal axis, and metabolism markers, flags for hypophysis-melanocortin receptors, are all up-regulated in island lizards under the RIS. This behavior, combined with the observed limited variation in the mitochondrial genome, agrees with the hypothesis that plasticity enables populations to persist in novel environmental conditions and that over time, natural selection will "fine-tune" the population to the environment by modifying the phenotype under selection. We believe that analysis of the transcriptome and the single gene expression, such that all the variations observed in the island populations, can be useful to shed light on evolutionary plasticity as a process affecting animals' populations in general.

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Reprogramming of cis-regulatory networks during skeletal muscle atrophy in male mice

A comprehensive atlas of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of cellular structure maintenance, metabolism, and responses to the environment. Here we show, using matched single-nucleus chromatin accessibility and RNA-sequencing from juvenile male C57BL6 mice, an atlas of accessible chromatin regions in both normal and denervated skeletal muscles. We identified cell-type-specific cis-regulatory networks, highlighting the dynamic regulatory circuits mediating transitions between myonuclear types. Through comparison of normal and perturbed muscle, we delineated the reprogramming of cis-regulatory networks in response to denervation, described the interplay of promoters/enhancers and target genes. We further unveil a hierarchical structure of transcription factors that delineate a regulatory network in atrophic muscle, identifying ELK4 as a key atrophy-related transcription factor that instigates muscle atrophy through TGF-β1 regulation. This study furnishes a rich genomic resource, essential for decoding the regulatory dynamics of skeletal muscle in both physiological and pathological states.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

iHerd: an integrative hierarchical graph representation learning framework to quantify network changes and prioritize risk genes in disease

Different genes form complex networks within cells to carry out critical cellular functions, while network alterations in this process can potentially introduce downstream transcriptome perturbations and phenotypic variations. Therefore, developing efficient and interpretable methods to quantify network changes and pinpoint driver genes across conditions is crucial. We propose a hierarchical graph representation learning method, called iHerd. Given a set of networks, iHerd first hierarchically generates a series of coarsened sub-graphs in a data-driven manner, representing network modules at different resolutions (e.g., the level of signaling pathways). Then, it sequentially learns low-dimensional node representations at all hierarchical levels via efficient graph embedding. Lastly, iHerd projects separate gene embeddings onto the same latent space in its graph alignment module to calculate a rewiring index for driver gene prioritization. To demonstrate its effectiveness, we applied iHerd on a tumor-to-normal GRN rewiring analysis and cell-type-specific GCN analysis using single-cell multiome data of the brain. We showed that iHerd can effectively pinpoint novel and well-known risk genes in different diseases. Distinct from existing models, iHerd's graph coarsening for hierarchical learning allows us to successfully classify network driver genes into early and late divergent genes (EDGs and LDGs), emphasizing genes with extensive network changes across and within signaling pathway levels. This unique approach for driver gene classification can provide us with deeper molecular insights. The code is freely available at https://github.com/aicb-ZhangLabs/iHerd. All other relevant data are within the manuscript and supporting information files.

Multi-modular structure of the gene regulatory network for specification and commitment of murine T cells

T cells develop from multipotent progenitors by a gradual process dependent on intrathymic Notch signaling and coupled with extensive proliferation. The stages leading them to T-cell lineage commitment are well characterized by single-cell and bulk RNA analyses of sorted populations and by direct measurements of precursor-product relationships. This process depends not only on Notch signaling but also on multiple transcription factors, some associated with stemness and multipotency, some with alternative lineages, and others associated with T-cell fate. These factors interact in opposing or semi-independent T cell gene regulatory network (GRN) subcircuits that are increasingly well defined. A newly comprehensive picture of this network has emerged. Importantly, because key factors in the GRN can bind to markedly different genomic sites at one stage than they do at other stages, the genes they significantly regulate are also stage-specific. Global transcriptome analyses of perturbations have revealed an underlying modular structure to the T-cell commitment GRN, separating decisions to lose "stem-ness" from decisions to block alternative fates. Finally, the updated network sheds light on the intimate relationship between the T-cell program, which depends on the thymus, and the innate lymphoid cell (ILC) program, which does not.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

Minimal frustration underlies the usefulness of incomplete regulatory network models in biology

Regulatory networks as large and complex as those implicated in cell-fate choice are expected to exhibit intricate, very high-dimensional dynamics. Cell-fate choice, however, is a macroscopically simple process. Additionally, regulatory network models are almost always incomplete and/or inexact, and do not incorporate all the regulators and interactions that may be involved in cell-fate regulation. In spite of these issues, regulatory network models have proven to be incredibly effective tools for understanding cell-fate choice across contexts and for making useful predictions. Here, we show that minimal frustration-a feature of biological networks across contexts but not of random networks-can compel simple, low-dimensional steady-state behavior even in large and complex networks. Moreover, the steady-state behavior of minimally frustrated networks can be recapitulated by simpler networks such as those lacking many of the nodes and edges and those that treat multiple regulators as one. The present study provides a theoretical explanation for the success of network models in biology and for the challenges in network inference.

The Spatial Organization of Bacterial Transcriptional Regulatory Networks

The transcriptional regulatory network (TRN) is the central pivot of a prokaryotic organism to receive, process and respond to internal and external environmental information. However, little is known about its spatial organization so far. In recent years, chromatin interaction data of bacteria such as Escherichia coli and Bacillus subtilis have been published, making it possible to study the spatial organization of bacterial transcriptional regulatory networks. By combining TRNs and chromatin interaction data of E. coli and B. subtilis, we explored the spatial organization characteristics of bacterial TRNs in many aspects such as regulation directions (positive and negative), central nodes (hubs, bottlenecks), hierarchical levels (top, middle, bottom) and network motifs (feed-forward loops and single input modules) of the TRNs and found that the bacterial TRNs have a variety of stable spatial organization features under different physiological conditions that may be closely related with biological functions. Our findings provided new insights into the connection between transcriptional regulation and the spatial organization of chromosome in bacteria and might serve as a factual foundation for trying spatial-distance-based gene circuit design in synthetic biology.

The SmMYB36-SmERF6/SmERF115 module regulates the biosynthesis of tanshinones and phenolic acids in salvia miltiorrhiza hairy roots

Tanshinone and phenolic acids are the most important active substances of Salvia miltiorrhiza, and the insight into their transcriptional regulatory mechanisms is an essential process to increase their content in vivo. SmMYB36 has been found to have important regulatory functions in the synthesis of tanshinone and phenolic acid; paradoxically, its mechanism of action in S. miltiorrhiza is not clear. Here, we demonstrated that SmMYB36 functions as a promoter of tanshinones accumulation and a suppressor of phenolic acids through the generation of SmMYB36 overexpressed and chimeric SmMYB36-SRDX (EAR repressive domain) repressor hairy roots in combination with transcriptomic-metabolomic analysis. SmMYB36 directly down-regulate the key enzyme gene of primary metabolism, SmGAPC, up-regulate the tanshinones biosynthesis branch genes SmDXS2, SmGGPPS1, SmCPS1 and down-regulate the phenolic acids biosynthesis branch enzyme gene, SmRAS. Meanwhile, SmERF6, a positive regulator of tanshinone synthesis activating SmCPS1, was up-regulated and SmERF115, a positive regulator of phenolic acid biosynthesis activating SmRAS, was down-regulated. Furthermore, the seven acidic amino acids at the C-terminus of SmMYB36 are required for both self-activating domain and activation of target gene expression. As a consequence, this study contributes to reveal the potential relevance of transcription factors synergistically regulating the biosynthesis of tanshinone and phenolic acid.

Transcription factor network analysis identifies REST/NRSF as an intrinsic regulator of CNS regeneration in mice

The inability of neurons to regenerate long axons within the CNS is a major impediment to improving outcome after spinal cord injury, stroke, and other CNS insults. Recent advances have uncovered an intrinsic program that involves coordinate regulation by multiple transcription factors that can be manipulated to enhance growth in the peripheral nervous system. Here, we use a systems genomics approach to characterize regulatory relationships of regeneration-associated transcription factors, identifying RE1-Silencing Transcription Factor (REST; Neuron-Restrictive Silencer Factor, NRSF) as a predicted upstream suppressor of a pro-regenerative gene program associated with axon regeneration in the CNS. We validate our predictions using multiple paradigms, showing that mature mice bearing cell type-specific deletions of REST or expressing dominant-negative mutant REST show improved regeneration of the corticospinal tract and optic nerve after spinal cord injury and optic nerve crush, which is accompanied by upregulation of regeneration-associated genes in cortical motor neurons and retinal ganglion cells, respectively. These analyses identify a role for REST as an upstream suppressor of the intrinsic regenerative program in the CNS and demonstrate the utility of a systems biology approach involving integrative genomics and bio-informatics to prioritize hypotheses relevant to CNS repair.

Single-cell network biology characterizes cell type gene regulation for drug repurposing and phenotype prediction in Alzheimer's disease

Dysregulation of gene expression in Alzheimer's disease (AD) remains elusive, especially at the cell type level. Gene regulatory network, a key molecular mechanism linking transcription factors (TFs) and regulatory elements to govern gene expression, can change across cell types in the human brain and thus serve as a model for studying gene dysregulation in AD. However, AD-induced regulatory changes across brain cell types remains uncharted. To address this, we integrated single-cell multi-omics datasets to predict the gene regulatory networks of four major cell types, excitatory and inhibitory neurons, microglia and oligodendrocytes, in control and AD brains. Importantly, we analyzed and compared the structural and topological features of networks across cell types and examined changes in AD. Our analysis shows that hub TFs are largely common across cell types and AD-related changes are relatively more prominent in some cell types (e.g., microglia). The regulatory logics of enriched network motifs (e.g., feed-forward loops) further uncover cell type-specific TF-TF cooperativities in gene regulation. The cell type networks are also highly modular and several network modules with cell-type-specific expression changes in AD pathology are enriched with AD-risk genes. The further disease-module-drug association analysis suggests cell-type candidate drugs and their potential target genes. Finally, our network-based machine learning analysis systematically prioritized cell type risk genes likely involved in AD. Our strategy is validated using an independent dataset which showed that top ranked genes can predict clinical phenotypes (e.g., cognitive impairment) of AD with reasonable accuracy. Overall, this single-cell network biology analysis provides a comprehensive map linking genes, regulatory networks, cell types and drug targets and reveals cell-type gene dysregulation in AD.

A novel approach GRNTSTE to reconstruct gene regulatory interactions applied to a case study for rat pineal rhythm gene

Accurate inference and prediction of gene regulatory network are very important for understanding dynamic cellular processes. The large-scale time series genomics data are helpful to reveal the molecular dynamics and dynamic biological processes of complex biological systems. Firstly, we collected the time series data of the rat pineal gland tissue in the natural state according to a fixed sampling rate, and performed whole-genome sequencing. The large-scale time-series sequencing data set of rat pineal gland was constructed, which includes 480 time points, the time interval between adjacent time points is 3 min, and the sampling period is 24 h. Then, we proposed a new method of constructing gene expression regulatory network, named the gene regulatory network based on time series data and entropy transfer (GRNTSTE) method. The method is based on transfer entropy and large-scale time-series gene expression data to infer the causal regulatory relationship between genes in a data-driven mode. The comparative experiments prove that GRNTSTE has better performance than dynamical gene network inference with ensemble of trees (dynGENIE3) and SCRIBE, and has similar performance to TENET. Meanwhile, we proved that the performance of GRNTSTE is slightly lower than that of SINCERITIES method and better than other gene regulatory network construction methods in BEELINE framework, which is based on the BEELINE data set. Finally, the rat pineal rhythm gene expression regulatory network was constructed by us based on the GRNTSTE method, which provides an important reference for the study of the pineal rhythm mechanism, and is of great significance to the study of the pineal rhythm mechanism.

Refactoring transcription factors for metabolic engineering

Due to the ability to regulate target metabolic pathways globally and dynamically, metabolic regulation systems composed of transcription factors have been widely used in metabolic engineering and synthetic biology. This review introduced the categories, action principles, prediction strategies, and related databases of transcription factors. Then, the application of global transcription machinery engineering technology and the transcription factor-based biosensors and quorum sensing systems are overviewed. In addition, strategies for optimizing the transcriptional regulatory tools' performance by refactoring transcription factors are summarized. Finally, the current limitations and prospects of constructing various regulatory tools based on transcription factors are discussed. This review will provide theoretical guidance for the rational design and construction of transcription factor-based metabolic regulation systems.

Hierarchical transcription factor and regulatory network for drought response in Betula platyphylla

Although many genes and biological processes involved in abiotic stress response have been identified, how they are regulated remains largely unclear. Here, to study the regulatory mechanism of birch (Betula platyphylla) responding to drought induced by polyethylene glycol (PEG) 6000 (20%, w/v), a partial correlation coefficient-based algorithm for constructing gene regulatory network (GRN) was proposed, and a three-layer hierarchical GRN was constructed, including 68 transcription factors (TFs), and 252 structural genes. Totally, 1448 predicted regulatory relationships are included, and most of them are novel. The reliability of GRN was verified by ChIP-PCR and qRT-PCR based on transient transformation. About 55% of genes in the bottom layer of GRN could confer drought tolerance. We selected the two TFs, BpMADS11 and BpNAC090, from the up layer and characterized their function in drought tolerance. Overexpression of BpMADS11 and BpNAC090 both reduces electrolyte leakage, ROS and MDA contents, displaying increased drought tolerance than wild-type birch. According to this GRN, the important biological processes involved in drought were identified, including "signaling hormone pathways", "water transport", "regulation of stomatal movement" and "response to oxidative stress". This work indicated that BpERF017, BpAGL61 and BpNAC090 are the key upstream regulators in birch drought tolerance. Our data clearly revealed the upstream regulators and TF-DNA interaction regulate different biological processes to adapt drought stress.

Gene regulatory network inference in long-lived C. elegans reveals modular properties that are predictive of novel aging genes

We design a “wisdom-of-the-crowds” GRN inference pipeline and couple it to complex network analysis to understand the organizational principles governing gene regulation in long-lived glp-1/Notch Caenorhabditis elegans. The GRN has three layers (input, core, and output) and is topologically equivalent to bow-tie/hourglass structures prevalent among metabolic networks. To assess the functional importance of structural layers, we screened 80% of regulators and discovered 50 new aging genes, 86% with human orthologues. Genes essential for longevity—including ones involved in insulin-like signaling (ILS)—are at the core, indicating that GRN's structure is predictive of functionality. We used in vivo reporters and a novel functional network covering 5,497 genetic interactions to make mechanistic predictions. We used genetic epistasis to test some of these predictions, uncovering a novel transcriptional regulator, sup-37, that works alongside DAF-16/FOXO. We present a framework with predictive power that can accelerate discovery in C. elegans and potentially humans.

•

Gene-regulatory inference provides global network of long-lived animals

•

The large-scale topology of the network has an hourglass structure

•

Membership to the core of the hourglass is a good predictor of functionality

•

Discovered 50 novel aging genes, including sup-37, a DAF-16 dependent gene

The economy of chromosomal distances in bacterial gene regulation

In the transcriptional regulatory network (TRN) of a bacterium, the nodes are genes and a directed edge represents the action of a transcription factor (TF), encoded by the source gene, on the target gene. It is a condensed representation of a large number of biological observations and facts. Nonrandom features of the network are structural evidence of requirements for a reliable systemic function. For the bacterium Escherichia coli we here investigate the (Euclidean) distances covered by the edges in the TRN when its nodes are embedded in the real space of the circular chromosome. Our work is motivated by 'wiring economy' research in Computational Neuroscience and starts from two contradictory hypotheses: (1) TFs are predominantly employed for long-distance regulation, while local regulation is exerted by chromosomal structure, locally coordinated by the action of structural proteins. Hence long distances should often occur. (2) A large distance between the regulator gene and its target requires a higher expression level of the regulator gene due to longer reaching times and ensuing increased degradation (proteolysis) of the TF and hence will be evolutionarily reduced. Our analysis supports the latter hypothesis.

Chromatin-state barriers enforce an irreversible mammalian cell fate decision

Stem and progenitor cells have the capacity to balance self-renewal and differentiation. Hematopoietic myeloid progenitors replenish more than 25 billion terminally differentiated neutrophils every day under homeostatic conditions and can increase this output in response to stress or infection. At what point along the spectrum of maturation do progenitors lose capacity for self-renewal and become irreversibly committed to differentiation? Using a system of conditional myeloid development that can be toggled between self-renewal and differentiation, we interrogate determinants of this "point of no return" in differentiation commitment. Irreversible commitment is due primarily to loss of open regulatory site access and disruption of a positive feedback transcription factor activation loop. Restoration of the transcription factor feedback loop extends the window of cell plasticity and alters the point of no return. These findings demonstrate how the chromatin state enforces and perpetuates cell fate and identify potential avenues for manipulating cell identity.

Paths through the Yeast Regulatory Network in Different Physiological States

We analyse paths through the regulatory networks that control gene-expression patterns in Yeast, in five different physiological states: cell cycle, DNA damage, stress response, diauxic shift, and sporulation. The network in each state is specified as a directed graph, containing different sets of edges connecting pairs selected from a combined set of 1475 nodes. Each network contains some nodes that have no parents, and others that have no children. We call these, respectively, 'source' and 'sink' nodes. For each network we enumerate paths between source and sink nodes. In a previous paper [1], we defined, extracted and compared the neighbourhoods of each transcription factor in different physiological states, and how the system reconfigures itself. Here we compare the usage of nodes and edges by different networks, and how they are assembled into paths. The picture that emerges is that the networks are not disjoint but show substantial sharing of nodes and edges; however, they assemble these materials into different sets of paths. Four of the networks, other than the cell-cycle network, contain paths between only a small fraction (< 13%) of possible source-sink pairs. Although the cell-cycle network is not an outlier in terms of total number of nodes and edges, and number of sink nodes, it is very much an outlier in having a greater proportion of source-to-sink paths than the other networks.

Neighbourhoods in the yeast regulatory network in different physiological states

MOTIVATION

The gene expression regulatory network in yeast controls the selective implementation of the information contained in the genome sequence. We seek to understand how, in different physiological states, the network reconfigures itself to produce a different proteome.

RESULTS

This article analyses this reconfiguration, focussing on changes in the local structure of the network. In particular, we define, extract and compare the 1-neighbourhoods of each transcription factor, where a 1-neighbourhood of a node in a network is the minimal subgraph of the network containing all nodes connected to the central node by an edge. We report the similarities and differences in the topologies and connectivities of these neighbourhoods in five physiological states for which data are available: cell cycle, DNA damage, stress response, diauxic shift and sporulation. Based on our analysis, it seems apt to regard the components of the regulatory network as 'software', and the responses to changes in state, 'reprogramming'.

Decoupling transcription factor expression and activity enables dimmer switch gene regulation

Gene-regulatory networks achieve complex mappings of inputs to outputs through mechanisms that are poorly understood. We found that in the galactose-responsive pathway in Saccharomyces cerevisiae, the decision to activate the transcription of genes encoding pathway components is controlled independently from the expression level, resulting in behavior resembling that of a mechanical dimmer switch. This was not a direct result of chromatin regulation or combinatorial control at galactose-responsive promoters; rather, this behavior was achieved by hierarchical regulation of the expression and activity of a single transcription factor. Hierarchical regulation is ubiquitous, and thus dimmer switch regulation is likely a key feature of many biological systems. Dimmer switch gene regulation may allow cells to fine-tune their responses to multi-input environments on both physiological and evolutionary time scales.

Scaling of intrinsic noise in an autocratic reaction network

Biochemical reactions in living cells often produce stochastic trajectories due to the fluctuations of the finite number of the macromolecular species present inside the cell. A significant number of computational and theoretical studies have previously investigated stochasticity in small regulatory networks to understand its origin and regulation. At the systems level regulatory networks have been determined to be hierarchical resembling social networks. In order to determine the stochasticity in networks with hierarchical architecture, here we computationally investigated intrinsic noise in an autocratic reaction network in which only the upstream regulators regulate the downstream regulators. We studied the effects of the qualitative and quantitative nature of regulatory interactions on the stochasticity in the network. We established an unconventional scaling of noise with average abundance in which the noise passes through a minimum indicating that the network can be noisy both in the low and high abundance regimes. We determined that the bursty kinetics of the trajectories are responsible for such scaling. The scaling of noise remains intact for a mixed network that includes democratic subnetworks within the autocratic network.

In silico identification of key regulators instigating the pre-meiotic phase during respiration in Saccharomyces cerevisiae

Like higher eukaryotes, diploid MATa/MATα budding yeasts can undergo both mitosis and meiosis. Although the potential reason for their phase switching is elucidated by two consecutive processes, i.e. transition from fermentation (mitotic growth) to respiration in glucose-deficient media and then complete shift to meiotic phase in combined nitrogen- and glucose-starved media, the genomic interactions and regulatory cascade operating this drive remain elusive. Here, we aim to explore the regulatory cross-talk that mediates the phase transition. We have hypothesized that pre-growth in glucose-starved condition (yeast extract-peptone-acetate media) not only causes switch from fermentation to respiration but also prepares them for meiosis via a myriad of signaling events regulated by transcription factors (TFs). We have identified 23 putative TFs from integrated protein-protein interaction and gene regulatory network that were reconstructed from predicted and experimentally validated data. A total of six TFs (Xbp1p, Abf1p, Cbf1p, Ste12p, Reb1p and Gcn4p) are found to be highly connected in the network and involved in the cross-talk between respiration and cellular preparation for meiosis. We have identified Abf1p and Adr1p as the master regulators of the integrated network. This study in yeast will help to decipher the pre-meiotic initiation that occurs in higher eukaryotes.

Carbon Catabolite Repression Governs Diverse Physiological Processes and Development in Aspergillus nidulans

Carbon catabolite repression (CCR) is a common phenomenon of microorganisms that enable efficient utilization of carbon nutrients, critical for the fitness of microorganisms in the wild and for pathogenic species to cause infection. In most filamentous fungal species, the conserved transcription factor CreA/Cre1 mediates CCR. Previous studies demonstrated a primary function for CreA/Cre1 in carbon metabolism; however, the phenotype of creA/cre1 mutants indicated broader roles. The global function and regulatory mechanism of this wide-domain transcription factor has remained elusive. Here, we applied two powerful genomics methods (transcriptome sequencing and chromatin immunoprecipitation sequencing) to delineate the direct and indirect roles of Aspergillus nidulans CreA across diverse physiological processes, including secondary metabolism, iron homeostasis, oxidative stress response, development, N-glycan biosynthesis, unfolded protein response, and nutrient and ion transport. The results indicate intricate connections between the regulation of carbon metabolism and diverse cellular functions. Moreover, our work also provides key mechanistic insights into CreA regulation and identifies CreA as a master regulator controlling many transcription factors of different regulatory networks. The discoveries for this highly conserved transcriptional regulator in a model fungus have important implications for CCR in related pathogenic and industrial species. IMPORTANCE The ability to scavenge and use a wide range of nutrients for growth is crucial for microorganisms' survival in the wild. Carbon catabolite repression (CCR) is a transcriptional regulatory phenomenon of both bacteria and fungi to coordinate the expression of genes required for preferential utilization of carbon sources. Since carbon metabolism is essential for growth, CCR is central to the fitness of microorganisms. In filamentous fungi, CCR is mediated by the conserved transcription factor CreA/Cre1, whose function in carbon metabolism has been well established. However, the global roles and regulatory mechanism of CreA/Cre1 are poorly defined. This study uncovers the direct and indirect functions of CreA in the model organism Aspergillus nidulans over diverse physiological processes and development and provides mechanistic insights into how CreA controls different regulatory networks. The work also reveals an interesting functional divergence between filamentous fungal and yeast CreA/Cre1 orthologues.

How higher goals are constructed and collapse under stress: A hierarchical Bayesian control systems perspective

In this paper, we show that organisms can be modeled as hierarchical Bayesian control systems with small world and information bottleneck (bow-tie) network structure. Such systems combine hierarchical perception with hierarchical goal setting and hierarchical action control. We argue that hierarchical Bayesian control systems produce deep hierarchies of goal states, from which it follows that organisms must have some form of 'highest goals'. For all organisms, these involve internal (self) models, external (social) models and overarching (normative) models. We show that goal hierarchies tend to decompose in a top-down manner under severe and prolonged levels of stress. This produces behavior that favors short-term and self-referential goals over long term, social and/or normative goals. The collapse of goal hierarchies is universally accompanied by an increase in entropy (disorder) in control systems that can serve as an early warning sign for tipping points (disease or death of the organism). In humans, learning goal hierarchies corresponds to personality development (maturation). The failure of goal hierarchies to mature properly corresponds to personality deficits. A top-down collapse of such hierarchies under stress is identified as a common factor in all forms of episodic mental disorders (psychopathology). The paper concludes by discussing ways of testing these hypotheses empirically.

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.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

How directed is a directed network?

The trophic levels of nodes in directed networks can reveal their functional properties. Moreover, the trophic coherence of a network, defined in terms of trophic levels, is related to properties such as cycle structure, stability and percolation. The standard definition of trophic levels, however, borrowed from ecology, suffers from drawbacks such as requiring basal nodes, which limit its applicability. Here we propose simple improved definitions of trophic levels and coherence that can be computed on any directed network. We demonstrate how the method can identify node function in examples including ecosystems, supply chain networks, gene expression and global language networks. We also explore how trophic levels and coherence relate to other topological properties, such as non-normality and cycle structure, and show that our method reveals the extent to which the edges in a directed network are aligned in a global direction.

Biological Networks Regulating Cell Fate Choice Are Minimally Frustrated

Characterization of the differences between biological and random networks can reveal the design principles that enable the robust realization of crucial biological functions including the establishment of different cell types. Previous studies, focusing on identifying topological features that are present in biological networks but not in random networks, have, however, provided few functional insights. We use a Boolean modeling framework and ideas from the spin glass literature to identify functional differences between five real biological networks and random networks with similar topological features. We show that minimal frustration is a fundamental property that allows biological networks to robustly establish cell types and regulate cell fate choice, and that this property can emerge in complex networks via Darwinian evolution. The study also provides clues regarding how the regulation of cell fate choice can go awry in a disease like cancer and lead to the emergence of aberrant cell types.

Distinctions between sex and time in patterns of DNA methylation across puberty

BACKGROUND

There are significant sex differences in human physiology and disease; the genomic sources of these differences, however, are not well understood. During puberty, a drastic neuroendocrine shift signals physical changes resulting in robust sex differences in human physiology. Here, we explore how shifting patterns of DNA methylation may inform these pathways of biological plasticity during the pubertal transition. In this study we analyzed DNA methylation (DNAm) in saliva at two time points across the pubertal transition within the same individuals. Our purpose was to compare two domains of DNAm patterns that may inform processes of sexual differentiation 1) sex related sites, which demonstrated differences between males from females and 2) time related sites in which DNAm shifted significantly between timepoints. We further explored the correlated network structure sex and time related DNAm networks and linked these patterns to pubertal stage, assays of salivary testosterone, a reliable diagnostic of free, unbound hormone that is available to act on target tissues, and overlap with androgen response elements.

RESULTS

Sites that differed by biological sex were largely independent of sites that underwent change across puberty. Time-related DNAm sites, but not sex-related sites, formed correlated networks that were associated with pubertal stage. Both time and sex DNAm networks reflected salivary testosterone levels that were enriched for androgen response elements, with sex-related DNAm networks being informative of testosterone levels above and beyond biological sex later in the pubertal transition.

CONCLUSIONS

These results inform our understanding of the distinction between sex- and time-related differences in DNAm during the critical period of puberty and highlight a novel linkage between correlated patterns of sex-related DNAm and levels of salivary testosterone.

Exploring disease comorbidity in a module-module interaction network

Understanding disease comorbidity contributes to improved quality of life in patients who are suffering from multiple diseases. Therefore, to better explore comorbid diseases, the clarification of associations between diseases based on biological functions is essential. In our study, we propose a method for identifying disease comorbidity in a module-based network, named the module-module interaction (MMI) network, which represents how biological functions influence each other. To construct the MMI network, we detected gene modules - sets of genes that have a higher probability of taking part in specific functions - and established a link between these modules. Subsequently, we constructed disease-related networks in the MMI network to understand inherent disease mechanisms and calculated comorbidity scores of disease pairs using Gene Ontology (GO) terms. Our results show that we can obtain further information on disease mechanisms by considering interactions between functional modules instead of between genes. In addition, we verified that predicted comorbid relationships of disease pairs based on the MMI network are more significant than those based on the protein-protein interaction (PPI) network. This study can be useful to elucidate the mechanisms underlying comorbidities for further study, which will provide a broader insight into the pathogenesis of diseases.

Flux-based hierarchical organization of Escherichia coli's metabolic network

Biological networks across scales exhibit hierarchical organization that may constrain network function. Yet, understanding how these hierarchies arise due to the operational constraint of the networks and whether they impose limits to molecular phenotypes remains elusive. Here we show that metabolic networks include a hierarchy of reactions based on a natural flux ordering that holds for every steady state. We find that the hierarchy of reactions is reflected in experimental measurements of transcript, protein and flux levels of Escherichia coli under various growth conditions as well as in the catalytic rate constants of the corresponding enzymes. Our findings point at resource partitioning and a fine-tuning of enzyme levels in E. coli to respect the constraints imposed by the network structure at steady state. Since reactions in upper layers of the hierarchy impose an upper bound on the flux of the reactions downstream, the hierarchical organization of metabolism due to the flux ordering has direct applications in metabolic engineering.

Construction of a hierarchical gene regulatory network centered around a transcription factor

We have modified a multitude of transcription factors (TFs) in numerous plant species and some animal species, and obtained transgenic lines that exhibit phenotypic alterations. Whenever we observe phenotypic changes in a TF's transgenic lines, we are always eager to identify its target genes, collaborative regulators and even upstream high hierarchical regulators. This issue can be addressed by establishing a multilayered hierarchical gene regulatory network (ML-hGRN) centered around a given TF. In this article, a practical approach for constructing an ML-hGRN centered on a TF using a combined approach of top-down and bottom-up network construction methods is described. Strategies for constructing ML-hGRNs are vitally important, as these networks provide key information to advance our understanding of how biological processes are regulated.

Chromosomal origin of replication coordinates logically distinct types of bacterial genetic regulation

For a long time it has been hypothesized that bacterial gene regulation involves an intricate interplay of the transcriptional regulatory network (TRN) and the spatial organization of genes in the chromosome. Here we explore this hypothesis both on a structural and on a functional level. On the structural level, we study the TRN as a spatially embedded network. On the functional level, we analyze gene expression patterns from a network perspective (“digital control”), as well as from the perspective of the spatial organization of the chromosome (“analog control”). Our structural analysis reveals the outstanding relevance of the symmetry axis defined by the origin (Ori) and terminus (Ter) of replication for the network embedding and, thus, suggests the co-evolution of two regulatory infrastructures, namely the transcriptional regulatory network and the spatial arrangement of genes on the chromosome, to optimize the cross-talk between two fundamental biological processes: genomic expression and replication. This observation is confirmed by the functional analysis based on the differential gene expression patterns of more than 4000 pairs of microarray and RNA-Seq datasets for E. coli from the Colombos Database using complex network and machine learning methods. This large-scale analysis supports the notion that two logically distinct types of genetic control are cooperating to regulate gene expression in a complementary manner. Moreover, we find that the position of the gene relative to the Ori is a feature of very high predictive value for gene expression, indicating that the Ori–Ter symmetry axis coordinates the action of distinct genetic control mechanisms.

The hourglass organization of the Caenorhabditis elegans connectome

We approach the C. elegans connectome as an information processing network that receives input from about 90 sensory neurons, processes that information through a highly recurrent network of about 80 interneurons, and it produces a coordinated output from about 120 motor neurons that control the nematode's muscles. We focus on the feedforward flow of information from sensory neurons to motor neurons, and apply a recently developed network analysis framework referred to as the "hourglass effect". The analysis reveals that this feedforward flow traverses a small core ("hourglass waist") that consists of 10-15 interneurons. These are mostly the same interneurons that were previously shown (using a different analytical approach) to constitute the "rich-club" of the C. elegans connectome. This result is robust to the methodology that separates the feedforward from the feedback flow of information. The set of core interneurons remains mostly the same when we consider only chemical synapses or the combination of chemical synapses and gap junctions. The hourglass organization of the connectome suggests that C. elegans has some similarities with encoder-decoder artificial neural networks in which the input is first compressed and integrated in a low-dimensional latent space that encodes the given data in a more efficient manner, followed by a decoding network through which intermediate-level sub-functions are combined in different ways to compute the correlated outputs of the network. The core neurons at the hourglass waist represent the information bottleneck of the system, balancing the representation accuracy and compactness (complexity) of the given sensory information.

The negative cofactor 2 complex is a key regulator of drug resistance in Aspergillus fumigatus

The frequency of antifungal resistance, particularly to the azole class of ergosterol biosynthetic inhibitors, is a growing global health problem. Survival rates for those infected with resistant isolates are exceptionally low. Beyond modification of the drug target, our understanding of the molecular basis of azole resistance in the fungal pathogen Aspergillus fumigatus is limited. We reasoned that clinically relevant antifungal resistance could derive from transcriptional rewiring, promoting drug resistance without concomitant reductions in pathogenicity. Here we report a genome-wide annotation of transcriptional regulators in A. fumigatus and construction of a library of 484 transcription factor null mutants. We identify 12 regulators that have a demonstrable role in itraconazole susceptibility and show that loss of the negative cofactor 2 complex leads to resistance, not only to the azoles but also the salvage therapeutics amphotericin B and terbinafine without significantly affecting pathogenicity.

Resistance to primary treatments of invasive aspergillosis is growing. Here, the authors generate a knockout library for 484 transcription factors in Aspergillus fumigatus, and show that loss of the NCT complex leads to cross-resistance to all primary and some salvage therapeutics without affecting pathogenicity.

Determinants of target prioritization and regulatory hierarchy for the bacterial small RNA SgrS

Small RNA (sRNA) regulators promote efficient responses to stress, but the mechanisms for prioritizing target mRNA regulation remain poorly understood. This study examines mechanisms underlying hierarchical regulation by the sRNA SgrS, found in enteric bacteria and produced under conditions of metabolic stress. SgrS posttranscriptionally coordinates a nine-gene regulon to restore growth and homeostasis. An in vivo reporter system quantified SgrS-dependent regulation of target genes and established that SgrS exhibits a clear target preference. Regulation of some targets is efficient even at low SgrS levels, whereas higher SgrS concentrations are required to regulate other targets. In vivo and in vitro analyses revealed that RNA structure and the number and position of base pairing sites relative to the start of translation impact the efficiency of regulation of SgrS targets. The RNA chaperone Hfq uses distinct modes of binding to different SgrS mRNA targets, which differentially influences positive and negative regulation. The RNA degradosome plays a larger role in regulation of some SgrS targets compared to others. Collectively, our results suggest that sRNA selection of target mRNAs and regulatory hierarchy are influenced by several molecular features and that the combination of these features precisely tunes the efficiency of regulation of multi-target sRNA regulons.

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

A periodic table of cell types

Single cell biology is currently revolutionizing developmental and evolutionary biology, revealing new cell types and states in an impressive range of biological systems. With the accumulation of data, however, the field is grappling with a central unanswered question: what exactly is a cell type? This question is further complicated by the inherently dynamic nature of developmental processes. In this Hypothesis article, we propose that a 'periodic table of cell types' can be used as a framework for distinguishing cell types from cell states, in which the periods and groups correspond to developmental trajectories and stages along differentiation, respectively. The different states of the same cell type are further analogous to 'isotopes'. We also highlight how the concept of a periodic table of cell types could be useful for predicting new cell types and states, and for recognizing relationships between cell types throughout development and evolution.

Regulatory feedback loops bridge the human gene regulatory network and regulate carcinogenesis

The development of disease involves a systematic disturbance inside cells and is associated with changes in the interactions or regulations among genes forming biological networks. The bridges inside a network are critical in shortening the distances between nodes. We observed that, inside the human gene regulatory network, one strongly connected core bridged the whole network. Other regulations outside the core formed a weakly connected component surrounding the core like a peripheral structure. Furthermore, the regulatory feedback loops (FBLs) inside the core compose an interface-like structure between the core and periphery. We then denoted the regulatory FBLs as the interface core. Notably, both the cancer-associated and essential biomolecules and regulations were significantly overrepresented in the interface core. These results implied that the interface core is not only critical for the network structure but central in cellular systems. Furthermore, the enrichment of the cancer-associated and essential regulations in the interface core might be attributed to its bridgeness in the network. More importantly, we identified one regulatory FBL between HNF4A and NR2F2 that possesses the highest bridgeness in the interface core. Further investigation suggested that the disturbance of the HNF4A-NR2F2 FBL might protect tumor cells from apoptotic processes. Our results emphasize the relevance of the regulatory network properties to cellular systems and might reveal a critical role of the interface core in cancer.

Hierarchical cooperation of transcription factors from integration analysis of DNA sequences, ChIP-Seq and ChIA-PET data

Background

Chromosomal architecture, which is constituted by chromatin loops, plays an important role in cellular functions. Gene expression and cell identity can be regulated by the chromatin loop, which is formed by proximal or distal enhancers and promoters in linear DNA (1D). Enhancers and promoters are fundamental non-coding elements enriched with transcription factors (TFs) to form chromatin loops. However, the specific cooperation of TFs involved in forming chromatin loops is not fully understood.

Results

Here, we proposed a method for investigating the cooperation of TFs in four cell lines by the integrative analysis of DNA sequences, ChIP-Seq and ChIA-PET data. Results demonstrate that the interaction of enhancers and promoters is a hierarchical and dynamic complex process with cooperative interactions of different TFs synergistically regulating gene expression and chromatin structure. The TF cooperation involved in maintaining and regulating the chromatin loop of cells can be regulated by epigenetic factors, such as other TFs and DNA methylation.

Conclusions

Such cooperation among TFs provides the potential features that can affect chromatin’s 3D architecture in cells. The regulation of chromatin 3D organization and gene expression is a complex process associated with the hierarchical and dynamic prosperities of TFs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5535-2) contains supplementary material, which is available to authorized users.

Comparative Transcriptome Analysis Shows Conserved Metabolic Regulation during Production of Secondary Metabolites in Filamentous Fungi

Filamentous fungi possess great potential as sources of medicinal bioactive compounds, such as antibiotics, but efficient production is hampered by a limited understanding of how their metabolism is regulated. We investigated the metabolism of six secondary metabolite-producing fungi of the Penicillium genus during nutrient depletion in the stationary phase of batch fermentations and assessed conserved metabolic responses across species using genome-wide transcriptional profiling. A coexpression analysis revealed that expression of biosynthetic genes correlates with expression of genes associated with pathways responsible for the generation of precursor metabolites for secondary metabolism. Our results highlight the main metabolic routes for the supply of precursors for secondary metabolism and suggest that the regulation of fungal metabolism is tailored to meet the demands for secondary metabolite production. These findings can aid in identifying fungal species that are optimized for the production of specific secondary metabolites and in designing metabolic engineering strategies to develop high-yielding fungal cell factories for production of secondary metabolites. IMPORTANCE Secondary metabolites are a major source of pharmaceuticals, especially antibiotics. However, the development of efficient processes of production of secondary metabolites has proved troublesome due to a limited understanding of the metabolic regulations governing secondary metabolism. By analyzing the conservation in gene expression across secondary metabolite-producing fungal species, we identified a metabolic signature that links primary and secondary metabolism and that demonstrates that fungal metabolism is tailored for the efficient production of secondary metabolites. The insight that we provide can be used to develop high-yielding fungal cell factories that are optimized for the production of specific secondary metabolites of pharmaceutical interest.