- Research article
- Open Access
Comparative transcriptome analysis revealed key factors for differential cadmium transport and retention in roots of two contrasting peanut cultivars
BMC Genomics volume 19, Article number: 938 (2018)
Peanut is the world’s fourth largest oilseed crop that exhibits wide cultivar variations in cadmium (Cd) accumulation. To establish the mechanisms of Cd distribution and accumulation in peanut plants, eight cDNA libraries from the roots of two contrasting cultivars, Fenghua 1 (low-Cd cultivar) and Silihong (high-Cd cultivar), were constructed and sequenced by RNA-sequencing. The expression patterns of 16 candidate DEGs were validated by RT-qPCR analysis.
A total of 75,634 genes including 71,349 known genes and 4484 novel genes were identified in eight cDNA libraries, among which 6798 genes were found to be Cd-responsive DEGs and/or DEGs between these two cultivars. Interestingly, 183 DEGs encoding ion transport related proteins and 260 DEGs encoding cell wall related proteins were identified. Among these DEGs, nine metal transporter genes (PDR1, ABCC4 and ABCC15, IRT1, ZIP1, ZIP11, YSL7, DTX43 and MTP4) and nine cell wall related genes (PEs, PGIPs, GTs, XYT12 CYP450s, LACs, 4CL2, C4H and CASP5) showed higher expression in Fenghua 1 than in Silihong.
Both the metal transporters and cell wall modification might be responsible for the difference in Cd accumulation and translocation between Fenghua 1 and Silihong. These findings would be useful for further functional analysis, and reveal the molecular mechanism responsible for genotype difference in Cd accumulation.
Cadmium (Cd) is one of the most widespread and toxic heavy metal pollutants, whose contamination in arable soil from natural and anthropogenic sources has drawn wide public health concern . Although Cd is a non-essential element for plant growth, it is easily taken up by the root system and transferred to aboveground parts of plants . The excess amount of Cd in plants can trigger negative effects on many physiological processes such as inhibition of photosynthesis, interference of nutrient uptake and induction of oxidative stress . Moreover, Cd can accumulate in the human body via the food chain, posing various health problems including itai-itai, cancer, renal dysfunction, osteoporosis and cardiovascular diseases . It is very important to take effective measures to reduce Cd accumulation in edible parts of crops. For this purpose, screening and breeding low-Cd cultivars are expected to be one of the sustainable strategies .
Peanut (Arachis hypogaea L.; AABB, 2n = 4x = 40) is an important oilseed crop that is cultivated extensively worldwide. China is one of the major exporters of peanut products, sharing around 40% of the world trade in peanut products . However, higher Cd enrichment capacity in the seeds of peanuts has become a serious problem for the peanut industry [6,7,8,9]. Seed Cd contents of peanuts in some regions of China such as Shandong and Liaoning have been shown to exceed the maximum level of the contaminant (0.5 mg/kg for Cd) following Chinese standard (GB 2762–2012) . Moreover, the higher Cd accumulation in peanut seeds has also been reported in other countries, such as Australia [9, 11] and Argentina . Exposure of plants to 2 and 4 mg kg− 1 Cd did not inhibit shoot biomass, seed yield, and oil content for most of the peanut cultivars . Thus, cultivation of peanuts in Cd contaminated soil may pose severe health risks for human beings.
Peanuts exhibit a wide variation in Cd accumulation in vegetative organs and seeds among cultivars [6,7,8, 13, 14]. Seed Cd concentrations of peanuts are not associated with Cd concentrations in vegetative organs at the mature stage, but positively correlated with Cd concentrations in shoots of the seedlings [6, 8]. The results suggest that shoot Cd concentrations at the vegetative growth stage may be important for determining cultivar differences in Cd accumulation in the seed [6, 8]. Although peanut seeds develop in the soil, their Cd concentration predominantly comes from the uptake of Cd by the root system [7, 9, 15]. Cd immobilization in the peanut root system could inhibit Cd translocation in peanut plants, and reduce Cd accumulation in peanut seed [16, 17]. Root morphology has been shown to significantly correlate with Cd accumulation in peanuts . Cultivars with longer fine roots have a high capacity of Cd accumulation . It seems that the root plays an important role in determining variation in Cd accumulation among peanut cultivars, however, the underlying mechanism is still unclear. Understanding the molecular mechanism of cultivar difference in Cd accumulation is essential for breeding low Cd cultivars of peanuts.
The uptake and translocation of Cd in plant is controlled by some metal transporter proteins such as natural resistance associated macrophage proteins (Nramp) [18, 19], zinc-regulated transporters, iron-regulated transporter-like protein (ZIP) , ATP-binding cassette transporters (ABC) , and P1B-ATPases [22, 23]. Furthermore, the translocation of Cd from roots to shoots was also determined by the capacity for retaining Cd in the roots, in which Cd binding to the cell wall and vacuolar Cd sequestration are involved [24, 25]. Therefore, we hypothesize that cultivar differences in Cd accumulation in peanut may be resulted from the differential expression of genes encoding metal transporters and cell wall modification in the roots.
In this study, a comparative transcriptome analysis was carried out on the roots of two peanut cultivars differing in seed Cd accumulation, Silihong (high Cd-accumulating cultivar) and Fenghua 1 (low Cd-accumulating cultivar) , under Cd-free and Cd-treated conditions. The main aims are: (i) to reveal the gene expression patterns in the roots of Fenghua 1 and Silihong in response to Cd exposure; (ii) to identify the key genes involved in Cd uptake, distribution and translocation in peanut roots; and (iii) to elucidate the gene regulatory network that is responsible for the cultivar differences of peanuts. The results presented here would be important for further illustrating the mechanism of Cd uptake and translocation as well as for breeding low Cd cultivars of peanuts.
Cd accumulation and translocation in plants
When exposed to 2 μM Cd for seven days, the concentration of Cd in roots and shoots was 224.7 and 18.3 mg kg− 1 dry weight for Fenghua 1, and 255.8 and 29.9 mg kg− 1 dry weight for Silihong, respectively (Fig. 1). Shoot Cd concentration was significantly higher in Silihong than in Fenghua 1, while the root Cd concentration was similar between the two cultivars. Similar to the shoot Cd concentration, the translocation factor of Cd from roots to shoots was significantly high in Silihong (0.12) than in Fenghua 1 (0.08) (Fig. 1). The current results, in agreement with previous findings , indicate that Fenghua 1 differed in the properties of Cd translocation and accumulation from Silihong.
RNA sequencing analysis of eight cDNA libraries
Eight cDNA libraries, namely, FCK_1, FCK_2, FCd_1, FCd_2, SCK_1, SCK_2, SCd_1 and SCd_2, were constructed from Cd-free (CK) and Cd-treated (Cd) peanut roots of both cultivars and sequenced using a BGISEQ-500 platform. A total of 164.87, 167.18, 164.59 and 167.03 Mb raw reads were produced from two biological replicate libraries for FCK, FCd, SCK and SCd, respectively. After removing the reads with adapters, low-quality and high content of unknown base (N), more than 63 million 100 bp paired-end clean reads with Q20 percentage over 97% remained in each cDNA library (Additional file 1: Table S1).
The clean reads from each cDNA library were aligned to reference genome sequences of the two progenitors of the cultivated peanut (A. duranensis-AA and A. ipaensis-BB). The average genome mapping rate is 90.16%, and the average gene mapping rate is 86.05% (Additional file 2: Table S2).
Analysis of gene expression levels
A total of 75,634 genes including 71,349 known genes and 4484 novel genes were identified from A and B genome (Additional file 3: Table S3 and Additional file 4: Table S4). Furthermore, 24,024 novel isoforms and 5264 noncoding transcripts were identified in eight cDNA libraries (Additional file 3: Table S3). The expression levels of these genes from the mapped libraries were normalized as FPKM (Additional file 5: Figure S1). The gene expression distribution (FPKM) was similar among eight libraries. For instance, the low expression level genes with FPKM interval 0~10 account for 63.50, 64.62, 64.66, 63.41, 64.04, 63.36, 59.24 and 63.56% in FCd_1, FCd_2, FCK_1, FCK_2, SCd_1, SCd_2, SCK_1 and SCK_2 respectively.
Identification of differentially expressed genes (DEGs)
Pairwise comparison analysis for each gene were performed between Fenghua 1 and Silihong (FCK/SCK and FCd/SCd) or between Cd-free and Cd-treated samples in each cultivar (FCd/FCK and SCd/SCK). DEGs were identified by the threshold of |log2 fold-change| ≥ 1 and Padj-value ≤ 0.05. A total of 6798 genes were differentially regulated in the four comparisons (Additional file 6: Table S5). Under control condition, 5297 DEGs were identified between two cultivars, while this value reached 2793 DEGs under Cd-treated condition (Fig. 2a). Of which, 1764 genes were common ones (Fig. 2c). Compared to Cd-free library, 331 genes were differentially regulated in Fenghua 1 in Cd exposure, including 278 up-regulated and 53 down-regulated genes. Whereas, 1302 Cd-treated DEGs were identified in Silihong, including 132 up-regulated and 1170 down-regulated genes (Fig. 2b). Among these genes, only 50 DEGs were common Cd-responsive genes in two cultivars (Fig. 2d), including 49 genes showed similar expression patterns and one gene showed opposing expression pattern (Additional file 6: Table S5).
Gene ontology (GO) analysis of DEGs
GO assignments were used to classify the functions of DEGs, a total of 2860 genes including 157 (FCd/FCK), 610 (SCd/SCK), 2214 (FCK/SCK) and 941 (FCd/SCd) genes, were assigned into 46 GO terms consisting of 22 biological process, 13 cellular component and 11 molecular function subcategories (Additional file 7: Figure S2a-b). Among these GO terms, the top two abundant categories of biological process (‘metabolic process’ and ‘cellular process’), molecular function (‘catalytic activity’ and ‘binding’), and cellular component (‘membrane’ and ‘membrane part’) were similar in the four comparisons.
Furthermore, the results of significantly enriched GO terms (corrected p-value < 0.01) are listed in Additional file 8: Table S6. For Cd-responsive DEGs, four GO terms such as response to inorganic substance, response to symbiont and response to metal ion in the biological process category and 16 GO terms such as zinc ion transmembrane transporter activity, divalent inorganic cation transmembrane transporter activity and transporter activity in the molecular function category were uniquely enriched in Fenghua 1. In Silihong, 11 GO terms such as transcription, DNA-templated, RNA biosynthetic process and nucleic acid-templated transcription in the biological process category and ten GO terms such as DNA binding transcription factor activity, transcription regulator activity and calcium ion binding in the molecular function category were uniquely enriched (Additional file 8: Table S6). For DEGs between Fenghua 1 and Silihong, three GO terms including lignin metabolic process, phenylpropanoid metabolic process and lignin catabolic process in the biological process category and the two GO terms including hydroquinone: oxygen oxidoreductase activity and oxidoreductase activity, acting on diphenols and related substances as donors, oxygen as acceptor (Additional file 8: Table S6) were significantly enriched under Cd-treated condition.
KEEG metabolic pathway analysis of DEGs
Here, DEGs were further annotated with Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.ad.jp/kegg/) to elucidate the molecular interactions among the genes. A total of 266 (FCd/FCK), 1073 (SCd/SCK), 4164 (FCK/SCK) and 2115 (FCd/SCd) DEGs were respectively assigned to 66, 116, 133, and 132 pathways (Additional file 9: Table S7), which were grouped into five groups including metabolism, genetic information processing, environmental information processing, cellular processes and organismal systems (Fig. 3). Among these pathways, the top two categories of them including metabolic pathway (59, 22.18%; 254, 23.67%; 960, 23.05%; 535, 25.30%) and biosynthesis of secondary metabolites (46, 17.29%; 193, 17.99%; 617, 14.82%; 337, 15.93%) were similar in the four comparisons (Additional file 9: Table S7). Followed by MAPK signaling pathway-plant, phenylpropanoid biosynthesis and ABC transporters in FCd/FCK; plant-pathogen interaction, MAPK signaling pathway-plant and flavonoid biosynthesis in SCd/SCK; Plant-pathogen interaction, MAPK signaling pathway-plant and plant hormone signal transduction in FCK/SCK; plant-pathogen interaction, RNA transport, protein processing in endoplasmic reticulum and phenylpropanoid biosynthesis in FCd/SCd (Additional file 9: Table S7).
In addition, one (FCd/FCK), six (SCd/SCK), four (FCK/SCK) and two (FCd/SCd) pathways were identified as significantly enriched (Q ≤ 0.01) (Additional file 10: Table S8). For DEGs between Fenghua 1 and Silihong, circadian rhythm-plant pathway was shared enriched under both control and Cd conditions (Additional file 10: Table S8). Besides, the pathway of cyanoamino acid metabolism was specifically significantly enriched for DEGs between two cultivars under Cd-treated. Further, Cd-responsive DEGs in Fenghua 1 were specifically significantly enriched in ABC transporters. The pathways of alpha-Linolenic acid metabolism, plant-pathogen interaction, flavonoid biosynthesis, circadian rhythm-plant, zeatin biosynthesis and MAPK signaling pathway-plant were significantly enriched for Cd-responsive DEGs in Silihong.
High expression level of DEGs between Fenghua 1 and Silihong under Cd exposure
To obtain a general statistical overview of transcriptional gene regulation the difference of shoot Cd accumulation in two peanut cultivars, the high expression level of DEGs (an FPKM value > 10 in at least one of four cDNA libraries and the threshold of |log2 (FCd/SCd)| ≥ 3) including 109 up- and 96 down-regulated were identified in FCd/SCd (Additional file 11: Table S9). Based on functional analysis, some of them were involved in cell wall metabolism (i.e., polygalacturonase inhibitory PGI, chitinase-like protein and GDSL esterase/lipase), ion transport (i.e., ABC transporter D family member 1 ABCD1, V-type proton ATPase), transcription regualtion (i.e., transcription factor bHLH83, bZIP transcription factor 11 and auxin-induced protein AUX28) and other metabolic processes (i.e., serine/threonine-protein kinase, metallothionein-like protein). The results suggested that these DEGs might play crucial regulatory roles in the the difference of shoot Cd accumulation in two peanut cultivars. Here, the DEGs were involved in the cell wall metabolism and ion transport as our main further investigation.
DEGs involved in heavy metal transport
According to GO functional annotation, 183 DEGs were identified to have high similarity with diverse transporters such as ATP-binding cassette (ABC) transporters (e.g. ABCAs, ABCBs, ABCCs and ABCGs), natural resistance-associated macrophage proteins (e.g. Nramp5 and Nramp6), ZIP transporters (e.g. IRT1, IRT3, ZIP1, ZIP4, ZIP5 and ZIP11), yellow stripe-like transporters (e.g. YSL3 and YSL7), MATE efflux family proteins (e.g. protein DETOXIFICATIONs DTXs) and metal tolerance proteins (e.g. MTP4 and MTP9) (Additional file 12: Table S10). These transporters were classified into seven categories based on their expression patterns (Additional file 13: Table S11). The classes from first to fourth were Cd-responsive DEGs. Of these, two genes in the first class were up-regulated, and one gene in the second class was down-regulated in both Fenghua 1 and Silihong after Cd exposure; the 22 genes in the third class and five genes in the fourth class were up-regulated in Fenghua 1 and Silihong roots under Cd exposure, respectively. Furthermore, the fifth and sixth classes were DEGs between two cultivars. The fifth class contained 58 genes whose expressions were higher in Fenghua 1 than those in Silihong, whereas, the sixth class had 61 genes expressed lower in Fenghua 1 than those in Silihong. The 34 genes in seventh class were Cd-responsive DEGs that is also differentially expressed between Fenghua 1 and Silihong.
We also found that the DEGs encoding heavy metal-mediated transporters, such as ABCs, IRTs, ZIPs, YSL3, Nramp5, MTP9, pleiotropic drug resistance proteins (PDRs) and copper-transporting ATPase HMA5 (HMA5), were all up-regulated by Cd in Fenghua 1 (Additional file 12: Table S10). In SCd/SCK, besides the DEGs encoding transporters including ABCB15, ABCG20, HMA5, Nramp6, DTX40 and two pore calcium channel protein 1 (TPCN1), the remaining 22 DEGs encoding transporters such as heavy metal-associated isoprenylated plant proteins (HIPPs), calcium-transporting ATPase 2 (Ca2+-ATPase2), MATE efflux family protein (MATE), and ABCG11, were down-regulated after Cd exposure (Additional file 12: Table S10). For the FCd/SCd, genes encoding ABC transporters (e.g. ABCB28, ABCC4/8/15, ABCD1 and ABCG15), oligopeptide transporter 7 (OPT7), DTX43/49, IRT1, PDR1, YSL7 and MTP4, were higher in Fenghua 1 than in Silihong. However, ABCA1, ABCB1/11, ABCC10, ABCG7, HIPP35, ZIP1, Nramp2, Ca2+-ATPase12, cation/H(+) antiporter (CHX18/20), DTX16/43/48/51 and MATE expressed higher in Silihong than in Fenghua 1 (Additional file 12: Table S10).
DEGs related to cell wall modification and vacuolar sequestration
Plant cell wall has been considered as the critical site for retention of Cd, especially in roots. Our study suggested that a total of 260 DEGs were homologous with cell wall metabolism-related genes, including 108 involved in cell wall degradation (75, 23 and five were involved in pectin-, hemicellulose- and cellulose-degrading process, respectively; five were other cell wall degrading-related enzymes), 143 involved in cell wall synthesis (53 in cellulose biosynthesis, 49 in suberin biosynthesis, 24 in lignin biosynthesis, 11 in hemicellulose, four in pectin, and two in Casparian strip formation). Furthermore, three encoding expansin proteins, five encoding COBRA proteins and one encoding extension 2 like protein were involved in cell wall organization (Additional file 14: Table S12).
For the 108 cell wall degradation DEGs, six were Cd-responsive DEGs including one upregulated in Fenghua 1 and five downregulated in Silihong, and 87 were DEGs between Fenghua 1 and Silihong (Additional file 14: Table S12). The remaining 15 genes were Cd responsive DEGs, which also expressed differentially between two cultivars, and 13 of them were down-regulated by Cd in Silihong (Additional file 14: Table S12). Furthermore, among 108 DEGs, only three upregulated were found in FCd/FCK. Remarkably, 32 of 108 DEGs were expressed higher in Fenghua 1 than in Silihong. Regarding the 143 cell wall synthesis DEGs, 19 were Cd-responsive DEGs and 97 were DEGs between Fenghua 1 and Silihong. For the 19 Cd responsive DEGs, 15 were up-regulated by Cd in Fenghua 1, and four were down-regulated by Cd in Silihong. The remaining 27 genes were Cd responsive DEGs, which also expressed differentially between two cultivars. 46 of 143 DEGs were expressed higher in Fenghua 1.
In addition, the formation of Cd complexes with glutathione (GSH), metallothionein (MT), phytochelatin (PC) and nicotianamine (NA) is another way used by plants to sequestrate Cd in vacuolar. In this study, two (107459596, 107612042) and one (107465178) nicotianamine synthase (NAS) genes were up-regulated by Cd in Fenghua 1 and Silihong, respectively. Moreover, NAS (107612042) and MT2 (107639467) were expressed higher in Fenghua 1 than in Silihong (Additional file 14: Table S12). However, no significant difference in GSH and PC biosynthesis-related genes expression were observed in each comparison.
To verify the differential expression patterns of DEGs identified in our transcriptome data, a total of ten candidate transport- and six cell wall-related genes were selected for RT-qPCR analysis. The results indicated that the expression patterns of 15 out of 16 DEGs showed good agreement with the RNA-Seq-based gene expression patterns, and only one gene (ABCC36) was not well consistent with the results of sequencing (Fig. 4; Additional file 6: Table S5). For example, several Cd-responsive up-regulated DEGs also indicated relatively high expression level in Fenghua 1 (ZIP4, Nramp5/6, and ABCC3) and Silihong (Nramp6) under Cd exposure in RT-qPCR analysis (Fig. 4). For DEGs between Fenghua 1 and Silihong, the expression level of ABCC4, PDR1, DTX43, β-Gal8, XYL1, EG10 genes were higher in Fenghua 1, whereas LAC11 was highly expressed in Silihong under Cd exposure (Fig. 4). These results suggested that the data detected by RNA-Seq analysis of this study provide reliable inference.
Comparison of gene expression in roots between Fenghua 1 and Silihong
There are clear physiological differences in shoot Cd concentrations between Fenghua 1 and Silihong following Cd exposure . However, the molecular mechanism underlying the difference in Cd accumulation among peanut cultivars has not been fully elucidated.
RNA sequencing (RNA-seq) analysis has been used for revealing the molecular mechanisms of differential Cd uptake, accumulation and translocation in many plants [26,27,28]. However, the comprehensive identification of the DEGs involved in root uptake, accumulation and transport of Cd in peanut remained unexplored. Based on peanut genome , a total of 6798 genes including 331 (FCK/FCd), 1302 (SCK/SCd), 5297 (FCK/SCK) and 2793 (FCd/SCd) genes, were identified significant differentially expression in pairwise comparisons (Additional file 6: Table S5). Of them, 183 DEGs encoding transporter proteins and 260 DEGs encoding cell wall metabolism-related proteins were identified. Most Cd-responsive DEGs (278/331, 83.99%) in Fenghua 1 were up-regulated by Cd, while the majority of Cd-responsive DEGs (1170/1302, 89.86%) in Silihong were down-regulated by Cd (Additional file 6: Table S5). The results suggested that the two cultivars differed in the molecular mechanisms in response to Cd exposure. To our knowledge, this is the first research to identify and dissect extensive Cd-responsive genes in peanuts, and insight into molecular mechanisms responsible for the differential shoot Cd accumulation in two cultivars.
Heavy metal transporters involved in differential root uptake and translocation of Cd between Fenghua 1 and Silihong
It has been widely proven that Zn2+, Ca2+, Fe2+, Mn2+, Cu2+ and Mg2+ transporters could transport Cd into plant root cells, and affect root-to-shoot Cd translocation . The physiological data indicated that Fenghua 1 showed lower shoot Cd concentrations than Silihong, for which the root-to-shoot translocation of Cd might be responsible . In this study, we identified 183 DEGs encoding ion transporter proteins, among which 30 were Cd-responsive DEGs. We also found that 34 Cd-responsive DEGs encoding ion transporter proteins were expressed differentially between the two cultivars (Additional file 12: Table S10). Compared with the control, Cd up-regulated the Cd-responsive DEGs in Fenghua 1, whereas most of DEGs in Silihong were down-regulated by Cd. Only four DEGs, 107471962 and 107623484 (HMA5), 107477738 (DTX40) and 107491040 (HIPP2), were induced by Cd treatment in both cultivars. These results indicate that the two peanut cultivars differed in the expression of ion transporter that may contribute to the difference in Cd uptake and translocation.
We found that six DEGs encoding metal transporters related to Cd influx, including IRT1, ZIP4 and Nramp5/6 (Table 1). AtIRT1 , NcZNT1/AtZIP4 , SaNramp6  and OsNramp5  have been shown to mediate cellular metals (Zn, Fe, Mn and Cd) uptake across the plasma membrane. IRT1, a Fe-regulated transporter, is expressed in the root epidermis . Overexpressing AtIRT1 in Arabidopsis leds to increased accumulation of Cd and Zn in the transgenic lines under Fe starvation , which indicates that AtIRT1 can affect the Cd uptake capacity. Therefore, under Cd exposure, the higher expression of IRT1 in the roots of Fenghua 1 compared with Silihong may be involved in the difference of root Cd absorption between Fenghua 1 and Silihong (Figs. 5 and 6). Overexpression of AtZIP4 and SaNramp6 in Arabidopsis increased accumulation of Cd in transgenic plants [20, 32]. AtNramp6 is also thought to be a major role in the intracellular distribution of Cd . Our results showed that ZIP4, Nramp5 and IRT1 were up-regulated by Cd in Fenghua 1, while in Silihong, they were unchanged (Fig. 5; Table 1). Nramp5 (107606033, 107470537 and 107460699) shows a high degree of homology to the AtNramp6 from Arabidopsis. However, Nramp6 (107618045) was upregualted in Silihong, while in Fenghua 1, it was unchanged (Fig. 5; Table 1). ZIP4 and Nramp5/6 were not DEGs as judged by our criteria between the two cultivars. These results suggested that ZIP4 and Nramp5/6 may be related to Cd-response in roots of two peanut cultivars. However, the functions of these genes in peanuts need further investigated.
Furthermore, several DEGs encoding metal transporter proteins related to sequestration of Cd into vacuolar (ABCC3/4 and MTP4) and the endomembrane systems (ZIP11) are also observed (Table 1). The tonoplast-localized AtABCC3/AtMRP3 has been reported to regulated Cd detoxification via sequestration of PC-Cd complexes into vacuolar . AtMRP4/AtABCC4 mediate the transport of folate into vacuolar . Here, although most genes encoding ABCC3 were up-regulated by Cd in Fenghua 1 (Fig. 5), they were unchanged between Fenghua 1 and Silihong under Cd exposure (Table 1). This result suggests that ABCC3 may not be involved in the difference of shoot Cd accumulation between the two cultivars. Two ABCC15 (107481980, 107632015) genes were close homologue of AtMRP3 in Arabidopsis (Additional file 12: Table S10). Importantly, ABCC4 (107466097, 107630594), ABCC15 (107481980, 107632015) and MTP4 showed higher expression in Fenghua 1 roots than in Silihong roots under Cd exposure (Fig. 6; Table 1). CsMTP4 was shown to function in Zn and Cd transport and tolerance by sequestering Cd into the vacuolar compartments in cucumber . Xu et al.  reported that ZIP11 is located in the endomembrane, which could compartmentalize Zn/Cd into vesicles and prevent metal translocation from root to shoot. Here, two ZIP11 (107462023, 107617536) genes were up-regulated by Cd in Fenghua 1, while in Silihong, they were unchanged (Fig. 5; Table 1). FPKM value analysis indicated that ZIP11 (107462023) showed higher expression in Fenghua 1 roots than in Silihong roots (Fig. 6). These results suggest that Cd-inducible ABCC4, ABCC15, MTP4 and ZIP11 may enhance root vacuoles or vesicles sequestration of Cd and consequently, reducing root-to-shoot Cd translocation in Fenghua 1. Whether these genes associate with sequestration of Cd in peanut plants or not require further studies.
We also identified some DEGs encoding plasma membrane-localized transporter proteins related to Cd efflux, including ABCG36/PDR8, PDR1, YSL3/7, MTP9, DTX43 and ZIP1/5 (Table 1). These genes have been reported to regulate in Cd accumulation or detoxification in other plants [38,39,40,41]. Migocka et al.  reported that CsMTP9 in cucumber was a plasma membrane H+-coupled antiporter involved in the Mn2+ and Cd2+ efflux from root cells. Here, MTP9 was up-regulated by Cd in Fenghua 1, but in Silihong, it was unchanged, and it showed no significant differential expression between the two cultivars under control and Cd-treated (Fig. 5; Table 1). Peanut PDR 1 and DTX43 were homology with AtPDR12 and AtFRD3 in Arabidopsis, respectively (Additional file 12: Table S10). Two PDR 1 (107647851, 107491711) were upregulated in Fenghua 1, but in Silihong, it was unchanged. Besides, PDR1, DTX43 (107462988) and YSL7 transporters showed higher expression in Fenghua 1 than in Silihong under Cd exposure (Fig. 6; Table 1). This observation is also consistent with previous reports that the expression of YSL7 and MATE efflux family protein indicated higher expression in QLQ (low-shoot-Cd) than in T308 (high-shoot-Cd) . In Arabidopsis, two ABC genes, AtABCG36/PDR8 and AtABCG40/PDR12, are induced by Cd as an efflux pump of Cd2+ conferring detoxification of Cd [40,41,42]. In contrast, ABCG36/PDR8 showed higher expression in Silihong than in Fenghua 1 under control condition, while it was unchanged under Cd-treated (Fig. 6; Table 1). Further, plasma membranes-localized YSL transporter family can involve in transporting of NA-metal complexes, and BjYSL7  and SnYSL3  have been reported to participate in the efflux of NA-Cd complexes. YSL3 was upregulated by Cd in Fenghua 1 (Fig. 5; Table 1), but no difference was observed between two cultivars. FRD3, a MATE efflux family protein, is thought to be involved in Cd tolerance via the removing Cd from the cytoplasm [44, 45]. These results indicated that PDR1, DTX43 and YSL7 may be the critical factors in determining the contrasting shoot-Cd accumulation in two cultivars. In addition to genes mentioned above, four genes belonging to ZIP family, ZIP1 (107461527, 107614121 and 107482454) and ZIP 5 (107638316) were upregulated by Cd in Fenghua 1, but in Silihong, it was unchanged (Fig. 5; Table 1). ZIP1 (107614121) and ZIP1 (107482454) of peanut exhibit high sequence similarity to Arabidopsis AtZIP2 and peanut ZIP5, respectively (Additional file 12: Table S10). Milner et al.  reported that AtZIP1 was localized at the vacuolar membrane, while AtZIP2 was localized to plasma membrane of root stele cells, and both transporters play a role in Mn and Zn translocation from roots to shoots. ZIP1 (107458727) shows lower expression in Fenghua 1 than in Silihong under Cd exposure (Fig. 6; Table 1). Therefore, the lower expression of ZIP1 in the roots of Fenghua 1 may reduce Cd translocation from roots to shoots, and the exact role of this transporter needs further investigation.
Retention of Cd in the root may contribute to the contrasting shoot-Cd accumulation in Fenghua 1 and Silihong
Cell wall is the first barrier to Cd uptake into the cytosol . It is also the major site for Cd2+ binding and retention in plant roots . The Cd binding capacity of root cell walls could influence the cultivar difference in Cd translocation to the shoot . In this study, a total of 260 DEGs (FCd/FCK: 27; SCd/SCK: 40; FCd/SCd: 93; FCK/SCK: 191) were identified to show high similarity with cell wall-related proteins (Additional file 14: Table S12).
Twenty two genes encoding enzymes involved in the biosynthesis of cellulose, hemicellulose, lignin and suberin were up-regulated by Cd in Fenghua 1, including 13 UDP-glycosyltransferase family genes (e.g., UGT71K1, UGT72D2 and UGT83A1), four cytochrome P450 family genes (CYP450s), three 4-coumarate: CoA ligase-like proteins (4CL9), cellulose synthase-like protein G2 (CSLG2) and caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) (Additional file 14: Table S12). In contrast, only one gene encoding suberin synthesis enzymes (CYP450 704C1, 107475575) was up-regulated by Cd in Silihong, while the remained 21 DEGs (95.24%) encoding cellulose, lignin and suberin biosynthesis enzymes were down-regulated (Additional file 14: Table S12). These results suggested that the cell wall biosynthesis was more active in Fenghua 1 than in Silihong, and this may lead to different binding capacity of Cd to cell walls.
It has been reported that pectin and hemicellulose components of the cell wall can bind Cd2+ by the negatively charged carboxyl and hydroxyl groups . Thus, increased pectin and hemicellulose content might prevent the translocation of Cd from the roots to shoots . PME catalyzes pectin demethylation, which enhances the capacity of pectin to bind metal ions . Overexpression OsPME14 in transgenic rice increased Al retention in the root tip cell wall . Polygalacturonase inhibitor proteins (PGIPs) can inhibit polygalacturonase (PG) activity and pectin-degrading. Moreover, xyloglucan hemicellulose component has been reported to be a much more effective binding of Al than pectin in Arabidopsis . Here, most of DEGs (59/93, 63.44%) were expressed higher in Fenghua 1 than that in Silihong under Cd-treated conditions, including 24 DEGs involved in pectin catabolic process, eight DEGs involved in suberin synthesis and cellulose catabolic process, respectively, seven DEGs involved in hemicellulose catabolic process, five DEGs involved in lignin synthesis, one DEG involved in the formation of Casparian strip and the other three involved in cell wall modification (Additional file 14: Table S12). Among of them, the higher expression of pectinesterase (PEs, 107475981, 107630347), PGIPs (i.e., 107481713, 107481714 and 107631575), xyloglucan glycosyltransferase 12 (XYT12) and glycosyltransferases (GTs, 107638833, 107459345) in Fenghua 1 might increase the retention of Cd in root cell wall and consequently, reducing Cd influx into the cytosol.
Furthermore, the intensive wall thickenings were observed in the inner tangential walls of endodermal cells after exposure to Cd , which could cause by the formation of Casparian strips and suberin lamellae, increase of suberin and lignin contents. These changes can serve as a barrier limiting apoplasmic flow of Cd to xylem and its translocation to shoots [53, 54]. In this study, the higher expression of casparian strip membrane protein 5 (CASP5, 107634327), CYP450s (i.e., 107475575, 107478900, and 107482119), 4-coumarate--CoA ligase 2 (4CL2), laccases (LACs, 107639475, 107496989 and 107625654) and cinnamate 4-monooxygenase (C4H) in Fenghua 1 may contribute to cell wall modification in the roots and consequently leading to lower root-to-shoot Cd translocation.
The retention of Cd in the roots can be mediated by metal-chelating molecules such as MTs, GSHs, PCs and NAs [25, 55]. In this study, NAS (107612042) was upregulated by Cd in Fenghua 1, while in Silihong, it was unchanged. Moreover, NAS (107612042) and MT2 (107639467) were expressed higher in Fenghua 1 than in Silihong under Cd exposure (Additional file 14: Table S12), suggesting that these two genes may be involved in the difference of root Cd sequestration between Fenghua 1 and Silihong. However, no significant difference in GSH and PC biosynthesis-related genes expression were observed in each comparison, suggesting that GSHs and PCs may not be the key factors regulating the cultivar difference in Cd retention in peanut roots.
In this study, 183 DEGs encoding ion transport related proteins (i.e., ZIPs, ABCs, Nramps, YSLs and MTPs) were identified. Importantly, one Cd influx gene (IRT1), four Cd efflux genes (PDR12, YSL7, ZIP1 and FRD3) and four tonoplast- and endomembrane-localized transport genes (ZIP11, MTP4, ABCC4 and ABCC15), showed higher expression in Fenghua 1 than in Silihong, might be responsible for the cultivar differences in shoot Cd accumulation in peanut. Furthermore, 260 DEGs encoding cell wall-related proteins were also identified. Cell wall biosynthesis may be more activated in Fenghua 1. The higher expression of PEs, PGIPs, GTs, XYT12 CYP450s, LACs, 4CL2, C4H and CASP5 genes in Fenghua1 might contribute to reducing movement of Cd to xylem and its translocation to shoots, thus these genes might be involved in cultivar difference in shoot Cd accumulation of peanut. These results could provide a new insight into the molecular mechanisms about cultivar difference in Cd accumulation, which would be essential for breeding low-Cd cultivars of peanuts.
Plant materials and experiment design
According to our previous study , two peanut cultivars differing in seed Cd concentrations, Fenghua 1 (F, low-Cd cultivar) and Silihong (S, high-Cd cultivar) were selected for this experiment. Seeds of the two peanut cultivars were obtained commercially from the Peanut Institute of Shandong Province, Qingdao, China. Seeds were surface sterilized in 5% sodium hypochlorite for 1 min, then soaked in distilled water for 24 h, and sown into pots with washed sand. One week after sowing, the seedlings with uniform size were transferred to polyethylene pot (six seedlings per pot) irrigated with 3.5 L of nutrient solution (pH 5.8) for two weeks. Then, the seedlings for each cultivar were treated with 0 (Cd-free) and 2 μM CdCl2 (Cd-treated), respectively. The experiment was arranged in a randomized complete design with three replications (pots). Peanut plants were cultured according to the conditions used by Su et al. . After one week of treatment, root samples for RNA-seq and RT-qPCR analysis were collected separately from each Cd-free and Cd-treated seedling. A total of five independent biological replicates were sampled from a pool of six seedlings for RNA-Seq (two biological replicates) and RT-qPCR validation (three biological replicates). All samples were immediately frozen in liquid nitrogen and stored at − 80 °C. The Cd concentration of roots and shoots was determined by atomic absorption spectrometry as previously described .
cDNA library construction and RNA sequencing
Total RNA was isolated from Cd-free and Cd-treated peanut roots of both cultivars using Trizol® Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. Then the RNA integrity and purity were assessed using Agilent 2100 Bioanalyzer (Agilent, USA) and NanoDrop™ spectrophotometer (Thermo Fisher Scientific, USA), respectively. Eight cDNA libraries named FCK_1, FCK_2, FCd_1, FCd_2, SCK_1, SCK_2, SCd_1 and SCd_2 were constructed with the BGISEQ-500 platform combining the DNA nanoball-based nanoarrays and stepwise sequencing using combinational probe-anchor synthesis sequencing method (BGI, Shenzhen, China). Briefly, Poly (A) mRNAs were captured from total RNA with poly-T oligo-attached magnetic beads (Thermo Fisher Scientific, USA), and then were fragmented into small pieces using divalent cations and heat treatments. Next, the first-strand cDNA was synthesized using reverse transcriptase and random primers with small fragments of mRNA as templates; and the second-strand cDNA was synthesized via DNA polymerase I and RNase H. Thereafter, the double-stranded cDNA ends were further repaired by a single nucleotide ‘A’ base addition, and subsequent ligation of the adapter. After PCR amplification, the PCR yields were quantified by Qubit (Invitrogen) and pooled samples together to make a single strand DNA circle (ssDNA circle), which gave the final cDNA library.
DNA nanoballs (DNBs) were generated with the ssDNA circle by rolling circle replication (RCR) to enlarge the fluorescent signals at the sequencing process. The DNBs were loaded into the patterned nanoarrays and pair-end reads of 100 bp were read through on the BGISEQ-500 platform for the following data analysis study.
Data filtering and mapping
Using SOAPnuke software (BGI, Shenzhen, China), the raw reads were initially processed by removing reads containing adaptors, those with more than 5% unknown N bases and low-quality reads, and the clean reads in each library were obtained as the basis of further data analysis. Then, the clean reads were aligned to the reference genome sequences of peanut (A. duranensis-AA and A. ipaensis-BB) (Aradu1.1, https://www.ncbi.nlm.nih.gov/assembly/GCF_000817695.2; Araip1.0, https://www.ncbi.nlm.nih.gov/assembly/GCA_000816755.1)  using HISAT2 (v2.2.4, http://www.ccb.jhu.edu/software/hisat) with default parameters. The mapping reads for each library were reconstructed using StringTie software (v1.0.4, http://ccb.jhu.edu/software/stringtie); and then Cuffmerge was used to integrate all the reconstructed transcripts from the eight libraries. Novel transcripts were obtained by comparing reconstructed transcripts to reference annotation using Cuffcompare (Cufflinks v2.2.1, http://cole-trapnell-lab.github.io/cufflinks); and then, the coding potential of novel transcripts were predicted using CPC (v0.9-r2, http://cpc.cbi.pku.edu.cn); and a complete reference was generated via merging the novel coding transcripts with the reference transcripts. Finally, the clean reads were aligned to the complete reference using Bowtie2 program (v2.2.5, http://bowtie-bio.sourceforge.net/bowtie2/index.shtml).
Gene expression analysis
Expected count tables for each gene were obtained by RSEM (v1.2.12, http://deweylab.biostat.wisc.edu/RSEM) package. The expression levels for each of the gene in each library were normalized to FPKM (fragments per kilobase of transcript per million mapped reads). Based on the negative binomial distribution, read counts were used to determine the differential expression analysis among four groups (FCK, FCd, SCK and SCd, two biological replicates per group) by the DEseq 2 package as described by Love et al. . The |log2 fold-change| ≥ 1 and Padj-value ≤ 0.05 were used as thresholds to determine the significance of the difference in gene expression.
Gene ontology (GO) and pathway analysis of DEGs
With GO and pathway annotation result, the GO and KEGG functional classification for the DEGs were performed according to official classification, respectively. Then, GO and KEGG pathway enrichment analysis was implemented using phyper in R package with a corrected p-value ≤ 0.01 to identify significantly enriched GO terms and biological pathways.
Total RNA was isolated from three biological replicates of root samples of both peanut cultivars exposed to 0 and 2 μM CdCl2 treatments, respectively. Then the first strand cDNA was synthesized by Prime Script® RT reagent Kit (Takara, Dalian, China). The primers of 16 DEGs were designed using Beacon Designer 7.0 software (Premier Biosoft International, USA), and Actin gene was used as the internal control (Additional file 15: Table S13). The primer specificity was evaluated by BLASTing primer sequences against the NCBI database. RT-qPCR was performed using a SYBR Premix EX Taq Kit (Takara) in a 20 μl reaction mixture on an ABI7300 (Applied Biosystems, Foster City, CA, USA) following the method described by Yu et al. . All reactions were performed in three technical replicates. The equation ratio = 2−ΔΔCτ was employed to calculate the relative expression level of the selected genes. Data statistical analysis with IBM SPSS Statistics Version 25 software (Armonk, NY, USA) was conducted using Duncan’s multiple range test at the P < 0.05 level of significance.
4-coumarate: CoA ligase-like protein
ATP-binding cassette transporter protein
Casparian strip membrane protein
Cation diffusion facilitator
Cellulose synthase-like protein
Differentially expressed gene
Cadmium/zinc/copper-transporting ATPase (heavy metal-associated domain)
Kyoto encyclopedia of genes and genomes
Metal tolerance protein
Natural resistance-associated macrophage protein
Pleiotropic drug resistance-type ABC transporter protein
Polygalacturonase inhibitor protein
Reverse-transcription quantitative PCR
Yellow stripe-like transporter
Zinc-regulated transporter/Iron-regulated transporter-like Protein
Tóth G, Hermann T, Da Silva MR, Montanarella L. Heavy metals in agricultural soils of the European Union with implications for food safety. Environ Int. 2016;88:299–309.
Clemens S, Aarts MGM, Thomine S, Verbruggen N. Plant science: the key to preventing slow cadmium poisoning. Trends Plant Sci. 2013;18(2):92–9.
Benavides MP, Gallego SM, Tomaro ML. Cadmium toxicity in plants. Braz J Plant Physiol. 2005;17:21–34.
Grant CA, Clarke JM, Duguid S, Chaney RL. Selection and breeding of plant cultivars to minimize cadmium accumulation. Sci Total Environ. 2008;390(2–3):301.
Dai X, Bai Y, Jiang J, Chen X, Zhou H, Yin N, Chen L, Ding X, Li P. Cadmium in Chinese postharvest peanuts and dietary exposure assessment in associated population. J Agr Food Chem. 2016;64(41):7849–55.
Shi G, Su G, Lu Z, Liu C, Wang X. Relationship between biomass, seed components and seed Cd concentration in various peanut (Arachis hypogaea L.) cultivars grown on Cd-contaminated soils. Ecotoxicol Environ Saf. 2014;110:174–81.
Su Y, Wang X, Liu C, Shi G. Variation in cadmium accumulation and translocation among peanut cultivars as affected by iron deficiency. Plant Soil. 2013;363(1–2):201–13.
Liu C, Yu R, Shi G. Effects of drought on the accumulation and redistribution of cadmium in peanuts at different developmental stages. Arch Agr Soil Sci. 2016;63(8):1049–57.
Mclaughlin MJ, Bell MJ, Wright GC, Cozens GD. Uptake and partitioning of cadmium by cultivars of peanut (Arachis hypogaea L.). Plant Soil. 2000;222(1–2):51–8.
Chen Z, Song W, Li P, Ding X, Wang M, Qian Y. Survey and dietary risk assessment of cadmium in peanut produced in China. J Agro-Environ Sci. 2012;31(2):237–44.
Mclaughlin MJ, Bell MJ, Wright GC, Cruickshank A. Inter- and intra-specific variation in accumulation of cadmium by peanut, soybean, and navybean. Aust J Agric Res. 1997;48(8):1151–60.
Bianucci E, Fabra A, Castro S. Cadmium accumulation and tolerance in Bradyrhizobium spp. (peanut microsymbionts). Curr Microbiol. 2011;62(1):96–100.
Su G, Li F, Lin J, Liu C, Shi G. Peanut as a potential crop for bioenergy production via Cd-phytoextraction: A life-cycle pot experiment. Plant Soil. 2013;365(1):337–45.
Lu Z, Zhang Z, Su Y, Liu C, Shi G. Cultivar variation in morphological response of peanut roots to cadmium stress and its relation to cadmium accumulation. Ecotoxic Environ Saf. 2013;91(4):147–55.
Popelka JC, Schubert S, Schulz R, Hansen AP. Cadmium uptake and translocation during reproductive development of peanut (Arachis hypogaea L.). Angew Bot. 1996;70(3):140–3.
Bianucci E, Sobrino-Plata J, Carpena-Ruiz RO, del Carmen Tordable M, Fabra A, Hernández LE, Castro S. Contribution of phytochelatins to cadmium tolerance in peanut plants. Metallomics. 2012;4(10):1119–24.
Su Y, Liu J, Lu Z, Wang X, Zhang Z, Shi G. Effects of iron deficiency on subcellular distribution and chemical forms of cadmium in peanut roots in relation to its translocation. Environ Exp Bot. 2014;97(1):40–8.
Takahashi R, Ishimaru Y, Senoura T, Shimo H, Ishikawa S, Arao T, Nakanishi H, Nishizawa NK. The OsNRAMP1 iron transporter is involved in cd accumulation in rice. J Exp Bot. 2011;62(14):4843–50.
Ishimaru Y, Takahashi R, Bashir K, Shimo H, Senoura T, Sugimoto K, Ono K, Yano M, Ishikawa S, Arao T, Nakanishi H, Nishizawa NK. Characterizing the role of rice NRAMP5 in manganese, iron and cadmium transport. Sci Rep. 2012;2(6071):286.
Lin YF, Hassan Z, Talukdar S, Schat H, Aarts MGM. Expression of the ZNT1 zinc transporter from the metal hyperaccumulator Noccaea caerulescens confers enhanced zinc and cadmium tolerance and accumulation to Arabidopsis thaliana. PLoS One. 11(3):e0149750.
Brunetti P, Zanella L, De Paolis A, Di Litta D, Cecchetti V, Falasca G, Barbieri M, Altamura MM, Costantino P, Cardarelli M. Cadmium-inducible expression of the ABC-type transporter AtABCC3 increases phytochelatin-mediated cadmium tolerance in Arabidopsis. J Exp Bot. 2015;66(13):3815–29.
Takahashi R, Ishimaru Y, Shimo H, Ogo Y, Senoura T, Nishizawa NK, Nakanishi H. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice. Plant Cell Environ. 2012;35(11):1948–57.
Verret F, Gravot A, Auroy P, Leonhardt N, David P, Nussaume L, Vavasseur A, Richaud P. Overexpression of AtHMA4 enhances root-to-shoot translocation of zinc and cadmium and plant metal tolerance. FEBS Lett. 2004;576(3):306–12.
Fernández R, Fernández-Fuego D, Bertrand A, González A. Strategies for Cd accumulation in Dittrichia viscosa (L.) Greuter: role of the cell wall, non-protein thiols and organic acids. Plant Physiol Biochem. 2014;78:63–70.
Nocito FF, Lancilli C, Dendena B, Lucchini G, Sacchi GA. Cadmium retention in rice roots is influenced by cadmium availability, chelation and translocation. Plant Cell Environ. 2011;34(6):994–1008.
Yu R, Li D, Du X, Xia S, Liu C, Shi G. Comparative transcriptome analysis reveals key cadmium transport-related genes in roots of two pak choi (Brassica rapa L. ssp. chinensis) cultivars. BMC Genomics. 2017;18(1):587.
Huang YY, Shen C, Chen JX, He CT, Zhou Q, Tan X, Yuan JG, Yang ZY. Comparative transcriptome analysis of two Ipomoea aquatica Forsk. Cultivars targeted to explore possible mechanism of genotype-dependent accumulation of cadmium. J Agr Food Chem. 2016;64(25):5241–50.
Xu J, Sun J, Du L, Liu X. Comparative transcriptome analysis of cadmium responses in Solanum nigrum and Solanum torvum. New Phytol. 2012;196(1):110–24.
Bertioli DJ, Cannon SB, Froenicke L, Huang G, Farmer AD, Cannon EKS, Liu X, Gao D, Clevenger J, Dash S, Ren L, Moretzsohn MC, Shirasawa K, Huang W, Vidigal B, Abernathy B, Chu Y, Niederhuth CE, Umale P, Araújo ACG, Kozik A, Kim KD, Burow MD, Varshney RK, Wang X, Zhang X, Barkley N, Guimarães PM, Isobe S, Guo B, Liao B, Stalker HT, Schmitz RJ, Scheffler BE, Leal-Bertioli SCM, Xun X, Jackson SA, Michelmore R, Ozias-Akins P. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat Genet. 2016;48(4):438–46.
Koźmińska A, Wiszniewska A, Hanus-Fajerska E, Muszyńska E. Recent strategies of increasing metal tolerance and phytoremediation potential using genetic transformation of plants. Plant Biotechnol Rep. 2018;12(1). https://doi.org/10.1007/s11816-11017-10467-11812.
Connolly EL, Fett JP, Guerinot ML. Expression of the IRT1 metal transporter is controlled by metals at the levels of transcript and protein accumulation. Plant Cell. 2002;14(6):1347–57.
Chen S, Han X, Fang J, Lu Z, Qiu W, Liu M, Sang J, Jiang J, Zhuo R. Sedum alfredii SaNramp6 metal transporter contributes to cadmium accumulation in transgenic Arabidopsis thaliana. Sci Rep. 2017;7(1):13318.
Sasaki A, Yamaji N, Yokosho K, Ma JF. Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. Plant Cell. 2012;24(5):2155–67.
Vert G, Grotz N, Dédaldéchamp F, Gaymard F, Guerinot ML, Briat JF, Curie C. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002;14(6):1223–33.
Cailliatte R, Lapeyre B, Briat JF, Mari S, Curie C. The NRAMP6 metal transporter contributes to cadmium toxicity. Biochem J. 2009;422(2):217–28.
Klein M, Geisler M, Suh SJ, Kolukisaoglu HÜ, Azevedo L, Plaza S, Curtis MD, Richter A, Weder B, Schulz B, Martinoia E. Disruption of AtMRP4, a guard cell plasma membrane ABCC-type ABC transporter, leads to deregulation of stomatal opening and increased drought susceptibility. Plant J. 2004;39(2):219–36.
Migocka M, Kosieradzka A, Papierniak A, Maciaszczyk-Dziubinska E, Posyniak E, Garbiec A, Filleur S. Two metal-tolerance proteins, MTP1 and MTP4, are involved in Zn homeostasis and Cd sequestration in cucumber cells. J Exp Bot. 2015;66(3):1001–15.
Migocka M, Papierniak A, Kosieradzka A, Posyniak E, Maciaszczyk-Dziubinska E, Biskup R, Garbiec A, Marchewka T. Cucumber metal transport protein CsMTP9 is a plasma membrane H+-coupled antiporter involved in the Mn2+ and Cd2+ efflux from root cells. Plant J. 2015;84(6):1045–58.
Wang J, Li Y, Zhang Y, Chai T. Molecular cloning and characterization of a Brassica juncea yellow stripe-like gene, BjYSL7, whose overexpression increases heavy metal tolerance of tobacco. Plant Cell Rep. 2013;32(5):651–62.
Kim DY, Bovet L, Maeshima M, Martinoia E, Lee Y. The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2007;50(2):207–18.
Migocka M, Papierniak A, Rajsz A. Cucumber PDR8/ABCG36 and PDR12/ABCG40 plasma membrane proteins and their up-regulation under abiotic stresses. Biol Plant. 2017;61(1):115–26.
Lee M, Lee K, Lee J, Noh EW, Lee Y. AtPDR12 contributes to lead resistance in Arabidopsis. Plant Physiol. 2005;138(2):827–36.
Feng S, Tan J, Zhang Y, Liang S, Xiang S, Wang H, Chai T. Isolation and characterization of a novel cadmium-regulated yellow stripe-like transporter (SnYSL3) in Solanum nigrum. Plant Cell Rep. 2017;36(2):281–96.
Durrett TP, Gassmann W, Rogers EE. The FRD3-mediated efflux of citrate into the root vasculature is necessary for efficient iron translocation. Plant Physiol. 2007;144(1):197–205.
Yamaguchi H, Fukuoka H, Arao T, Ohyama A, Nunome T, Miyatake K, Negoro S. Gene expression analysis in cadmium-stressed roots of a low cadmium-accumulating solanaceous plant, Solanum torvum. J Exp Bot. 2010;61(2):423–37.
Milner MJ, Seamon J, Craft E, Kochian LV. Transport properties of members of the ZIP family in plants and their role in Zn and Mn homeostasis. J Exp Bot. 2013;64(1):369–81.
Krzesłowska M. The cell wall in plant cell response to trace metals: polysaccharide remodeling and its role in defense strategy. Acta Physiol Plant. 2011;33(1):35–51.
Xiong J, An LY, Lu H, Zhu C. Exogenous nitric oxide enhances cadmium tolerance of rice by increasing pectin and hemicellulose contents in root cell wall. Planta. 2009;230(4):755–65.
Douchiche O, Driouich A, Morvan C. Spatial regulation of cell-wall structure in response to heavy metal stress: cadmium-induced alteration of the methyl-esterification pattern of homogalacturonans. Ann Bot. 2010;105(3):481–91.
Yang XY, Zeng ZH, Yan JY, Fan W, Bian HW, Zhu MY, Yang JL, Zheng SJ. Association of specific pectin methylesterases with Al-induced root elongation inhibition in rice. Physiol Plant. 2013;148(4):502–11.
Zhu XF, Shi YZ, Lei GJ, Fry SC, Zhang BC, Zhou YH, Braam J, Jiang T, Xu XY, Mao CZ. XTH31, encoding an in vitro XEH/XET-active enzyme, regulates aluminum sensitivity by modulating in vivo XET action, cell wall xyloglucan content, and aluminum binding capacity in Arabidopsis. Plant Cell. 2012;24(11):4731–47.
Degenhardt B, Gimmler H. Cell wall adaptations to multiple environmental stresses in maize roots. J Exp Bot. 2000;51(344):595–603.
Kováčik J, Klejdus B. Dynamics of phenolic acids and lignin accumulation in metal-treated Matricaria chamomilla roots. Plant Cell Rep. 2008;27(3):605–15.
Martinka M, Dolan L, Pernas M, Abe J, Lux A. Endodermal cell–cell contact is required for the spatial control of Casparian band development in Arabidopsis thaliana. Ann Bot. 2012;110(2):361–71.
Koen E, Besson-Bard A, Duc C, Astier J, Gravot A, Richaud P, Lamotte O, Boucherez J, Gaymard F, Wendehenne D. Arabidopsis thaliana nicotianamine synthase 4 is required for proper response to iron deficiency and to cadmium exposure. Plant Sci. 2013;209:1–11.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
This work was in part supported by grants from the Natural Science Foundation of China (Nos. 31671599 and 31370515) and Natural Science Foundation of Anhui Province of China (1708085MC55).
Availability of data and materials
The RNA sequence dataset supporting the results of this article is available in the repository of NCBI Sequence Read Archive (SRA) with the GenBank accession No.: (SRX3808415–3808422).
Ethics approval and consent to participate
Consent for publication
All the authors have signed the consent form.
The authors declare that they have no competing interests for this research.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Overview of raw and clean reads in Cd-free (CK) and Cd-treated (Cd) for two peanut cultivars. (XLSX 10 kb)
Table S2. Mapping results of clean reads against the peanut reference sequence. (XLSX 11 kb)
Table S3. Summary of novel transcripts. (XLSX 9 kb)
Table S4. Gene expression analysis for all sample libraries. (XLSX 5503 kb)
Figure S1. The gene expression distribution (FPKM) in eight samples. X axis represents the sample name. Y axis represents the gene amount. The dark color means the moderate and high expression level which FPKM value ≥10, while the light color means the low expression level which FPKM value < 10. (DOCX 282 kb)
Table S5. List of differentially expressed genes in four groups based on pairwise comparisons. (XLSX 27764 kb)
Figure S2. Gene ontology classification of DEGs identified. The enriched biological process, cellular component and molecular function GO terms of Cd-responsive DEGs in two peanut cultivars (a) and DEGs between Fenghua and Silihong under CK and Cd-treated conditions (b). The x-axis represents the GO term; the y-axis denotes the number of genes. (DOCX 314 kb)
Table S6. GO enrichment analysis of all DEGs. The GO terms with corrected p-value ≤0.01 were considered significantly enriched among the DEGs. (XLSX 15 kb)
Table S7. Overview of all 133 KEGG pathways for DEGs in four comparisons. (XLSX 33 kb)
Table S8. List of enriched pathways for DEGs in four comparisons. Pathways with a corrected P-value ≤ 0.01 were regard as significantly enriched pathways in DEGs. (XLSX 12 kb)
Table S9. The analysis of DEGs in FCd/SCd with an FPKM value ≥10 in at least one of four cDNA libraries and |log2 (FCd/SCd)| ≥ 3. (XLSX 41 kb)
Table S10. The DEGs involved in transport activities in two peanut cultivars under Cd-free and Cd-treated. Items highlighted in red mean up-regulated, while in green mean down-regulated. (XLSX 59 kb)
Table S11. The seven expression patterns of transport-related DEGs in two peanut cultivars under Cd-free and Cd-treated. Items highlighted in red mean up-regulated, while in green mean down-regulated. (XLSX 28 kb)
Table S12. The DEGs involved in cell wall metabolism and Cd chelating in two peanut cultivars. Items highlighted in red mean up-regulated, while in green mean down-regulated. (XLSX 71 kb)
Table S13. The primers used in RT-qPCR analysis. (XLSX 13 kb)
About this article
Cite this article
Yu, R., Ma, Y., Li, Y. et al. Comparative transcriptome analysis revealed key factors for differential cadmium transport and retention in roots of two contrasting peanut cultivars. BMC Genomics 19, 938 (2018). https://doi.org/10.1186/s12864-018-5304-7
- Gene expression