RNA-Seq Analysis of Rice Roots Reveals the Involvement of Post-Transcriptional Regulation in Response to Cadmium Stress

Athyrium yokoscense, a cadmium-hypertolerant fern, exhibits two cadmium stress mitigation strategies in its roots and aerial parts

Athyrium yokoscense is hypertolerant to cadmium (Cd) and can grow normally under a high Cd concentration despite Cd being a highly toxic heavy metal. To mitigate Cd stress in general plant species, Cd is promptly chelated with a thiol compound and is isolated into vacuoles. Generated active oxygen species (ROS) in the cytoplasm are removed by reduced glutathione. However, we found many differences in the countermeasures in A. yokoscense. Thiol compounds accumulated in the stele of the roots, although a long-term Cd exposure induced Cd accumulation in the aerial parts. Synchrotron radiation-based X-ray fluorescence (SR-XRF) analysis indicated that a large amount of Cd was localized in the cell walls of the roots. Overexpression of AyNramp5a, encoding a representative Fe and Mn transporter of A. yokoscense, increased both Cd uptake and Fe and Mn uptake in rice calli under the Cd exposure conditions. Organic acids are known to play a key role in reducing Cd availability to the plants by forming chelation and preventing its entry in free form into the roots. In A. yokoscense roots, Organic acids were abundantly detected. Investigating the chemical forms of the Cd molecules by X-ray absorption fine structure (XAFS) analysis detected many compounds with Cd-oxygen (Cd-O) binding in A. yokoscense roots, whereas in the aerial parts, the ratio of the compounds with Cd-sulfur (Cd-S) binding was increased. Together, our results imply that the strong Cd tolerance of A. yokoscense is an attribute of the following two mechanisms: Cd-O compound formation in the cell wall is a barrier to reduce Cd uptake into aerial parts. Thiol compounds in the region of root stele are involved in detoxication of Cd by formation of Cd-S compounds.

Integrated Review of Transcriptomic and Proteomic Studies to Understand Molecular Mechanisms of Rice's Response to Environmental Stresses

Rice (Oryza sativa L.) is grown nearly worldwide and is a staple food for more than half of the world's population. With the rise in extreme weather and climate events, there is an urgent need to decode the complex mechanisms of rice's response to environmental stress and to breed high-yield, high-quality and stress-resistant varieties. Over the past few decades, significant advancements in molecular biology have led to the widespread use of several omics methodologies to study all aspects of plant growth, development and environmental adaptation. Transcriptomics and proteomics have become the most popular techniques used to investigate plants' stress-responsive mechanisms despite the complexity of the underlying molecular landscapes. This review offers a comprehensive and current summary of how transcriptomics and proteomics together reveal the molecular details of rice's response to environmental stresses. It also provides a catalog of the current applications of omics in comprehending this imperative crop in relation to stress tolerance improvement and breeding. The evaluation of recent advances in CRISPR/Cas-based genome editing and the application of synthetic biology technologies highlights the possibility of expediting the development of rice cultivars that are resistant to stress and suited to various agroecological environments.

Integrative study of transcriptome and microbiome to reveal the response of Rhododendron decorum to cadmium stress

The anomalies of cadmium (Cd) in karst region pose a severe threat to plant growth and development. In this study, the responses of Rhododendron decorum to Cd stress were investigated at physiological, molecular, and endophytic microbial levels, and the potential correlation among these responses was assessed. The Cd stress impeded R. decorum growth and led to an increase in malondialdehyde (MDA) and hydrogen peroxide (H2O2) levels, as well as enhanced superoxide dismutase (SOD) and catalase (CAT) activities. Meanwhile, Cd stress increased the Cd (up to 80 times compared to the control), sodium (Na), aluminum (Al), and zinc (Zn) contents, while decreased the magnesium (Mg) and manganese (Mn) contents in R. decorum leaves. Transcriptome suggested that Cd significantly regulated the pathways including "protein repair", "hormone-mediated signaling pathway", and "ATP-binding cassette (ABC) transporters". Additionally, q-PCR analysis showed that Cd stress significantly up-regulated the expressions of ABCB19-like and pleiotropic drug resistance, while down-regulated the expressions of indole-3-acetic acid-amido synthetase and cytokinin dehydrogenase. The Cd stress influenced the composition of endophytic microbial communities in R. decorum leaves and enhanced the interspecific bacterial associations. Furthermore, the bacterial genera Achromobacter, Aureimonas and fungal genus Vishniacozyma exhibited a high degree of connectivity with other nodes in networks constructed by the metal element contents, differentially expressed genes (DEGs), and microbial communities, respectively. These findings provide a comprehensive insight into the response of R. decorum to Cd-induced stress, which might facilitate the breeding of the Cd-tolerant R. decorum.

Underground communication: Long non-coding RNA signaling in the plant rhizosphere

Long non-coding RNAs (lncRNAs) have emerged as integral gene-expression regulators underlying plant growth, development, and adaptation. To adapt to the heterogeneous and dynamic rhizosphere, plants use interconnected regulatory mechanisms to optimally fine-tune gene-expression-governing interactions with soil biota, as well as nutrient acquisition and heavy metal tolerance. Recently, high-throughput sequencing has enabled the identification of plant lncRNAs responsive to rhizosphere biotic and abiotic cues. Here, we examine lncRNA biogenesis, classification, and mode of action, highlighting the functions of lncRNAs in mediating plant adaptation to diverse rhizosphere factors. We then discuss studies that reveal the significance and target genes of lncRNAs during developmental plasticity and stress responses at the rhizobium interface. A comprehensive understanding of specific lncRNAs, their regulatory targets, and the intricacies of their functional interaction networks will provide crucial insights into how these transcriptomic switches fine-tune responses to shifting rhizosphere signals. Looking ahead, we foresee that single-cell dissection of cell-type-specific lncRNA regulatory dynamics will enhance our understanding of the precise developmental modulation mechanisms that enable plant rhizosphere adaptation. Overcoming future challenges through multi-omics and genetic approaches will more fully reveal the integral roles of lncRNAs in governing plant adaptation to the belowground environment.

Exploration the interaction of cadmium and copper toxic effects in pakchoi (Brassica chinensis L) roots through combinatorial transcriptomic and weighted gene co-expression network analysis

The interaction between cadmium(Cd) and copper(Cu) during combined pollution can lead to more complex toxic effects on humans and plants.However, there is still a lack of sufficient understanding regarding the types of interactions at the plant molecular level and the response strategies of plants to combined pollution. To assess this, we investigated the phenotypic and transcriptomic patterns of pakchoi (Brassica chinensis L) roots in response to individual and combined pollution of Cd and Cu. The results showed that compared to single addition, the translocation factor of heavy metals in roots significantly decreased (p < 0.05) under the combined addition, resulting in higher accumulation of Cd and Cu in the roots. Transcriptomic analysis of pakchoi roots revealed that compared to single pollution, there were 312 and 1926 differentially expressed genes (DEGs) specifically regulated in the Cd2Cu20 and Cd2Cu100 combined treatments, respectively. By comparing the expression of these DEGs among different treatments, we found that the combined pollution of Cd and Cu mainly affected the transcriptome of the roots in an antagonistic manner. Enrichment analysis indicated that pakchoi roots upregulated the expression of genes involved in glucosetransferase activity, phospholipid homeostasis, proton transport, and the biosynthesis of phenylpropanoids and flavonoids to resist Cd and Cu combined pollution. Using weighted gene co-expression network analysis (WGCNA), we identified hub genes related to the accumulation of Cd and Cu in the roots, which mainly belonged to the LBD, thaumatin-like protein, ERF, MYB, WRKY, and TCP transcription factor families. This may reflect a transcription factor-driven trade-off strategy between heavy metal accumulation and growth in pakchoi roots. Additionally, compared to single metal pollution, the expression of genes related to Nramp, cation/H+ antiporters, and some belonging to the ABC transporter family in the pakchoi roots was significantly upregulated under combined pollution. This could lead to increased accumulation of Cd and Cu in the roots. These findings provide new insights into the interactions and toxic mechanisms of multiple metal combined pollution at the molecular level in plants.

The C2H2 family protein ZAT17 engages in the cadmium stress response by interacting with PRL1 in Arabidopsis

Cadmium (Cd) is a heavy metal and a toxic substance. Soil Cd pollution has emerged as a significant environmental issue that jeopardizes both the safety of agricultural products and human health. PLEIOTROPIC REGULATORY LOCUS 1 (PRL1) has been identified as a crucial factor in Cd stress and a series of defence mechanisms. However, the mechanism through which PRL1 mediates its downstream signalling has remained poorly understood. Here, we discovered a prl1-2 suppressor (sup8) for prl1-2 that complemented the defective development phenotype of prl1-2 under Cd stress. Gene cloning revealed a mutation in the C2H2 transcription factor ZAT17 as the basis for the sup8 phenotype. Genetic and biochemical studies indicated that ZAT17 acts as a negative regulator of Cd tolerance. Transcriptome analysis revealed that ZAT17 influences the alternative splicing (AS) process of multiple Cd-responsive genes by interacting with members of the MAC splicing complex, including PRL1 and CDC5. In conclusion, the identification of the novel gene ZAT17 enriches the understanding of the Cd stress response pathway and provides a valuable candidate locus for breeding Cd-resistant plant varieties.

Stable Isotope Analyses Reveal Impact of Fe and Zn on Cd Uptake and Translocation by Theobroma cacao

High concentrations of toxic cadmium (Cd) in soils are problematic as the element accumulates in food crops such as rice and cacao. A mitigation strategy to minimise Cd accumulation is to enhance the competitive uptake of plant-essential metals. Theobroma cacao seedlings were grown hydroponically with added Cd. Eight different treatments were used, which included/excluded hydroponic or foliar zinc (Zn) and/or iron (Fe) for the final growth period. Analyses of Cd concentrations and natural stable isotope compositions by multiple collector ICP-MS were conducted. Cadmium uptake and translocation decreased when Fe was removed from the hydroponic solutions, while the application of foliar Zn-EDTA may enhance Cd translocation. No significant differences in isotope fractionation during uptake were found between treatments. Data from all treatments fit a single Cd isotope fractionation model associated with sequestration (seq) of isotopically light Cd in roots and unidirectional mobilisation (mob) of isotopically heavier Cd to the leaves (ε114Cdseq-mob = -0.13‱). This result is in excellent agreement with data from an investigation of 19 genetically diverse cacao clones. The different Cd dynamics exhibited by the clones and seen in response to different Fe availability may be linked to similar physiological processes, such as the regulation of specific transporter proteins.

Zinc accumulation in Atriplex lentiformis is driven by plant genes and the soil microbiome

Successful phytoremediation of acidic metal-contaminated mine tailings requires amendments to condition tailings properties prior to plant establishment. This conditioning process is complex and includes multiple changes in tailings bio-physico-chemical properties. The objective of this project is to identify relationships between tailings properties, the soil microbiome, and plant stress response genes during growth of Atriplex lentiformis in compost-amended (10 %, 15 %, 20 % w/w) mine tailings. Analyses include RNA-Seq for plant root gene expression, 16S rRNA amplicon sequencing for bacterial/archaeal communities, metal concentrations in both tailings and plant organs, and phenotypic measures of plant stress. Zn accumulation in A. lentiformis leaves varied with compost levels and was the highest in the intermediate treatment (15 %, TC15). Microbial analysis identified Alicyclobacillus, Hydrotalea, and Pseudolabrys taxa with the highest relative abundance in TC15, and these taxa were strongly associated with Zn accumulation. Furthermore, we identified 190 root genes with significant gene expression changes. These root genes were associated with different pathways including, abscisic acid and auxin signaling, defense responses, ion channels, metal ion binding, oxidative stress, transcription regulation, and transmembrane transport. However, root gene expression changes were not driven by the increasing levels of compost. For example, there were 15 genes that were up-regulated in TC15, whereas 106 genes were down-regulated in TC15. The variables analyzed explained 86 % of the variance in Zn accumulation in A. lentiformis leaves. Importantly, Zn accumulation was driven by Zn shoot concentrations, leaf stress symptoms, plant root genes, and microbial taxa. Therefore, our results suggest there are strong plant-microbiome associations that drive Zn accumulation in A. lentiformis and different plant gene pathways are involved in alleviating varying levels of metal stress. Future work is needed to gain a mechanistic understanding of these plant-microbiome interactions to optimize phytoremediation strategies as they will govern the success or failure of the revegetation process.

Transcriptome-based strategies for identifying aluminum tolerance genes in popcorn (Zea mays L. var. everta)

Aluminum (Al) toxicity limits crop production worldwide. Although studies have identified genes associated with Al tolerance in crops, a large amount of data remains unexplored using other strategies. Here, we searched for single substitutions and InDels across differentially expressed genes (DEGs), linked DEGs to Al-tolerance QTLs reported in the literature for common maize, and investigated the alternative splicing regulated by Al3+ toxicity. We found 929 substitutions between DEGs in Al-tolerant and 464 in Al-sensitive inbred lines, of which 165 and 80 were non-synonymous, respectively. Only 12 NS variants had deleterious predicted effect on protein function in Al-tolerant and 13 in Al-sensitive. Moreover, 378 DEGs were mapped in Al-QTL regions for the Al-tolerant and 213 for the Al-sensitive. Furthermore, Al stress is primarily regulated at the transcriptional level in popcorn. Important genes identified, such as HDT1, SWEET4a, GSTs, SAD9, PIP2-2, CASP-like 5, and AGP, may benefit molecular assisted popcorn breeding or be useful in biotechnological approaches. These findings offer insights into the mechanisms of Al tolerance in popcorn and provide a 'hypothesis-free' strategy for identifying and prioritizing candidate genes that could be used to develop molecular markers or cultivars resilient to acidic soils.

Molecular-Assisted Breeding of Cadmium Pollution-Safe Cultivars

Cadmium (Cd) contamination in edible agricultural products, especially in crops intended for consumption, has raised worldwide concerns regarding food safety. Breeding of Cd pollution-safe cultivars (Cd-PSCs) is an effective solution to preventing the entry of Cd into the food chain from contaminated agricultural soil. Molecular-assisted breeding methods, based on molecular mechanisms for cultivar-dependent Cd accumulation and bioinformatic tools, have been developed to accelerate and facilitate the breeding of Cd-PSCs. This review summarizes the recent progress in the research of the low Cd accumulation traits of Cd-PSCs in different crops. Furthermore, the application of molecular-assisted breeding methods, including transgenic approaches, genome editing, marker-assisted selection, whole genome-wide association analysis, and transcriptome, has been highlighted to outline the breeding of Cd-PSCs by identifying critical genes and molecular biomarkers. This review provides a comprehensive overview of the development of Cd-PSCs and the potential future for breeding Cd-PSC using modern molecular technologies.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

Evaluation of nanoplastics toxicity in the soil nematode Caenorhabditis elegans by iTRAQ-based quantitative proteomics

Plastic pollution is recognized as a major threat to ecosystems in the 21st century. Large plastic objects undergo biotic and abiotic degradation to generate micro- and nano-sized plastic pieces. Despite tremendous efforts to evaluate the adverse effects of microplastics, a comprehensive understanding of the toxicity of nanoplastics remains elusive, especially at the protein level. To this end, we used isobaric-tag-for-relative-and-absolute-quantitation-based quantitative proteomics to investigate the proteome dynamics of the soil nematode Caenorhabditis elegans in response to exposure to 100 nm polystyrene nanoplastics (PS-NPs). After 48 h of exposure to 0.1, 1, or 10 mg/L PS-NPs, 136 out of 1684 proteins were differentially expressed and 108 of these proteins were upregulated. These proteins were related to ribosome biogenesis, translation, proteolysis, kinases, protein processing in the endoplasmic reticulum, and energy metabolism. Remarkably, changes in proteome dynamics in response to exposure to PS-NPs were consistent with the phenotypic defects of C. elegans. Collectively, our findings demonstrate that disruption of proteome homeostasis is a biological consequence of PS-NPs accumulation in C. elegans, which provides insights into the molecular mechanisms underlying the toxicology of nanoplastics.

IlAP2, an AP2/ERF Superfamily Gene, Mediates Cadmium Tolerance by Interacting with IlMT2a in Iris lactea var. chinensis

Cadmium (Cd) stress has a major impact on ecosystems, so it is important to find suitable Cd-tolerant plants while elucidating the responsible molecular mechanism for phytoremediation to manage Cd soil contamination. Iris lactea var. chinensis is an ornamental perennial groundcover plant with strong tolerance to Cd. Previous studies found that IlAP2, an AP2/ERF superfamily gene, may be an interacting partner of the metallothionein gene IlMT2a, which plays a key role in Cd tolerance. To study the role of IlAP2 in regulating Cd tolerance in I. lactea, we analyzed its regulation function and mechanism based on a yeast two-hybrid assay, a bimolecular fluorescence complementation test, quantitative real-time PCR, transgenics and transcriptome sequencing. The results showed that IlAP2 interacts with IlMT2a and may cooperate with other transcription factors to regulate genes involved in signal transduction and plant hormones, leading to reduced Cd toxicity by hindering Cd transport. These findings provide insights into the mechanism of IlAP2-mediated stress responses to Cd and important gene resources for improving plant stress tolerance in phytoremediation.

Genetic Regulation Mechanism of Cadmium Accumulation and Its Utilization in Rice Breeding

Cadmium (Cd) is a heavy metal whose pollution in rice fields leads to varying degrees of Cd accumulation in rice. Furthermore, the long-term consumption of Cd-contaminated rice is harmful to human health. Therefore, it is of great theoretical significance and application value to clarify the genetic regulation mechanism of Cd accumulation in rice and cultivate rice varieties with low Cd accumulation for the safe use of Cd-contaminated soils. This review summarizes the effects of Cd on rice growth, yield, and quality; the physiological and molecular mechanisms of Cd absorption in the roots, loading, and transport of Cd in the xylem, the distribution of Cd in nodes, redistribution of Cd in leaves, and accumulation of Cd in the grains; the regulation mechanism of the Cd stress response; and the breeding of rice with low Cd accumulation. Future directions on the genetic regulation of Cd in rice and application are also discussed. This review provides a theoretical basis for studies exploring the genetic regulation of Cd stress in rice. It also offers a basis for formulating effective strategies to reduce the Cd content in rice.

Post-transcriptional regulation of 2-acetyl-1-pyrroline (2-AP) biosynthesis pathway, silicon, and heavy metal transporters in response to Zn in fragrant rice

Fragrant rice (Oryza sativa L.) has a high economic and nutritional value, and the application of micronutrients regulates 2-acetyl-1-pyrroline (2-AP) production, which is responsible for aroma in fragrant rice. Alternative splicing (AS) is an important post-transcriptional regulatory mechanism to generate transcript variability and proteome diversity in plants. However, no systematic investigation of AS events in response to micronutrients (Zn) has been performed in fragrant rice. Furthermore, the post-transcriptional regulation of genes involved in 2-AP biosynthesis is also not known. In this study, a comprehensive analysis of AS events under two gradients of Zn treatment in two different fragrant rice cultivars (Meixiangzhan-2 and Xiangyaxiangzhan) was performed based on RNA-seq analysis. A total of 386 and 598 significant AS events were found in Meixiangzhan-2 treated with low and high doses of Zn, respectively. In Xiangyaxiangzhan, a total of 449 and 598 significant AS events were found in low and high doses of Zn, respectively. Go analysis indicated that these genes were highly enriched in physiological processes, metabolism, and cellular processes in both cultivars. However, genotype and dose-dependent AS events were also detected in both cultivars. By comparing differential AS (DAS) events with differentially expressed genes (DEGs), we found a weak overlap among DAS and DEGs in both fragrant rice cultivars indicating that only a few genes are post-transcriptionally regulated in response to Zn treatment. We further report that Zn differentially regulates the expression of 2-AP biosynthesis-related genes in both cultivars and Zn treatment altered the editing frequency of single nucleotide polymorphism (SNPs) in the genes involved in 2-AP biosynthesis. Finally, we showed that epigenetic modifications associated with active gene transcription are generally enriched over 2-AP biosynthesis-related genes. Similar to the 2-AP pathway, we found that heavy metal transporters (genes related to silicon, iron, Zn and other metal transport) are also regulated at transcriptional and post-transcriptional levels in response to Zn in fragrant rice. Taken together, our results provide evidence of the post-transcriptional gene regulation in fragrant rice in response to Zn treatment and highlight that the 2-AP biosynthesis pathway and heavy metal transporters may also be regulated through epigenetic modifications. These findings will serve as a cornerstone for further investigation to understand the molecular mechanisms of 2-AP biosynthesis and regulation of heavy metal transporters in fragrant rice.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Single-molecule long-read sequencing reveals the potential impact of posttranscriptional regulation on gene dosage effects on the avian Z chromosome

Background

Mammalian sex chromosomes provide dosage compensation, but avian lack a global mechanism of dose compensation. Herein, we employed nanopore sequencing to investigate the genetic basis of gene expression and gene dosage effects in avian Z chromosomes at the posttranscriptional level.

Results

In this study, the gonad and head skin of female and male duck samples (n = 4) were collected at 16 weeks of age for Oxford nanopore sequencing. Our results revealed a dosage effect and local regulation of duck Z chromosome gene expression. Additionally, AS and APA achieve tissue-specific gene expression, and male-biased lncRNA regulates its Z-linked target genes, with a positive regulatory role for gene dosage effects on the duck Z chromosome. In addition, GO enrichment and KEGG pathway analysis showed that the dosage effects of Z-linked genes were mainly associated with the cellular response to hormone stimulus, melanin biosynthetic, metabolic pathways, and melanogenesis, resulting in sex differences.

Conclusions

Our data suggested that post transcriptional regulation (AS, APA and lncRNA) has a potential impact on the gene expression effects of avian Z chromosomes. Our study provides a new view of gene regulation underlying the dose effects in avian Z chromosomes at the RNA post transcriptional level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08360-8.

A Comprehensive Review on the Heavy Metal Toxicity and Sequestration in Plants

Heavy metal (HM) toxicity has become a global concern in recent years and is imposing a severe threat to the environment and human health. In the case of plants, a higher concentration of HMs, above a threshold, adversely affects cellular metabolism because of the generation of reactive oxygen species (ROS) which target the key biological molecules. Moreover, some of the HMs such as mercury and arsenic, among others, can directly alter the protein/enzyme activities by targeting their -SH group to further impede the cellular metabolism. Particularly, inhibition of photosynthesis has been reported under HM toxicity because HMs trigger the degradation of chlorophyll molecules by enhancing the chlorophyllase activity and by replacing the central Mg ion in the porphyrin ring which affects overall plant growth and yield. Consequently, plants utilize various strategies to mitigate the negative impact of HM toxicity by limiting the uptake of these HMs and their sequestration into the vacuoles with the help of various molecules including proteins such as phytochelatins, metallothionein, compatible solutes, and secondary metabolites. In this comprehensive review, we provided insights towards a wider aspect of HM toxicity, ranging from their negative impact on plant growth to the mechanisms employed by the plants to alleviate the HM toxicity and presented the molecular mechanism of HMs toxicity and sequestration in plants.

Comparative physiological and transcriptomic analyses illuminate common mechanisms by which silicon alleviates cadmium and arsenic toxicity in rice seedlings

The inessential heavy metal/loids cadmium (Cd) and arsenic (As), which often co-occur in polluted paddy soils, are toxic to rice. Silicon (Si) treatment is known to reduce Cd and As toxicity in rice plants. To better understand the shared mechanisms by which Si alleviates Cd and As stress, rice seedlings were hydroponically exposed to Cd or As, then treated with Si. The addition of Si significantly ameliorated the inhibitory effects of Cd and As on rice seedling growth. Si supplementation decreased Cd and As translocation from roots to shoots, and significantly reduced Cd- and As-induced reactive oxygen species generation in rice seedlings. Transcriptomics analyses were conducted to elucidate molecular mechanisms underlying the Si-mediated response to Cd or As stress in rice. The expression patterns of the differentially expressed genes in Cd- or As-stressed rice roots with and without Si application were compared. The transcriptomes of the Cd- and As-stressed rice roots were similarly and profoundly reshaped by Si application, suggesting that Si may play a fundamental, active role in plant defense against heavy metal/loid stresses by modulating whole genome expression. We also identified two novel genes, Os01g0524500 and Os06g0514800, encoding a myeloblastosis (MYB) transcription factor and a thionin, respectively, which may be candidate targets for Si to alleviate Cd and As stress in rice, as well as for the generation of Cd- and/or As-resistant plants. This study provides valuable resources for further clarification of the shared molecular mechanisms underlying the Si-mediated alleviation of Cd and As toxicity in rice.

A transcriptomic view of cadmium retention in roots of cadmium-safe rice line (Oryza sativa L.)

A better understanding of the mechanisms controlling cadmium (Cd) accumulation in rice will benefit the development of strategies to minimize Cd accumulation in grains. A Cd-safe rice line designated D62B accumulated less than 0.2 mg Cd kg^-1 in brown rice due to its strong capacity for Cd retention in roots. Here transcriptomic was used to clarify the underlying mechanisms of Cd response in roots of D62B compared with a high Cd-accumulating line (Wujin4B). There were 777, 1058 differentially expressed genes (DEGs) in D62B and Wujin4B, respectively, when exposed to Cd. The functions of DEGs were clearly line-specific. Cell wall biosynthesis responded more intensively to Cd stress in D62B, facilitating Cd restriction. Meanwhile, more glutathione (GSH) and phytochelatins synthesized in D62B with the upregulation of sulphur and GSH metabolism. Besides, membrane proteins played critical roles in Cd response in D62B, whereas 18 terms involved in regulation were enriched in Wujin4B. Exogenous GSH further induced the expression of genes related to GSH metabolism and cell wall biosynthesis, leading to the retention of more Cd. Great responsiveness of cell wall biosynthesis and GSH metabolism could be considered the most important specific mechanisms for Cd retention in the roots of Cd-safe rice line.

Physiological, biochemical, and transcriptional regulation in a leguminous forage Trifolium pratense L. responding to silver ions

Trifolium pratense L. (red clover) is an important leguminous crop with great potential for Ag-contaminated environment remediation. Whereas, the molecular mechanisms of Ag tolerance in red clover are largely unknown. Red clover seedlings were used for physiological and transcriptomic investigation under 0, 20, 50, and 100 mg/L Ag⁺ stress in our research to reveal potential molecular resistance mechanism. Research showed that red clover possessed fairly strong Ag absorbance capacity, the Ag level reached 0.14 and 2.35 mg/g·FW in the leaves and roots under 100 mg/L AgNO₃ stress condition. Root fresh weight, root dry weight, root water content, and photosynthetic pigments contents were significantly decreased with elevating AgNO₃ concentration. Obvious withered plant tissue, microstructure disorder, and disrupted organelles were observed. In vitro evaluations (e.g., PI and DCFH-DA staining) represented that AgNO₃ at high concentration (100 mg/L) exhibited obvious inhibition on cell viability, which was due possibly to the induction of reactive oxygen species (ROS) accumulation. A total of 44643 differentially expressed genes (DEGs) were identified under Ag stress, covering 27155 upregulated and 17488 downregulated genes. 12 stress-responsive DEGs was authenticated utilizing real-time quantitative PCR (qRT-PCR). Gene ontology (GO) analysis revealed that the DEGs were mostly related to metal ion binding (molecular function), nucleus (cellular component), and defense response (biological process). Involved DEGs in sequence-specific DNA binding transcription factor activity, response to various hormones (e.g., abscisic acid, IAA/Auxin, salicylic acid, and etc), calcium signal transduction, and protein ubiquitination were concluded to play crucial roles in Ag tolerance of red clover. On the other hand, Kyoto Encyclopedia of Genes and Genomes (KEGG) database annotated several stress responsive pathways such as plant-pathogen interaction, phenylpropanoid biosynthesis, ubiquitin mediated proteolysis, hormone signal transduction, and autophagy. Several down-regulated genes (e.g., RSF2, RCD1, DOX1, and etc) were identified indicating possible metabolic disturbance. Besides, protein-protein interaction network (PPI) identified several pivotal genes such as ribosomal proteins, TIR, and ZAT.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.

A meta-analysis of transcriptomic profiles reveals molecular pathways response to cadmium stress of Gramineae

As a non-essential heavy metal, cadmium (Cd) is toxic to plants. In the last 15 years, over 70 transcriptome studies have been published to decipher the molecular response mechanism against Cd stress in different plants. To extract generalization results from transcriptomic data across different plants and obtain some hub genes that respond to Cd stress, we carried out a meta-analysis of 32 published datasets. Cluster analysis revealed that plant species played a more decisive role than the media used and exposure time in the transcriptome patterns of plant roots response to Cd. The datasets from a Gramineae-like (GL) group were closer in clustering. 838 DEGs were commonly Cd-regulated in at least nine of 18 GL datasets. Gene ontology and KEGG pathway analyses revealed that oxidative stress-related terms and lignin synthesis-related terms were significantly enriched. Mapman analysis revealed that these common DEGs were mainly involved in regulation, cellular response, secondary metabolism, transport, cell wall and lipid metabolism. In Oryza sativa, 15 DEGs were up-regulated in at least four of five HM (As, Cr, Cd, Hg and Pb) groups, such as Os10g0517500 (methionine gamma-lyase) and Os01g0159800 (bHLH107). Moreover, our datasets can be used to retrieve log2FC value of specific genes across 29 studies (48 datasets), which provides data reference for the subsequent selection of HM-related genes. Our results provide the basis for further understanding of Cd tolerance mechanisms in plants.

From Trash to Luxury: The Potential Role of Plant LncRNA in DNA Methylation During Abiotic Stress

Remarkable progress has been made in elucidating important roles of plant non-coding RNAs. Among these RNAs, long noncoding RNAs (lncRNAs) have gained widespread attention, especially their role in plant environmental stress responses. LncRNAs act at different levels of gene expression regulation, and one of these mechanisms is by recruitment of DNA methyltransferases or demethylases to regulate the target gene transcription. In this mini-review, we highlight the function of lncRNAs, including their potential role in RNA-directed DNA Methylation (RdDM) silencing pathway and their potential function under abiotic stresses conditions. Moreover, we also present and discuss studies of lncRNAs in crops. Finally, we propose a path outlook for future research that may be important for plant breeding.

Comparative transcriptomic analysis reveals key genes and pathways in two different cadmium tolerance kenaf (Hibiscus cannabinus L.) cultivars

Soil cadmium (Cd) contamination has become a massive environmental problem. Kenaf is an industrial fiber crop with high tolerance to heavy metals and could be potentially used for soil phytoremediation. However, the molecular mechanism of Cd in kenaf tolerance remains largely unknown. In the present study, using two contrasting Cd sensitive kenaf (GH and YJ), the key factors accounting for differential Cd tolerance were investigated. GH has a stronger Cd transport and accumulation ability than YJ. In addition, physiological index investigation on malondialdehyde (MDA) contents and antioxidant enzyme (SOD, POD, and CAT) activities showed GH has a stronger detoxification capacity than YJ. Furthermore, the cell ultrastructure of GH is more stable than that of YJ under Cd stress. Transcriptome analysis revealed 2221 (689 up and 1532 down) and 3321 (2451 up and 870 down) genes were differentially expressed in GH and YJ, respectively. More DEGs (differentially expressed genes) were characterized as up-regulated in GH, indicating GH is inclined to activate gene expression to cope with cadmium stress. GO and KEGG analyses indicate that DEGs were assigned and enriched in different pathways. Plenty of critical Cd-induced DEGs such as SOD2, PODs, MT1, DTXs, NRT1, ABCs, CES, AP2/ERF, MYBs, NACs, and WRKYs were identified. The DEGs involved pathways, including antioxidant, heavy metal transport or detoxification, substance transport, plant hormone and calcium signals, ultrastructural component, and a wide range of transcription factors were suggested to play crucial roles in kenaf Cd tolerance, and accounting for the difference in Cd stress sensitivities.

Chemical Speciation, Plant Uptake, and Toxicity of Heavy Metals in Agricultural Soils

Heavy metals in agricultural soils exist in diverse dissolved (free cations and complexed species of positive, neutral, or negative charges), particulate (sorbed, structural, and coprecipitated), and colloidal (micro- and nanometer-sized particles) species. The fate of different heavy metal species is controlled by the master variables: pH (solubility), ionic strength (activity and charge-shielding), and dissolved organic carbon (complexation). In the rhizosphere, chemical speciation controls toxicokinetics (uptake and transport of metals by plants) while toxicodynamics (interaction between the plant and absorbed species) drives the toxicity outcome. Based on the critical review, the authors recommend omics and data mining techniques to link discrete knowledge bases from the speciation dynamics, soil microbiome, and plant transporter/gene expression relevant to homeostasis conditions of modern agriculture. Such efforts could offer a disruptive application tool to improve and sustain plant tolerance, food safety, and environmental quality.

Rice Biofortification With Zinc and Selenium: A Transcriptomic Approach to Understand Mineral Accumulation in Flag Leaves

Human malnutrition due to micronutrient deficiencies, particularly with regards to Zinc (Zn) and Selenium (Se), affects millions of people around the world, and the enrichment of staple foods through biofortification has been successfully used to fight hidden hunger. Rice (Oryza sativa L.) is one of the staple foods most consumed in countries with high levels of malnutrition. However, it is poor in micronutrients, which are often removed during grain processing. In this study, we have analyzed the transcriptome of rice flag leaves biofortified with Zn (900 g ha-1), Se (500 g ha-1), and Zn-Se. Flag leaves play an important role in plant photosynthesis and provide sources of metal remobilization for developing grains. A total of 3170 differentially expressed genes (DEGs) were identified. The expression patterns and gene ontology of DEGs varied among the three sets of biofortified plants and were limited to specific metabolic pathways related to micronutrient mobilization and to the specific functions of Zn (i.e., its enzymatic co-factor/coenzyme function in the biosynthesis of nitrogenous compounds, carboxylic acids, organic acids, and amino acids) and Se (vitamin biosynthesis and ion homeostasis). The success of this approach should be followed in future studies to understand how landraces and other cultivars respond to biofortification.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

In Silico Identification of QTL-Based Polymorphic Genes as Salt-Responsive Potential Candidates through Mapping with Two Reference Genomes in Rice

Recent advances in next generation sequencing have created opportunities to directly identify genetic loci and candidate genes for abiotic stress responses in plants. With the objective of identifying candidate genes within the previously identified QTL-hotspots, the whole genomes of two divergent cultivars for salt responses, namely At 354 and Bg 352, were re-sequenced using Illumina Hiseq 2500 100PE platform and mapped to Nipponbare and R498 genomes. The sequencing results revealed approximately 2.4 million SNPs and 0.2 million InDels with reference to Nipponbare while 1.3 million and 0.07 million with reference to R498 in two parents. In total, 32,914 genes were reported across all rice chromosomes of this study. Gene mining within QTL hotspots revealed 1236 genes, out of which 106 genes were related to abiotic stress. In addition, 27 abiotic stress-related genes were identified in non-QTL regions. Altogether, 32 genes were identified as potential genes containing polymorphic non-synonymous SNPs or InDels between two parents. Out of 10 genes detected with InDels, tolerant haplotypes of Os01g0581400, Os10g0107000, Os11g0655900, Os12g0622500, and Os12g0624200 were found in the known salinity tolerant donor varieties. Our findings on different haplotypes would be useful in developing resilient rice varieties for abiotic stress by haplotype-based breeding studies.

Full-Length Transcriptome Assembly of Italian Ryegrass Root Integrated with RNA-Seq to Identify Genes in Response to Plant Cadmium Stress

Cadmium (Cd) is a toxic heavy metal element. It is relatively easily absorbed by plants and enters the food chain, resulting in human exposure to Cd. Italian ryegrass (Lolium multiflorum Lam.), an important forage cultivated widely in temperate regions worldwide, has the potential to be used in phytoremediation. However, genes regulating Cd translocation and accumulation in this species are not fully understood. Here, we optimized PacBio ISO-seq and integrated it with RNA-seq to construct a de novo full-length transcriptomic database for an un-sequenced autotetraploid species. With the database, we identified 2367 differentially expressed genes (DEGs) and profiled the molecular regulatory pathways of Italian ryegrass with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis in response to Cd stress. Overexpression of a DEG LmAUX1 in Arabidopsis thaliana significantly enhanced plant Cd concentration. We also unveiled the complexity of alternative splicing (AS) with a genome-free strategy. We reconstructed full-length UniTransModels using the reference transcriptome, and 29.76% of full-length models had more than one isoform. Taken together, the results enhanced our understanding of the genetic diversity and complexity of Italian ryegrass under Cd stress and provided valuable genetic resources for its gene identification and molecular breeding.

Long Non-coding RNA in Plants in the Era of Reference Sequences

The discovery of non-coding RNAs (ncRNAs), and the subsequent elucidation of their functional roles, was largely delayed due to the misidentification of non-protein-coding parts of DNA as "junk DNA," which forced ncRNAs into the shadows of their protein-coding counterparts. However, over the past decade, insight into the important regulatory roles of ncRNAs has led to rapid progress in their identification and characterization. Of the different types of ncRNAs, long non-coding RNAs (lncRNAs), has attracted considerable attention due to their mRNA-like structures and gene regulatory functions in plant stress responses. While RNA sequencing has been commonly used for mining lncRNAs, a lack of widespread conservation at the sequence level in addition to relatively low and highly tissue-specific expression patterns challenges high-throughput in silico identification approaches. The complex folding characteristics of lncRNA molecules also complicate target predictions, as the knowledge about the interaction interfaces between lncRNAs and potential targets is insufficient. Progress in characterizing lncRNAs and their targets from different species may hold the key to efficient identification of this class of ncRNAs from transcriptomic and potentially genomic resources. In wheat and barley, two of the most important crops, the knowledge about lncRNAs is very limited. However, recently published high-quality genomes of these crops are considered as promising resources for the identification of not only lncRNAs, but any class of molecules. Considering the increasing demand for food, these resources should be used efficiently to discover molecular mechanisms lying behind development and a/biotic stress responses. As our understanding of lncRNAs expands, interactions among ncRNA classes, as well as interactions with the coding sequences, will likely define novel functional networks that may be modulated for crop improvement.

Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales)

The limitations of RNA sequencing make it difficult to accurately predict alternative splicing (AS) and alternative polyadenylation (APA) events and long non-coding RNAs (lncRNAs), all of which reveal transcriptomic diversity and the complexity of gene regulation. Gnetum, a genus with ambiguous phylogenetic placement in seed plants, has a distinct stomatal structure and photosynthetic characteristics. In this study, a full-length transcriptome of Gnetum luofuense leaves at different developmental stages was sequenced with the latest PacBio Sequel platform. After correction by short reads generated by Illumina RNA-Seq, 80,496 full-length transcripts were obtained, of which 5269 reads were identified as isoforms of novel genes. Additionally, 1660 lncRNAs and 12,998 AS events were detected. In total, 5647 genes in the G. luofuense leaves had APA featured by at least one poly(A) site. Moreover, 67 and 30 genes from the bHLH gene family, which play an important role in stomatal development and photosynthesis, were identified from the G. luofuense genome and leaf transcripts, respectively. This leaf transcriptome supplements the reference genome of G. luofuense, and the AS events and lncRNAs detected provide valuable resources for future studies of investigating low photosynthetic capacity of Gnetum.

Genome-wide alternative polyadenylation dynamics in response to biotic and abiotic stresses in rice

Alternative polyadenylation (APA) is an important way to regulate gene expression at the post-transcriptional level, and is extensively involved in plant stress responses. However, the systematic roles of APA regulation in response to abiotic and biotic stresses in rice at the genome scale remain unknown. To take advantage of available RNA-seq datasets, using a novel tool APAtrap, we identified thousands of genes with significantly differential usage of polyadenylation [poly(A)] sites in response to the abiotic stress (drought, heat shock, and cadmium) and biotic stress [bacterial blight (BB), rice blast, and rice stripe virus (RSV)]. Genes with stress-responsive APA dynamics commonly exhibited higher expression levels when their isoforms with short 3' untranslated region (3' UTR) were more abundant. The stress-responsive APA events were widely involved in crucial stress-responsive genes and pathways: e.g. APA acted as a negative regulator in heat stress tolerance; APA events were involved in DNA repair and cell wall formation under Cd stress; APA regulated chlorophyll metabolism, being associated with the pathogenesis of leaf diseases under RSV and BB challenges. Furthermore, APA events were found to be involved in glutathione metabolism and MAPK signaling pathways, mediating a crosstalk among the abiotic and biotic stress-responsive regulatory networks in rice. Analysis of large-scale datasets revealed that APA may regulate abiotic and biotic stress-responsive processes in rice. Such post-transcriptome diversities contribute to rice adaption to various environmental challenges. Our study would supply useful resource for further molecular assisted breeding of multiple stress-tolerant cultivars for rice.

Gene identification and transcriptome analysis of low cadmium accumulation rice mutant (lcd1) in response to cadmium stress using MutMap and RNA-seq

BACKGROUND

Cadmium (Cd) is a widespread toxic heavy metal pollutant in agricultural soil, and Cd accumulation in rice grains is a major intake source of Cd for Asian populations that adversely affect human health. However, the molecular mechanism underlying Cd uptake, translocation and accumulation has not been fully understood in rice plants.

RESULTS

In this study, a mutant displaying extremely low Cd accumulation (lcd1) in rice plant and grain was generated by EMS mutagenesis from indica rice cultivar 9311 seeds. The candidate SNPs associated with low Cd accumulation phenotype in the lcd1 mutant were identified by MutMap and the transcriptome changes between lcd1 and WT under Cd exposure were analyzed by RNA-seq. The lcd1 mutant had lower Cd uptake and accumulation in rice root and shoot, as well as less growth inhibition compared with WT in the presence of 5 μM Cd. Genetic analysis showed that lcd1 was a single locus recessive mutation. The SNP responsible for low Cd accumulation in the lcd1 mutant located at position 8,887,787 on chromosome 7, corresponding to the seventh exon of OsNRAMP5. This SNP led to a Pro236Leu amino acid substitution in the highly conserved region of OsNRAMP5 in the lcd1 mutant. A total of 1208 genes were differentially expressed between lcd1 and WT roots under Cd exposure, and DEGs were enriched in transmembrane transport process GO term. Increased OsHMA3 expression probably adds to the effect of OsNRAMP5 mutation to account for the significant decreases in Cd accumulation in rice plant and grain of the lcd1 mutant.

CONCLUSIONS

An extremely low Cd mutant lcd1 was isolated and identified using MutMap and RNA-seq. A Pro236Leu amino acid substitution in the highly conserved region of OsNRAMP5 is likely responsible for low Cd accumulation in the lcd1 mutant. This work provides more insight into the mechanism of Cd uptake and accumulation in rice, and will be helpful for developing low Cd accumulation rice by marker-assisted breeding.

A transcriptomic (RNA-seq) analysis of genes responsive to both cadmium and arsenic stress in rice root

Cadmium (Cd) and arsenic (As) are nonessential and toxic elements in rice that often occur together in contaminated paddy field soils. To understand whether rice has a common molecular response mechanism against Cd and As toxicity, 30-day seedlings (Oryza sativa L. indica) were exposed separately to Cd and As³⁺ in hydroponic cultures for up to 7 days. Root transcriptomic analysis of plants exposed to Cd and As for 3 days revealed that a total of 2224 genes in rice roots responded to Cd stress, while 1503 genes responded to As stress. Of these, 841 genes responded to both stressors. The genes in common to Cd and As stress were associated with redox control, stress response, transcriptional regulation, transmembrane transport, signal transduction, as well as biosynthesis and metabolism of macromolecules and sulfur compounds. In plants exposed to Cd and As separately or in combination for 3 and 7 days, qRT-PCR verification revealed that the glutathione metabolism associated gene Os09g0367700 was significantly up-regulated with respect to unexposed controls and had a positive synergistic effect under combined Cd and As stress. In addition, the redox control related genes Os06g0216000, Os07g0638300 and Os01g0294500, the glutathione metabolism related gene Os01g0530900, the cell wall biogenesis related genes Os05g0247800, Os11g0592000 and Os03g0416200, the expression regulation related genes Os07g0597200 and Os02g0168200, and the transmembrane transport related genes Os04g0524500, also varied significantly with respect to an unexposed control and displayed synergistic effects after 7 days of simultaneous exposure to Cd and As. Our identification of a novel set of genes in rice which responded to both Cd and As³⁺ stress may be of value in mitigating the toxicity of co-contaminated soils. These results also provide a deeper understanding of the molecular mechanisms involved in response to multi-metal/loids stress, and may be used in the genetic improvement of rice varieties.

Root Hair Single Cell Type Specific Profiles of Gene Expression and Alternative Polyadenylation Under Cadmium Stress

Transcriptional networks are tightly controlled in plant development and stress responses. Alternative polyadenylation (APA) has been found to regulate gene expression under abiotic stress by increasing the heterogeneity at mRNA 3'-ends. Heavy metals like cadmium pollute water and soil due to mining and industry applications. Understanding how plants cope with heavy metal stress remains an interesting question. The Arabidopsis root hair was chosen as a single cell model to investigate the functional role of APA in cadmium stress response. Primary root growth inhibition and defective root hair morphotypes were observed. Poly(A) tag (PAT) libraries from single cell types, i.e., root hair cells, non-hair epidermal cells, and whole root tip under cadmium stress were prepared and sequenced. Interestingly, a root hair cell type-specific gene expression under short term cadmium exposure, but not related to the prolonged treatment, was detected. Differentially expressed poly(A) sites were identified, which largely contributed to altered gene expression, and enriched in pentose and glucuronate interconversion pathways as well as phenylpropanoid biosynthesis pathways. Numerous genes with poly(A) site switching were found, particularly for functions in cell wall modification, root epidermal differentiation, and root hair tip growth. Our findings suggest that APA plays a functional role as a potential stress modulator in root hair cells under cadmium treatment.

Transcriptome analysis reveals that barnyard grass exudates increase the allelopathic potential of allelopathic and non-allelopathic rice (Oryza sativa) accessions

BACKGROUND

Allelopathy in rice (Oryza sativa) is a chemically induced response that is elevated by the exogenous application of chemical compounds and barnyard grass root exudates. An in-depth understanding of the response mechanisms of rice to chemical induction is necessary for the identification of target genes for increasing the allelopathic potential of rice. However, no previous studies have evaluated the transcriptomic changes associated with allelopathy in rice in response to barnyard grass exudates treatment. Thus, the aim of the present study was to reveal differentially expressed genes (DEGs) in allelopathic and non-allelopathic rice seedlings treated with barnyard grass exudates to identify target allelopathy genes.

RESULTS

The inhibitory effect of the culture solutions on the allelopathic rice accession PI312777 (PI) and the non-allelopathic rice accession Lemont (LE) significantly increased (P < 0.05) after treatment with barnyard grass root exudates. The RNA sequencing results revealed that 14,891 genes in PI(+B) vs. LE(+B), 12,505 genes in PI(+B) vs. PI(-B), and 5857 genes in LE(+B) vs. LE(-B) were differentially expressed following root exudates treatment. These DEGs were classified into three categories and 32 functional groups, i.e., 12 groups in the biological process category, 12 groups in the cellular component category, and eight groups in the molecular function category. There were 5857 and 2846 upregulated genes and 135 and 50 upregulated Gene Ontology terms (P < 0.05) in the biological process category in PI(+B) vs. PI(-B) and LE(+B) vs. LE(-B), respectively. These results indicated that the allelopathic accession PI is more sensitive than the non-allelopathic accession LE to exogenous root exudates treatment. Genes related to rice allelochemical-related biosynthesis pathways, particularly the shikimic acid and acetic acid pathways, were significantly differentially expressed in both rice accessions. These findings suggested that phenolic acids, fatty acids, and flavonoids, which constitute the downstream metabolites of the shikimic acid and acetic acid pathways, are significantly expressed in response to root exudates of barnyard grass.

CONCLUSIONS

The allelopathic potential of both rice accessions could be significantly enhanced by barnyard grass root exudates application. Furthermore, genes related to the biosynthesis pathways of reported rice allelochemicals were significantly differentially expressed in both accessions. Phenylalanine ammonia lyase was determined to be a potential target for the regulation of chemical induction.

Bioinformatic Exploration of the Targets of Xylem Sap miRNAs in Maize under Cadmium Stress

Cadmium (Cd) has the potential to be chronically toxic to humans through contaminated crop products. MicroRNAs (miRNAs) can move systemically in plants. To investigate the roles of long-distance moving xylem miRNAs in regulating maize response to Cd stress, three xylem sap small RNA (sRNA) libraries were constructed for high-throughput sequencing to identify potential mobile miRNAs in Cd-stressed maize seedlings and their putative targets in maize transcriptomes. In total, about 199 miRNAs (20–22 nucleotides) were identified in xylem sap from maize seedlings, including 97 newly discovered miRNAs and 102 known miRNAs. Among them, 10 miRNAs showed differential expression in xylem sap after 1 h of Cd treatment. Two miRNAs target prediction tools, psRNAtarget (reporting the inhibition pattern of cleavage) and DPMIND (discovering Plant MiRNA-Target Interaction with degradome evidence), were used in combination to identify, via bioinformatics, the targets of 199 significantly expressed miRNAs in maize xylem sap. The integrative results of these two bioinformatic tools suggested that 27 xylem sap miRNAs inhibit 34 genes through cleavage with degradome evidence. Moreover, nearly 300 other genes were also the potential miRNAs cleavable targets without available degradome data support, and the majority of them were enriched in abiotic stress response, cell signaling, transcription regulation, as well as metal handling. These approaches and results not only enhanced our understanding of the Cd-responsive long-distance transported miRNAs from the view of xylem sap, but also provided novel insights for predicting the molecular genetic mechanisms mediated by miRNAs.

Cold-Dependent Expression and Alternative Splicing of Arabidopsis Long Non-coding RNAs

Plants re-program their gene expression when responding to changing environmental conditions. Besides differential gene expression, extensive alternative splicing (AS) of pre-mRNAs and changes in expression of long non-coding RNAs (lncRNAs) are associated with stress responses. RNA-sequencing of a diel time-series of the initial response of Arabidopsis thaliana rosettes to low temperature showed massive and rapid waves of both transcriptional and AS activity in protein-coding genes. We exploited the high diversity of transcript isoforms in AtRTD2 to examine regulation and post-transcriptional regulation of lncRNA gene expression in response to cold stress. We identified 135 lncRNA genes with cold-dependent differential expression (DE) and/or differential alternative splicing (DAS) of lncRNAs including natural antisense RNAs, sORF lncRNAs, and precursors of microRNAs (miRNAs) and trans-acting small-interfering RNAs (tasiRNAs). The high resolution (HR) of the time-series allowed the dynamics of changes in transcription and AS to be determined and identified early and adaptive transcriptional and AS changes in the cold response. Some lncRNA genes were regulated only at the level of AS and using plants grown at different temperatures and a HR time-course of the first 3 h of temperature reduction, we demonstrated that the AS of some lncRNAs is highly sensitive to small temperature changes suggesting tight regulation of expression. In particular, a splicing event in TAS1a which removed an intron that contained the miR173 processing and phased siRNAs generation sites was differentially alternatively spliced in response to cold. The cold-induced reduction of the spliced form of TAS1a and of the tasiRNAs suggests that splicing may enhance production of the siRNAs. Our results identify candidate lncRNAs that may contribute to the regulation of expression that determines the physiological processes essential for acclimation and freezing tolerance.

Activation of a gene network in durum wheat roots exposed to cadmium

BACKGROUND

Among cereals, durum wheat (Triticum turgidum L. subsp. durum) accumulates cadmium (Cd) at higher concentration if grown in Cd-polluted soils. Since cadmium accumulation is a risk for human health, the international trade organizations have limited the acceptable concentration of Cd in edible crops. Therefore, durum wheat cultivars accumulating low cadmium in grains should be preferred by farmers and consumers. To identify the response of durum wheat to the presence of Cd, the transcriptomes of roots and shoots of Creso and Svevo cultivars were sequenced after a 50-day exposure to 0.5 μM Cd in hydroponic solution.

RESULTS

No phytotoxic effects or biomass reduction was observed in Creso and Svevo plants at this Cd concentration. Despite this null effect, cadmium was accumulated in root tissues, in shoots and in grains suggesting a good cadmium translocation rate among tissues. The mRNA sequencing revealed a general transcriptome rearrangement after Cd treatment and more than 7000 genes were found differentially expressed in root and shoot tissues. Among these, the up-regulated genes in roots showed a clear correlation with cadmium uptake and detoxification. In particular, about three hundred genes were commonly up-regulated in Creso and Svevo roots suggesting a well defined molecular strategy characterized by the transcriptomic activation of several transcription factors mainly belonging to bHLH and WRKY families. bHLHs are probably the activators of the strong up-regulation of three NAS genes, responsible for the synthesis of the phytosiderophore nicotianamine (NA). Moreover, we found the overall up-regulation of the methionine salvage pathway that is tightly connected with NA synthesis and supply the S-adenosyl methionine necessary for NA biosynthesis. Finally, several vacuolar NA chelating heavy metal transporters were vigorously activated.

CONCLUSIONS

In conclusion, the exposure of durum wheat to cadmium activates in roots a complex gene network involved in cadmium translocation and detoxification from heavy metals. These findings are confident with a role of nicotianamine and methionine salvage pathway in the accumulation of cadmium in durum wheat.

Activation of rice WRKY transcription factors: an army of stress fighting soldiers?

Rice WRKYs comprise a large family of transcription factors and present remarkable structure features and a unique DNA binding site. Their importance in plants goes beyond the response to stressful stimuli, since they participate in hormonal pathways and developmental processes. Indeed, the majority of WRKYs present an independent activation since they are able to perform self-transcriptional regulation. However, some WRKY activation depends on epigenetic and transcript regulation by micro RNAs. Their protein function depends, almost always, on the posttranslational changes. Taking to account its properties of auto-activation, all these regulators process are extremely important for complete WRKY regulation. In this sense, here we provide an overview of transcriptional activation and posttranscriptional and posttranslational regulation of rice WRKY genes under stresses.

Comparative analysis of Cd-responsive maize and rice transcriptomes highlights Cd co-modulated orthologs

Background

Metal tolerance is often an integrative result of metal uptake and distribution, which are fine-tuned by a network of signaling cascades and metal transporters. Thus, with the goal of advancing the molecular understanding of such metal homeostatic mechanisms, comparative RNAseq-based transcriptome analysis was conducted to dissect differentially expressed genes (DEGs) in maize roots exposed to cadmium (Cd) stress.

Results

To unveil conserved Cd-responsive genes in cereal plants, the obtained 5166 maize DEGs were compared with 2567 Cd-regulated orthologs in rice roots, and this comparison generated 880 universal Cd-responsive orthologs groups composed of 1074 maize DEGs and 981 rice counterparts. More importantly, most of the orthologous DEGs showed coordinated expression pattern between Cd-treated maize and rice, and these include one large orthologs group of pleiotropic drug resistance (PDR)-type ABC transporters, two clusters of amino acid transporters, and 3 blocks of multidrug and toxic compound extrusion (MATE) efflux family transporters, and 3 clusters of heavy metal-associated domain (HMAD) isoprenylated plant proteins (HIPPs), as well as all 4 groups of zinc/iron regulated transporter protein (ZIPs). Additionally, several blocks of tandem maize paralogs, such as germin-like proteins (GLPs), phenylalanine ammonia-lyases (PALs) and several enzymes involved in JA biosynthesis, displayed consistent co-expression pattern under Cd stress.

Out of the 1074 maize DEGs, approximately 30 maize Cd-responsive genes such as ZmHIPP27, stress-responsive NAC transcription factor (ZmSNAC1) and 9-cis-epoxycarotenoid dioxygenase (NCED, vp14) were also common stress-responsive genes reported to be uniformly regulated by multiple abiotic stresses. Moreover, the aforementioned three promising Cd-upregulated genes with rice counterparts were identified to be novel Cd-responsive genes in maize.

Meanwhile, one maize glutamate decarboxylase (ZmGAD1) with Cd co-modulated rice ortholog was selected for further analysis of Cd tolerance via heterologous expression, and the results suggest that ZmGAD1 can confer Cd tolerance in yeast and tobacco leaves.

Conclusions

These novel findings revealed the conserved function of Cd-responsive orthologs and paralogs, which would be valuable for elucidating the genetic basis of the plant response to Cd stress and unraveling Cd tolerance genes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5109-8) contains supplementary material, which is available to authorized users.

Extensive Differential Splicing Underlies Phenotypically Plastic Aphid Morphs

Phenotypic plasticity results in a diversity of phenotypes from a single genotype in response to environmental cues. To understand the molecular basis of phenotypic plasticity, studies have focused on differential gene expression levels between environmentally determined phenotypes. The extent of alternative splicing differences among environmentally determined phenotypes has largely been understudied. Here, we study alternative splicing differences among plastically produced morphs of the pea aphid using RNA-sequence data. Pea aphids express two separate polyphenisms (plasticity with discrete phenotypes): a wing polyphenism consisting of winged and wingless females and a reproduction polyphenism consisting of asexual and sexual females. We find that pea aphids alternatively splice 34% of their genes, a high percentage for invertebrates. We also find that there is extensive use of differential spliced events between genetically identical, polyphenic females. These differentially spliced events are enriched for exon skipping and mutually exclusive exon events that maintain the open reading frame, suggesting that polyphenic morphs use alternative splicing to produce phenotype-biased proteins. Many genes that are differentially spliced between polyphenic morphs have putative functions associated with their respective phenotypes. We find that the majority of differentially spliced genes is not differentially expressed genes. Our results provide a rich candidate gene list for future functional studies that would not have been previously considered based solely on gene expression studies, such as ensconsin in the reproductive polyphenism, and CAKI in the wing polyphenism. Overall, this study suggests an important role for alternative splicing in the expression of environmentally determined phenotypes.

A comprehensive analysis of the Korean fir (Abies koreana) genes expressed under heat stress using transcriptome analysis

Korean fir (Abies koreana), a rare species endemic to South Korea, is sensitive to climate change. Here, we used next-generation massively parallel sequencing technology and de novo transcriptome assembly to gain a comprehensive overview of the Korean fir transcriptome under heat stress. Sequencing control and heat-treated samples of Korean fir, we obtained more than 194,872,650 clean reads from each sample. After de novo assembly and quantitative assessment, 42,056 unigenes were generated with an average length of 908 bp. In total, 6,401 differentially expressed genes were detected, of which 2,958 were up-regulated and 3,443 down-regulated, between the heat-treated and control samples. A gene ontology analysis of these unigenes revealed heat-stress-related terms, such as "response to stimulus". Further, in depth analysis revealed 204 transcription factors and 189 Hsps as differentially expressed. Finally, 12 regulated candidate genes associated with heat stress were examined using quantitative real-time PCR (qRT-PCR). In this study, we present the first comprehensive characterisation of Korean fir subjected to heat stress using transcriptome analysis. It provides an important resource for future studies of Korean fir with the objective of identifying heat stress tolerant lines.

Methylmercury exposure alters RNA splicing in human neuroblastoma SK-N-SH cells: Implications from proteomic and post-transcriptional responses

The neurotoxic effects of methylmercury (MeHg) have been intensively studied. However, the molecular mechanisms responsible for the neurotoxicity of MeHg are not fully understood. To decipher these mechanisms, proteomic and high-throughput mRNA sequencing (RNA-seq) technique were utilized, comprehensively evaluating the cellular responses of human neuroblastoma SK-N-SH cells to MeHg exposure. Proteomic results revealed that MeHg exposure interfered with RNA splicing via splicesome, along with the known molecular mechanisms of mercury-related neurotoxicity (e.g. oxidative stress, protein folding, immune system processes, and cytoskeletal organization). The effects of MeHg on RNA splicing were further verified using RNA-seq. Compared to control, a total of 658 aberrant RNA alternative splicing (AS) events were observed after MeHg exposure. Proteomics and RNA-seq results also demonstrated that mercury chloride (HgCl₂) influenced the expression levels of several RNA splicing related proteins and 676 AS events compared to control. These results suggested that RNA splicing could be a new molecular mechanism involved in MeHg and HgCl₂ neurotoxicity.

Genome-wide analysis of long non-coding RNAs affecting roots development at an early stage in the rice response to cadmium stress

Background

Long non-coding RNAs (lncRNAs) have been found to play a vital role in several gene regulatory networks involved in the various biological processes in plants related to stress response. However, systematic analyses of lncRNAs expressed in rice Cadmium (Cd) stress are seldom studied. Thus, we presented the characterization and expression of lncRNAs in rice root development at an early stage in response to Cd stress.

Results

The lncRNA deep sequencing revealed differentially expressed lncRNAs among Cd stress and normal condition. In the Cd stress group, 69 lncRNAs were up-regulated and 75 lncRNAs were down-regulated. Furthermore, 386 matched lncRNA-mRNA pairs were detected for 120 differentially expressed lncRNAs and 362 differentially expressed genes in cis, and target gene-related pathway analyses exhibited significant variations in cysteine and methionine metabolism pathway-related genes. For the genes in trans, overall, 28,276 interaction relationships for 144 lncRNAs and differentially expressed protein-coding genes were detected. The pathway analyses found that secondary metabolites, such as phenylpropanoids and phenylalanine, and photosynthesis pathway-related genes were significantly altered by Cd stress. All of these results indicate that lncRNAs may regulate genes of cysteine-rich peptide metabolism in cis, as well as secondary metabolites and photosynthesis in trans, to activate various physiological and biochemical reactions to respond to excessive Cd.

Conclusion

The present study could provide a valuable resource for lncRNA studies in response to Cd treatment in rice. It also expands our knowledge about lncRNA biological function and contributes to the annotation of the rice genome.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4807-6) contains supplementary material, which is available to authorized users.

Long Non-Coding RNAs Responsive to Salt and Boron Stress in the Hyper-Arid Lluteño Maize from Atacama Desert

Long non-coding RNAs (lncRNAs) have been defined as transcripts longer than 200 nucleotides, which lack significant protein coding potential and possess critical roles in diverse cellular processes. Long non-coding RNAs have recently been functionally characterized in plant stress-response mechanisms. In the present study, we perform a comprehensive identification of lncRNAs in response to combined stress induced by salinity and excess of boron in the Lluteño maize, a tolerant maize landrace from Atacama Desert, Chile. We use deep RNA sequencing to identify a set of 48,345 different lncRNAs, of which 28,012 (58.1%) are conserved with other maize (B73, Mo17 or Palomero), with the remaining 41.9% belonging to potentially Lluteño exclusive lncRNA transcripts. According to B73 maize reference genome sequence, most Lluteño lncRNAs correspond to intergenic transcripts. Interestingly, Lluteño lncRNAs presents an unusual overall higher expression compared to protein coding genes under exposure to stressed conditions. In total, we identified 1710 putatively responsive to the combined stressed conditions of salt and boron exposure. We also identified a set of 848 stress responsive potential trans natural antisense transcripts (trans-NAT) lncRNAs, which seems to be regulating genes associated with regulation of transcription, response to stress, response to abiotic stimulus and participating of the nicotianamine metabolic process. Reverse transcription-quantitative PCR (RT-qPCR) experiments were performed in a subset of lncRNAs, validating their existence and expression patterns. Our results suggest that a diverse set of maize lncRNAs from leaves and roots is responsive to combined salt and boron stress, being the first effort to identify lncRNAs from a maize landrace adapted to extreme conditions such as the Atacama Desert. The information generated is a starting point to understand the genomic adaptabilities suffered by this maize to surpass this extremely stressed environment.

Oxidative Stress and Heavy Metals in Plants

Oxidative stress is a pathological process related to not only animal kingdom but also plants. Regarding oxidative stress in plants, heavy metals are frequently discussed as causative stimuli with relevance to ecology. Because heavy metals have broad technological importance, they can easily contaminate the environment. Much of previous effort regarding the harmful impact of the heavy metals was given to their toxicology in the animals and humans. Their implication in plant pathogeneses is less known and remains underestimated.The current paper summarizes basic facts about heavy metals, their distribution in soil, mobility, accumulation by plants, and initiation of oxidative stress including the decline in basal metabolism. The both actual and frontier studies in the field are summarized and discussed. The major pathophysiological pathways are introduced as well and link between heavy metals toxicity and their ability to initiate an oxidative damage is provided. Mobility and bioaccessibility of the metals is also considered as key factors in their impact on oxidative stress development in the plant. The metals like lead, mercury, copper, cadmium, iron, zinc, nickel, vanadium are depicted in the text.Heavy metals appear to be significant contributors to pathological processes in the plants and oxidative stress is probably an important contributor to the effect. The most sensitive plant species are enlisted and discussed in this review. The facts presented here outline next effort to investigate pathological processes in the plants.