Mapping and Identifying Candidate Genes Enabling Cadmium Accumulation in Brassica napus Revealed by Combined BSA-Seq and RNA-Seq Analysis

• Published in: International journal of molecular sciences 2023-06, Vol.24 (12), p.10163
• Main Authors: Wang, Huadong, Liu, Jiajia, Huang, Juan, Xiao, Qing, Hayward, Alice, Li, Fuyan, Gong, Yingying, Liu, Qian, Ma, Miao, Fu, Donghui, Xiao, Meili
• Format: Article
• Language: English
• Subjects: Accumulation; Biosynthesis; Brassica napus; Brassica napus - genetics; bulk segregation analysis; Cadmium; cadmium accumulation; Chromosomes; Enrichment; Flowers & plants; Gene mapping; Genes; Genomes; Glutathione transferase; Heavy metals; Humans; Inductively coupled plasma mass spectrometry; Kinases; Leaves; Mapping; Mass spectrometry; Mass spectroscopy; Metabolism; Molecular modelling; Morphology; Phosphoglycerate kinase; Plant Breeding; Population genetics; Quantitative Trait Loci; Rape plants; Rapeseed; Senescence; Sequence Analysis, RNA; Single-nucleotide polymorphism; Soil contamination; Soil pollution; Soil remediation; transcriptome; whole genome resequencing
• ISSN: 1422-0067; 1661-6596; 1422-0067
• DOI: 10.3390/ijms241210163

Rapeseed has the ability to absorb cadmium in the roots and transfer it to aboveground organs, making it a potential species for remediating soil cadmium (Cd) pollution. However, the genetic and molecular mechanisms underlying this phenomenon in rapeseed are still unclear. In this study, a 'cadmium-enriched' parent, 'P1', with high cadmium transport and accumulation in the shoot (cadmium root: shoot transfer ratio of 153.75%), and a low-cadmium-accumulation parent, 'P2', (with a cadmium transfer ratio of 48.72%) were assessed for Cd concentration using inductively coupled plasma mass spectrometry (ICP-MS). An F genetic population was constructed by crossing 'P1' with 'P2' to map QTL intervals and underlying genes associated with cadmium enrichment. Fifty extremely cadmium-enriched F individuals and fifty extremely low-accumulation F individuals were selected based on cadmium content and cadmium transfer ratio and used for bulk segregant analysis (BSA) in combination with whole genome resequencing. This generated a total of 3,660,999 SNPs and 787,034 InDels between these two segregated phenotypic groups. Based on the delta SNP index (the difference in SNP frequency between the two bulked pools), nine candidate Quantitative trait loci (QTLs) from five chromosomes were identified, and four intervals were validated. RNA sequencing of 'P1' and 'P2' in response to cadmium was also performed and identified 3502 differentially expressed genes (DEGs) between 'P1' and 'P2' under Cd treatment. Finally, 32 candidate DEGs were identified within 9 significant mapping intervals, including genes encoding a glutathione S-transferase (GST), a molecular chaperone (DnaJ), and a phosphoglycerate kinase (PGK), among others. These genes are strong candidates for playing an active role in helping rapeseed cope with cadmium stress. Therefore, this study not only sheds new light on the molecular mechanisms of Cd accumulation in rapeseed but could also be useful for rapeseed breeding programs targeting this trait.