Comparative analyses of stress-responsive genes in Arabidopsis thaliana: insight from genomic data mining, functional enrichment, pathway analysis and phenomics

Transcriptomics, proteomics, and metabolomics interventions prompt crop improvement against metal(loid) toxicity

The escalating challenges posed by metal(loid) toxicity in agricultural ecosystems, exacerbated by rapid climate change and anthropogenic pressures, demand urgent attention. Soil contamination is a critical issue because it significantly impacts crop productivity. The widespread threat of metal(loid) toxicity can jeopardize global food security due to contaminated food supplies and pose environmental risks, contributing to soil and water pollution and thus impacting the whole ecosystem. In this context, plants have evolved complex mechanisms to combat metal(loid) stress. Amid the array of innovative approaches, omics, notably transcriptomics, proteomics, and metabolomics, have emerged as transformative tools, shedding light on the genes, proteins, and key metabolites involved in metal(loid) stress responses and tolerance mechanisms. These identified candidates hold promise for developing high-yielding crops with desirable agronomic traits. Computational biology tools like bioinformatics, biological databases, and analytical pipelines support these omics approaches by harnessing diverse information and facilitating the mapping of genotype-to-phenotype relationships under stress conditions. This review explores: (1) the multifaceted strategies that plants use to adapt to metal(loid) toxicity in their environment; (2) the latest findings in metal(loid)-mediated transcriptomics, proteomics, and metabolomics studies across various plant species; (3) the integration of omics data with artificial intelligence and high-throughput phenotyping; (4) the latest bioinformatics databases, tools and pipelines for single and/or multi-omics data integration; (5) the latest insights into stress adaptations and tolerance mechanisms for future outlooks; and (6) the capacity of omics advances for creating sustainable and resilient crop plants that can thrive in metal(loid)-contaminated environments.

Computational analysis of potential candidate genes involved in the cold stress response of ten Rosaceae members

Background

Plant species from Rosaceae family are economically important. One of the major environmental factors impacting those species is cold stress. Although several Rosaceae plant genomes have recently been sequenced, there have been very few research conducted on cold upregulated genes and their promoter binding sites. In this study, we used computational approaches to identify and analyse potential cold stress response genes across ten Rosaceae family members.

Results

Cold stress upregulated gene data from apple and strawberry were used to identify syntelogs in other Rosaceae species. Gene duplication analysis was carried out to better understand the distribution of these syntelog genes in different Rosaceae members. A total of 11,145 popular abiotic stress transcription factor-binding sites were identified in the upstream region of these potential cold-responsive genes, which were subsequently categorised into distinct transcription factor (TF) classes. MYB classes of transcription factor binding site (TFBS) were abundant, followed by bHLH, WRKY, and AP2/ERF. TFBS patterns in the promoter regions were compared among these species and gene families, found to be quite different even amongst functionally related syntelogs. A case study on important cold stress responsive transcription factor family, AP2/ERF showed less conservation in TFBS patterns in the promoter regions. This indicates that syntelogs from the same group may be comparable at the gene level but not at the level of cis-regulatory elements. Therefore, for such genes from the same family, different repertoire of TFs could be recruited for regulation and expression. Duplication events must have played a significant role in the similarity of TFBS patterns amongst few syntelogs of closely related species.

Conclusions

Our study overall suggests that, despite being from the same gene family, different combinations of TFs may play a role in their regulation and expression. The findings of this study will provide information about potential genes involved in the cold stress response, which will aid future functional research of these gene families involved in many important biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08751-x.

Advances in "Omics" Approaches for Improving Toxic Metals/Metalloids Tolerance in Plants

Food safety has emerged as a high-urgency matter for sustainable agricultural production. Toxic metal contamination of soil and water significantly affects agricultural productivity, which is further aggravated by extreme anthropogenic activities and modern agricultural practices, leaving food safety and human health at risk. In addition to reducing crop production, increased metals/metalloids toxicity also disturbs plants' demand and supply equilibrium. Counterbalancing toxic metals/metalloids toxicity demands a better understanding of the complex mechanisms at physiological, biochemical, molecular, cellular, and plant level that may result in increased crop productivity. Consequently, plants have established different internal defense mechanisms to cope with the adverse effects of toxic metals/metalloids. Nevertheless, these internal defense mechanisms are not adequate to overwhelm the metals/metalloids toxicity. Plants produce several secondary messengers to trigger cell signaling, activating the numerous transcriptional responses correlated with plant defense. Therefore, the recent advances in omics approaches such as genomics, transcriptomics, proteomics, metabolomics, ionomics, miRNAomics, and phenomics have enabled the characterization of molecular regulators associated with toxic metal tolerance, which can be deployed for developing toxic metal tolerant plants. This review highlights various response strategies adopted by plants to tolerate toxic metals/metalloids toxicity, including physiological, biochemical, and molecular responses. A seven-(omics)-based design is summarized with scientific clues to reveal the stress-responsive genes, proteins, metabolites, miRNAs, trace elements, stress-inducible phenotypes, and metabolic pathways that could potentially help plants to cope up with metals/metalloids toxicity in the face of fluctuating environmental conditions. Finally, some bottlenecks and future directions have also been highlighted, which could enable sustainable agricultural production.

Functional characterization of a heterologously expressed Brassica napus WRKY41-1 transcription factor in regulating anthocyanin biosynthesis in Arabidopsis thaliana

Previous studies have shown that a plant WRKY transcription factor, WRKY41, has multiple functions, and regulates seed dormancy, hormone signaling pathways, and both biotic and abiotic stress responses. However, it is not known about the roles of AtWRKY41 from the model plant, Arabidopsis thaliana, and its ortholog, BnWRKY41, from the closely related and important oil-producing crop, Brassica napus, in the regulation of anthocyanin biosynthesis. Here, we found that the wrky41 mutation in A. thaliana resulted in a significant increase in anthocyanin levels in rosette leaves, indicating that AtWRKY41 acts as repressor of anthocyanin biosynthesis. RNA sequencing and quantitative real-time PCR analysis revealed increased expression of three regulatory genes AtMYB75, AtMYB111, and AtMYBD, and two structural genes, AT1G68440 and AtGSTF12, all of which contribute to anthocyanin biosynthesis, in the sixth rosette leaves of wrky41-2 plants at 20 days after germination. We cloned the full length complementary DNA of BnWRKY41-1 from the C2 subgenome of the B. napus genotype Westar and observed that, when overexpressed in tobacco leaves as a fusion protein with green fluorescent protein, BnWRKY41-1 is localized to the nucleus. We further showed that overexpression of BnWRKY41-1 in the A. thaliana wrky41-2 mutant rescued the higher anthocyanin content phenotype in rosette leaves of the mutant. Moreover, the elevated expression levels in wrky41-2 rosette leaves of several important regulatory and structural genes regulating anthocyanin biosynthesis were not observed in the BnWRKY41-1 overexpressing lines. These results reveal that BnWRKY41-1 has a similar role with AtWRKY41 in regulating anthocyanin biosynthesis when overexpressed in A. thaliana. This gene represents a promising target for genetically manipulating B. napus to increase the amounts of anthocyanins in rosette leaves.

Nataraja K. An Approach to Function Annotation for Proteins of Unknown Function (PUFs) in the Transcriptome of Indian Mulberry

The modern sequencing technologies are generating large volumes of information at the transcriptome and genome level. Translation of this information into a biological meaning is far behind the race due to which a significant portion of proteins discovered remain as proteins of unknown function (PUFs). Attempts to uncover the functional significance of PUFs are limited due to lack of easy and high throughput functional annotation tools. Here, we report an approach to assign putative functions to PUFs, identified in the transcriptome of mulberry, a perennial tree commonly cultivated as host of silkworm. We utilized the mulberry PUFs generated from leaf tissues exposed to drought stress at whole plant level. A sequence and structure based computational analysis predicted the probable function of the PUFs. For rapid and easy annotation of PUFs, we developed an automated pipeline by integrating diverse bioinformatics tools, designated as PUFs Annotation Server (PUFAS), which also provides a web service API (Application Programming Interface) for a large-scale analysis up to a genome. The expression analysis of three selected PUFs annotated by the pipeline revealed abiotic stress responsiveness of the genes, and hence their potential role in stress acclimation pathways. The automated pipeline developed here could be extended to assign functions to PUFs from any organism in general. PUFAS web server is available at http://caps.ncbs.res.in/pufas/ and the web service is accessible at http://capservices.ncbs.res.in/help/pufas.

Transcriptional regulatory networks in Arabidopsis thaliana during single and combined stresses

Differentially evolved responses to various stress conditions in plants are controlled by complex regulatory circuits of transcriptional activators, and repressors, such as transcription factors (TFs). To understand the general and condition-specific activities of the TFs and their regulatory relationships with the target genes (TGs), we have used a homogeneous stress gene expression dataset generated on ten natural ecotypes of the model plant Arabidopsis thaliana, during five single and six combined stress conditions. Knowledge-based profiles of binding sites for 25 stress-responsive TF families (187 TFs) were generated and tested for their enrichment in the regulatory regions of the associated TGs. Condition-dependent regulatory sub-networks have shed light on the differential utilization of the underlying network topology, by stress-specific regulators and multifunctional regulators. The multifunctional regulators maintain the core stress response processes while the transient regulators confer the specificity to certain conditions. Clustering patterns of transcription factor binding sites (TFBS) have reflected the combinatorial nature of transcriptional regulation, and suggested the putative role of the homotypic clusters of TFBS towards maintaining transcriptional robustness against cis-regulatory mutations to facilitate the preservation of stress response processes. The Gene Ontology enrichment analysis of the TGs reflected sequential regulation of stress response mechanisms in plants.

Plant phenomics and the need for physiological phenotyping across scales to narrow the genotype-to-phenotype knowledge gap

Plants are affected by complex genome×environment×management interactions which determine phenotypic plasticity as a result of the variability of genetic components. Whereas great advances have been made in the cost-efficient and high-throughput analyses of genetic information and non-invasive phenotyping, the large-scale analyses of the underlying physiological mechanisms lag behind. The external phenotype is determined by the sum of the complex interactions of metabolic pathways and intracellular regulatory networks that is reflected in an internal, physiological, and biochemical phenotype. These various scales of dynamic physiological responses need to be considered, and genotyping and external phenotyping should be linked to the physiology at the cellular and tissue level. A high-dimensional physiological phenotyping across scales is needed that integrates the precise characterization of the internal phenotype into high-throughput phenotyping of whole plants and canopies. By this means, complex traits can be broken down into individual components of physiological traits. Since the higher resolution of physiological phenotyping by 'wet chemistry' is inherently limited in throughput, high-throughput non-invasive phenotyping needs to be validated and verified across scales to be used as proxy for the underlying processes. Armed with this interdisciplinary and multidimensional phenomics approach, plant physiology, non-invasive phenotyping, and functional genomics will complement each other, ultimately enabling the in silico assessment of responses under defined environments with advanced crop models. This will allow generation of robust physiological predictors also for complex traits to bridge the knowledge gap between genotype and phenotype for applications in breeding, precision farming, and basic research.

Genome-wide transcriptomic comparison of cotton (Gossypium herbaceum) leaf and root under drought stress

In this study, the 454 pyrosequencing platform was used for analyzing the comparative transcriptomic profiles of leaf and root tissues of 1-month-old cotton (Gossypium herbaceum) plants under drought stress. A total of 56,354 and 49,308 reads were generated from leaf and root tissues, respectively, and clustered into 6,313 and 5,858 unigenes. The differentially expressed unigenes that showed up-regulation (≥2-fold) or down-regulation (2≤-fold) were considered for further analysis. A total of 3,517 unigenes were differentially expressed in both tissues. The 1,528 genes specific to leaves and 1,128 specific to roots were obtained. The 28 biological pathways in two tissues were found to respond significantly to drought stress. A total of 289 in leaf and 277 in root unknown (novel) unigenes were found to be remarkably regulated by drought stress. Some key regulatory genes involved in abiotic stress such as WRKY, ERF, AP2 EREBP, MYB, and LEA were highly expressed in leaves. The genes RHD3, LBD, and transcription factor WRKY75, known for root development under various stress conditions, were expressed specifically in root. The genes related to chlorophyll a/b binding protein and photosystem-related proteins showed significant higher expression in roots and as compared to leaves. It can be concluded that cotton leaves are distinct from roots in terms of molecular mechanisms for responses to drought stress.

Genes Upregulated in Winter Wheat (Triticum aestivum L.) during Mild Freezing and Subsequent Thawing Suggest Sequential Activation of Multiple Response Mechanisms

Exposing fully cold-acclimated wheat plants to a mild freeze-thaw cycle of -3 °C for 24h followed by +3 °C for 24 or 48 h results in dramatically improved tolerance of subsequent exposure to sub-freezing temperatures. Gene enrichment analysis of crown tissue from plants collected before or after the -3 °C freeze or after thawing at +3 °C for 24 or 48 h revealed that many biological processes and molecular functions were activated during the freeze-thaw cycle in an increasing cascade of responses such that over 150 processes or functions were significantly enhanced by the end of the 48 h, post-freeze thaw. Nearly 2,000 individual genes were upregulated more than 2-fold over the 72 h course of freezing and thawing, but more than 70% of these genes were upregulated during only one of the time periods examined, suggesting a series of genes and gene functions were involved in activation of the processes that led to enhanced freezing tolerance. This series of functions appeared to include extensive cell signaling, activation of stress response mechanisms and the phenylpropanoid biosynthetic pathway, extensive modification of secondary metabolites, and physical restructuring of cell membranes. By identifying plant lines that are especially able to activate these multiple mechanisms it may be possible to develop lines with enhanced winterhardiness.

POEAS: Automated Plant Phenomic Analysis Using Plant Ontology

Biological enrichment analysis using gene ontology (GO) provides a global overview of the functional role of genes or proteins identified from large-scale genomic or proteomic experiments. Phenomic enrichment analysis of gene lists can provide an important layer of information as well as cellular components, molecular functions, and biological processes associated with gene lists. Plant phenomic enrichment analysis will be useful for performing new experiments to better understand plant systems and for the interpretation of gene or proteins identified from high-throughput experiments. Plant ontology (PO) is a compendium of terms to define the diverse phenotypic characteristics of plant species, including plant anatomy, morphology, and development stages. Adoption of this highly useful ontology is limited, when compared to GO, because of the lack of user-friendly tools that enable the use of PO for statistical enrichment analysis. To address this challenge, we introduce Plant Ontology Enrichment Analysis Server (POEAS) in the public domain. POEAS uses a simple list of genes as input data and performs enrichment analysis using Ontologizer 2.0 to provide results in two levels, enrichment results and visualization utilities, to generate ontological graphs that are of publication quality. POEAS also offers interactive options to identify user-defined background population sets, various multiple-testing correction methods, different enrichment calculation methods, and resampling tests to improve statistical significance. The availability of such a tool to perform phenomic enrichment analyses using plant genes as a complementary resource will permit the adoption of PO-based phenomic analysis as part of analytical workflows. POEAS can be accessed using the URL http://caps.ncbs.res.in/poeas.