Transcriptome analysis of rice (Oryza sativa L.) shoots responsive to cadmium stress

Identification of candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq

BACKGROUND

The cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement.

RESULTS

Transcriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to 'oxidation-reduction process' and 'response to stress' differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to 'response to stress') decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to 'oxidation-reduction process', 'response to stress' and 'protein phosphorylation' interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type.

CONCLUSIONS

In general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.

Transcription factors involved in plant responses to cadmium-induced oxidative stress

Cadmium (Cd) is a heavy metal highly toxic to living organisms. Cd pollution of soils has become a serious problem worldwide, posing a severe threat to crop production and human health. When plants are poisoned by Cd, their growth and development are inhibited, chloroplasts are severely damaged, and respiration and photosynthesis are negatively affected. Therefore, elucidating the molecular mechanisms that underlie Cd tolerance in plants is important. Transcription factors can bind to specific plant cis-acting genes. Transcription factors are frequently reported to be involved in various signaling pathways involved in plant growth and development. Their role in the resistance to environmental stress factors, particularly Cd, should not be underestimated. The roles of several transcription factor families in the regulation of plant resistance to Cd stress have been widely demonstrated. In this review, we summarize the mechanisms of five major transcription factor families-WRKY, ERF, MYB, bHLH, and bZIP-in plant resistance to Cd stress to provide useful information for using molecular techniques to solve Cd pollution problems in the future.

Mapping of Candidate Genes for Nitrogen Uptake and Utilization in Japonica Rice at Seedling Stage

Nitrogen (N) is one of the essential nutrients for the growth and development of crops. The adequate application of N not only increases the yield of crops but also improves the quality of agricultural products, but the excessive application of N can cause many adverse effects on ecology and the environment. In this study, genome-wide association analysis (GWAS) was performed under low- and high-N conditions based on 788,396 SNPs and phenotypic traits relevant to N uptake and utilization (N content and N accumulation). A total of 75 QTLs were obtained using GWAS, which contained 811 genes. Of 811 genes, 281 genes showed different haplotypes, and 40 genes had significant phenotypic differences among different haplotypes. Of these 40 genes, 5 differentially expressed genes (Os01g0159250, Os02g0618200, Os02g0618400, Os02g0630300, and Os06g0619000) were finally identified as the more valuable candidate genes based on the transcriptome data sequenced from Longjing31 (low-N-tolerant variety) and Songjing 10 (low-N-sensitive variety) under low- and high-N treatments. These new findings enrich the genetic resources for N uptake and utilization in rice, as well as lay a theoretical foundation for improving the efficiency of N uptake and utilization in rice.

More effective than direct contact: Nano hydroxyapatite pre-treatment regulates the growth and Cd uptake of rice (Oryza sativa L.) seedlings

Cd contamination in rice urgently needs to be addressed. Nano hydroxyapatite (n-HAP) is an eco-friendly material with excellent Cd fixation ability. However, due to its own high reactivity, innovative application of n-HAP in the treatment of Cd contamination in rice is needed. In this study, we proposed a new application, namely n-HAP pre-treatment, which can effectively reduce Cd accumulation in rice and alleviate Cd stress. The results showed that 80 mg/L n-HAP pre-treatment significantly reduced Cd content in rice shoot by 35.1%. Biochemical and combined transcriptomic-proteomic analysis revealed the possible molecular mechanisms by which n-HAP pre-treatment promoted rice growth and reduced Cd accumulation. (1) n-HAP pre-treatment regulated gibberellin and jasmonic acid synthesis-related pathways, increased gibberellin content and decreased jasmonic acid content in rice root, which promoted rice growth; (2) n-HAP pre-treatment up-regulated gene CATA1 expression and down-regulated gene OsGpx1 expression, which increased rice CAT activity and GSH content; (3) n-HAP pre-treatment up-regulated gene OsZIP1 expression and down-regulated gene OsNramp1 expression, which reduced Cd uptake, increased Cd efflux from rice root cells.

Transcriptome Analysis Reveals the Stress Tolerance Mechanisms of Cadmium in Zoysia japonica

Cadmium (Cd) is a severe heavy metal pollutant globally. Zoysia japonica is an important perennial warm-season turf grass that potentially plays a role in phytoremediation in Cd-polluted soil areas; however, the molecular mechanisms underlying its Cd stress response are unknown. To further investigate the early gene response pattern in Z. japonica under Cd stress, plant leaves were harvested 0, 6, 12, and 24 h after Cd stress (400 μM CdCl2) treatment and used for a time-course RNA-sequencing analysis. Twelve cDNA libraries were constructed and sequenced, and high-quality data were obtained, whose mapped rates were all higher than 94%, and more than 601 million bp of sequence were generated. A total of 5321, 6526, and 4016 differentially expressed genes were identified 6, 12, and 24 h after Cd stress treatment, respectively. A total of 1660 genes were differentially expressed at the three time points, and their gene expression profiles over time were elucidated. Based on the analysis of these genes, the important mechanisms for the Cd stress response in Z. japonica were identified. Specific genes participating in glutathione metabolism, plant hormone signal and transduction, members of protein processing in the endoplasmic reticulum, transporter proteins, transcription factors, and carbohydrate metabolism pathways were further analyzed in detail. These genes may contribute to the improvement of Cd tolerance in Z. japonica. In addition, some candidate genes were highlighted for future studies on Cd stress resistance in Z. japonica and other plants. Our results illustrate the early gene expression response of Z. japonica leaves to Cd and provide some new understanding of the molecular mechanisms of Cd stress in Zosia and Gramineae species.

The Effect of Cadmium on Plants in Terms of the Response of Gene Expression Level and Activity

Cadmium (Cd) is a heavy metal that can cause damage to living organisms at different levels. Even at low concentrations, Cd can be toxic to plants, causing harm at multiple levels. As they are unable to move away from areas contaminated by Cd, plants have developed various defence mechanisms to protect themselves. Hyperaccumulators, which can accumulate and detoxify heavy metals more efficiently, are highly valued by scientists studying plant accumulation and detoxification mechanisms, as they provide a promising source of genes for developing plants suitable for phytoremediation techniques. So far, several genes have been identified as being upregulated when plants are exposed to Cd. These genes include genes encoding transcription factors such as iron-regulated transporter-like protein (ZIP), natural resistance associated macrophage protein (NRAMP) gene family, genes encoding phytochelatin synthases (PCs), superoxide dismutase (SOD) genes, heavy metal ATPase (HMA), cation diffusion facilitator gene family (CDF), Cd resistance gene family (PCR), ATP-binding cassette transporter gene family (ABC), the precursor 1-aminocyclopropane-1-carboxylic acid synthase (ACS) and precursor 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) multigene family are also influenced. Thanks to advances in omics sciences and transcriptome analysis, we are gaining more insights into the genes involved in Cd stress response. Recent studies have also shown that Cd can affect the expression of genes related to antioxidant enzymes, hormonal pathways, and energy metabolism.

The endo-beta mannase MAN7 contributes to cadmium tolerance by modulating root cell wall binding capacity in Arabidopsis thaliana

The heavy metal cadmium (Cd) is detrimental to crop growth and threatens human health through the food chain. To cope with Cd toxicity, plants employ multiple strategies to decrease Cd uptake and its root-to-shoot translocation. However, genes that participate in the Cd-induced transcriptional regulatory network, including those encoding transcription factors, remain largely unidentified. In this study, we demonstrate that ENDO-BETA-MANNASE 7 (MAN7) is necessary for the response of Arabidopsis thaliana to toxic Cd levels. We show that MAN7 is responsible for mannase activity and modulates mannose content in the cell wall, which plays a role in Cd compartmentalization in the cell wall under Cd toxicity conditions. Additionally, the repression of root growth by Cd was partially reversed via exogenous application of mannose, suggesting that MAN7-mediated cell wall Cd redistribution depends on the mannose pathway. Notably, we identified a basic leucine zipper (bZIP) transcription factor, bZIP44, that acts upstream of MAN7 in response to Cd toxicity. Transient dual-luciferase assays indicated that bZIP44 directly binds to the MAN7 promoter region and activates its transcription. Loss of bZIP44 function was associated with greater sensitivity to Cd treatment and higher accumulation of the heavy metal in roots and shoots. Moreover, MAN7 overexpression relieved the inhibition of root elongation seen in the bzip44 mutant under Cd toxicity conditions. This study thus reveals a pathway showing that MAN7-associated Cd tolerance in Arabidopsis is controlled by bZIP44 upon Cd exposure.

Transcriptome Analysis Reveals Differentially Expressed Genes Involved in Cadmium and Arsenic Accumulation in Tea Plant (Camellia sinensis)

Tea (Camellia sinensis) is the second most consumed drink in the world. Rapid industrialization has caused various impacts on nature and increased pollution by heavy metals. However, the molecular mechanisms of cadmium (Cd) and arsenic (As) tolerance and accumulation in tea plants are poorly understood. The present study focused on the effects of heavy metals Cd and As on tea plants. Transcriptomic regulation of tea roots after Cd and As exposure was analyzed to explore the candidate genes involved in Cd and As tolerance and accumulation. In total, 2087, 1029, 1707, and 366 differentially expressed genes (DEGs) were obtained in Cd1 (with Cd treatment for 10 days) vs. CK (without Cd treatment), Cd2 (with Cd treatment for 15 days) vs. CK, As1 (with As treatment for 10 days) vs. CK (without Cd treatment), and As2 (with As treatment for 15 days) vs. CK, respectively. Analysis of DEGs showed that a total of 45 DEGs with the same expression patterns were identified in four pairwise comparison groups. One ERF transcription factor (CSS0000647) and six structural genes (CSS0033791, CSS0050491, CSS0001107, CSS0019367, CSS0006162, and CSS0035212) were only increased at 15 d of Cd and As treatments. Using weighted gene co-expression network analysis (WGCNA) revealed that the transcription factor (CSS0000647) was positively correlated with five structural genes (CSS0001107, CSS0019367, CSS0006162, CSS0033791, and CSS0035212). Moreover, one gene (CSS0004428) was significantly upregulated in both Cd and As treatments, suggesting that these genes might play important roles in enhancing the tolerance to Cd and As stresses. These results provide candidate genes to enhance multi-metal tolerance through the genetic engineering technology.

Temporal Comparative Transcriptome Analysis on Wheat Response to Acute Cd Toxicity at the Seedling Stage

Cadmium (Cd) is a non-essential and toxic metal that accumulates in plant's tissues and diminishes plant growth and productivity. In the present study, differential root transcriptomic analysis was carried out to identify Cd stress-responsive gene networks and functional annotation under Cd stress in wheat seedlings. For this purpose, the Yannong 0428 wheat cultivar was incubated with 40 µm/L of CdCl₂·2.5H₂O for 6 h at three different seedling growth days. After the quality screening, using the Illumina Hiseq 2000 platform, more than 2482 million clean reads were retrieved. Following this, 84.8% to 89.3% of the clean reads at three time points under normal conditions and 86.5% to 89.1% of the reads from the Cd stress condition were mapped onto the wheat reference genome. In contrast, at three separate seedling growth days, the data analysis revealed a total of 6221 differentially expressed genes (DEGs), including 1543 (24.8%) up-regulated genes and 4678 (75.8%) down-regulated genes. In total, 120 DEGs were co-expressed throughout all the growth days, whereas 1096, 1088, and 2265 DEGs were found to be selectively up-/down-regulated at 7d, 14d, and 30d, respectively. However, the clustering of DEGs, through utilizing the Kyoto Encyclopedia of Genes and Genomes (KEGG), revealed that the DEGs in the metabolic category were frequently annotated for phenylpropanoid biosynthesis. In comparison, a considerable number of DEGs were linked to protein processing in the endoplasmic reticulum under the process of genetic information processing. Similarly, in categories in organismal systems and cellular processes, DEGs were found in plant hormone signal transduction pathways, and DEGs were identified in the plant-pathogen interaction pathway, respectively. However, DEGs in "endocytosis pathways" were enriched in environmental information processing. In addition, in-depth annotations of roughly specific heavy metal stress-response genes and pathways were also mined, and the expression patterns of eight DEGs were studied using quantitative real-time PCR. The results were congruent with the findings of RNA sequencing regarding transcript abundance in the studied wheat cultivar.

Transcriptome alterations of radish shoots exposed to cadmium can be interpreted in the context of leaf senescence

Till now few transcriptome studies have described shoot responses of heavy metal (HM)-sensitive plants to excess Cd and still a unifying model of Cd action is lacking. Using RNA-seq technique, the transcriptome responses of radish (Raphanus sativus L.) leaves to Cd stress were investigated in plants raised hydroponically under control and 5.0 mg L^-1 Cd. The element was mainly accumulated in roots and led to declined biomass and photosynthetic pigments, increased H₂O₂ and lipid peroxidation, and the accumulation of sugars, protein thiols, and phytochelatins. Out of 524 differentially expressed genes (DEGs), 244 and 280 upregulated and downregulated ones were assigned to 82 and 115 GO terms, respectively. The upregulated DEGs were involved in osmotic regulation, protein metabolism, chelators, and carbohydrate metabolisms, whereas downregulated DEGs were related to photosynthesis, response to oxidative stress, glucosinolate, and secondary metabolite biosynthesis. Our transcriptome data suggest that Cd triggers ROS production and photosynthesis decline associated with increased proteolysis through ubiquitin-proteasome system (UPS)- and chloroplast-proteases and in this way brings about re-mobilization of N and C stores into amino acids and sugars. Meanwhile, declined glucosinolate metabolism in favor of chelator synthesis and upregulation of dehydrins as inferred from transcriptome analysis confers shoots some tolerance to the HM-derived ionic/osmotic imbalances. Thus, the induction of leaf senescence might be a major long-term response of HM-sensitive plants to Cd toxicity.

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.

Gene Expression in the Developing Seed of Wild and Domesticated Rice

The composition and nutritional properties of rice are the product of the expression of genes in the developing seed. RNA-Seq was used to investigate the level of gene expression at different stages of seed development in domesticated rice (Oryza sativa ssp. japonica var. Nipponbare) and two Australian wild taxa from the primary gene pool of rice (Oryza meridionalis and Oryza rufipogon type taxa). Transcriptome profiling of all coding sequences in the genome revealed that genes were significantly differentially expressed at different stages of seed development in both wild and domesticated rice. Differentially expressed genes were associated with metabolism, transcriptional regulation, nucleic acid processing, and signal transduction with the highest number of being linked to protein synthesis and starch/sucrose metabolism. The level of gene expression associated with domestication traits, starch and sucrose metabolism, and seed storage proteins were highest at the early stage (5 days post anthesis (DPA)) to the middle stage (15 DPA) and declined late in seed development in both wild and domesticated rice. However, in contrast, black hull colour (Bh4) gene was significantly expressed throughout seed development. A substantial number of novel transcripts (38) corresponding to domestication genes, starch and sucrose metabolism, and seed storage proteins were identified. The patterns of gene expression revealed in this study define the timing of metabolic processes associated with seed development and may be used to explain differences in rice grain quality and nutritional value.

Transcriptional memory of gene expression across generations participates in transgenerational plasticity of field pennycress in response to cadmium stress

Transgenerational plasticity (TGP) occurs when maternal environments influence the expression of traits in offspring, and in some cases may increase fitness of offspring and have evolutionary significance. However, little is known about the extent of maternal environment influence on gene expression of offspring, and its relationship with trait variations across generations. In this study, we examined TGP in the traits and gene expression of field pennycress (Thlaspi arvense) in response to cadmium (Cd) stress. In the first generation, along with the increase of soil Cd concentration, the total biomass, individual height, and number of seeds significantly decreased, whereas time to flowering, superoxide dismutase (SOD) activity, and content of reduced glutathione significantly increased. Among these traits, only SOD activity showed a significant effect of TGP; the offspring of Cd-treated individuals maintained high SOD activity in the absence of Cd stress. According to the results of RNA sequencing and bioinformatic analysis, 10,028 transcripts were identified as Cd-responsive genes. Among them, only 401 were identified as transcriptional memory genes (TMGs) that maintained the same expression pattern under normal conditions in the second generation as in Cd-treated parents in the first generation. These genes mainly participated in Cd tolerance-related processes such as response to oxidative stress, cell wall biogenesis, and the abscisic acid signaling pathways. The results of weighted correlation network analysis showed that modules correlated with SOD activity recruited more TMGs than modules correlated with other traits. The SOD-coding gene CSD2 was found in one of the modules correlated with SOD activity. Furthermore, several TMGs co-expressed with CSD2 were hub genes that were highly connected to other nodes and critical to the network's topology; therefore, recruitment of TMGs in offspring was potentially related to TGP. These findings indicated that, across generations, transcriptional memory of gene expression played an important role in TGP. Moreover, these results provided new insights into the trait evolution processes mediated by phenotypic plasticity.

Transcriptome Analysis and SNP Identification Reveal That Heterologous Overexpression of Two Uncharacterized Genes Enhances the Tolerance of Magnaporthe oryzae to Manganese Toxicity

Manganese is a crucial trace element that constitutes the cofactors of many enzymes. However, excessive Mn2+ can be toxic for both prokaryotes and eukaryotes. The mechanism of fungal genetics and metabolism in response to Mn2+ stress remains understudied, warranting further studies. Magnaporthe oryzae is well-established as the most destructive pathogen of rice. A field strain, YN2046, more sensitive to Mn2+ toxicity than other strains, was obtained from a previous study. Herein, we explored the genetic mechanisms of Mn2+ sensitivity in YN2046 through comparative transcriptomic analyses. We found that many genes previously reported to participate in Mn2+ stress were not regulated in YN2046. These non-responsive genes might cause Mn2+ sensitivity in YN2046. Weight gene correlation network analysis (WGCNA) was performed to characterize the expression profile in YN2046. Some overexpressed genes were only found in the Mn2+ tolerant isolate YN125. Among these, many single nucleotide polymorphism (SNP) were identified between YN125 and YN2046, which might disrupt the expression levels of Mn responsive genes. We cloned two uncharacterized genes, MGG_13347 and MGG_16609, from YN125 and transformed them to YN2046 with a strong promoter. Our results showed that the heterologous overexpression of two genes in YN2046 restored its sensitivity. Transcriptomic and biochemical analyses were performed to understand Mn tolerance mechanisms mediated by the two heterologous overexpressed genes. Our results showed that heterologous overexpression of these two genes activated downstream gene expression and metabolite production to restore M. oryzae sensitivity to Mn, implying that SNPs in responsive genes account for different phenotypes of the two strains under Mn stress. IMPORTANCE Heavy metals are used for fungicides as they target phytopathogen in multiple ways. Magnaporthe oryzae is the most destructive rice pathogen and is threatening global rice production. In the eukaryotes, the regulation mechanisms of Mn homeostasis often focus on the posttranslation, there were a few results about regulation at transcript level. The comparative transcriptome analysis showed that fewer genes were regulated in the Mn-sensitive strain. WGCNA and SNP analyses found that mutations in promoter and coding sequence regions might disrupt the expression of genes involved in Mn detoxification in the sensitive strain. We transferred two unannotated genes that were cloned from the Mn-tolerant strain into a sensitive strain with strong promoters, and the transformants exhibited an enhanced tolerance to Mn2+ toxicity. Transcriptome and biochemistry results indicated that heterologous overexpression of the two genes enhanced the tolerance to Mn toxicity by reactivation of downstream genes in M. oryzae.

Stress signaling convergence and nutrient crosstalk determine zinc-mediated amelioration against cadmium toxicity in rice

Consumption of rice (Oryza sativa L.) is one of the major pathways for heavy metal bioaccumulation in humans over time. Understanding the molecular responses of rice to heavy metal contamination in agriculture is useful for eco-toxicological assessment of cadmium (Cd) and its interaction with zinc (Zn). In certain crops, the impacts of Cd stress or Zn nutrition on the biophysical chemistry and gene expression have been widely investigated, but their molecular interactions at transcriptomic level, particularly in rice roots, are still elusive. Here, hydroponic investigations were carried out with two rice genotypes (Yinni-801 and Heizhan-43), varying in Cd contents in plant tissues to determine their transcriptomic responses upon Cd15 (15 µM) and Cd15+Zn50 (50 µM) treatments. High throughput RNA-sequencing analysis confirmed that 496 and 2407 DEGs were significantly affected by Cd15 and Cd15+Zn50, respectively, among which 1016 DEGs were commonly induced in both genotypes. Multitude of DEGs fell under the category of protein kinases, such as calmodulin (CaM) and calcineurin B-like protein-interacting protein kinases (CBL), indicating a dynamic shift in hormonal signal transduction and Ca2+ involvement with the onset of treatments. Both genotypes expressed a mutual regulation of transcription factors (TFs) such as WRKY, MYB, NAM, AP2, bHLH and ZFP families under both treatments, whereas genes econding ABC transporters (ABCs), high affinity K+ transporters (HAKs) and Glutathione-S-transferases (GSTs), were highly up-regulated under Cd15+Zn50 in both genotypes. Zinc addition triggered more signaling cascades and detoxification related genes in regulation of immunity along with the suppression of Cd-induced DEGs and restriction of Cd uptake. Conclusively, the effective integration of breeding techniques with candidate genes identified in this study as well as economically and technologically viable methods, such as Zn nutrient management, could pave the way for selecting cultivars with promising agronomic qualities and reduced Cd for sustainable rice production.

Cross-Kingdom Comparative Transcriptomics Reveals Conserved Genetic Modules in Response to Cadmium Stress

It is known that organisms have developed various mechanisms to cope with cadmium (Cd) stress, while we still lack a system-level understanding of the functional isomorphy among them. In the present study, a cross-kingdom comparison was conducted among Escherichia coli, Saccharomyces cerevisiae, and Chlamydomonas reinhardtii, through toxicological tests, comparative transcriptomics, as well as conventional functional genomics. An equivalent level of Cd stress was determined via inhibition tests. Through transcriptome comparison, the three organisms exhibited differential gene expression under the same Cd stress relative to the corresponding no-treatment control. Results from functional enrichment analysis of differentially expressed genes (DEGs) showed that four metabolic pathways responsible for combating Cd stress were commonly regulated in the three organisms, including antioxidant reactions, sulfur metabolism, cell wall remodeling, and metal transport. In vivo expression patterns of 43 DEGs from the four pathways were further examined using quantitative PCR and resulted in a relatively comparable dynamic of gene expression patterns with transcriptome sequencing (RNA-seq). Cross-kingdom comparison of typical Cd stress-responding proteins resulted in the detection of 12 groups of homologous proteins in the three species. A class of potential metal transporters were subjected to cross-transformation to test their functional complementation. An ABC transporter gene in E. coli, possibly homologous to the yeast ycf1, was heterologously expressed in S. cerevisiae, resulting in enhanced Cd tolerance. Overall, our findings indicated that conserved genetic modules against Cd toxicity were commonly regulated among distantly related microbial species, which will be helpful for utilizing them in modifying microbial traits for bioremediation.

IMPORTANCE Research is establishing a systems biology view of biological response to Cd stress. It is meaningful to explore whether there is regulatory isomorphy among distantly related organisms. A transcriptomic comparison was done among model microbes, leading to the identification of a conserved cellular model pinpointing the generic strategies utilized by microbes for combating Cd stress. A novel E. coli transporter gene substantially increased yeast’s Cd tolerance. Knowledge on systems understanding of the cellular response to metals provides the basis for developing bioengineering remediation technology.

Short-term effects of cadmium on leaf growth and nutrient transport in rice plants

Consumption of rice grains contaminated with high concentrations of cadmium (Cd) can cause serious long-term health problems. Moreover, even low Cd concentrations present in the soil can result in the abatement of plant performance, leading to lower grain yield. Studies examining the molecular basis of plant defense against Cd-induced oxidative stress could pave the way in creating superior rice varieties that display an optimal antioxidative defense system to cope with Cd toxicity. In this study, we showed that after one day of Cd exposure, hydroponically grown rice plants exhibited adverse shoot biomass and leaf growth effects. Cadmium accumulates especially in the roots and the leaf meristematic region, leading to a disturbance of manganese homeostasis in both the roots and leaves. The leaf growth zone showed an increased amount of lipid peroxidation indicating that Cd exposure disturbed the oxidative balance. We propose that an increased expression of genes related to the glutathione metabolism such as glutathione synthetase 2, glutathione reductase and phytochelatin synthase 2, rather than genes encoding for antioxidant enzymes, is important in combating early Cd toxicity within the leaves of rice plants. Furthermore, the upregulation of two RESPIRATORY BURST OXIDASE HOMOLOG genes together with a Cd concentration-dependent increase of abscisic acid might cause stomatal closure or cell wall modification, potentially leading to the observed leaf growth reduction. Whereas abscisic acid was also elevated at long term exposure, a decrease of the growth hormone auxin might further contribute to growth inhibition and concomitantly, an increase in salicylic acid might stimulate the activity of antioxidative enzymes after a longer period of Cd exposure. In conclusion, a clear interplay between phytohormones and the oxidative challenge affect plant growth and acclimation during exposure to Cd stress.

Toward Integrated Multi-Omics Intervention: Rice Trait Improvement and Stress Management

Rice (Oryza sativa) is an imperative staple crop for nearly half of the world's population. Challenging environmental conditions encompassing abiotic and biotic stresses negatively impact the quality and yield of rice. To assure food supply for the unprecedented ever-growing world population, the improvement of rice as a crop is of utmost importance. In this era, "omics" techniques have been comprehensively utilized to decipher the regulatory mechanisms and cellular intricacies in rice. Advancements in omics technologies have provided a strong platform for the reliable exploration of genetic resources involved in rice trait development. Omics disciplines like genomics, transcriptomics, proteomics, and metabolomics have significantly contributed toward the achievement of desired improvements in rice under optimal and stressful environments. The present review recapitulates the basic and applied multi-omics technologies in providing new orchestration toward the improvement of rice desirable traits. The article also provides a catalog of current scenario of omics applications in comprehending this imperative crop in relation to yield enhancement and various environmental stresses. Further, the appropriate databases in the field of data science to analyze big data, and retrieve relevant information vis-à-vis rice trait improvement and stress management are described.

Comparative root transcriptome analysis of two soybean cultivars with different cadmium sensitivities reveals the underlying tolerance mechanisms

Soybean can provide rich protein and fat and has great economic value worldwide. Cadmium (Cd) is a toxic heavy metal to organisms. It can accumulate in plants and be transmitted to the human body via the food chain. Cd is a serious threat to soybean development, particularly root growth. Some soybean cultivars present tolerant symptoms under Cd stress; however, the potential mechanisms are not fully understood. Here, we optimized RNA-seq to identify the differentially expressed genes (DEGs) in Cd-sensitive (KUAI) and Cd-tolerant (KAIYU) soybean roots and compared the DEGs between KAIYU and KUAI. A total of 1506 and 1870 DEGs were identified in the roots of KUAI and KAIYU, respectively. Through Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and gene function analyses, we found that genes related to antioxidants and sequestration were responsible for Cd tolerance in KAIYU. In addition, overexpression of Glyma11g02661, which encodes a heavy metal-transporting ATPase, significantly improved Cd tolerance in transgenic hairy roots. These results provide a preliminary understanding of the tolerance mechanisms in response to Cd stress in soybean root development and are of great importance in developing Cd-resistant soybean cultivars by using the identified DEGs through genetic modification.

An update on redox signals in plant responses to biotic and abiotic stress crosstalk: insights from cadmium and fungal pathogen interactions

Complex signalling pathways are involved in plant protection against single and combined stresses. Plants are able to coordinate genome-wide transcriptional reprogramming and display a unique programme of transcriptional responses to a combination of stresses that differs from the response to single stresses. However, a significant overlap between pathways and some defence genes in the form of shared and general stress-responsive genes appears to be commonly involved in responses to multiple biotic and abiotic stresses. Reactive oxygen and nitrogen species, as well as redox signals, are key molecules involved at the crossroads of the perception of different stress factors and the regulation of both specific and general plant responses to biotic and abiotic stresses. In this review, we focus on crosstalk between plant responses to biotic and abiotic stresses, in addition to possible plant protection against pathogens caused by previous abiotic stress. Bioinformatic analyses of transcriptome data from cadmium- and fungal pathogen-treated plants focusing on redox gene ontology categories were carried out to gain a better understanding of common plant responses to abiotic and biotic stresses. The role of reactive oxygen and nitrogen species in the complex network involved in plant responses to changes in their environment is also discussed.

Metabolome and transcriptome analysis reveals the molecular profiles underlying the ginseng response to rusty root symptoms

BACKGROUND

Ginseng rusty root symptoms (GRS) is one of the primary diseases of ginseng. This disease leads to a severe decline in the quality of ginseng. It has been shown that the occurrence of GRS is associated with soil environmental degradation, which may involve changes in soil microbiology and physicochemical properties.

RESULTS

In this study, GRS and healthy ginseng (HG) samples were used as experimental materials for comparative analysis of transcriptome and metabolome. Compared with those in HG samples, 949 metabolites and 9451 genes were significantly changed at the metabolic and transcriptional levels in diseased samples. The diseased tissues' metabolic patterns changed, and the accumulation of various organic acids, alkaloids, alcohols and phenols in diseased tissues increased significantly. There were significant differences in the expression of genes involved in plant hormone signal transduction, phenylpropanoid biosynthesis, the peroxidase pathway, and the plant-pathogen interaction pathway.

CONCLUSION

The current study involved a comparative metabolome and transcriptome analysis of GRS and HG samples. Based on the findings at the transcriptional and metabolic levels, a mechanism model of the ginseng response to GRS was established. Our results provide new insights into ginseng's response to GRS, which will reveal the potential molecular mechanisms of this disease in ginseng.

Transcriptome analysis and association mapping reveal the genetic regulatory network response to cadmium stress in Populus tomentosa

Long non-coding RNAs (lncRNAs) play essential roles in plant abiotic stress responses, but the response of lncRNA-mediated genetic networks to cadmium (Cd) treatment remain elusive in trees, the promising candidates for phytoremediation of Cd contamination. We identified 172 Cd-responsive lncRNAs and 295 differentially expressed target genes in the leaves of Cd-treated Populus tomentosa. Functional annotation revealed that these lncRNAs were involved in various processes, including photosynthesis, hormone regulation, and phenylalanine metabolism. Association studies identified 78 significant associations, representing 14 Cd-responsive lncRNAs and 28 target genes for photosynthetic and leaf physiological traits. Epistasis uncovered 83 pairwise interactions among these traits, revealing Cd-responsive lncRNA-mediated genetic networks for photosynthesis and leaf physiology in P. tomentosa. We focused on the roles of two Cd-responsive lncRNA-gene pairs, MSTRG.22608.1-PtoMYB73 and MSTRG.5634.1-PtoMYB27, in Cd tolerance of Populus, and detected insertions/deletions within lncRNAs as polymorphisms driving target gene expression. Genotype analysis of lncRNAs and heterologous overexpression of PtoMYB73 and PtoMYB27 in Arabidopsis indicated their effects on enhancing Cd tolerance, photosynthetic rate, and leaf growth, and the potential interaction mechanisms of PtoMYB73 with abiotic stresses. Our study identifies the genetic basis for the response of Populus to Cd treatment, facilitating genetic improvement of Cd tolerance in trees.

Auxins, the hidden player in chloroplast development

Throughout decades of plant research, the plant hormones known as auxins have been found to be of vital importance in most plant development processes. Indole-3-acetic acid (IAA) represents the most common auxin in plants and can be synthesized from its tryptophan precursor, which is synthesized in the chloroplast. The chloroplast constitutes an organelle of great relevance to plants since the photosynthesis process by which plants get most of their energy is carried out there. The role of auxins in photosynthesis has been studied for at least 50 years, and in this time, it has been shown that auxins have an effect on several of the essential components and structure of the chloroplast. In recent decades, a high number of genes have been reported to be expressed in the chloroplast and some of their mutants have been shown to alter different auxin-mediated pathways. Genes in signaling pathways such as IAA/AUX, ARF, GH.3, SAUR and TIR, biosynthesis-related genes such as YUCCA and transport-related genes such as PIN have been identified among the most regulated genes in mutants related to alterations in the chloroplast. This review aims to provide a complete and updated summary of the relationship between auxins and several processes that involve the chloroplast, including chloroplast development, plant albinism, redox regulation and pigment synthesis.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

Genome-wide analysis of PHD finger gene family and identification of potential miRNA and their PHD finger gene specific targets in Oryza sativa indica

Rice (Oryza sativa L.) is one of the most important cereal crops for one third of the world population. However, the grain quality as well as yield of rice is severely affected by various abiotic stresses. Environmental stresses affect the expression of various microRNAs (miRNAs) which in turn negatively regulate gene expression at the post-transcriptional level either by degrading the target mRNA genes or suppressing translation in plants. Plant homeo-domain (PHD) finger proteins are known to be involved in the plant response to salinity stress. In the present study, we identified 44 putative OsPHD finger genes in Oryza sativa Indica, using Ensembl Plants Database. Using computational approach, potential miRNAs that target OsPHD finger genes were identified. Out of the 44 OsPHD finger genes only three OsPHD finger genes i.e., OsPHD2, OsPHD35 and OsPHD11, were found to be targeted by five newly identified putative miRNAs i.e., ath-miRf10010-akr, ath-miRf10110-akr, osa-miR1857–3p, osa-miRf10863-akr, and osa-miRf11806-akr. This is the first report of these five identified miRNAs on targeting PHD finger in Oryza sativa Indica. Further, expression analysis of 44 PHD finger genes under salinity was also performed using quantitative Real-Time PCR. The expression profile of 8 genes were found to be differentially regulated, among them two genes were significantly up regulated i.e., OsPHD6 and OsPHD12. In silico protein-protein interaction analysis using STRING database showed interaction of the OsPHD finger proteins with other protein partners that are directly or indirectly involved in development and abiotic stress tolerance.

•

Total of 44 Plant homeo-domain (PHD) finger proteins were identified & classified into 10 groups in Oryza sativa Indica.

•

This is the first report showing 5 newly identified putative miRNAs targeting three OsPHD genes i.e., OsPHD2, 11 and 35.

•

Expression analysis of PHD finger genes showed up-regulation of the 2 genes OsPHD 6 & 12 under salinity stress treatment.

•

Protein-protein network of OsPHDs showed protein partners that are involved in plant growth and abiotic stress tolerance.

In silico identification and validation of miRNA and their DIR specific targets in Oryza sativa Indica under abiotic stress

Several biotic (bacterial and viral pathogenesis) and abiotic stress factors like salt, drought, cold, and extreme temperatures significantly reduce crop productivity and grain quality throughout the world. MicroRNAs (miRNAs) are small (~22 nucleotides) non-coding endogenous RNA molecules which negatively regulate gene expression at the post-transcriptional level either by degrading the target protein-coding mRNA genes or suppressing translation in plants. Dirigent (DIR) gene protein plays a crucial role as they are involved to dictate the stereochemistry of a compound synthesized by other enzymes as well as in lignifications against biotic and abiotic stress. In plants, several miRNAs, as well as their targets, are known to regulate stress response but systematic identification of the same is limited. The present work has been designed for in silico identification of miRNAs against a total of sixty-one DIR genes in Oryza sativa Indica followed by target prediction of identified miRNAs through the computational approach and thereafter validation of potential miRNAs in rice genotypes. We systematically identified 3 miRNA and their respective DIR specific target gene in Oryza sativa Indica. The expression of these three miRNAs and their respective DIR specific targets were validated in rice seedlings subjected to five different abiotic stress conditions (heavy metal, high temperature, low temperature, salinity and drought) by quantitative Real-Time PCR (qRT-PCR). Expression analysis indicated that miRNA under stress conditions regulates the gene expression of the DIR gene in rice. To the best of our knowledge this is this is the first report in any organism showing the expression of ath-miRf10317-akr, and osamiRf10761-akr miRNAs in response to various abiotic stresses.

•

Total 61 DIR proteins were identified & classified into 6 groups based on phylogeny analysis in Oryza sativa Indica.

•

Three miRNAs ath-miRf10317-akr, cre-miR910 and osa-miRf10761-akr were identified via computational approach.

•

These 3 miRNAs in response to abiotic stresses showed inverse expression pattern in the respective target genes.

•

This is the first report on expression of ath-miRf10317-akr, and osa-miRf10761-akr miRNAs in response to abiotic stresses.

Regulation of Carbohydrate Metabolism by Trehalose-6-Phosphate Synthase 3 in the Brown Planthopper, Nilaparvata lugens

Nilaparvata lugens (Stål) (Hemiptera: Delphacidae) is one of the pests that harm rice. In this paper, a new trehalose-6-phosphate synthase gene, TPS3, was identified by transcriptome sequencing and gene cloning. To explore its role in the energy metabolism of N. lugens we examined the carbohydrate contents at different stages of development, the tissue expression of TPS, and some physiological and biochemical indicators by injecting dsTPS3 and dsTPSs (a proportional mixture of dsTPS1, dsTPS2, and dsTPS3). The glucose content at the fifth instar was significantly higher than that in the fourth instar and the adult stages. The trehalose and glycogen contents before molting were higher than those after molting. TPS1, TPS2, and TPS3 were expressed in the head, leg, wing bud, and cuticle, with the highest expression in the wing bud. In addition, compared with the control group, the glucose content increased significantly at 48 h after RNA interference, and the trehalose content decreased significantly after 72 h. qRT-PCR showed that the expression level of UGPase decreased significantly at 48 h after injection, whereas GS expression increased significantly at 48 h after injecting dsTPS3. After dsTPS injection, the expression levels of PPGM2, UGPase, GP, and GS increased significantly at 72 h. After interfering with the expression of TPS3 gene alone, UGPase expression decreased significantly at 48 h, and GS expression increased significantly at 72 h. Finally, combined with the digital gene expression and pathway analysis, 1439 and 1346 genes were upregulated, and 2127 and 1927 genes were downregulated in the dsTPS3 and dsTPSs groups, respectively. The function of most differential genes was concentrated in sugar metabolism, lipid metabolism, and amino acid metabolism. The results indicated that TPS3 plays a key role in the energy metabolism of N. lugens and confirmed that TPS3 is a feasible target gene for RNA interference in N. lugens. Simultaneously, they provide a theoretical basis for the development and utilization of TPS3 to control pests.

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.

Transcriptome analysis reveals the roles of stem nodes in cadmium transport to rice grain

BACKGROUND

Node is the central organ of transferring nutrients and ions in plants. Cadmium (Cd) induced crop pollution threatens the food safety. Breeding of low Cd accumulation cultivar is a chance to resolve this universal problem. This study was performed to identify tissue specific genes involved in Cd accumulation in different rice stem nodes. Panicle node and the first node under panicle (node I) were sampled in two rice cultivars: Xiangwanxian No. 12 (low Cd accumulation cultivar) and Yuzhenxiang (high Cd accumulation cultivar). RNA-seq analysis was performed to identify differentially expressed genes (DEGs) and microRNAs.

RESULTS

Xiangwanxian No. 12 had lower Cd concentration in panicle node, node I and grain compared with Yuzhenxiang, and node I had the highest Cd concentration in the two cultivars. RNA seq analysis identified 4535 DEGs and 70 miRNAs between the two cultivars. Most genesrelated to the "transporter activity", such as OsIRT1, OsNramp5, OsVIT2, OsNRT1.5A, and OsABCC1, play roles in blocking the upward transport of Cd. Among the genes related to "response to stimulus", we identified OsHSP70 and OsHSFA2d/B2c in Xiangwanxian No. 12, but not in Yuzhenxiang, were all down-regulated by Cd stimulus. The up-regulation of miRNAs (osa-miR528 and osa-miR408) in Xiangwanxian No. 12 played a potent role in lowering Cd accumulation via down regulating the expression of candidate genes, such as bZIP, ERF, MYB, SnRK1 and HSPs.

CONCLUSIONS

Both panicle node and node I of Xiangwanxian No. 12 played a key role in blocking the upward transportation of Cd, while node I played a critical role in Yuzhenxiang. Distinct expression patterns of various transporter genes such as OsNRT1.5A, OsNramp5, OsIRT1, OsVIT2 and OsABCC1 resulted in differential Cd accumulation in different nodes. Likewise, distinct expression patterns of these transporter genes are likely responsible for the low Cd accumulation in Xiangwanxian No. 12 cultivar. MiRNAs drove multiple transcription factors, such as OsbZIPs, OsERFs, OsMYBs, to play a role in Cd stress response.