Dataset Information

Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection.

ABSTRACT:

Background

Cadmium (Cd) accumulation in rice followed by transfer to the food chain causes severe health problems in humans. Breeding of low Cd accumulation varieties is one of the most economical ways to solve the problem. However, information on the identity of rice germplasm with low Cd accumulation is limited, particularly in indica, and the genetic basis of Cd accumulation in rice is not well understood.

Results

Screening of 312 diverse rice accessions revealed that the grain Cd concentrations of these rice accessions ranged from 0.12 to 1.23 mg/kg, with 24 accessions less than 0.20 mg/kg. Three of the 24 accessions belong to indica. Japonica accumulated significantly less Cd than indica (p < 0.001), while tropical japonica accumulated significantly less Cd than temperate japonica (p < 0.01). GWAS in all accessions identified 14 QTLs for Cd accumulation, with 7 identified in indica and 7 identified in japonica subpopulations. No common QTL was identified between indica and japonica. The previously identified genes (OsHMA3, OsNRAMP1, and OsNRAMP5) from japonica were colocalized with QTLs identified in japonica instead of indica. Expression analysis of OsNRAMP2, the candidate gene of the novel QTL (qCd3-2) identified in the present study, demonstrated that OsNRAMP2 was mainly induced in the shoots of high Cd accumulation accessions after Cd treatment. Four amino acid differences were found in the open reading frame of OsNRAMP2 between high and low Cd accumulation accessions. The allele from low Cd accumulation accessions significantly increased the Cd sensitivity and accumulation in yeast. Subcellular localization analysis demonstrated OsNRAMP2 expressed in the tonoplast of rice protoplast.

Conclusion

The results suggest that grain Cd concentrations are significantly different among subgroups, with Cd concentrations decreasing from indica to temperate japonica to tropical japonica. However, considerable variations exist within subgroups. The fact that no common QTL was identified between indica and japonica implies that there is a different genetic basis for determining Cd accumulation between indica and japonica, or that some QTLs for Cd accumulation in rice are subspecies-specific. Through further integrated analysis, it is speculated that OsNRAMP2 could be a novel functional gene associated with Cd accumulation in rice.

SUBMITTER: Zhao J

PROVIDER: S-EPMC6249348 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Root system architecture plays a crucial role in nutrient and water absorption during rice production. Genetic improvement of the rice root system requires elucidating its genetic control. Genome-wide association studies (GWASs) have identified genomic regions responsible for rice root phenotypes. However, candidate gene prioritization around the peak region often suffers from low statistical power and resolution. Transcriptomics enables other statistical mappings, such as transcriptome-wide association study (TWAS) and expression GWAS (eGWAS), which improve candidate gene identification by leveraging the natural variation of the expression profiles. To explore the genes responsible for root phenotypes, we conducted GWAS, TWAS, and eGWAS for 12 root phenotypes in 57 rice accessions using 427,751 single nucleotide polymorphisms (SNPs) and the expression profiles of 16,901 genes expressed in the roots. The GWAS identified three significant peaks, of which the most significant peak responsible for seven root phenotypes (crown root length, crown root surface area, number of crown root tips, lateral root length, lateral root surface area, lateral root volume, and number of lateral root tips) was detected at 6,199,732 bp on chromosome 8. In the most significant GWAS peak region, OsENT1 was prioritized as the most plausible candidate gene because its expression profile was strongly negatively correlated with the seven root phenotypes. In addition to OsENT1, OsEXPA31, OsSPL14, OsDEP1, and OsDEC1 were identified as candidate genes responsible for root phenotypes using TWAS. Furthermore, a cis-eGWAS peak SNP was detected for OsDjA6, which showed the eighth strongest association with lateral root volume in the TWAS. The cis-eGWAS peak SNP for OsDjA6 was in strong linkage disequilibrium (LD) with a GWAS peak SNP on the same chromosome for lateral root volume and in perfect LD with another SNP variant in a putative cis-element at the 518 bp upstream of the gene. These candidate genes provide new insights into the molecular breeding of root system architecture.