Identification of QTLs for Cadmium Tolerance During Seedling Stage and Validation of qCDSL1 in Rice

Ding Shilin, Liu Chaolei, Shang Lianguang, Yang Shenglong, Zhang Anpeng, Jiang Hongzhen, Ruan Banpu, Fang Guonan, Tian Biao, Ye Guoyou, Guo Longbiao, Qian Qian, Gao Zhenyu

Research Paper

Identification of QTLs for Cadmium Tolerance During Seedling Stage and Validation ofin Rice

Ding Shilin1, 2, #, Liu Chaolei2, #, Shang Lianguang3, Yang Shenglong2, Zhang Anpeng2, Jiang Hongzhen2, Ruan Banpu2, Fang Guonan2, Tian Biao2, Ye Guoyou4, Guo Longbiao2, Qian Qian2, Gao Zhenyu2

(Root Biology Center, College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou 350002, China; State Key Laboratory of Rice Biology, China National Rice Research Institute, Hangzhou 310006, China;Agricultural Genomics Institute, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China; Strategic Innovation Platform, International Rice Research Institute, DAPO Box 7777, Metro Manila, the Philippines; These authors contributed equally to this work)

Cadmium (Cd) is a non-essential toxic metal that is harmful to plants. To investigate the genetic mechanism of Cd tolerance in rice, quantitative trait loci (QTLs) associated with Cd tolerance at the seedling stage were analyzed using a recombinant inbred line (RIL) population derived from a cross between PA64s and 93-11. A total of 36 QTLs associated with shoot length, root length, shoot dry weight, root dry weight and total dry weight were detected in Hangzhou and Lingshui of China. Among them, 15 QTLs were identified under the control condition and 15 QTLs were identified under the Cd stress condition, and 6 QTLs for Cd tolerant coefficient were detected on chromosomes 1, 3, 7 and 9. Theandwere identified in Hangzhou and Lingshui, respectively, and had overlapping intervals on chromosome 1. To further confirm the effects ofand, we developed a chromosome segment substitution line (CSSL), CSSL, in 93-11 background harboring/from PA64s. Compared to 93-11, CSSLhad increased shoot length under the Cd stress condition. These results pave the way for further isolation of those genes controlling Cd tolerance in rice and marker-assistant selection of rice elite varieties with Cd tolerance.

cadmium tolerance; recombinant inbred line;; rice; quantitative trait locus

Cadmium (Cd) is a highly toxic heavy metal element to both plants and humans. In recent years, due to unreasonable industrial discharge, poor treatment of solid waste, sewage irrigation and application of fertilizers containing heavy metal elements, Cd content in the soil has increased sharply and its pollution has become increasingly serious (Clemens et al, 2013). Cd enters human body mainly through the food chain and has a half-life period of 10?35 years in human body (Clemens et al, 2013). Cd is mainly accumulated in the kidneys after entering the human body, which will be harmful to human bones and respiratory system, thus causing a series of diseases (Jarup and Akesson, 2009; Nawrot et al, 2010).

Rice is one of the most important crops in China. Cd pollution not only deteriorates soil quality of paddy rice fields, but also damages rice growth and development. Meanwhile, it can be accumulated in rice grains through absorption and transportation in rice plants, which has become an important source for endangering human health (Yao et al, 2003). Genetic variation in Cd tolerance of rice indicates that it is possible to select cultivars with strong resistance to Cd and low Cd concentration in grains. Therefore, the study of genetic mechanism underlies Cd tolerance in rice is helpful to develop Cd-tolerant rice cultivars.

For complex traits controlled by multiple genes, quantitative trait locus (QTL) mapping has become a powerful means to identify the number, location and effects of genetic factors, and also an important step to understand molecular genetic mechanism for complicated phenotype. So far, many QTLs and genes related to Cd accumulation have been identified and isolated in rice, including(Agrawal et al, 2002),(Agrawal et al, 2003),(Agrawal et al, 2003),(Mukhopadhyay et al, 2004),(Koike et al, 2004),(Nakanishi et al, 2006),(Nakanishi et al, 2006),(Lee et al, 2007),(Kuramata et al, 2009),(Shim et al, 2009),(Shimo et al, 2011),(Oda et al, 2011),(Miyadate et al, 2011),(Uraguchi et al, 2011),(Takahashi et al, 2011),(Ishimaru et al, 2011),(Takahashi et al, 2012),(Yuan et al, 2012),(Ishikawa et al, 2012),(Ramegowda et al, 2013),(Lim et al, 2014),(Yu et al, 2015),(Wang et al, 2015),(Wang et al, 2016),(Das et al, 2017),(Tan et al, 2017),(Luo et al, 2018) and(Yan et al, 2019). Although genetic studies were mainly focused on Cd accumulation in rice, genetic analysis of QTLs related to Cd tolerance in shoots and roots of rice was relatively fewer.

To identify QTLs for Cd toxicity tolerance in rice at the seedling stage, we treated 118 recombinant inbred lines (RILs) derived from a cross between parents 93-11 and PA64s with 50 μmol/L CdCl2under controlled conditions. With a high-density SNP linkage map, a total of 36 QTLs associated with 5 important traits were detected in Hangzhou and Lingshui of China. A novel QTL,for shoot length under Cd stress was validated by a chromosome segment substitution line CSSLin 93-11 background harboringfrom PA64s, and candidate genes were predicted, which will be helpful for further cloning of the QTL and understanding of the genetic basis for rice Cd tolerance.

Table 1. Phenotypic values of parents under different concentrations of Cd.

Data are Mean ± SD (= 6). * and ** represent 5% and 1% significant differences, respectively, according to the Student’s-test.

RESULTS

Phenotypic variance in RIL population of 93-11 and PA64s in Hangzhou and Lingshui

The traits of root length, shoot length, root dry weight, shoot dry weight and total dry weight were used to evaluate the Cd stress tolerance of the parents 93-11 and PA64s, and their RIL population in Hangzhou and Lingshui. To determine the suitable Cd concentration for treatment, gradient Cd concentrations (0, 10, 20, 30, 40, 50, 100 and 150 μmol/L) were applied on the two parents (93-11 and PA64s) in 2018, and 50 μmol/L Cd concentration was selected for Cd treatment (Table 1). The seeds of parents (93-11 and PA64s) and RILs harvested in 2019 were treated with 50 μmol/L CdCl2. Except for root length, root dry weight in Hangzhou and Lingshui, and total dry weight in Lingshui of PA64s, the values of traits were reduced in Cd stress compared to the control for two parents (Table 2). Cd tolerance coefficients of all the measured traits were higher in PA64s than 93-11 both in Hangzhou and Lingshui, indicating that PA64s is relatively tolerant to Cd compared with 93-11.

The mean, range, skewness and kurtosis of each examined trait in RIL population were summarized in Table 2. Except for root dry weight, the phenotypic values of these traits in both locations under the Cd stress were lower than those under the control condition without Cd. And the segregation of these traits was continuously distributed among all the RILs (Figs. S1 and S2). Additionally, transgressive segregation in both directions was observed for all the traits, suggesting that both parents transmitted favorable alleles for each trait. All the traits of 118 RIL lines showed approximately normal distribution, which reflected polygenic segregation and satisfied the request of QTL analysis. In addition, all the tested traits exhibited a significantly positive correlation under the control and Cd stress conditions in Hangzhou and Lingshui, respectively (Table S1).

Table 2. Phenotypic values of recombinant inbred line (RIL) population and their parents under control and Cd stress.

CTC, Cd tolerance coefficient.Data are Mean ± SD (= 6). * and ** indicate 5% and 1% significant differences between the two parents, respectively, according to the Student’s-test.

QTL analysis of five traits associated with cadmium tolerance

To identify QTLs for the five traits associated with cadmium tolerance under the control and Cd stress, a total of 36 putative QTLs were detected on rice chromosomes 1, 2, 3, 4, 5, 7, 9, 10 and 11 (Table 3 and Fig. 1). Among them, 15 QTLs were identified under the control condition with each QTL accounted for 7.8%?18.1% of phenotypic variation, 15 QTLs were detected under the Cd stress condition with each QTL accounted for 0.9%?12.8% of phenotypic variation, and 6 QTLs for Cd tolerance coefficient were also mapped on chromosomes 1, 3, 7 and 9 with each QTL accounted for 6.2%?15.4% of phenotypic variation.

Under the control and Cd stress conditions in Hangzhou and Lingshui, most of QTLs had positive allele coming from PA64s, indicating that most alleles from PA64s increase Cd tolerance. Total of seven clusters of QTLs were found, including two QTLs(and) mapped on the same region of chromosome 1, three QTLs(,and) on the same or overlapping region of chromosome 1, four QTLs (,,and) on the same or overlapping region of chromosome 2, three QTLs(,and) on the same or overlapping region of chromosome 4, three QTLs(,and) on the same or overlapping region of chromosome 5, four QTLs (,,and) and two QTLs(and) on the same or overlapping region of chromosome 11, respectively. These clusters of QTLs were all involved in dry matter accumulation under the control and Cd stress conditions. Onlyandwere detected both in Lingshui and Hangzhou, indicating the QTLs are genetically stable and independent on the two different environments (Table 3 and Fig. 1).

Table 3. Putative QTLs with LOD > 2.0 detected in rice recombinant inbred line population.

Individual QTL is designated with the italicized abbreviation of the character and the chromosome number. When more than one QTL affecting a character is identified on the same chromosome, they are distinguished using decimal numbers.Maximum-likelihood LOD score for the QTL calculated by MultiQTL package.The positive or negative value indicates that allele from 93-11 or PA64s increases the trait score, respectively. CTC, Cd tolerance coefficient; RL, Root length; SL, Shoot length; RDW, Root dry weight; SDW, Shoot dry weight; TDW, Total dry weight; Chr, Chromosome; Add, Additive effect; Var, Variation explained by the putative QTL.

Fig. 1. Chromosomal locations of all QTLs for Cd tolerance in rice recombinant inbred line population at seedling stage.

The genetic distance of marker (cM) is annotated on the left of each chromosome. Chr, Chromosome; CTC, Cd tolerance coefficient; RL, Root length; SL, Shoot length; RDW, Root dry weight; SDW, Shoot dry weight; TDW, Total dry weight.

Validation of qCDSL1 for shoot length under Cd stress and determination of candidate genes

To validate the physiological role ofin Cd stress response, we developed CSSL(Fig. 2-A). This line harbored the PA64s-derivedallele in 93-11 genetic background. To investigate whetherregulates the response to Cd stress, we exposed the 7-day-old plants of 93-11 and CSSLto 0 and 50 μmol/L Cd for 10 d in a hydroponic experiment. The shoot length of CSSLwas significantly shorter than that of 93-11 under the control condition and significantly longer than that of 93-11 under the Cd stress condition (Fig. 2-B and -C). There were 127 annotated genes in the 873.7 kb region for. Besides 13 transposon protein genes, 16 retrotransposon protein genes, 7 hypothetical protein genes and 25 expressed protein genes, 66 functionally annotated genes were included in the region (Table S2). Previous studies showed that thegene plays an important role in response to Cd stress and it is also located in theregion(Yu et al, 2015). Sequencing of thegene revealed six SNPs in the promoter ofbetween 93-11 and PA64s (Fig. 2-D). The expression level ofin CSSLwas significantly lower than that of 93-11 under the control condition and significantly higher than that of 93-11 under the Cd stress condition (Fig. 2-E), which was consisted with previous observation thatexpression increased5-fold in roots when treated with Cd (Yu et al, 2015).

DISCUSSION

Rice tolerance to Cd is a quantitative trait with complex genetic basis (Xue et al, 2009). To study the genetic basis of rice tolerance to Cd stress, it is necessary to select appropriate specific traits. Rice seedlings and seed germination are sensitive to heavy metal stress. It was found that with the increase of Cd concentration, the seedling length, root length, fresh weight and dry weight were all decreased (Yang et al, 2017). Different concentrations of Cd showed inhibition on seedling growth with different degrees (Li et al, 2019). Thus, it is assumed that the suitable Cd level should be selected in order to detect QTLs with large effects on Cd tolerance. Here, when the concentration of Cd was 50 μmol/L, the differences of shoot length, shoot dry weight, root dry weight and total dry weight between the two parents were the most significant (Table 1). Therefore, 50 μmol/L of Cd was chosen to treat RIL seedlings.

The detection of QTLs is affected by population size, threshold value, marker number, heritability, genetic background and experimental conditions. In previous studies, Xue et al (2009) detected 22 QTLs for 6 traits under the control and Cd stress conditions, including shoot height, root length, shoot dry weight, root dry weight, total dry weight and chlorophyll content, and 6 QTLs related to Cd tolerance coefficient at the seedling stage. Among them,for total dry weight were mapped in the region of G249?G164, where thelocated for relative root dry weight in our study. Li et al (2019) recently detected six QTLs for Cd tolerance coefficient on chromosomes 1, 4, 7, 8 and 10, among which thefor tolerance index of root fresh weight was also located in the region overlapped withhere, suggesting a major QTL controlling Cd tolerance on chromosome 7 in rice seedlings. A total of 34 new QTLs detected in the study indicated that Cd stress has different effects on various organs of rice seedlings and the genetic complexity of Cd tolerance in rice.

Up to now, a few studies have been performed on mapping of QTLs for Cd stress tolerance at the rice seedling stage. However, no reports have been published on cloning of QTLs/genes for Cd tolerance at the rice seedling stage. After searching for genes related to Cd metabolism,was found to be located inorregion,was found to be located inandregion, indicating thatandare likely to control the expression of QTLs for these traits under the Cd condition. After annotation analysis of genes within theregion, thegene was found to be associated with Cd stress responses (Yu et al, 2015). Sequencing and expression analyses suggestedwas most probably responsible for regulating tolerance to Cd stress. Further fine mapping and functional confirmation will be conducted with F2population by crossing the CSSL with 93-11. Pyramiding of major QTLs is a powerful strategy in rice breeding. Therefore, QTLs related to high Cd resistance identified here and QTLs for low Cd accumulation in grains can also be used to develop novel elite rice varieties by QTL pyramid.

METHODS

Rice materials

93-11 is anvariety, while PA64s is an-like variety with maternal origin of. The RIL population was generated from an advanced self-fertilization population of 93-11/PA64s F1plants. A total of 118 lines were obtained, and the seeds were harvested at the China National Rice Research Institute in Hangzhou and Lingshui of China with normal conventional field cultivation. The chromosomal segment substitution line, CSSL, was selected from the advanced backcross population (BC4F2) derived from a cross of the recurrent parent 93-11 and the donor parent PA64s (Zhang et al, 2019).

Fig. 2. Validation of.

A, Schematic graph of chromosome 1 of CSSL, the parents 93-11 and PA64s. The white and black bars represent 93-11 and PA64s alleles, respectively. B, Comparison of seedling growth morphology of 93-11 and CSSLunder the control and Cd stress conditions. Scale bar, 2 cm.C, Shoot lengths of 93-11 and CSSLunder the control and Cd stress conditions. Data are Mean ± SD (= 6). * and ** indicate 5% and 1% significant levels compared to 93-11 under the control and Cd stress conditions, respectively, according to the Student’s-test.D, Gene structure and sequence differences of() between PA64s and 93-11.E, Relative expression level ofin shoots of 93-11 and CSSLunder the control and Cd stress conditions. Data are Mean ± SD (= 3). ** indicates 1% significant level compared to 93-11 under the control and Cd stress conditions, respectively, according to the Student’s-test.

Culture conditions

The uniform seeds of the two parents and RILs were surface sterilized in 3% H2O2solution for 10 min, and then rinsed five times with deionized water. The seeds were then soaked in deionized water in the dark at 28 oC for 2 d, and then transferred to a net floating on deionized water for further 5 d in a controlled chamber with a photoperiod of 16 h light/8 h dark. After 5 d, rice seedlings with similar size were selected and transplanted into 40 L plastic containers containing nutrient solution (96 seedlings per container) and cultivated in the hydroponic form. The light/dark temperatures were set at 32 oC/ 28 oC, and relative humidity was kept at 80%. The seedlings were cultured in a half-strength Kimura B nutrient solution (pH 5.4) with the following composition: 90 μmol/L KH2PO4, 270 μmol/L MgSO4, 180 μmol/L (NH4)2SO4, 90 μmol/L KNO3, 180 μmol/L Ca(NO3)2, 3 μmol/L H3BO3, 0.5 μmol/L MnCl2, 1 μmol/L (NH4)6Mo7O24, 0.4 μmol/L ZnSO4and 20 μmol/L Fe3+-EDTA. The solution was renewed every 2 d. At 3 d after transplanting to the basic solution, two treatments were established: control, in which nine plants grew in the nutrient solution without Cd addition; and cadmium stress, CdCl2was added to the solution in equal increments every day, and after 5 d, the final Cd concentration was 50 μmol/L.

Trait measurements

At the 10 d after Cd treatment, shoot length, root length, shoot dry weight, root dry weight and total dry weight of each treatment were measured. Shoot length was measured from the coleoptile node to the tip of the longest leaf, and root length was measured from the coleoptile node to the tip of the longest root. The mean values were calculated from measured 4–6 plants. The sampled plants were separated into roots and shoots, dried at 70 oC for 3 d in an oven, and then weighted root and shoot dry weights. Total dry weight was calculated according to shoot and root dry weights. Cd tolerance coefficient (CTC) of relative root length, relative shoot length, relative shoot dry weight, relative root dry weight and relative total dry weight were calculated using the following formula: CTC = The value in Cd stress treatment / The value in the control.

Data analysis and QTL mapping

Statistical analysis was conducted with the SAS software (version 9.0). By resequencing parents 93-11 and PA64s and 118 RILs, we obtained 2 622 single nucleotide polymorphism (SNP) markers with high quality and polymorphism, and constructed a high-density SNP linkage map. The SNP linkage map covered a total of 1 381.9 cM of the rice genome, with an average linkage map spacing of 0.392 cM (Gao et al, 2013). QTL analysis was performed with the MultiQTL package (www.mutiqtl.com) using the maximum likelihood interval mapping approach for the RILs. QTL was determined with threshold< 0.005 and the threshold of LOD > 2.0 was chosen for claiming a putative QTL. The genetic parameters, additive effects and explained variation of each QTL were also estimated. We followed the rules for the QTL nomenclature by McCouch et al (1997).

Expression analysis at transcript level

To analyze the transcriptional expression of candidate genes, the 7-day-old plants of 93-11 and CSSLwere exposed to 0 or 50 μmol/L CdCl2for 10 d, and then the shoots were excised for RNA extraction. Total RNA was extracted using the RNA Extraction kit (Axygen, New York, America). DNase I- treated RNA (1 μg) was used to synthesize the first-strand cDNAs by using a ReverTra?Ace qPCR-RT kit (TOYOBO, Osaka, Japan). The cDNA products and the 2× SYBR Green PCR Master Mix (TOYOBO, Osaka, Japan) were used for real-time PCR analysis.was used as an internal control. Data were collected in accordance with the ABI PRISM 7900HT Sequence Detector system. Primers forandare listed in Table S3.

ACKNOWLEDGEMENTS

This research was supported by the National Natural Science Foundation of China (Grant No. 31671761) and the Agricultural Science and Technology Innovation Program, Shenzhen Science and Technology Program (Grant No. 2017050414212249).

SUPPLEMENTAL DATA

The following materials are available in the online version of this article at http://www.sciencedirect.com/science/journal/ 16726308; http://www.ricescience.org.

Fig. S1. Frequency distribution of Cd tolerance at the seedling stage in RIL population with seeds from Hangzhou under control and Cd stress conditions.

Fig. S2.Frequency distribution of Cd tolerance at the seedling stage in RIL population with seeds from Lingshui under control and Cd stress conditions.

Table S1. Correlation coefficients between five traits under control and Cd stress conditions.

Table S2. Functionally annotated genes in the 873.7 kb region for.

Table S3. Primers used for qRT-PCR in the study.

Agrawal G K, Rakwal R, Iwahashi H. 2002. Isolation of novel rice (L.) multiple stress responsive MAP kinase gene,, whose mRNA accumulates rapidly in response to environmental cues., 294(5): 1009–1016.

Agrawal G K, Agrawal S K, Shibato J, Iwahashi H, Rakwal R, 2003. Novel rice MAP kinasesandinvolved in encountering diverse environmental stresses and developmental regulation., 300(3): 775–783.

Clemens S, Aarts M G M, Thomine S, Verbruggen N. 2013. Plant science: The key to preventing slow cadmium poisoning., 18(2): 92–99.

Das N, Bhattacharya S, Bhattacharyya S, Maiti M K. 2017. Identification of alternatively spliced transcripts of rice phytochelatin synthase 2 gene, involved in mitigation of cadmium and arsenic stresses., 94: 167–183.

Gao Z Y, Zhao S C, He W M, Guo L B, Peng Y L, Wang J J, Guo X S, Zhang X M, Rao Y C, Zhang C, Dong G J, Zhang F Y, Lu C X, Hu J, Zhou Q, Liu H J, Wu H Y, Xu J, Ni P X, Zeng D L, Liu D H, Tian P, Gong L H, Ye C, Zhang G H, Wang J, Tian F K, Xue D W, Liao Y, Zhu L, Cheng M S, Li J Y, Cheng S H, Zhang G Y, Wang J, Qian Q. 2013. Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences., 110: 14492–14497.

Ishikawa S, Ishimaru Y, Igura M, Kuramata M, Abe T, Senoura T, Hase Y, Arao T, Nishizawa N K, Nakanishi H. 2012. Ion-beam irradiation, gene identification, and marker-assisted breeding in the development of low-cadmium rice., 109: 19166–19171.

Ishimaru Y, Kakei Y, Shimo H, Bashir K, Sato Y, Sato Y, Uozumi N, Nakanishi H, Nishizawa N K. 2011. A rice phenolic efux transporter is essential for solubilizing precipitated apoplasmic iron in the plant stele., 286: 24649–24655.

Jarup L, Akesson A. 2009. Current status of cadmium as an environ- mental health problem., 238(3): 201–208.

Koike S, Inoue H, Mizuno D, Takahashi M, Nakanishi H, Mori S, Nishizawa N K. 2004. OsYSL2 is a rice metal-nicotianamine transporter that is regulated by iron and expressed in the phloem., 39(3): 415–424.

Kuramata M, Masuya S, Takahashi Y, Kitagawa E, Inoue C, Ishikawa S, Youssefan S, Kusano T. 2009. Novel cysteine-rich peptides fromandenhance tolerance to cadmium by limiting its cellular accumulation., 50(1): 106–117.

Lee S, Kim Y Y, Lee Y, An G. 2007. Rice P1B-type heavy-metal ATPase, OsHMA9, is a metal efflux protein., 145(3): 831–842.

Li W X, Ou Yang L J, Wen W, Xiong Y Y, Xu W Q, Peng X S, Chen X R, He X P, Fu J R, Bian J M, Xu J, Zhou D H, He H H, Sun X T, Zhu C L. 2019. Identification of QTL for cadmium tolerance at seedling stage of rice (L.)., 41(1): 19–24. (in Chinese with English abstract)

Lim S D, Hwang J G, Han A R, Park Y C, Lee C, Ok Y S, Jang C S. 2014. Positive regulation of rice RING E3 ligase OsHIR1 in arsenic and cadmium uptakes., 85: 365–379.

Luo J S, Huang J, Zeng D L, Peng J S, Zhang G B, Ma H L, Guan Y, Yi H Y, Fu Y L, Han B, Lin H X, Qian Q, Gong J M. 2018. A defensin-like protein drives cadmium efflux and allocation in rice., 9(1): 645.

McCouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Morishima H, Kinoshita T. 1997. Rice: Report on QTL nomenclature., 14: 11–13.

Miyadate H, Adachi S, Hiraizumi A, Tezuka K, Nakazawa N, Kawamoto T, Katou K, Kodama I, Sakurai K, Takahashi H, Satoh-Nagasawa N, Watanabe A, Fujimura T, Akagi H. 2011. OsHMA3, a P1B-type of ATPase affects root-to-shoot cadmium translocation in rice by mediating efflux into vacuoles., 189(1): 190–199.

Mukhopadhyay A, Vij S, Tyagi A K. 2004. Overexpression of a zinc-finger protein gene from rice confers tolerance to cold, dehydration, and salt stress in transgenic tobacco., 101(16): 6309–6314.

Nakanishi H, Ogawa I, Ishimaru Y, Mori S, Nishizawa N K. 2006. Iron deficiency enhances cadmium uptake and translocation mediated by the Fe transporters OsIRT1 and OSIRT2 in rice., 52(4): 464–469.

Nawrot T S, Staessen J A, Roels H A, Munters E, Cuypers A, Richart T, Ruttens A, Smeets K, Clijsters H, Vangronsveld J. 2010. Cadmium exposure in the population: From health risks to strategies of prevention., 23(5): 769–782.

Oda K, Otani M, Uraguchi S, Akihiro T, Fujiwara T. 2011. Riceis Cd inducible and confers Cd tolerance on yeast., 75(6): 1211–1213.

Ramegowda Y, Venkategowda R, Jagadish P, Govind G, Hanuman- thareddy R R, Makarla U, Guligowda S. 2013. Expression of a rice Zn transporter, OsZIP1, increases Zn concentration in tobacco and finger millet transgenic plants., 7(3): 309–319.

Shim D, Hwang J U, Lee J, Lee S, Choi Y, An G, Martinoia E, Lee Y. 2009. Orthologs of the class A4 heat shock transcription factor HsfA4a confer cadmium tolerance in wheat and rice., 21(12): 4031–4043.

Shimo H, Ishimaru Y, An G, Yamakawa T, Nakanishi H, Nishizawa N K. 2011.(), a novel gene related to cadmium tolerance and accumulation in rice., 62(15): 5727–5734.

Takahashi R, Ishimaru Y, Senoura T, Shimo H, Ishikawa S, Arao T, Nakanishi H, Nishizawa N K. 2011. The OsNRAMP1 iron transporter is involved in Cd accumulation in rice., 62(14): 4843–4850.

Takahashi R, Ishimaru Y, Shimo H, Ogo Y, Senoura T, Nishizawa N K, Nakanishi H. 2012. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice., 35(11): 1948–1957.

Tan M P, Cheng D, Yang Y N, Zhang G Q, Qin M J, Chen J, Chen Y H, Jiang M Y. 2017. Co-expression network analysis of the transcriptomes of rice roots exposed to various cadmium stresses reveals universal cadmium responsive genes., 17(1): 194.

Uraguchi S, Kamiya T, Sakamoto T, Kassai K, Sato Y, NagamuraY, Yoshida A, Kyozuka J, Ishikawa S, Fujiwara T. 2011. Low-affinity cation transporter (OsLCT1) regulates cadmium transport into rice grains., 108: 20959–20964.

Wang C H, Guo W L, Shan Y, Wei P C, David W O. 2016. Reduction of Cd in rice through expression of-like gene fragments., 9(2): 301–304.

Wang F J, Wang M, Liu Z P, Shi Y, Han T Q, Ye Y Y, Gong N, Sun J W, Zhu C. 2015. Different responses of low grain-Cd- accumulating and high grain-Cd-accumulating rice cultivars to Cd stress., 96: 261–269.

Xue D W, Chen M C, Zhang G P. 2009. Mapping of QTLs associated with cadmium tolerance and accumulation during seedling stage in rice (L.)., 165(3): 587–596.

Yan H L, Xu W X, Xie J Y, Gao Y W, Wu L L, Sun L, Feng L, Chen X, Zhang T, Dai C H, Li T, Lin X N, Zhang Z Y, Wang X Q, Li F M, Zhu X Y, Li J J, Li Z C, Chen C Y, Ma M, Zhang H L, He Z Y. 2019. Variation of a major facilitator superfamily gene contributes to differential cadmium accumulation between rice subspecies., 10: 2562.

Yang M, Chen L, Xu Q G, Sun Y L. 2017. Effects of cadmium stress on seed germination and growth characteristic of different rice cultivars., 31(6): 659–663.

Yao H Y, Xu J M, Huang C Y. 2003. Substrate utilization pattern, biomass and activity of microbial communities in a sequence of heavy metal-polluted paddy soils., 115(1/2): 139–148.

Yu C L, Sun C D, Shen C J, Wang S K, Liu F, Liu Y, Chen Y L, Li C Y, Qian Q, Aryal B, Geisler M, Jiang D A, Qi Y H. 2015. The auxin transporter, OsAUX1, is involved in primary root and root hair elongation and in Cd stress responses in rice (L)., 83(5): 818–830.

Yuan L Y, Yang S G, Liu B X, Zhang M, Wu K Q. 2012. Molecular characterization of a rice metal tolerance protein, OsMTP1., 31(1): 67–79.

Zhang B, Shang L G, Ruan B P, Zhang A P, Yang S L, Jiang H Z, Liu C L, Hong K, Lin H, Gao Z Y, Hu J, Zeng D L, Guo L B, Qian Q. 2019. Development of three sets of high throughput genotyped rice chromosome segment substitution lines and QTL mapping for eleven traits., 12(1): 12–33.

Copyright ? 2021, China National Rice Research Institute. Hosting by Elsevier B V

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Peer review under responsibility of China National Rice Research Institute

http://dx.doi.org/10.1016/j.rsci.2020.11.009

12 December 2019;

30 March 2020

s:Gao Zhenyu (gaozhenyu@caas.cn); Qian Qian (qianqian188@hotmail.com)

(Managing Editor: Wu Yawen)