Dynamic modeling of transcriptional gene regulatory network uncovers distinct pathways during the onset of Arabidopsis leaf senescence

Integrative systems biology framework discovers common gene regulatory signatures in multiple mechanistically distinct inflammatory skin diseases

More than 20% of the population across the world is affected by non-communicable inflammatory skin diseases including psoriasis, atopic dermatitis, hidradenitis suppurativa, rosacea, etc. Many of these chronic diseases are painful and debilitating with limited effective therapeutic interventions. However, recent advances in psoriasis treatment have improved the effectiveness and provide better management of the disease. This study aims to identify common regulatory pathways and master regulators that regulate molecular pathogenesis. We designed an integrative systems biology framework to identify the significant regulators across several inflammatory skin diseases. With conventional transcriptome analysis, we identified 55 shared genes, which are enriched in several immune-associated pathways in eight inflammatory skin diseases. Next, we exploited the gene co-expression-, and protein-protein interaction-based networks to identify shared genes and protein components in different diseases with relevant functional implications. Additionally, the network analytics unravels 55 high-value proteins as significant regulators in molecular pathogenesis. We believe that these significant regulators should be explored with critical experimental approaches to identify the putative drug targets for more effective treatments. As an example, we identified IKZF1 as a shared significant master regulator in three inflammatory skin diseases, which can serve as a putative drug target with known disease-derived molecules for developing efficacious combinatorial treatments for hidradenitis suppurativa, atopic dermatitis, and rosacea. The proposed framework is very modular, which can indicate a significant path of molecular mechanism-based drug development from complex transcriptomics data and other multi-omics data.

Integrated Systems Biology Pipeline to Compare Co-Expression Networks in Plants and Elucidate Differential Regulators

To identify sets of genes that exhibit similar expression characteristics, co-expression networks were constructed from transcriptome datasets that were obtained from plant samples at various stages of growth and development or treated with diverse biotic, abiotic, and other environmental stresses. In addition, co-expression network analysis can provide deeper insights into gene regulation when combined with transcriptomics. The coordination and integration of all these complex networks to deduce gene regulation are major challenges for plant biologists. Python and R have emerged as major tools for managing complex scientific data over the past decade. In this study, we describe a reproducible protocol POTFUL (pant co-expression transcription factor regulators), implemented in Python 3, for integrating co-expression and transcription factor target protein networks to infer gene regulation.

Transcriptome meta-analysis of abiotic stresses-responsive genes and identification of candidate transcription factors for broad stress tolerance in wheat

Under field conditions, wheat is subjected to single or multiple stress conditions. The elucidation of the molecular mechanism of stress response is a key step to identify candidate genes for stress resistance in plants. In this study, RNA-seq data analysis identified 17.324, 10.562, 5.510, and 8.653 differentially expressed genes (DEGs) under salt, drought, heat, and cold stress, respectively. Moreover, the comparison of DEGs from each stress revealed 2374 shared genes from which 40% showed highly conserved expression patterns. Moreover, co-expression network analysis and GO enrichment revealed co-expression modules enriched with genes involved in transcription regulation, protein kinase pathway, and genes responding to phytohormones or modulating hormone levels. The expression of 15 selected transcription factor encoding genes was analyzed under abiotic stresses and ABA treatment in durum wheat. The identified transcription factor genes are excellent candidates for genetic engineering of stress tolerance in wheat.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Cytotrap: An Innovative Approach for Protein-Protein Interaction Studies for Cytoplasmic Proteins

Protein-protein interaction mapping has gained immense importance in understanding protein functions in diverse biological pathways. There are various in vivo and in vitro techniques associated with the protein-protein interaction studies but generally, the focus is confined to understanding the protein interaction in the nucleus of the cell, and thus it limits the availability to explore protein interactions that are happening in the cytoplasm of the cell. Since posttranslational modification is a crucial step in signaling pathways and cellular protein interactions harnessing the cytoplasmic protein and evaluating the interaction in the cytoplasm, this protocol will provide more information about studying these types of protein interactions. Cytotrap is a type of yeast-two-hybrid system that differs in its ability to anchor along the membrane, thus directing the protein of interest to anchor along the membrane through the myristoylation signaling unit. The vector containing the target protein contains the myristoylation unit, called the prey, and the bait unit contains the protein of interest as a fusion with the hSos protein. In an event of interaction between the target and the protein of interest, the hSos protein unit will be localized to the membrane and the GDP/GTP exchange unit will trigger the activation of the Ras pathway that leads to the survival of the temperature-sensitive yeast strain at a higher temperature.

Dynamic Enrichment for Evaluation of Protein Networks (DEEPN): A High Throughput Yeast Two-Hybrid (Y2H) Protocol to Evaluate Networks

Proteins are the building blocks of life, and a vast array of cellular processes is handled by protein-protein interactions (PPIs). The protein complexes formed via PPIs lead to tangled networks that, with their continuous remodeling, build up systematic functional units. Over the years, PPIs have become an area of interest for many researchers, leading to the development of multiple in vitro and in vivo methods to reveal these interactions. The yeast-two-hybrid (Y2H) system is a potent genetic way to map PPIs in both a micro- and high-throughput manner. Y2H is a technique that involves using modified yeast cells to identify protein-protein interactions. For Y2H, the yeast cells are engineered only to grow when there is a significant interaction between a specific protein with its interacting partner. PPIs are identified in the Y2H system by stimulating reporter genes in response to a restored transcription factor. However, Y2H results may be constrained by stringency requirements, as the limited number of colony screenings through this technique could result in the possible elimination of numerous genuine interactions. Therefore, DEEPN (dynamic enrichment for evaluation of protein networks) can be used, offering the potential to study the multiple static and transient protein interactions in a single Y2H experiment. DEEPN utilizes next-generation DNA sequencing (NGS) data in a high-throughput manner and subsequently applies computational analysis and statistical modeling to identify interacting partners. This protocol describes customized reagents and protocols through which DEEPN analysis can be utilized efficiently and cost-effectively.

Protein-Protein Interaction Network Exploration Using Cytoscape

As the protein-protein interaction (PPI) data increase exponentially, the development and usage of computational methods to analyze these datasets have become a new research horizon in systems biology. The PPI network analysis and visualization can help identify functional modules of the network, pathway genes involved in common cellular functions, and functional annotations of novel genes. Currently, a variety of tools are available for network graph visualization and analysis. Cytoscape, an open-source software tool, is one of them. It provides an interactive visualization interface along with other core features to import, navigate, filter, cluster, search, and export networks. It comes with hundreds of in-built Apps in App Manager to resolve research questions related to network visualization and integration. This chapter aims to illustrate the Cytoscape application to visualize and analyze the PPI network using Arabidopsis interactome-1 main (AI-1MAIN) PPI network dataset from Plant Interactome Database.

Building Protein-Protein Interaction Graph Database Using Neo4j

A cell's various components interact with each other in a coordinated manner to respond to environmental cues and intracellular signals. Compared to the other biological networks, the protein-protein interaction (PPI) is mostly responsible for maintaining signaling pathways. Increasing numbers of experimentally verified and predicted PPIs in plants demand a scalable platform to deal with large and complex datasets. Network/graph data can be organized and analyzed using different tools. This chapter uses Neo4j, a graph database management system, to store and analyze plant PPI networks. To make the graph database and analyze network centrality, we used Arabidopsis interactome-1 main (AI-1MAIN) PPI network.

Protein-Protein Interaction Network Analysis Using NetworkX

In recent years, extracting information from biological data has become a particularly valuable way of gaining knowledge. Molecular interaction networks provide a framework for visualizing cellular processes, but their complexity frequently makes their interpretation difficult. Proteins are one of the primary determinants of biological function. Indeed, most biological activities in the living cells are functionally regulated by protein-protein interactions (PPIs). Thus, studying protein interactions is critical for understanding their roles within the cell. Exploring the PPI networks can open new avenues for future experimental studies and offer interspecies predictions for effective interaction mapping. In this chapter we will demonstrate how to construct, visualize, and analyze a protein-protein interaction network using NetworkX.

Preparation and Utilization of a Versatile GFP-Protein Trap-Like System for Protein Complex Immunoprecipitation in Plants

Protein complex immunoprecipitation (co-IP) is an in vitro technique used to study protein-protein interaction between two or more proteins. This method relies on affinity purification of recombinant epitope-tagged proteins followed by western blotting detection using tag-specific antibodies for the confirmation of positive interaction. The traditional co-IP method relies on the use of porous beaded support with immobilized antibodies to precipitate protein complexes. However, this method is time-consuming, labor-intensive, and provides lower reproducibility and yield of protein complexes. Here, we describe the implementation of magnetic beads and high-affinity anti-green fluorescent protein (GFP) antibodies to develop an in vitro GFP-protein trap-like system. This highly reproducible system utilizes a combination of small sample size, versatile lysis buffer, and lower amounts of magnetic beads to obtain protein complexes and aggregates that are compatible with functional assays, Western blotting, and mass spectrometry. In addition to protein-protein interactions, this versatile method can be employed to study protein-nucleic acid interactions. This protocol also highlights troubleshooting and includes recommendations to optimize its application.

A pipeline of integrating transcriptome and interactome to elucidate central nodes in host-pathogens interactions

Investigating the complexity of host-pathogen interactions is challenging. Here, we outline a pipeline to identify important proteins and signaling molecules in human-viral interactomes. Firstly, we curate a comprehensive human interactome. Subsequently, we infer viral targets and transcriptome-specific human interactomes (VTTSHI) for papillomavirus and herpes viruses by integrating viral targets and transcriptome data. Finally, we reveal the common and shared nodes and pathways in viral pathogenesis following network topology and pathway enrichment analyses.

For complete details on the use and execution of this protocol, please refer to Kumar et al. (2020).

•

Protocol for integrative analysis of transcriptome and proteome network data

•

Network subgroup enrichment of two host-pathogen interaction networks

•

Preprocessing of data and heterogeneous network integration

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

A genome-wide comparative evolutionary analysis of zinc finger-BED transcription factor genes in land plants

Zinc finger (Zf)-BED proteins are a novel superfamily of transcription factors that controls numerous activities in plants including growth, development, and cellular responses to biotic and abiotic stresses. Despite their important roles in gene regulation, little is known about the specific functions of Zf-BEDs in land plants. The current study identified a total of 750 Zf-BED-encoding genes in 35 land plant species including mosses, bryophytes, lycophytes, gymnosperms, and angiosperms. The gene family size was somewhat proportional to genome size. All identified genes were categorized into 22 classes based on their specific domain architectures. Of these, class I (Zf-BED_DUF-domain_Dimer_Tnp_hAT) was the most common in the majority of the land plants. However, some classes were family-specific, while the others were species-specific, demonstrating diversity at different classification levels. In addition, several novel functional domains were also predicated including WRKY and nucleotide-binding site (NBS). Comparative genomics, transcriptomics, and proteomics provided insights into the evolutionary history, duplication, divergence, gene gain and loss, species relationship, expression profiling, and structural diversity of Zf-BEDs in land plants. The comprehensive study of Zf-BEDs in Gossypium sp., (cotton) also demonstrated a clear footprint of polyploidization. Overall, this comprehensive evolutionary study of Zf-BEDs in land plants highlighted significant diversity among plant species.

A Rice Protein Interaction Network Reveals High Centrality Nodes and Candidate Pathogen Effector Targets

Network science identifies key players in diverse biological systems including host-pathogen interactions. We demonstrated a scale-free network property for a comprehensive rice protein–protein interactome (RicePPInets) that exhibits nodes with increased centrality indices. While weighted k-shell decomposition was shown efficacious to predict pathogen effector targets in Arabidopsis, we improved its computational code for a broader implementation on large-scale networks including RicePPInets. We determined that nodes residing within the internal layers of RicePPInets are poised to be the most influential, central, and effective information spreaders. To identify central players and modules through network topology analyses, we integrated RicePPInets and co-expression networks representing susceptible and resistant responses to strains of the bacterial pathogens Xanthomonas oryzae pv. oryzae and X. oryzae pv. oryzicola (Xoc) and generated a RIce-Xanthomonas INteractome (RIXIN). This revealed that previously identified candidate targets of pathogen transcription activator-like (TAL) effectors are enriched in nodes with enhanced connectivity, bottlenecks, and information spreaders that are located in the inner layers of the network, and these nodes are involved in several important biological processes. Overall, our integrative multi-omics network-based platform provides a potentially useful approach to prioritizing candidate pathogen effector targets for functional validation, suggesting that this computational framework can be broadly translatable to other complex pathosystems.

Specificity of H2O2 signaling in leaf senescence: is the ratio of H2O2 contents in different cellular compartments sensed in Arabidopsis plants?

Leaf senescence is an integral part of plant development and is driven by endogenous cues such as leaf or plant age. Developmental senescence aims to maximize the usage of carbon, nitrogen and mineral resources for growth and/or for the sake of the next generation. This requires efficient reallocation of the resources out of the senescing tissue into developing parts of the plant such as new leaves, fruits and seeds. However, premature senescence can be induced by severe and long-lasting biotic or abiotic stress conditions. It serves as an exit strategy to guarantee offspring in an unfavorable environment but is often combined with a trade-off in seed number and quality. In order to coordinate the very complex process of developmental senescence with environmental signals, highly organized networks and regulatory cues have to be in place. Reactive oxygen species, especially hydrogen peroxide (H2O2), are involved in senescence as well as in stress signaling. Here, we want to summarize the role of H2O2 as a signaling molecule in leaf senescence and shed more light on how specificity in signaling might be achieved. Altered hydrogen peroxide contents in specific compartments revealed a differential impact of H2O2 produced in different compartments. Arabidopsis lines with lower H2O2 levels in chloroplasts and cytoplasm point to the possibility that not the actual contents but the ratio between the two different compartments is sensed by the plant cells.

Network biology to uncover functional and structural properties of the plant immune system

In the last two decades, advances in network science have facilitated the discovery of important systems' entities in diverse biological networks. This graph-based technique has revealed numerous emergent properties of a system that enable us to understand several complex biological processes including plant immune systems. With the accumulation of multiomics data sets, the comprehensive understanding of plant-pathogen interactions can be achieved through the analyses and efficacious integration of multidimensional qualitative and quantitative relationships among the components of hosts and their microbes. This review highlights comparative network topology analyses in plant-pathogen co-expression networks and interactomes, outlines dynamic network modeling for cell-specific immune regulatory networks, and discusses the new frontiers of single-cell sequencing as well as multiomics data integration that are necessary for unraveling the intricacies of plant immune systems.

Transcriptional circuitry atlas of genetic diverse unstimulated murine and human macrophages define disparity in population-wide innate immunity

Macrophages are ubiquitous custodians of tissues, which play decisive role in maintaining cellular homeostasis through regulatory immune responses. Within tissues, macrophage exhibit extremely heterogeneous population with varying functions orchestrated through regulatory response, which can be further exacerbated in diverse genetic backgrounds. Gene regulatory networks (GRNs) offer comprehensive understanding of cellular regulatory behavior by unfolding the transcription factors (TFs) and regulated target genes. RNA-Seq coupled with ATAC-Seq has revolutionized the regulome landscape influenced by gene expression modeling. Here, we employ an integrative multi-omics systems biology-based analysis and generated GRNs derived from the unstimulated bone marrow-derived macrophages of five inbred genetically defined murine strains, which are reported to be linked with most of the population-wide human genetic variants. Our probabilistic modeling of a basal hemostasis pan regulatory repertoire in diverse macrophages discovered 96 TFs targeting 6279 genes representing 468,291 interactions across five inbred murine strains. Subsequently, we identify core and distinctive GRN sub-networks in unstimulated macrophages to describe the system-wide conservation and dissimilarities, respectively across five murine strains. Our study concludes that discrepancies in unstimulated macrophage-specific regulatory networks not only drives the basal functional plasticity within genetic backgrounds, additionally aid in understanding the complexity of racial disparity among the human population during stress.

Inference of Gene Regulatory Network from Single-Cell Transcriptomic Data Using pySCENIC

With the advent of recent next-generation sequencing (NGS) technologies in genomics, transcriptomics, and epigenomics, profiling single-cell sequencing became possible. The single-cell RNA sequencing (scRNA-seq) is widely used to characterize diverse cell populations and ascertain cell type-specific regulatory mechanisms. The gene regulatory network (GRN) mainly consists of genes and their regulators-transcription factors (TF). Here, we describe the lightning-fast Python implementation of the SCENIC (Single-Cell reEgulatory Network Inference and Clustering) pipeline called pySCENIC. Using single-cell RNA-seq data, it maps TFs onto gene regulatory networks and integrates various cell types to infer cell-specific GRNs. There are two fast and efficient GRN inference algorithms, GRNBoost2 and GENIE3, optionally available with pySCENIC. The pipeline has three steps: (1) identification of potential TF targets based on co-expression; (2) TF-motif enrichment analysis to identify the direct targets (regulons); and (3) scoring the activity of regulons (or other gene sets) on single cell types.

Dynamic Regulatory Event Mining by iDREM in Large-Scale Multi-omics Datasets During Biotic and Abiotic Stress in Plants

The system-wide complexity of genome regulation encoding the organism phenotypic diversity is well understood. However, a major challenge persists about the appropriate method to describe the systematic dynamic genome regulation event utilizing enormous multi-omics datasets. Here, we describe Interactive Dynamic Regulatory Events Miner (iDREM) which reconstructs gene-regulatory networks from temporal transcriptome, proteome, and epigenome datasets during stress to envisage "master" regulators by simulating cascades of temporal transcription-regulatory and interactome events. The iDREM is a Java-based software that integrates static and time-series transcriptomics and proteomics datasets, transcription factor (TF)-target interactions, microRNA (miRNA)-target interaction, and protein-protein interactions to reconstruct temporal regulatory network and identify significant regulators in an unsupervised manner. The hidden Markov model detects specialized manipulated pathways as well as genes to recognize statistically significant regulators (TFs/miRNAs) that diverge in temporal activity. This method can be translated to any biotic or abiotic stress in plants and animals to predict the master regulators from condition-specific multi-omics datasets including host-pathogen interactions for comprehensive understanding of manipulated biological pathways.

Dark-Induced Senescence Causes Localized Changes in DNA Methylation

Senescence occurs in a programmed manner to dismantle the vegetative tissues and redirect nutrients towards metabolic pathways supporting reproductive success. External factors can trigger the senescence program as an adaptive strategy, indicating that this terminal program is controlled at different levels. It has been proposed that epigenetic factors accompany the reprogramming of the senescent genome; however, the mechanism and extent of this reprogramming remain unknown. Using bisulphite conversion followed by sequencing, we assessed changes in the methylome of senescent Arabidopsis (Arabidopsis thaliana) leaves induced by darkness and monitored their effect on gene and transposable element (TE) expression with transcriptome sequencing. Upon dark-induced senescence, genes controlling chromatin silencing were collectively down-regulated. As a consequence, the silencing of TEs was impaired, causing in particular young TEs to become preferentially reactivated. In parallel, heterochromatin at chromocenters was decondensed. Despite the disruption of the chromatin maintenance network, the global DNA methylation landscape remained highly stable, with localized changes mainly restricted to CHH methylation. Together, our data show that the terminal stage of plant life is accompanied by global changes in chromatin structure but only localized changes in DNA methylation, adding another example of the dynamics of DNA methylation during plant development.