Transcriptomic changes during maize roots development responsive to Cadmium (Cd) pollution using comparative RNAseq-based approach

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.

Comparative transcriptome analysis reveals key genes and coordinated mechanisms in two rice cultivars differing in cadmium accumulation

Although Cd accumulation varies among rice varieties is recognized, the underlying mechanisms are not well clarified. In this study, comparative transcriptome analysis were performed by hydroponic culture system with two rice varieties, Y1540 (high Cd accumulator) and Y15 (low Cd accumulator) under 20 μM Cd stress. Results revealed 17,320 differentially expressed genes (DEGs) in roots of Y15 (7,655 upregulated and 9,665 downregulated) and 17,386 DEGs in roots of Y1540 (8,823 upregulated and 8,563 downregulated) expose to 20 μM Cd stress. Gene ontology (GO) analysis enriched 24 and 26 terms in Y15 and Y1540 respectively, including 23 common terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment showed 27 and 28 significant pathways in Y15 and Y1540 respectively, with 19 common pathways. Different responses to Cd stress between cultivars were not only reflected in differently enriched GO terms and KEGG pathways but also in different DEGs of 23 common GO terms and significant sequences represented by p-values of 19 common KEGG pathways. Both cultivars resist Cd through common processes with different weights; hence glutathione metabolism, mineral absorption, biosynthesis of secondary metabolites, and degradation of aromatic compounds could be playing a more important role in Y1540, whereas ribosome biogenesis in eukaryotes, mismatch repair, aminoacyl-tRNA biosynthesis, and the cell cycle maybe playing a more important role in Y15. Weighted gene co-expression network analysis (WGCNA) showed that five and three modules were clustered in Y15 and Y1540, respectively, with yellow and brown modules in Y15 and brown modules in Y1540 being significantly related to Cd stress. Further analysis showed that most of hub genes in Y15 were related to signal transduction or transcription factors, while most of hub genes in Y1540 were related to binding, metabolic, and secondary metabolic processes, which demonstrated their different response patterns at transcriptomic level to Cd stress.

Integrated transcriptome and metabonomic analysis of key metabolic pathways in response to cadmium stress in novel buckwheat and cultivated species

INTRODUCTION

Buckwheat (Fagopyrum tataricum), an important food crop, also has medicinal uses. It is widely planted in Southwest China, overlapping with planting areas remarkably polluted by cadmium (Cd). Therefore, it is of great significance to study the response mechanism of buckwheat under Cd stress and further develop varieties with excellent Cd tolerance.

METHODS

In this study, two critical periods of Cd stress treatment (days 7 and 14 after Cd treatment) of cultivated buckwheat (Pinku-1, named K33) and perennial species (F. tatari-cymosum Q.F. Chen) (duoku, named DK19) were analyzed using transcriptome and metabolomics.

RESULTS

The results showed that Cd stress led to changes in reactive oxygen species (ROS) and the chlorophyll system. Moreover, Cd-response genes related to stress response, amino acid metabolism, and ROS scavenging were enriched or activated in DK19. Transcriptome and metabolomic analyses highlighted the important role of galactose, lipid (glycerophosphatide metabolism and glycerophosphatide metabolism), and glutathione metabolism in response to Cd stress in buckwheat, which are significantly enriched at the gene and metabolic levels in DK19.

DISCUSSION

The results of the present study provide valuable information for a better understanding of the molecular mechanisms underlying Cd tolerance in buckwheat and useful clues for the genetic improvement of drought tolerance in buckwheat.

WGCNA Analysis Revealed the Hub Genes Related to Soil Cadmium Stress in Maize Kernel (Zea mays L.)

Soil contamination by heavy metals has become a prevalent topic due to their widespread release from industry, agriculture, and other human activities. Great progress has been made in elucidating the uptake and translocation of cadmium (Cd) accumulation in rice. However, there is still little known about corresponding progress in maize. In the current study, we performed a comparative RNA-Seq-based approach to identify differentially expressed genes (DEGs) of maize immature kernel related to Cd stress. In total, 55, 92, 22, and 542 DEGs responsive to high cadmium concentration soil were identified between XNY22-CHS-8 vs. XNY22-YA-8, XNY22-CHS-24 vs. XNY22-YA-24, XNY27-CHS-8 vs. XNY27-YA-8, and XNY27-CHS-24 vs. XNY27-YA-24, respectively. The weighted gene co-expression network analysis (WGCNA) categorized the 9599 Cd stress-responsive hub genes into 37 different gene network modules. Combining the hub genes and DEGs, we obtained 71 candidate genes. Gene Ontology (GO) enrichment analysis of genes in the greenyellow module in XNY27-YA-24 and connectivity genes of these 71 candidate hub genes showed that the responses to metal ion, inorganic substance, abiotic stimulus, hydrogen peroxide, oxidative stress, stimulus, and other processes were enrichment. Moreover, five candidate genes that were responsive to Cd stress in maize kernel were detected. These results provided the putative key genes and pathways to response to Cd stress in maize kernel, and a useful dataset for unraveling the underlying mechanism of Cd accumulation in maize kernel.

Metabolomic and transcriptomic analyses reveal that sucrose synthase regulates maize pollen viability under heat and drought stress

Maize pollen is highly sensitive to heat and drought, but few studies have investigated the combined effects of heat and drought on pollen viability. In this study, pollen's structural and physiological characteristics were determined after heat, drought, and combined stressors. Furthermore, integrated metabolomic and transcriptomic analyses of maize pollen were conducted to identify potential mechanisms of stress responses. Tassel growth and spikelet development were considerably suppressed, pollen viability was negatively impacted, and pollen starch granules were depleted during anthesis under stress. The inhibitory effects were more significant due to combined stresses than to heat or drought individually. The metabolic analysis identified 71 important metabolites in the combined stress compared to the other treatments, including sugars and their derivatives related to pollen viability. Transcriptomics also revealed that carbohydrate metabolism was significantly altered under stress. Moreover, a comprehensive metabolome-transcriptome analysis identified a central mechanism in the biosynthesis of UDP-glucose involved in reducing the activity of sucrose synthase SH-1 (shrunken 1) and sus1 (sucrose synthase 1) that suppressed sucrose transfer to UDP-glucose, leading to pollen viability exhaustion under stress. In conclusion, the lower pollen viability after heat and drought stress was associated with poor sucrose synthase activity due to the stress treatments.

Genome-wide analysis of metallothionein gene family in maize to reveal its role in development and stress resistance to heavy metal

BACKGROUND

Maize (Zea mays L.) is a widely cultivated cereal and has been used as an optimum heavy metal phytoremediation crop. Metallothionein (MT) proteins are small, cysteine-rich, proteins that play important roles in plant growth and development, and the regulation of stress response to heavy metals. However, the MT genes for maize have not been fully analyzed so far.

METHODS

The putative ZmMT genes were identified by HMMER.The heat map of ZmMT genes spatial expression analysis was generated by using R with the log2 (FPKM + 1).The expression profiles of ZmMT genes under three kinds of heavy metal stresses were quantified by using qRT-PCR. The metallothionein proteins was aligned using MAFFT and phylogenetic analysis were constructed by ClustalX 2.1. The protein theoretical molecular weight and pI, subcellular localization, TFs binding sites, were predicted using ProtParam, PSORT, PlantTFDB, respectively.

RESULTS

A total of 9 ZmMT genes were identified in the whole genome of maize. The results showed that eight of the nine ZmMT proteins contained one highly conserved metallothio_2 domain, while ZmMT4 contained a Metallothio_PEC domain. All the ZmMT proteins could be classified into three major groups and located on five chromosomes. The ZmMT promoters contain a large number of hormone regulatory elements and hormone-related transcription factor binding sites. The ZmMT genes exhibited spatiotemporal specific expression patterns in 23 tissues of maize development stages and showed the different expression patterns in response to Cu, Cd, and Pb heavy metal stresses.

CONCLUSIONS

We identified the 9 ZmMT genes, and explored their conserved motif, tissue expression patterns, evolutionary relationship. The expression profiles of ZmMT genes under three kinds of heavy metal stresses (Cu, Cd, Pb) were analyzed. In summary, the expression of ZmMTs have poteintial to be regulated by hormones. The specific expression of ZmMTs in different tissues of maize and the response to different heavy metal stresses are revealed that the role of MT in plant growth and development, and stress resistance to heavy metals.

Weighted gene co-expression network analysis unveils gene networks regulating folate biosynthesis in maize endosperm

Folates are essential elements for human growth and development, and their deficiency can lead to serious disorders. Waxy maize is a rich source of folates; however, the regulatory mechanism underlying folate biosynthesis in the endosperm remains unclear. Here, we examined changes in the folate content of maize endosperm collected at 15, 18, 21, 24, and 27 days after pollination (DAP) using liquid chromatograph-mass spectrometry and identified genes related to folate biosynthesis using transcriptome sequencing data. The results showed that 5-methyl-tetrahydrofolate and 5,10-methylene tetrahydrofolate were the main storage forms of folates in the endosperm, and their contents were relatively high at 21-24 days. We also identified 569, 3183, 4365, and 5513 differentially expressed genes (DEGs) in different days around milk stage. Functional annotation revealed 518 transcription factors (TFs) belonging to 33 families exhibiting specific expression in at least one sampling time. The key hub genes involved in folate biosynthesis were identified by weighted gene co-expression network analysis. In total, 24,976 genes were used to construct a co-expression network with 29 co-expression modules, among which the brown and purple modules were highly related to folate biosynthesis. Further, 187 transcription factors in the brown and purple modules were considered potential transcription factors related to endosperm folate biosynthesis. These results may improve the understanding of the molecular mechanism underlying folate biosynthesis in waxy maize and lead to the development of nutritionally fortified varieties.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02974-7.

Genetic dissection of maize seedling traits in an IBM Syn10 DH population under the combined stress of lead and cadmium

The heavy metals lead and cadmium have become important pollutants in the environment, which exert negative effects on plant morphology, growth and photosynthesis. It is particularly significant to uncover the genetic loci and the causal genes for lead and cadmium tolerance in plants. This study used an IBM Syn10 DH population to identify the quantitative trait loci (QTL) controlling maize seedling tolerance to lead and cadmium by linkage mapping. The broad-sense heritability of these seedling traits ranged from 65.8-97.3% and 32.0-98.8% under control (CK) and treatment (T) conditions, respectively. A total of 53 and 64 QTL were detected under CK and T conditions, respectively. Moreover, 42 QTL were identified using lead and cadmium tolerance coefficient (LCTC). Among these QTL, five and two major QTL that explained > 10% of phenotypic variation were identified under T condition and using LCTC, respectively. Furthermore, eight QTL were simultaneously identified by T and LCTC, explaining 5.23% to 9.21% of the phenotypic variations. Within these major and common QTL responsible for the combined heavy metal tolerance, four candidate genes (Zm00001d048759, Zm00001d004689, Zm00001d004843, Zm00001d033527) were previously reported to correlate with heavy metal transport and tolerance. These findings will contribute to functional gene identification and molecular marker-assisted breeding for improving heavy metal tolerance in maize.

Main Molecular Pathways Associated with Copper Tolerance Response in Imperata cylindrica by de novo Transcriptome Assembly

The metallophyte Imperata cylindrica inhabits copper (Cu) polluted soils in large areas from Central Chile. Here, we subjected clonal vegetative plantlets to 300 mg Cu kg-1 of substrate for 21 days to identify the main molecular pathways involved in the response to Cu stress. Transcriptomic analyses were performed for shoots and roots, with and without Cu supply. RNA-Seq and de novo transcriptome assembly were performed to identify the gene response associated with molecular mechanisms of Cu tolerance in I. cylindrica. De novo transcriptome revealed a total of 200,521 transcripts (1777 bp) comprising ~91% complete ultra-conserved genes in the eukaryote and Plantae database. The differentially expressed genes (DEGs) in roots were 7386, with 3558 of them being up-regulated and the other 3828 down-regulated. The transcriptome response in shoots was significantly less, showing only 13 up-regulated and 23 down-regulated genes. Interestingly, DEGs mainly related with actin and cytoskeleton formation, and to a minor degree, some DEGs associated with metal transporters and superoxide dismutase activity in root tissues were found. These transcriptomic results suggest that cytoskeleton could be acting as a mechanism of Cu-binding in the root, resulting in a high Cu tolerance response in this metallophyte, which deserve to be analyzed ultra-structurally. Our study contributes to reinforcing the potential of I. cylindrica as a candidate plant species to be used as a phytoremediation agent in Cu-contaminated environments.

Quantitative proteomic analysis reveals complex regulatory and metabolic response of Iris lactea Pall. var. chinensis to cadmium toxicity

Cadmium pollution has become a serious environmental problem. Iris lactea var. chinensis showed strong Cd tolerance and accumulation ability, which has significant potential to be applied for the phytoremediation of Cd-contaminated soil. However, the lack of molecular information on the mechanism of I. lactea response to Cd limited the improvement of phytoremediation efficiency. In this study, label-free proteomics analysis of Cd response in I. lactea showed that there were 163 and 196 differentially expressed proteins (DEPs) in the shoots and roots, respectively. Bioinformatics analysis indicated the DEPs responding to Cd stress mainly involved in signal transduction, ion transport, redox etc., and participate in the pathway of amino acid biosynthesis, lignin biosynthesis, glycerolipid metabolism and glutathione metabolism. Besides, differential expression of seven DEPs was validated via gene expression analysis. Finally, we found that a Cd-induced mannose-specific lectin (IlMSL) from I. lactea enhanced the Cd sensitivity and increased Cd accumulation in yeast. The results of this study will enhance our understanding of the molecular mechanism of Cd tolerance and accumulation in I. lactea and ultimately provide valuable resources for using Cd tolerant genes for developing efficient strategies for phytoremediation of Cd-contaminated soils or limiting Cd accumulation in food crops.

A transcriptomic analysis reveals the adaptability of the growth and physiology of immature tassel to long-term soil water deficit in Zea mays L

Drought is a key threat to maize growth and yield. Understanding the mechanism of immature tassel (IT) response to long term drought is of paramount importance. Here, the maize inbred line PH6WC was tested under well-watered (CK) and two water deficit treatments (WD1 and WD2). The final IT length in the WD1 and WD2 treatments decreased by nearly 6.2% and 21.2% compared to the CK, respectively, and the average accumulation rate IT dry matter was 1.5-fold and 1.8-fold slower, respectively. Furthermore, RNA sequencing analysis was conducted on the IT sampled at 30 days after the WD treatments. In total, the cellular component in gene ontology (GO) analysis suggested that the differentially expressed genes were significantly enriched in three common terms (apoplast, plant-type cell wall, and anchored component of membrane) among the CK vs WD1, CK vs WD2, and WD1 vs WD2 comparisons. Next, a co-expression network analysis identified 44 modules that contained global expression genes. Finally, by combining the GO analysis with modules, nine genes involved in carbohydrate metabolism and the antioxidant system were screened out, and the six corresponding physiological parameters were all significantly increased under the WD treatments. These results showed that, although the IT length and dry matter decreased, the IT enhanced the adaptation to drought by regulating their own genetic and physiological changes.

Physiological and transcriptomic analyses of cadmium stress response in Dendrobium officinale seedling

Dendrobium officinale is an economically important Chinese herb with ornamental and medicinal values. However, the mechanisms by which D. officinale adapts to cadmium (Cd) stress is unknown. Here, physiological changes in D. officinale roots and leaves exposed to increasing levels of Cd stress (CdSO₄ concentration of 2, 5, 9, 14 mg L^-1) were analyzed at 7, 15, 30, and 45 days after treatment. The Cd stress of 14 mg L^-1 significantly increased the levels of antioxidants and induced malondialdehyde and proline accumulation (P < 0.05). Cd subcellular distribution showed that Cd sequestration into soluble fraction is the major detoxification mechanism in D. officinale roots. Subsequently, the transcriptome profile of D. officinale roots treated with 14 mg L^-1 Cd for 15 and 30 days was analyzed. Compared to control, 2,469 differentially expressed genes (DEGs) were identified, comprising 1,486 up-regulated genes and 983 down-regulated genes. The DEGs associated with metabolic pathways for Cd uptake, transportation and detoxification were analyzed. Several processes such as metal transporter, sulfate glutathione metabolism, cell wall metabolism, phenylpropanoid metabolism were identified to be important for Cd stress adaptation. More genes were expressed at 15 days after treatment compared to 30 days. WRKY, Trihelix, NF-YC, MYB, bZIP and bHLH transcription factors were over-expressed at both time points. Furthermore, candidate genes from the glutathione metabolism pathway were identified, and qRT-PCR analysis of ten DEGs indicated a high coorelation with RNA-seq expression profiles. Our findings provide significant information for further research of Cd stress responsive genes functions in D. officinale, especially the genes from the glutathione metabolism pathway.

Transcriptome -wide modulation combined with morpho-physiological analyses of Typha orientalis roots in response to lead challenge

Lead (Pb) is a common pollutant in many environments, including in the soil, water, and/or air. Typha orientalis Presl, a large emergent aquatic plant, has been reported to function as a Pb-tolerant and Pb-accumulating plant; however, very little molecular information regarding the tolerance of T. orientalis towards Pb is known. In this study, Pb accumulation and key factors involved in the Pb stress response at different Pb concentrations were investigated. Pb was primarily accumulated in the roots and was mainly located in the cell wall and membrane systems. Differentially expressed genes (DEGs) were identified in T. orientalis roots after Pb exposure via RNA-seq analyses. In the 0.10 mM and 0.25 mM Pb²⁺-treated groups, a total of 3275 DEGs were detected relative to the control. Many of these genes were associated with oxidation-reduction processes, metal transport, protein kinase/phosphorylation, and DNA binding transcription factors, which were shown to be Pb-responsive DEGs. Mapping Kyoto Encyclopedia of Genes and Genomes (KEGG) database, "phenylpropanoid biosynthesis" was analyzed as the major pathway of the important modules of overlapping DEGs of 0.10 mM and 0.25 mM Pb²⁺ treatments. Furthermore, a lead response gene named ToLR1 with unknown function was of particular interest. The full-length of ToLR1 sequence was cloned using rapid amplification of cDNA ends (RACE) and overexpressed in Arabidopsis thaliana, which resulted in enhanced resistance to Pb stress. This is the first report providing genomic information detailing Pb responsive genes in T. orientalis. Moreover, this study provides novel insights into the molecular mechanisms underlying the response of T. orientalis and other accumulators towards Pb stress. The key genes identified in this study may serve as potential targets for genetic engineering targeting phytoremediation.

The Endophytic Fungus Chaetomium cupreum Regulates Expression of Genes Involved in the Tolerance to Metals and Plant Growth Promotion in Eucalyptus globulus Roots

The endophytic strain Chaetomium cupreum isolated from metal-contaminated soil was inoculated in Eucalyptus globulus roots to identify genes involved in metal stress response and plant growth promotion. We analyzed the transcriptome of E. globulus roots inoculated with C. cupreum. De novo sequencing, assembly, and analysis were performed to identify molecular mechanisms involved in metal stress tolerance and plant growth promotion. A total of 393,371,743 paired-end reads were assembled into 135,155 putative transcripts. It was found that 663 genes significantly changed their expression in the presence of treatment, of which 369 were up-regulated and 294 were down-regulated. We found differentially expressed genes (DEGs) encoding metal transporters, transcription factors, stress and defense response proteins, as well as DEGs involved in auxin biosynthesis and metabolism. Our results showed that the inoculation of C. cupreum enhanced tolerance to metals and growth promotion on E. globulus. This study provides new information to understand molecular mechanisms involved in plant–microbe interactions under metals stress.

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.

Immunosuppression, oxidative stress, and glycometabolism disorder caused by cadmium in common carp (Cyprinus carpio L.): Application of transcriptome analysis in risk assessment of environmental contaminant cadmium

Cadmium (Cd), a hazardous environmental contaminant with irreversible toxicity to fish, has been detected in aquatic environment of many countries. The common carp is one of the most widely distributed fish in the world, so we used common carp to assess environmental contaminant risk. In present study, we investigated effects of Cd on immune function, oxidative defense, and glycometabolism in the spleens of common carp by transcriptome analysis. Obtained 3794 differentially expressed genes (including 1848 up-regulated and 1946 down-regulated genes) were enriched using databases of Kyoto Encyclopedia of Genes and Genomes, and Gene Ontology in David bioinformatics software (version 6.8). The pathways and gene functions of immune, oxidative defense, and glycometabolism were obtained and identified. Some relative genes were validated using qRT-PCR and gene expression of IL-1β, INF-γ, IL-6, Cxcl18b, HO-1a, CAT, GPx1, GCK, and FBA decreased; and gene expression of B4GALT1, GPAT3, and CYP26B1 increased. Our results indicated that Cd exposure led to immunosuppression, oxidative stress, and glycometabolism disorder in the common carp spleens. The present study gives a novel insight and method on environmental risk assessment.

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.

Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus

BACKGROUND

Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus associated with betasatellite (Tobacco curly shoot betasatellite, TbCSB), which causes serious leaf curl disease on tomato and tobacco in China. It is interesting that TbCSV induced severe upward leaf curling in Nicotiana benthamiana, but in the presence of TbCSB, symptoms changed to be downward leaf curling. However, the mechanism of interactions between viral pathogenicity, host defense, viral-betasatellite interactions and virus-host interactions remains unclear.

METHODS

In this study, RNA-seq was used to analyze differentially expressed genes (DEGs) in N. benthamiana plants infected by TbCSV (Y35A) and TbCSV together with TbCSB (Y35AB) respectively.

RESULTS

Through mapping to N. benthamiana reference genome, 59,814 unigenes were identified. Transcriptome analysis revealed that a total of 4081 and 3196 DEGs were identified in Y35AB vs CK (control check) and Y35A vs CK, respectively. Both GO and KEGG analyses were conducted to classify the DEGs. Ten of the top 15 GO terms were enriched in both DEGs of Y35AB vs CK and Y35A vs CK, and these enriched GO terms mainly classified into three categories including biological process, cellular component and molecular function. KEGG pathway analysis indicated that 118 and 111 pathways were identified in Y35AB vs CK and Y35A vs CK, respectively, of which nine and six pathways were significantly enriched. Three major pathways in Y35AB vs CK involved in metabolic pathways, carbon metabolism and photosynthesis, while those in Y35A vs CK were related to Ribosome, Glyoxylate and dicarboxylate metabolism and DNA replication. We observed that 8 PR genes were significantly up-regulated and 44 LRR-RLK genes were significantly differentially expressed in Y35A treatment or in Y35AB treatment. In addition, 7 and 13 genes were identified to be significantly changed in biosynthesis and signal transduction pathway of brassinosteroid (BR) and jasmonic acid (JA) respectively.

CONCLUSIONS

These results presented here would be particularly useful to further elucidate the response of the host plant against virus infection.

Transcriptome analysis for identifying possible gene regulations during maize root emergence and formation at the initial growth stage

The root plays an important role during plant development and growth, i.e., the plant body maintenance, nutrient storage, absorption of water, oxygen and nutrient from the soil, and storage of water and carbohydrates, etc. The objective of this study was attempted to determine root-specific genes at the initial developmental stages of maize by using network-based transcriptome analysis. The raw data obtained using RNA-seq were filtered for quality control of the reads with the FASTQC tool, and the filtered reads were pre-proceed using the TRIMMOMATIC tool. The enriched BINs of the DEGs were detected using PageMan analysis with the ORA_FISHER statistical test, and genes were assigned to metabolic pathways by using the MapMan tool, which was also used for detecting transcription factors (TFs). For reconstruction of the co-expression network, we used the algorithm for the reconstruction of accurate cellular networks (ARACNE) in the R package, and then the reconstructed co-expression network was visualized using the Cytoscape tool. RNA-seq. was performed using maize shoots and roots at different developmental stages of root emergence (6-10 days after planting, VE) and 1 week after plant emergence (V2). A total of 1286 differentially expressed genes (DEGs) were detected in both tissues. Many DEGs involved in metabolic pathways exhibited altered mRNA levels between VE and V2. In addition, we observed gene expression changes for 113 transcription factors and found five enriched cis-regulatory elements in the 1-kb upstream regions of both DEGs. The network-based transcriptome analysis showed two modules as co-expressed gene clusters differentially expressed between the shoots and roots during plant development. The DEGs of one module exhibited gene expressional coherence in the maize root tips, suggesting that their functional relationships are associated with the initial developmental stage of the maize root. Finally, we confirmed reliable mRNA levels of the hub genes in the potential sub-network related to initial root development at the different developmental stages of VE, V2, and 2 weeks after plant emergence.

Co-expression network analysis of the transcriptomes of rice roots exposed to various cadmium stresses reveals universal cadmium-responsive genes

BACKGROUND

The migration of cadmium (Cd) from contaminated soil to rice is a cause for concern. However, the molecular mechanism underlying the response of rice roots to various Cd stresses remains to be clarified from the viewpoint of the co-expression network at a system-wide scale.

RESULTS

We employed a comparative RNAseq-based approach to identify early Cd-responsive differentially expressed genes (DEGs) in rice 'Nipponbare' seedling roots after 1 h of high-Cd treatment. A multiplicity of the identified 1772 DEGs were implicated in hormone signaling and transcriptional regulation, particularly NACs and WRKYs were all upregulated under Cd stress. All of the 6 Cd-upregulated ABC transporters were pleiotropic drug resistance proteins (PDRs), whereas all of the 6 ZRT/IRT-like proteins (ZIPs) were consistently downregulated by Cd treatment. To further confirm our results of this early transcriptomic response to Cd exposure, we then conducted weighted gene co-expression network analysis (WGCNA) to re-analyze our RNAseq data in combination with other 11 previously published RNAseq datasets for rice roots exposed to diverse concentrations of Cd for extended treatment periods. This integrative approach identified 271 transcripts as universal Cd-regulated DEGs that are key components of the Cd treatment coupled co-expression module. A global view of the 164 transcripts with annotated functions in pathway networks revealed several Cd-upregulated key functional genes, including transporter ABCG36/OsPDR9, 12-oxo-phytodienoic acid reductases (OPRs) for JA synthesis, and ZIM domain proteins JAZs in JA signaling, as well as OsWRKY10, NAC, and ZFP transcription factors. More importantly, 104 of these, including ABCG36/OsPDR9, OsNAC3, as well as several orthologs in group metalloendoproteinase, plastocyanin-like domain containing proteins and pectin methylesterase inhibitor, may respond specifically to various Cd pressures, after subtracting the 60 general stress-responsive genes reported to be commonly upregulated following multiple stresses.

CONCLUSION

An integrative approach was implemented to identify DEGs and co-expression network modules in response to various Cd pressures, and 104 of the 164 annotatable universal Cd-responsive DEGs may specifically respond to various Cd pressures. These results provide insight into the universal molecular mechanisms beneath the Cd response in rice roots, and suggest many promising targets for improving the rice acclimation process against Cd toxicity.

Physiological and transcriptome response to cadmium in cosmos (Cosmos bipinnatus Cav.) seedlings

To date, several species of Asteraceae have been considered as Cd-accumulators. However, little information on the Cd tolerance and associated mechanisms of Asteraceae species Cosmos bipinnatus, is known. Presently, several physiological indexes and transcriptome profiling under Cd stress were investigated. C. bipinnatus exhibited strong Cd tolerance and recommended as a Cd-accumulator, although the biomasses were reduced by Cd. Meanwhile, Cd stresses reduced Zn and Ca uptake, but increased Fe uptake. Subcellular distribution indicated that the vacuole sequestration in root mainly detoxified Cd under lower Cd stress. Whilst, cell wall binding and vacuole sequestration in root co-detoxified Cd under high Cd exposure. Meanwhile, 66,407 unigenes were assembled and 41,674 (62.75%) unigenes were annotated in at least one database. 2,658 DEGs including 1,292 up-regulated unigenes and 1,366 down-regulated unigenes were identified under 40 μmol/L Cd stress. Among of these DEGs, ZIPs, HMAs, NRAMPs and ABC transporters might participate in Cd uptake, translocation and accumulation. Many DEGs participating in several processes such as cell wall biosynthesis, GSH metabolism, TCA cycle and antioxidant system probably play critical roles in cell wall binding, vacuole sequestration and detoxification. These results provided a novel insight into the physiological and transcriptome response to Cd in C. bipinnatus seedlings.

RNAseq analysis reveals pathways and candidate genes associated with salinity tolerance in a spaceflight-induced wheat mutant

Salinity stress has become an increasing threat to food security worldwide and elucidation of the mechanism for salinity tolerance is of great significance. Induced mutation, especially spaceflight mutagenesis, is one important method for crop breeding. In this study, we show that a spaceflight-induced wheat mutant, named salinity tolerance 1 (st1), is a salinity-tolerant line. We report the characteristics of transcriptomic sequence variation induced by spaceflight, and show that mutations in genes associated with sodium ion transport may directly contribute to salinity tolerance in st1. Furthermore, GO and KEGG enrichment analysis of differentially expressed genes (DEGs) between salinity-treated st1 and wild type suggested that the homeostasis of oxidation-reduction process is important for salt tolerance in st1. Through KEGG pathway analysis, "Butanoate metabolism" was identified as a new pathway for salinity responses. Additionally, key genes for salinity tolerance, such as genes encoding arginine decarboxylase, polyamine oxidase, hormones-related, were not only salt-induced in st1 but also showed higher expression in salt-treated st1 compared with salt-treated WT, indicating that these genes may play important roles in salinity tolerance in st1. This study presents valuable genetic resources for studies on transcriptome variation caused by induced mutation and the identification of salt tolerance genes in crops.

RNAseq analysis reveals pathways and candidate genes associated with salinity tolerance in a spaceflight-induced wheat mutant

Salinity stress has become an increasing threat to food security worldwide and elucidation of the mechanism for salinity tolerance is of great significance. Induced mutation, especially spaceflight mutagenesis, is one important method for crop breeding. In this study, we show that a spaceflight-induced wheat mutant, named salinity tolerance 1 (st1), is a salinity-tolerant line. We report the characteristics of transcriptomic sequence variation induced by spaceflight, and show that mutations in genes associated with sodium ion transport may directly contribute to salinity tolerance in st1. Furthermore, GO and KEGG enrichment analysis of differentially expressed genes (DEGs) between salinity-treated st1 and wild type suggested that the homeostasis of oxidation-reduction process is important for salt tolerance in st1. Through KEGG pathway analysis, "Butanoate metabolism" was identified as a new pathway for salinity responses. Additionally, key genes for salinity tolerance, such as genes encoding arginine decarboxylase, polyamine oxidase, hormones-related, were not only salt-induced in st1 but also showed higher expression in salt-treated st1 compared with salt-treated WT, indicating that these genes may play important roles in salinity tolerance in st1. This study presents valuable genetic resources for studies on transcriptome variation caused by induced mutation and the identification of salt tolerance genes in crops.

The role of ZmWRKY4 in regulating maize antioxidant defense under cadmium stress

WRKY transcription factors act as positive regulators in abiotic stress responses by activation of the cellular antioxidant systems. However, there are few reports on the response of WRKY genes to cadmium (Cd) stress. In this study, the role of maize ZmWRKY4 in regulating antioxidant enzymes in Cd stress was investigated. The results indicated that Cd induced up-regulation of the expression and the activities of ZmWRKY4 and superoxide dismutase (SOD) and ascorbate peroxidase (APX). Transient expression and RNA interference (RNAi) silencing of ZmWRKY4 in maize mesophyll protoplasts further revealed that ZmWRKY4 was required for the abscisic acid (ABA)-induced increase in expression and activity of SOD and APX. Overexpression of ZmWRKY4 in protoplasts upregulated the expression and the activities of antioxidant enzymes, whereas ABA induced increases in the expression and the activities of antioxidant enzymes were blocked by the RNAi silencing of ZmWRKY4. Bioinformatic analysis indicated that ZmSOD4 and ZmcAPX both harbored two W-boxes, binding motif for WRKY transcription factors, in their promoter region. Intriguingly, ZmWRKY4 belongs to group I WRKYs with two WRKY domains. Moreover, the synchronized expression patterns indicate that ZmWRKY4 might play a critical role in either regulating the ZmSOD4 and ZmcAPX expression or cooperating with them in response to stress and phytohormone.

Transcriptome Analysis of Cadmium-Treated Roots in Maize (Zea mays L.)

Cadmium (Cd) is a heavy metal and is highly toxic to all plant species. However, the underlying molecular mechanism controlling the effects of auxin on the Cd stress response in maize is largely unknown. In this study, the transcriptome produced by maize 'Zheng 58' root responses to Cd stress was sequenced using Illumina sequencing technology. In our study, six RNA-seq libraries yielded a total of 244 million clean short reads and 30.37 Gb of sequence data. A total of 6342 differentially expressed genes (DEGs) were grouped into 908 Gene Ontology (GO) categories and 198 Kyoto Encyclopedia of Genes and Genomes terms. GO term enrichment analysis indicated that various auxin signaling pathway-related GO terms were significantly enriched in DEGs. Comparison of the transcript abundances for auxin biosynthesis, transport, and downstream response genes revealed a universal expression response under Cd treatment. Furthermore, our data showed that free indole-3-acetic acid (IAA) levels were significantly reduced; but IAA oxidase activity was up-regulated after Cd treatment in maize roots. The analysis of Cd activity in maize roots under different Cd and auxin conditions confirmed that auxin affected Cd accumulation in maize seedlings. These results will improve our understanding of the complex molecular mechanisms underlying the response to Cd stress in maize roots.

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

Widely-spread cadmium (Cd) pollution in the soil threatens both crop production and human health. How plants deal with the excess Cd are largely unknown. To evaluate the molecular mechanism by which plants respond to Cd stress, rice seedlings were treated with two concentrations of Cd and subjected to deep RNA sequencing. Comprehensive RNA-Seq analysis of rice roots under two gradients of Cd treatment revealed 1169 Cd toxicity-responsive genes. These genes were involved in the reactive oxygen species scavenging system, stress response, cell wall formation, ion transport, and signal transduction. Nine out of 93 predicted long non-coding RNAs (lncRNAs) were detected as Cd-responsive lncRNAs due to their high correlation with the Cd stress response. In addition, we analyzed alternative splicing (AS) events under different Cd concentrations. Four hundred and seventy-six differential alternatively spliced genes with 542 aberrant splicing events were identified. GO analysis indicated that these genes were highly enriched in oxidation reduction and cellular response to chemical stimulus. Real-time qRT-PCR validation analysis strengthened the reliability of our RNA-Seq results. The results suggest that post-transcriptional AS regulation may also be involved in plant responses to high Cd stress.