Genome-Wide Transcriptional Profiling and Metabolic Analysis Uncover Multiple Molecular Responses of the Grass Species Lolium perenne Under Low-Intensity Xenobiotic Stress

Deciphering Macromolecular Interactions Involved in Abiotic Stress Signaling: A Review of Bioinformatics Analysis

Plant functioning and responses to abiotic stresses largely involve regulations at the transcriptomic level via complex interactions of signal molecules, signaling cascades, and regulators. Nevertheless, all the signaling networks involved in responses to abiotic stresses have not yet been fully established. The in-depth analysis of transcriptomes in stressed plants has become a relevant state-of-the-art methodology to study these regulations and signaling pathways that allow plants to cope with or attempt to survive abiotic stresses. The plant science and molecular biology community has developed databases about genes, proteins, protein-protein interactions, protein-DNA interactions and ontologies, which are valuable sources of knowledge for deciphering such regulatory and signaling networks. The use of these data and the development of bioinformatics tools help to make sense of transcriptomic data in specific contexts, such as that of abiotic stress signaling, using functional biological approaches. The aim of this chapter is to present and assess some of the essential online tools and resources that will allow novices in bioinformatics to decipher transcriptomic data in order to characterize the cellular processes and functions involved in abiotic stress responses and signaling. The analysis of case studies further describes how these tools can be used to conceive signaling networks on the basis of transcriptomic data. In these case studies, particular attention was paid to the characterization of abiotic stress responses and signaling related to chemical and xenobiotic stressors.

Adaptive and maladaptive expression plasticity underlying herbicide resistance in an agricultural weed

Plastic phenotypic responses to environmental change are common, yet we lack a clear understanding of the fitness consequences of these plastic responses. Here, we use the evolution of herbicide resistance in the common morning glory (Ipomoea purpurea) as a model for understanding the relative importance of adaptive and maladaptive gene expression responses to herbicide. Specifically, we compare leaf gene expression changes caused by herbicide to the expression changes that evolve in response to artificial selection for herbicide resistance. We identify a number of genes that show plastic and evolved responses to herbicide and find that for the majority of genes with both plastic and evolved responses, plastic responses appear to be adaptive. We also find that selection for herbicide response increases gene expression plasticity. Overall, these results show the importance of adaptive plasticity for herbicide resistance in a common weed and that expression changes in response to strong environmental change can be adaptive.

Transcriptome analysis reveals gene responses to herbicide, tribenuron methyl, in Brassica napus L. during seed germination

Background

Tribenuron methyl (TBM) is an herbicide that inhibits sulfonylurea acetolactate synthase (ALS) and is one of the most widely used broad-leaved herbicides for crop production. However, soil residues or drifting of the herbicide spray might affect the germination and growth of rapeseed, Brassica napus, so it is imperative to understand the response mechanism of rape to TBM during germination. The aim of this study was to use transcriptome analysis to reveal the gene responses in herbicide-tolerant rapeseed to TBM stress during seed germination.

Results

2414, 2286, and 1068 differentially expressed genes (DEGs) were identified in TBM-treated resistant vs sensitive lines, treated vs. control sensitive lines, treated vs. control resistant lines, respectively. GO analysis showed that most DEGs were annotated to the oxidation-reduction pathways and catalytic activity. KEGG enrichment was mainly involved in plant-pathogen interactions, α-linolenic acid metabolism, glucosinolate biosynthesis, and phenylpropanoid biosynthesis. Based on GO and KEGG enrichment, a total of 137 target genes were identified, including genes involved in biotransferase activity, response to antioxidant stress and lipid metabolism. Biotransferase genes, CYP450, ABC and GST, detoxify herbicide molecules through physical or biochemical processes. Antioxidant genes, RBOH, WRKY, CDPK, MAPK, CAT, and POD regulate plant tolerance by transmitting ROS signals and triggering antioxidant enzyme expression. Lipid-related genes and hormone-related genes were also found, such as LOX3, ADH1, JAZ6, BIN2 and ERF, and they also played an important role in herbicide resistance.

Conclusions

This study provides insights for selecting TBM-tolerant rapeseed germplasm and exploring the molecular mechanism of TBM tolerance during germination.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07614-1.

A comparison of shared patterns of differential gene expression and gene ontologies in response to water-stress in roots and leaves of four diverse genotypes of Lolium and Festuca spp. temperate pasture grasses

Ryegrasses (Lolium spp.) and fescues (Festuca spp.) are closely related and widely cultivated perennial forage grasses. As such, resilience in the face of abiotic stresses is an important component of their traits. We have compared patterns of differentially expressed genes (DEGs) in roots and leaves of two perennial ryegrass genotypes and a single genotype of each of a festulolium (predominantly Italian ryegrass) and meadow fescue with the onset of water stress, focussing on overall patterns of DEGs and gene ontology terms (GOs) shared by all four genotypes. Plants were established in a growing medium of vermiculite watered with nutrient solution. Leaf and root material were sampled at 35% (saturation) and, as the medium dried, at 15%, 5% and 1% estimated water contents (EWCs) and RNA extracted. Differential gene expression was evaluated comparing the EWC sampling points from RNAseq data using a combination of analysis methods. For all genotypes, the greatest numbers of DEGs were identified in the 35/1 and 5/1 comparisons in both leaves and roots. In total, 566 leaf and 643 root DEGs were common to all 4 genotypes, though a third of these leaf DEGs were not regulated in the same up/down direction in all 4 genotypes. For roots, the equivalent figure was 1% of the DEGs. GO terms shared by all four genotypes were often enriched by both up- and down-regulated DEGs in the leaf, whereas generally, only by either up- or down-regulated DEGs in the root. Overall, up-regulated leaf DEGs tended to be more genotype-specific than down-regulated leaf DEGs or root DEGs and were also associated with fewer GOs. On average, only 5–15% of the DEGs enriching common GO terms were shared by all 4 genotypes, suggesting considerable variation in DEGs between related genotypes in enacting similar biological processes.

Local-scale dynamics of plant-pesticide interactions in a northern Brittany agricultural landscape

Soil pollution by anthropogenic chemicals is a major concern for sustainability of crop production and of ecosystem functions mediated by natural plant biodiversity. Understanding the complex effects of soil pollution requires multi-level and multi-scale approaches. Non-target and agri-environmental plant communities of field margins and vegetative filter strips are confronted with agricultural xenobiotics through soil contamination, drift, run-off and leaching events that result from chemical applications. Plant-pesticide dynamics in vegetative filter strips was studied at field scale in the agricultural landscape of a long-term ecological research network in northern Brittany (France). Vegetative filter strips effected significant pesticide abatement between the field and riparian compartments. However, comparison of pesticide usage modalities and soil chemical analysis revealed the extent and complexity of pesticide persistence in fields and vegetative filter strips, and suggested the contribution of multiple sources (yearly carry-over, interannual persistence, landscape-scale contamination). In order to determine the impact of such persistence, plant dynamics was followed in experimentally-designed vegetative filter strips of identical initial composition (Agrostis stolonifera, Anthemis tinctoria/Cota tinctoria, Centaurea cyanus, Fagopyrum esculentum, Festuca rubra, Lolium perenne, Lotus corniculatus, Phleum pratense, Trifolium pratense). After homogeneous vegetation establishment, experimental vegetative filter strips underwent rapid changes within the following two years, with Agrostis stolonifera, Festuca rubra, Lolium perenne and Phleum pratense becoming dominant and with the establishment of spontaneous vegetation. Co-inertia analysis showed that plant dynamics and soil residual pesticides could be significantly correlated, with the triazole fungicide epoxiconazole, the imidazole fungicide prochloraz and the neonicotinoid insecticide thiamethoxam as strong drivers of the correlation. However, the correlation was vegetative-filter-strip-specific, thus showing that correlation between plant dynamics and soil pesticides likely involved additional factors, such as threshold levels of residual pesticides. This situation of complex interactions between plants and soil contamination is further discussed in terms of agronomical, environmental and health issues.

Metabolomic and Transcriptomic Investigation of Metabolic Perturbations in Oryza sativa L. Triggered by Three Pesticides

Inappropriate application of pesticides often triggers molecular alterations in crops, which inadvertently disturbs metabolites and finally affects crop quality. Therefore, understanding the mechanism of action of pesticides on crops is essential for evaluating the potential environmental impact of pesticides. Our findings indicated that three typical pesticides, including herbicide butachlor, insecticide chlorpyrifos, and fungicide tricyclazole, induced the expression regulation of different key genes, exhibiting considerable distinction on metabolic responses in rice (Oryza sativa L.). Butachlor mainly affected five carbohydrate metabolism pathways (38.5%), and more than 48.0% of differentially expressed genes (DEGs) were involved in the starch and sucrose metabolism as well as photosynthesis, thereby disturbing the distribution of starch-sucrose. Chlorpyrifos dramatically affected six amino acid metabolism pathways (60.0%), and key DEGs mainly enriched in the aspartate and glutamate metabolism, inducing an increase in free amino acid contents (up to 29.02% of the control) and degradation of soluble proteins (down to 48.72% of the control). Tricyclazole remarkably affected six fatty acid metabolism pathways (53.9%) and significantly upregulated DEGs which primarily code oil-body membrane proteins, which resulted in the decline of saturated fatty acids (palmitic acid and stearic acid) and the increase of unsaturated fatty acids (linolenic acid and octadecadienoic acid). These findings provide a molecular-scale perspective on the response of crops to pesticides.

Species- and organ-specific responses of agri-environmental plants to residual agricultural pollutants

Soil pollution by anthropogenic chemicals is a major concern for sustainability of crop production and of ecosystem functions mediated by natural plant biodiversity. The complex effects on plants are however difficult to apprehend. Plant communities of field margins, vegetative filter strips or rotational fallows are confronted with agricultural pollutants through residual soil contamination and/or through drift, run-off and leaching events that result from chemical applications. Exposure to xenobiotics and heavy metals causes biochemical, physiological and developmental effects. However, the range of doses, modalities of exposure, metabolization of contaminants into derived xenobiotics, and combinations of contaminants result in variable and multi-level effects. Understanding these complex plant-pollutant interactions cannot directly rely on toxicological or agronomical approaches that focus on the effects of field-rate pesticide applications. It must take into account exposure at root level, sublethal concentrations of bioactive compounds and functional biodiversity of the plant species that are affected. The present study deals with agri-environmental plant species of field margins, vegetative filter strips or rotational fallows in European agricultural landscapes. Root and shoot physiological and growth responses were compared under controlled conditions that were optimally adjusted for each plant species. Contrasted responses of growth inhibition, no adverse effect or growth enhancement depended on species, organ and nature of contaminant. However, all of the agricultural contaminants under study (pesticides, pesticide metabolites, heavy metals, polycyclic aromatic hydrocarbons) had significant effects under conditions of sublethal exposure on at least some of the plant species. The fungicide tebuconazole and polycyclic aromatic hydrocarbon fluoranthene, which gave highest levels of responses, induced both activation or inhibition effects, in different plant species or in different organs of the same plant species. These complex effects are discussed in terms of dynamics of agri-environmental plants and of ecological consequences of differential root-shoot growth under conditions of soil contamination.

Non-target site-based resistance to tribenuron-methyl and essential involved genes in Myosoton aquaticum (L.)

BACKGROUND

Water chickweed (Myosoton aquaticum (L.)) is a dicot broadleaf weed that is widespread in winter fields in China, and has evolved serious resistance to acetolactate synthase (ALS) inhibiting herbicides.

RESULTS

We identified a M. aquaticum population exhibiting moderate (6.15-fold) resistance to tribenuron-methyl (TM). Target-site ALS gene sequencing revealed no known resistance mutations in these plants, and the in vitro ALS activity assays showed no differences in enzyme sensitivity between susceptible and resistant populations; however, resistance was reversed by pretreatment with the cytochrome P450 (CYP) monooxygenase inhibitor malathion. An RNA sequencing transcriptome analysis was performed to identify candidate genes involved in metabolic resistance, and the unigenes obtained by de novo transcriptome assembly were annotated across seven databases. In total, 34 differentially expressed genes selected by digital gene expression analysis were validated by quantitative real-time (qRT)-PCR. Ten consistently overexpressed contigs, including four for CYP, four for ATP-binding cassette (ABC) transporter, and two for peroxidase were further validated by qRT-PCR using additional plants from resistant and susceptible populations. Three CYP genes (with homology to CYP734A1, CYP76C1, and CYP86B1) and one ABC transporter gene (with homology to ABCC10) were highly expressed in all resistant plants.

CONCLUSION

The mechanism of TM resistance in M. aquaticum is controlled by NTSR rather than TSR. Four genes, CYP734A1, CYP76C1, CYP86B1, and ABCC10 could play essential role in metabolic resistance to TM and justify further functional studies. To our knowledge, this is the first large-scale transcriptome analysis of genes associated with NTSR in M. aquaticum using the Illumina platform. Our data provide resource for M. aquaticum biology, and will facilitate the study of herbicide resistance mechanism at the molecular level in this species as well as in other weeds.

Enhanced Herbicide Metabolism and Metabolic Resistance Genes Identified in Tribenuron-Methyl Resistant Myosoton aquaticum L

The evolved resistance of Myosoton aquaticum L. to acetolactate synthase (ALS) inhibitors is well established, but most research has focused on target-site resistance, while nontarget-site resistance remains neglected. Here, we investigated mechanisms of the latter. The pretreatment with the P450 inhibitor malathion significantly increased the sensitivity of resistant plants to tribenuron-methyl. The rapid P450-mediated tribenuron-methyl metabolism in resistant plants was confirmed by LC-MS/MS analysis. Besides, GST activity was higher among resistant than susceptible individuals. The next transcriptome analysis generated 544,102,236 clean reads from RNA sequencing libraries. De novo assembly yielded 102,529 unigenes with an average length of 866 bp, annotated across seven databases. Digital gene expression selected 25 differentially expressed genes, further validated with qRT-PCR. Three P450 genes, two GST genes, two glucosyltransferase genes, four ABC transporter genes, and four additional contigs were constitutively up-regulated in resistant individuals. Overall, our research confirmed that enhanced herbicide metabolism drives tribenuron-methyl resistance in M. aquaticum.

Low doses of triazine xenobiotics mobilize ABA and cytokinin regulations in a stress- and low-energy-dependent manner

The extent of residual contaminations of pesticides through drift, run-off and leaching is a potential threat to non-target plant communities. Arabidopsis thaliana responds to low doses of the herbicide atrazine, and of its degradation products, desethylatrazine and hydroxyatrazine, not only in the long term, but also under conditions of short-term exposure. In order to investigate underlying molecular mechanisms of low-dose responses and to decipher commonalities and specificities between different chemical treatments, parallel transcriptomic studies of the early effects of the atrazine-desethylatrazine-hydroxyatrazine chemical series were undertaken using whole-genome microarrays. All of the triazines under study produced coordinated and specific changes in gene expression. Hydroxyatrazine-responsive genes were mainly linked to root development, whereas atrazine and desethylatrazine mostly affected molecular signaling networks implicated in stress and hormone responses. Analysis of signaling-related genes, promoter sites and shared-function interaction networks highlighted the involvement of energy-, stress-, abscisic acid- and cytokinin-regulated processes, and emphasized the importance of cold-, heat- and drought-related signaling in the perception of low doses of triazines. These links between low-dose xenobiotic impacts and stress-hormone crosstalk pathways give novel insights into plant-pesticide interactions and plant-pollution interactions that are essential for toxicity evaluation in the context of environmental risk assessment.

Stress-induced evolution of herbicide resistance and related pleiotropic effects

Herbicide-resistant weeds, especially those with resistance to multiple herbicides, represent a growing worldwide threat to agriculture and food security. Natural selection for resistant genotypes may act on standing genetic variation, or on a genetic and physiological background that is fundamentally altered because of stress responses to sublethal herbicide exposure. Stress-induced changes include DNA mutations, epigenetic alterations, transcriptional remodeling, and protein modifications, all of which can lead to herbicide resistance and a wide range of pleiotropic effects. Resistance selected in this manner is termed systemic acquired herbicide resistance, and the associated pleiotropic effects are manifested as a suite of constitutive transcriptional and post-translational changes related to biotic and abiotic stress adaptation, representing the evolutionary signature of selection. This phenotype is being investigated in two multiple herbicide-resistant populations of the hexaploid, self-pollinating weedy monocot Avena fatua that display such changes as well as constitutive reductions in certain heat shock proteins and their transcripts, which are well known as global regulators of diverse stress adaptation pathways. Herbicide-resistant populations of most weedy plant species exhibit pleiotropic effects, and their association with resistance genes presents a fertile area of investigation. This review proposes that more detailed studies of resistant A. fatua and other species through the lens of plant evolution under stress will inform improved resistant weed prevention and management strategies. © 2018 Society of Chemical Industry.

Herbicide-related signaling in plants reveals novel insights for herbicide use strategies, environmental risk assessment and global change assessment challenges

Herbicide impact is usually assessed as the result of a unilinear mode of action on a specific biochemical target with a typical dose-response dynamics. Recent developments in plant molecular signaling and crosstalk between nutritional, hormonal and environmental stress cues are however revealing a more complex picture of inclusive toxicity. Herbicides induce large-scale metabolic and gene-expression effects that go far beyond the expected consequences of unilinear herbicide-target-damage mechanisms. Moreover, groundbreaking studies have revealed that herbicide action and responses strongly interact with hormone signaling pathways, with numerous regulatory protein-kinases and -phosphatases, with metabolic and circadian clock regulators and with oxidative stress signaling pathways. These interactions are likely to result in mechanisms of adjustment that can determine the level of sensitivity or tolerance to a given herbicide or to a mixture of herbicides depending on the environmental and developmental status of the plant. Such regulations can be described as rheostatic and their importance is discussed in relation with herbicide use strategies, environmental risk assessment and global change assessment challenges.

Eigenvector metabolite analysis reveals dietary effects on the association among metabolite correlation patterns, gene expression, and phenotypes

INTRODUCTION

'Multi-omics' datasets obtained from an organism of interest reared under different environmental treatments are increasingly common. Identifying the links among metabolites and transcripts can help to elucidate our understanding of the impact of environment at different levels within the organism. However, many methods for characterizing physiological connections cannot address unidentified metabolites.

OBJECTIVES

Here, we use Eigenvector Metabolite Analysis (EvMA) to examine links between metabolomic, transcriptomic, and phenotypic variation data and to assess the impact of environmental factors on these associations. Unlike other methods, EvMA can be used to analyze datasets that include unidentified metabolites and unannotated transcripts.

METHODS

To demonstrate the utility of EvMA, we analyzed metabolomic, transcriptomic, and phenotypic datasets produced from 20 Drosophila melanogaster genotypes reared on four dietary treatments. We used a hierarchical distance-based method to cluster the metabolites. The links between metabolite clusters, gene expression, and overt phenotypes were characterized using the eigenmetabolite (first principal component) of each cluster.

RESULTS

EvMA recovered chemically related groups of metabolites within the clusters. Using the eigenmetabolite, we identified genes and phenotypes that significantly correlated with each cluster. EvMA identifies new connections between the phenotypes, metabolites, and gene transcripts. Conclusion EvMA provides a simple method to identify correlations between metabolites, gene expression, and phenotypes, which can allow us to partition multivariate datasets into meaningful biological modules and identify under-studied metabolites and unannotated gene transcripts that may be central to important biological processes. This can be used to inform our understanding of the effect of environmental mechanisms underlying physiological states of interest.