An Arabidopsis expression predictor enables inference of transcriptional regulators for gene modules

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.

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.

Key Pathways and Genes of Arabidopsis thaliana and Arabidopsis halleri Roots under Cadmium Stress Responses: Differences and Similarities

Cadmium (Cd) is among the world's major health concerns, as it renders soils unsuitable and unsafe for food and feed production. Phytoremediation has the potential to remediate Cd-polluted soils, but efforts are still needed to develop a deep understanding of the processes underlying it. In this study, we performed a comprehensive analysis of the root response to Cd stress in A. thaliana, which can phytostabilize Cd, and in A. halleri, which is a Cd hyperaccumulator. Suitable RNA-seq data were analyzed by WGCNA to identify modules of co-expressed genes specifically associated with Cd presence. The results evidenced that the genes of the hyperaccumulator A. halleri mostly associated with the Cd presence are finely regulated (up- and downregulated) and related to a general response to chemical and other stimuli. Additionally, in the case of A. thaliana, which can phytostabilize metals, the genes upregulated during Cd stress are related to a general response to chemical and other stimuli, while downregulated genes are associated with functions which, affecting root growth and development, determine a deep modification of the organ both at the cellular and physiological levels. Furthermore, key genes of the Cd-associated modules were identified and confirmed by differentially expressed gene (DEG) detection and external knowledge. Together, key functions and genes shed light on differences and similarities among the strategies that the plants use to cope with Cd and may be considered as possible targets for future research.

Decoding transcriptional regulation via a human gene expression predictor

Transcription factors (TFs) regulate cellular activities by controlling gene expression, but a predictive model describing how TFs quantitatively modulate human transcriptomes is lacking. We construct a universal human gene expression predictor and utilize it to decode transcriptional regulation. Using the expression of 1613 TFs, the predictor reconstitutes highly accurate transcriptomes for samples derived from a wide range of tissues and conditions. The broad applicability of the predictor indicates that it recapitulates the quantitative relationships between TFs and target genes ubiquitous across tissues. Significant interacting TF-target gene pairs are extracted from the predictor and enable downstream inference of TF regulators for diverse pathways involved in development, immunity, metabolism, and stress response. A detailed analysis of the hematopoiesis process reveals an atlas of key TFs regulating the development of different hematopoietic cell lineages, and a portion of these TFs are conserved between humans and mice. The results demonstrate that our method is capable of delineating the TFs responsible for fate determination. Compared to other existing tools, our approach shows better performance in recovering the correct TF regulators. Thus, we present a novel approach that can be used to study human transcriptional regulation in general.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

Rice co-expression network analysis identifies gene modules associated with agronomic traits

Identifying trait-associated genes is critical for rice (Oryza sativa) improvement, which usually relies on map-based cloning, quantitative trait locus analysis, or genome-wide association studies. Here we show that trait-associated genes tend to form modules within rice gene co-expression networks, a feature that can be exploited to discover additional trait-associated genes using reverse genetics. We constructed a rice gene co-expression network based on the graphical Gaussian model using 8,456 RNA-seq transcriptomes, which assembled into 1,286 gene co-expression modules functioning in diverse pathways. A number of the modules were enriched with genes associated with agronomic traits, such as grain size, grain number, tiller number, grain quality, leaf angle, stem strength, and anthocyanin content, and these modules are considered to be trait-associated gene modules. These trait-associated gene modules can be used to dissect the genetic basis of rice agronomic traits and to facilitate the identification of trait genes. As an example, we identified a candidate gene, OCTOPUS-LIKE 1 (OsOPL1), a homolog of the Arabidopsis (Arabidopsis thaliana) OCTOPUS gene, from a grain size module and verified it as a regulator of grain size via functional studies. Thus, our network represents a valuable resource for studying trait-associated genes in rice.

Feature importance network reveals novel functional relationships between biological features in Arabidopsis thaliana

Understanding how the different cellular components are working together to form a living cell requires multidisciplinary approaches combining molecular and computational biology. Machine learning shows great potential in life sciences, as it can find novel relationships between biological features. Here, we constructed a dataset of 11,801 gene features for 31,522 Arabidopsis thaliana genes and developed a machine learning workflow to identify linked features. The detected linked features are visualised as a Feature Important Network (FIN), which can be mined to reveal a variety of novel biological insights pertaining to gene function. We demonstrate how FIN can be used to generate novel insights into gene function. To make this network easily accessible to the scientific community, we present the FINder database, available at finder.plant.tools.

EXPLICIT-Kinase: A gene expression predictor for dissecting the functions of the Arabidopsis kinome

Protein kinases regulate virtually all cellular processes, but it remains challenging to determine the functions of all protein kinases, collectively called the "kinome", in any species. We developed a computational approach called EXPLICIT-Kinase to predict the functions of the Arabidopsis kinome. Because the activities of many kinases can be regulated transcriptionally, their gene expression patterns provide clues to their functions. A universal gene expression predictor for Arabidopsis was constructed to predict the expression of 30,172 non-kinase genes based on the expression of 994 kinases. The model reconstituted highly accurate transcriptomes for diverse Arabidopsis samples. It identified the significant kinases as predictor kinases for predicting the expression of Arabidopsis genes and pathways. Strikingly, these predictor kinases were often regulators of related pathways, as exemplified by those involved in cytokinesis, tissue development, and stress responses. Comparative analyses revealed that portions of these predictor kinases are shared and conserved between Arabidopsis and maize. As an example, we identified a conserved predictor kinase, RAF6, from a stomatal movement module. We verified that RAF6 regulates stomatal closure. It can directly interact with SLAC1, a key anion channel for stomatal closure, and modulate its channel activity. Our approach enables a systematic dissection of the functions of the Arabidopsis kinome.

The Arabidopsis gene co-expression network

Identifying genes that interact to confer a biological function to an organism is one of the main goals of functional genomics. High-throughput technologies for assessment and quantification of genome-wide gene expression patterns have enabled systems-level analyses to infer pathways or networks of genes involved in different functions under many different conditions. Here, we leveraged the publicly available, information-rich RNA-Seq datasets of the model plant Arabidopsis thaliana to construct a gene co-expression network, which was partitioned into clusters or modules that harbor genes correlated by expression. Gene ontology and pathway enrichment analyses were performed to assess functional terms and pathways that were enriched within the different gene modules. By interrogating the co-expression network for genes in different modules that associate with a gene of interest, diverse functional roles of the gene can be deciphered. By mapping genes differentially expressing under a certain condition in Arabidopsis onto the co-expression network, we demonstrate the ability of the network to uncover novel genes that are likely transcriptionally active but prone to be missed by standard statistical approaches due to their falling outside of the confidence zone of detection. To our knowledge, this is the first A. thaliana co-expression network constructed using the entire mRNA-Seq datasets (>20,000) available at the NCBI SRA database. The developed network can serve as a useful resource for the Arabidopsis research community to interrogate specific genes of interest within the network, retrieve the respective interactomes, decipher gene modules that are transcriptionally altered under certain condition or stage, and gain understanding of gene functions.