ConnecTF: A platform to integrate transcription factor-gene interactions and validate regulatory networks

Fruit sugar hub: gene regulatory network associated with soluble solids content (SSC) in Prunus persica

Chilean peach growers have achieved worldwide recognition for their high-quality fruit products. Among the main factors influencing peach fruit quality, sweetness is pivotal for maintaining the market's competitiveness. Numerous studies have been conducted in different peach-segregating populations to unravel SSC regulation. However, different cultivars may also have distinct genetic conformation, and other factors, such as environmental conditions, can significantly impact SSC. Using a transcriptomic approach with a gene co-expression network analysis, we aimed to identify the regulatory mechanism that controls the sugar accumulation process in an 'O × N' peach population. This population was previously studied through genomic analysis, associating LG5 with the genetic control of the SSC trait. The results obtained in this study allowed us to identify 91 differentially expressed genes located on chromosome 5 of the peach genome as putative new regulators of sugar accumulation in peach, together with a regulatory network that involves genes directly associated with sugar transport (PpSWEET15), cellulose biosynthesis (PpCSLG2), flavonoid biosynthesis (PpPAL1), pectin modifications (PpPG, PpPL and PpPMEi), expansins (PpEXPA1 and PpEXPA8) and several transcription factors (PpC3H67, PpHB7, PpRVE1 and PpCBF4) involved with the SSC phenotype. These results contribute to a better understanding of the genetic control of the SSC trait for future breeding programs in peaches.

Voice from both sides: a molecular dialogue between transcriptional activators and repressors in seed-to-seedling transition and crop adaptation

High-quality seeds provide valuable nutrients to human society and ensure successful seedling establishment. During maturation, seeds accumulate storage compounds that are required to sustain seedling growth during germination. This review focuses on the epigenetic repression of the embryonic and seed maturation programs in seedlings. We begin with an extensive overview of mutants affecting these processes, illustrating the roles of core proteins and accessory components in the epigenetic machinery by comparing mutants at both phenotypic and molecular levels. We highlight how omics assays help uncover target-specific functional specialization and coordination among various epigenetic mechanisms. Furthermore, we provide an in-depth discussion on the Seed dormancy 4 (Sdr4) transcriptional corepressor family, comparing and contrasting their regulation of seed germination in the dicotyledonous species Arabidopsis and two monocotyledonous crops, rice and wheat. Finally, we compare the similarities in the activation and repression of the embryonic and seed maturation programs through a shared set of cis-regulatory elements and discuss the challenges in applying knowledge largely gained in model species to crops.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Bipartite networks represent causality better than simple networks: evidence, algorithms, and applications

A network, whose nodes are genes and whose directed edges represent positive or negative influences of a regulatory gene and its targets, is often used as a representation of causality. To infer a network, researchers often develop a machine learning model and then evaluate the model based on its match with experimentally verified "gold standard" edges. The desired result of such a model is a network that may extend the gold standard edges. Since networks are a form of visual representation, one can compare their utility with architectural or machine blueprints. Blueprints are clearly useful because they provide precise guidance to builders in construction. If the primary role of gene regulatory networks is to characterize causality, then such networks should be good tools of prediction because prediction is the actionable benefit of knowing causality. But are they? In this paper, we compare prediction quality based on "gold standard" regulatory edges from previous experimental work with non-linear models inferred from time series data across four different species. We show that the same non-linear machine learning models have better predictive performance, with improvements from 5.3% to 25.3% in terms of the reduction in the root mean square error (RMSE) compared with the same models based on the gold standard edges. Having established that networks fail to characterize causality properly, we suggest that causality research should focus on four goals: (i) predictive accuracy; (ii) a parsimonious enumeration of predictive regulatory genes for each target gene g; (iii) the identification of disjoint sets of predictive regulatory genes for each target g of roughly equal accuracy; and (iv) the construction of a bipartite network (whose node types are genes and models) representation of causality. We provide algorithms for all goals.

Identification of Novel Regulators of Leaf Senescence Using a Deep Learning Model

Deep learning has emerged as a powerful tool for investigating intricate biological processes in plants by harnessing the potential of large-scale data. Gene regulation is a complex process that transcription factors (TFs), cooperating with their target genes, participate in through various aspects of biological processes. Despite its significance, the study of gene regulation has primarily focused on a limited number of notable instances, leaving numerous aspects and interactions yet to be explored comprehensively. Here, we developed DEGRN (Deep learning on Expression for Gene Regulatory Network), an innovative deep learning model designed to decipher gene interactions by leveraging high-dimensional expression data obtained from bulk RNA-Seq and scRNA-Seq data in the model plant Arabidopsis. DEGRN exhibited a compared level of predictive power when applied to various datasets. Through the utilization of DEGRN, we successfully identified an extensive set of 3,053,363 high-quality interactions, encompassing 1430 TFs and 13,739 non-TF genes. Notably, DEGRN's predictive capabilities allowed us to uncover novel regulators involved in a range of complex biological processes, including development, metabolism, and stress responses. Using leaf senescence as an example, we revealed a complex network underpinning this process composed of diverse TF families, including bHLH, ERF, and MYB. We also identified a novel TF, named MAF5, whose expression showed a strong linear regression relation during the progression of senescence. The mutant maf5 showed early leaf decay compared to the wild type, indicating a potential role in the regulation of leaf senescence. This hypothesis was further supported by the expression patterns observed across four stages of leaf development, as well as transcriptomics analysis. Overall, the comprehensive coverage provided by DEGRN expands our understanding of gene regulatory networks and paves the way for further investigations into their functional implications.

Nitrogen sensing and regulatory networks: it's about time and space

A plant's response to external and internal nitrogen signals/status relies on sensing and signaling mechanisms that operate across spatial and temporal dimensions. From a comprehensive systems biology perspective, this involves integrating nitrogen responses in different cell types and over long distances to ensure organ coordination in real time and yield practical applications. In this prospective review, we focus on novel aspects of nitrogen (N) sensing/signaling uncovered using temporal and spatial systems biology approaches, largely in the model Arabidopsis. The temporal aspects span: transcriptional responses to N-dose mediated by Michaelis-Menten kinetics, the role of the master NLP7 transcription factor as a nitrate sensor, its nitrate-dependent TF nuclear retention, its "hit-and-run" mode of target gene regulation, and temporal transcriptional cascade identified by "network walking." Spatial aspects of N-sensing/signaling have been uncovered in cell type-specific studies in roots and in root-to-shoot communication. We explore new approaches using single-cell sequencing data, trajectory inference, and pseudotime analysis as well as machine learning and artificial intelligence approaches. Finally, unveiling the mechanisms underlying the spatial dynamics of nitrogen sensing/signaling networks across species from model to crop could pave the way for translational studies to improve nitrogen-use efficiency in crops. Such outcomes could potentially reduce the detrimental effects of excessive fertilizer usage on groundwater pollution and greenhouse gas emissions.

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana. A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

Insights into molecular links and transcription networks integrating drought stress and nitrogen signaling

Drought and the availability of nitrate, the predominant source of nitrogen (N) in agriculture, are major factors limiting plant growth and crop productivity. The dissection of the transcriptional networks' components integrating droght stress and nitrate responses provides valuable insights into how plants effectively balance stress response with growth programs. Recent evidence in Arabidopsis thaliana indicates that transcription factors (TFs) involved in abscisic acid (ABA) signaling affect N metabolism and nitrate responses, and reciprocally, components of nitrate signaling might affect ABA and drought gene responses. Advances in understanding regulatory circuits of nitrate and drought crosstalk in plant tissues empower targeted genetic modifications to enhance plant development and stress resistance, critical traits for optimizing crop yield and promoting sustainable agriculture.

KLF17 is an important regulatory component of the transcriptomic response of Atlantic salmon macrophages to Piscirickettsia salmonis infection

Piscirickettsia salmonis is the most important health problem facing Chilean Aquaculture. Previous reports suggest that P. salmonis can survive in salmonid macrophages by interfering with the host immune response. However, the relevant aspects of the molecular pathogenesis of P. salmonis have been poorly characterized. In this work, we evaluated the transcriptomic changes in macrophage-like cell line SHK-1 infected with P. salmonis at 24- and 48-hours post-infection (hpi) and generated network models of the macrophage response to the infection using co-expression analysis and regulatory transcription factor-target gene information. Transcriptomic analysis showed that 635 genes were differentially expressed after 24- and/or 48-hpi. The pattern of expression of these genes was analyzed by weighted co-expression network analysis (WGCNA), which classified genes into 4 modules of expression, comprising early responses to the bacterium. Induced genes included genes involved in metabolism and cell differentiation, intracellular transportation, and cytoskeleton reorganization, while repressed genes included genes involved in extracellular matrix organization and RNA metabolism. To understand how these expression changes are orchestrated and to pinpoint relevant transcription factors (TFs) controlling the response, we established a curated database of TF-target gene regulatory interactions in Salmo salar, SalSaDB. Using this resource, together with co-expression module data, we generated infection context-specific networks that were analyzed to determine highly connected TF nodes. We found that the most connected TF of the 24- and 48-hpi response networks is KLF17, an ortholog of the KLF4 TF involved in the polarization of macrophages to an M2-phenotype in mammals. Interestingly, while KLF17 is induced by P. salmonis infection, other TFs, such as NOTCH3 and NFATC1, whose orthologs in mammals are related to M1-like macrophages, are repressed. In sum, our results suggest the induction of early regulatory events associated with an M2-like phenotype of macrophages that drives effectors related to the lysosome, RNA metabolism, cytoskeleton organization, and extracellular matrix remodeling. Moreover, the M1-like response seems delayed in generating an effective response, suggesting a polarization towards M2-like macrophages that allows the survival of P. salmonis. This work also contributes to SalSaDB, a curated database of TF-target gene interactions that is freely available for the Atlantic salmon community.

Phloem-Mobile MYB44 Negatively Regulates Expression of PHOSPHATE TRANSPORTER 1 in Arabidopsis Roots

Phosphorus (P) is an essential plant macronutrient; however, its availability is often limited in soils. Plants have evolved complex mechanisms for efficient phosphate (Pi) absorption, which are responsive to changes in external and internal Pi concentration, and orchestrated through local and systemic responses. To explore these systemic Pi responses, here we identified AtMYB44 as a phloem-mobile mRNA, an Arabidopsis homolog of Cucumis sativus MYB44, that is responsive to the Pi-starvation stress. qRT-PCR assays revealed that AtMYB44 was up-regulated and expressed in both shoot and root in response to Pi-starvation stress. The atmyb44 mutant displayed higher shoot and root biomass compared to wild-type plants, under Pi-starvation conditions. Interestingly, the expression of PHOSPHATE TRANSPORTER1;2 (PHT1;2) and PHT1;4 was enhanced in atmyb44 in response to a Pi-starvation treatment. A split-root assay showed that AtMYB44 expression was systemically regulated under Pi-starvation conditions, and in atmyb44, systemic controls on PHT1;2 and PHT1;4 expression were moderately disrupted. Heterografting assays confirmed graft transmission of AtMYB44 transcripts, and PHT1;2 and PHT1;4 expression was decreased in heterografted atmyb44 rootstocks. Taken together, our findings support the hypothesis that mobile AtMYB44 mRNA serves as a long-distance Pi response signal, which negatively regulates Pi transport and utilization in Arabidopsis.

Gene regulatory network inference in soybean upon infection by Phytophthora sojae

Phytophthora sojae is a soil-borne oomycete and the causal agent of Phytophthora root and stem rot (PRR) in soybean (Glycine max [L.] Merrill). Yield losses attributed to P. sojae are devastating in disease-conducive environments, with global estimates surpassing 1.1 million tonnes annually. Historically, management of PRR has entailed host genetic resistance (both vertical and horizontal) complemented by disease-suppressive cultural practices (e.g., oomicide application). However, the vast expansion of complex and/or diverse P. sojae pathotypes necessitates developing novel technologies to attenuate PRR in field environments. Therefore, the objective of the present study was to couple high-throughput sequencing data and deep learning to elucidate molecular features in soybean following infection by P. sojae. In doing so, we generated transcriptomes to identify differentially expressed genes (DEGs) during compatible and incompatible interactions with P. sojae and a mock inoculation. The expression data were then used to select two defense-related transcription factors (TFs) belonging to WRKY and RAV families. DNA Affinity Purification and sequencing (DAP-seq) data were obtained for each TF, providing putative DNA binding sites in the soybean genome. These bound sites were used to train Deep Neural Networks with convolutional and recurrent layers to predict new target sites of WRKY and RAV family members in the DEG set. Moreover, we leveraged publicly available Arabidopsis (Arabidopsis thaliana) DAP-seq data for five TF families enriched in our transcriptome analysis to train similar models. These Arabidopsis data-based models were used for cross-species TF binding site prediction on soybean. Finally, we created a gene regulatory network depicting TF-target gene interactions that orchestrate an immune response against P. sojae. Information herein provides novel insight into molecular plant-pathogen interaction and may prove useful in developing soybean cultivars with more durable resistance to P. sojae.

A multiomics integrative analysis of color de-synchronization with softening of 'Hass' avocado fruit: A first insight into a complex physiological disorder

Exocarp color de-synchronization with softening of 'Hass' avocado is a relevant recurrent problem for the avocado supply chain. This study aimed to unravel the mechanisms driving this de-synchronization integrating omics datasets from avocado exocarp of different storage conditions and color phenotypes. In addition, we propose potential biomarkers to predict color synchronized/de-synchronized fruit. Integration of transcriptomics, proteomics and metabolomics and network analysis revealed eight transcription factors associated with differentially regulated genes between regular air (RA) and controlled atmosphere (CA) and twelve transcription factors related to avocado fruit color de-synchronization control in ready-to-eat stage. CA was positively correlated to auxins, ethylene, cytokinins and brassinosteroids-related genes, while RA was characterized by enrichment of cell wall remodeling and abscisic acid content associated genes. At ready-to-eat higher contents of flavonoids, abscisic acid and brassinosteroids were associated with color-softening synchronized avocados. In contrast, de-synchronized fruit revealed increases of jasmonic acid, salicylic acid and auxin levels.

A gene regulatory network in Arabidopsis roots reveals features and regulators of the plant response to elevated CO2

The elevation of CO₂ in the atmosphere increases plant biomass but decreases their mineral content. The genetic and molecular bases of these effects remain mostly unknown, in particular in the root system, which is responsible for plant nutrient uptake. To gain knowledge about the effect of elevated CO₂ on plant growth and physiology, and to identify its regulatory in the roots, we analyzed genome expression in Arabidopsis roots through a combinatorial design with contrasted levels of CO₂ , nitrate, and iron. We demonstrated that elevated CO₂ has a modest effect on root genome expression under nutrient sufficiency, but by contrast leads to massive expression changes under nitrate or iron deficiencies. We demonstrated that elevated CO₂ negatively targets nitrate and iron starvation modules at the transcriptional level, associated with a reduction in high-affinity nitrate uptake. Finally, we inferred a gene regulatory network governing the root response to elevated CO₂ . This network allowed us to identify candidate transcription factors including MYB15, WOX11, and EDF3 which we experimentally validated for their role in the stimulation of growth by elevated CO₂ . Our approach identified key features and regulators of the plant response to elevated CO₂ , with the objective of developing crops resilient to climate change.

Characterization of Gene Regulatory Networks in Plants Using New Methods and Data Types

A major question in plant biology is to understand how plant growth, development, and environmental responses are controlled and coordinated by the activities of regulatory factors. Gene regulatory network (GRN) analyses require integrated approaches that combine experimental approaches with computational analyses. A wide range of experimental approaches and tools are now available, such as targeted perturbation of gene activities, quantitative and cell-type specific measurements of dynamic gene activities, and systematic analysis of the molecular 'hard-wiring' of the systems. At the computational level, different tools and databases are available to study regulatory sequences, including intuitive visualizations to explore data-driven gene regulatory networks in different plant species. Furthermore, advanced data integration approaches have recently been developed to efficiently leverage complementary regulatory data types and learn context-specific networks.

AGENT for Exploring and Analyzing Gene Regulatory Networks from Arabidopsis

Gene regulatory networks (GRNs) are important for determining how an organism develops and how it responds to external stimuli. In the case of Arabidopsis thaliana, several GRNs have been identified covering many important biological processes. We present AGENT, the Arabidopsis GEne Network Tool, for exploring and analyzing published GRNs. Using tools in AGENT, regulatory motifs such as feed-forward loops can be easily identified. Nodes with high centrality-and hence importance-can likewise be identified. Gene expression data can also be overlaid onto GRNs to help discover subnetworks acting in specific tissues or under certain conditions.

Gene Regulatory Network Modeling Using Single-Cell Multi-Omics in Plants

Single-cell multi-omics technology can be applied to plant cells to characterize gene expression and open chromatin regions in individual cells. In this chapter, we describe a computational pipeline for the analysis of single-cell data to construct gene regulatory networks. The major steps of this pipeline include the following: (1) normalize and integrate scRNA-seq and scATAC-seq data (2) identify cluster maker genes (3) perform motif finding for selected marker genes, and (4) identify regulatory networks with machine learning. The pipeline has been tested using data from the model species Arabidopsis and is generally applicable to other plant and animal species to characterize regulatory networks using single-cell multi-omics data.

Building High-Confidence Gene Regulatory Networks by Integrating Validated TF-Target Gene Interactions Using ConnecTF

Many methods are now available to identify or predict the target genes of transcription factors (TFs) in plants. These include experimental approaches such as in vivo or in vitro TF-target gene-binding assays and various methods for identifying regulated targets in mutants, transgenics, or isolated plant cells. In addition, computational approaches are used to infer TF-target gene interactions from the regulatory elements or gene expression changes across treatments. While each of these approaches has now been applied to a large number of TFs from many species, each method has its own limitations which necessitates that multiple data types are integrated to build the most accurate representation of the gene regulatory networks operating in plants. To make the analyses of TF-target interaction datasets available to the broader research community, we have developed the ConnecTF web platform ( https://connectf.org/ ). In this chapter, we describe how ConnecTF can be used to integrate validated and predicted TF-target gene interactions in order to dissect the regulatory role of TFs in developmental and stress response pathways. Using as our examples KN1 and RA1, two well-characterized maize TFs involved in developing floral tissue, we demonstrate how ConnecTF can be used to (1) compare the target genes between TFs, (2) identify direct vs. indirect targets by combining TF-binding and TF-regulation datasets, (3) chart and visualize network paths between TFs and their downstream targets, and (4) prune inferred user networks for high-confidence predicted interactions using validated TF-target gene data. Finally, we provide instructions for setting up a private version of ConnecTF that enables research groups to store and analyze their own TF-target gene interaction datasets.

The TARGET System: Rapid Identification of Direct Targets of Transcription Factors by Gene Regulation in Plant Cells

The TARGET system allows for the rapid identification of direct regulated gene targets of transcription factors (TFs). It employs the transient transformation of plant protoplasts with inducible nuclear entry of the TF and subsequent transcriptomic and/or ChIP-seq analysis. The ability to separate direct TF-target gene regulatory interactions from indirect downstream responses and the significantly shorter amount of time required to perform the assay, compared to the generation of transgenics, make this plant cell-based approach a valuable tool for a higher throughput approach to identify the genome-wide targets of multiple TFs, to build validated transcriptional networks in plants. Here, we describe the use of the TARGET system in Arabidopsis seedling root protoplasts to map the gene regulatory network downstream of transcription factors-of-interest.

Transcriptome and Gene Regulatory Network Analyses Reveal New Transcription Factors in Mature Fruit Associated with Harvest Date in Prunus persica

Harvest date is a critical parameter for producers and consumers regarding agro-industrial performance. It involves a pleiotropic effect controlling the development of other fruit quality traits through finely controlling regulatory mechanisms. Fruit ripening is a process in which various signals and biological events co-occur and are regulated by hormone signaling that produces the accumulation/degradation of multiple compounds. However, the regulatory mechanisms that control the hormone signaling involved in fruit development and ripening are still unclear. To investigate the issue, we used individuals with early, middle and late harvest dates from a peach segregating population to identify regulatory candidate genes controlling fruit quality traits at the harvest stage and validate them in contrasting peach varieties for this trait. We identified 467 and 654 differentially expressed genes for early and late harvest through a transcriptomic approach. In addition, using the Arabidopsis DAP-seq database and network analysis, six transcription factors were selected. Our results suggest significant hormonal balance and cell wall composition/structure differences between early and late harvest samples. Thus, we propose that higher expression levels of the transcription factors HB7, ERF017 and WRKY70 in early harvest individuals would induce the expression of genes associated with the jasmonic acid pathway, photosynthesis and gibberellins inhibition. While on the other hand, the high expression levels of LHY, CDF3 and NAC083 in late harvest individuals would promote the induction of genes associated with abscisic acid biosynthesis, auxins and cell wall remodeling.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

The chromatin accessibility landscape of pistils and anthers in rice

Transcription activation is tightly associated with the openness of chromatin and allows direct contact between transcriptional regulators and their targeted DNA for gene expression. However, there are limited studies on the annotation of open chromatin regions (OCRs) in rice (Oryza sativa), especially those in reproductive organs. Here, we characterized OCRs in rice pistils and anthers with an assay for transposase-accessible chromatin using sequencing. Despite a large overlap, we found more OCRs in pistils than in anthers. These OCRs were enriched in gene transcription start sites (TSSs) and showed tight associations with gene expression. Transcription factor (TF) binding motifs were enriched at these OCRs as validated by TF chromatin immunoprecipitation followed by sequencing. Pistil-specific OCRs provided potential regulatory networks by binding directly to the targets, indicating that pistil-specific OCRs may be indicators of cis-regulatory elements in regulating pistil development, which are absent in anthers. We also found that open chromatin of pistils and anthers responded differently to low temperature (LT). These data offer a comprehensive overview of OCRs regulating reproductive organ development and LT responses in rice.

Validation of a high-confidence regulatory network for gene-to-NUE phenotype in field-grown rice

Nitrogen (N) and Water (W) - two resources critical for crop productivity - are becoming increasingly limited in soils globally. To address this issue, we aim to uncover the gene regulatory networks (GRNs) that regulate nitrogen use efficiency (NUE) - as a function of water availability - in Oryza sativa, a staple for 3.5 billion people. In this study, we infer and validate GRNs that correlate with rice NUE phenotypes affected by N-by-W availability in the field. We did this by exploiting RNA-seq and crop phenotype data from 19 rice varieties grown in a 2x2 N-by-W matrix in the field. First, to identify gene-to-NUE field phenotypes, we analyzed these datasets using weighted gene co-expression network analysis (WGCNA). This identified two network modules ("skyblue" & "grey60") highly correlated with NUE grain yield (NUEg). Next, we focused on 90 TFs contained in these two NUEg modules and predicted their genome-wide targets using the N-and/or-W response datasets using a random forest network inference approach (GENIE3). Next, to validate the GENIE3 TF→target gene predictions, we performed Precision/Recall Analysis (AUPR) using nine datasets for three TFs validated in planta. This analysis sets a precision threshold of 0.31, used to "prune" the GENIE3 network for high-confidence TF→target gene edges, comprising 88 TFs and 5,716 N-and/or-W response genes. Next, we ranked these 88 TFs based on their significant influence on NUEg target genes responsive to N and/or W signaling. This resulted in a list of 18 prioritized TFs that regulate 551 NUEg target genes responsive to N and/or W signals. We validated the direct regulated targets of two of these candidate NUEg TFs in a plant cell-based TF assay called TARGET, for which we also had in planta data for comparison. Gene ontology analysis revealed that 6/18 NUEg TFs - OsbZIP23 (LOC_Os02g52780), Oshox22 (LOC_Os04g45810), LOB39 (LOC_Os03g41330), Oshox13 (LOC_Os03g08960), LOC_Os11g38870, and LOC_Os06g14670 - regulate genes annotated for N and/or W signaling. Our results show that OsbZIP23 and Oshox22, known regulators of drought tolerance, also coordinate W-responses with NUEg. This validated network can aid in developing/breeding rice with improved yield on marginal, low N-input, drought-prone soils.

Transcription Factors Controlling the Rhizobium-Legume Symbiosis: Integrating Infection, Organogenesis and the Abiotic Environment

Legume roots engage in a symbiotic relationship with rhizobia, leading to the development of nitrogen-fixing nodules. Nodule development is a sophisticated process and is under the tight regulation of the plant. The symbiosis initiates with a signal exchange between the two partners, followed by the development of a new organ colonized by rhizobia. Over two decades of study have shed light on the transcriptional regulation of rhizobium-legume symbiosis. A large number of transcription factors (TFs) have been implicated in one or more stages of this symbiosis. Legumes must monitor nodule development amidst a dynamic physical environment. Some environmental factors are conducive to nodulation, whereas others are stressful. The modulation of rhizobium-legume symbiosis by the abiotic environment adds another layer of complexity and is also transcriptionally regulated. Several symbiotic TFs act as integrators between symbiosis and the response to the abiotic environment. In this review, we trace the role of various TFs involved in rhizobium-legume symbiosis along its developmental route and highlight the ones that also act as communicators between this symbiosis and the response to the abiotic environment. Finally, we discuss contemporary approaches to study TF-target interactions in plants and probe their potential utility in the field of rhizobium-legume symbiosis.

Plant-DTI: Extending the landscape of TF protein and DNA interaction in plants by a machine learning-based approach

As a sessile organism, plants hold elaborate transcriptional regulatory systems that allow them to adapt to variable surrounding environments. Current understanding of plant regulatory mechanisms is greatly constrained by limited knowledge of transcription factor (TF)-DNA interactions. To mitigate this problem, a Plant-DTI predictor (Plant DBD-TFBS Interaction) was developed here as the first machine-learning model that covered the largest experimental datasets of 30 plant TF families, including 7 plant-specific DNA binding domain (DBD) types, and their transcription factor binding sites (TFBSs). Plant-DTI introduced a novel TFBS feature construction, called TFBS base-preference, which enhanced the specificity of TFBS to DBD types. The proposed model showed better predictive performance with the TFBS base-preference than the simple binary representation. Plant-DTI was validated with 22 independent ChIP-seq datasets. It accurately predicted the measured DBD-TFBS pairs along with their TFBS motifs, and effectively predicted interactions of other TFs containing similar DBD types. Comparing to the existing state-of-art methods, Plant-DTI prediction showed a figure of merit in sensitivity and specificity with respect to the position weight matrix (PWM) and TSPTFBS methods. Finally, the proposed Plant-DTI model helped to fill the knowledge gap in the regulatory mechanisms of the cassava sucrose synthase 1 gene (MeSUS1). Plant-DTI predicted MeERF72 as a regulator of MeSUS1 in consistence with the yeast one-hybrid (Y1H) experiment. Taken together, Plant-DTI would help facilitate the prediction of TF-TFBS and TF-target gene (TG) interactions, thereby accelerating the study of transcriptional regulatory systems in plant species.

Single cell gene regulatory networks in plants: Opportunities for enhancing climate change stress resilience

Global warming poses major challenges for plant survival and agricultural productivity. Thus, efforts to enhance stress resilience in plants are key strategies for protecting food security. Gene regulatory networks (GRNs) are a critical mechanism conferring stress resilience. Until recently, predicting GRNs of the individual cells that make up plants and other multicellular organisms was impeded by aggregate population scale measurements of transcriptome and other genome-scale features. With the advancement of high-throughput single cell RNA-seq and other single cell assays, learning GRNs for individual cells is now possible, in principle. In this article, we report on recent advances in experimental and analytical methodologies for single cell sequencing assays especially as they have been applied to the study of plants. We highlight recent advances and ongoing challenges for scGRN prediction, and finally, we highlight the opportunity to use scGRN discovery for studying and ultimately enhancing abiotic stress resilience in plants.

Time-Based Systems Biology Approaches to Capture and Model Dynamic Gene Regulatory Networks

All aspects of transcription and its regulation involve dynamic events. However, capturing these dynamic events in gene regulatory networks (GRNs) offers both a promise and a challenge. The promise is that capturing and modeling the dynamic changes in GRNs will allow us to understand how organisms adapt to a changing environment. The ability to mount a rapid transcriptional response to environmental changes is especially important in nonmotile organisms such as plants. The challenge is to capture these dynamic, genome-wide events and model them in GRNs. In this review, we cover recent progress in capturing dynamic interactions of transcription factors with their targets-at both the local and genome-wide levels-and how they are used to learn how GRNs operate as a function of time. We also discuss recent advances that employ time-based machine learning approaches to forecast gene expression at future time points, a key goal of systems biology.

Inferring and analyzing gene regulatory networks from multi-factorial expression data: a complete and interactive suite

Background

High-throughput transcriptomic datasets are often examined to discover new actors and regulators of a biological response. To this end, graphical interfaces have been developed and allow a broad range of users to conduct standard analyses from RNA-seq data, even with little programming experience. Although existing solutions usually provide adequate procedures for normalization, exploration or differential expression, more advanced features, such as gene clustering or regulatory network inference, often miss or do not reflect current state of the art methodologies.

Results

We developed here a user interface called DIANE (Dashboard for the Inference and Analysis of Networks from Expression data) designed to harness the potential of multi-factorial expression datasets from any organisms through a precise set of methods. DIANE interactive workflow provides normalization, dimensionality reduction, differential expression and ontology enrichment. Gene clustering can be performed and explored via configurable Mixture Models, and Random Forests are used to infer gene regulatory networks. DIANE also includes a novel procedure to assess the statistical significance of regulator-target influence measures based on permutations for Random Forest importance metrics. All along the pipeline, session reports and results can be downloaded to ensure clear and reproducible analyses.

Conclusions

We demonstrate the value and the benefits of DIANE using a recently published data set describing the transcriptional response of Arabidopsis thaliana under the combination of temperature, drought and salinity perturbations. We show that DIANE can intuitively carry out informative exploration and statistical procedures with RNA-Seq data, perform model based gene expression profiles clustering and go further into gene network reconstruction, providing relevant candidate genes or signalling pathways to explore. DIANE is available as a web service (https://diane.bpmp.inrae.fr), or can be installed and locally launched as a complete R package.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07659-2).