Expanding the genetic variation of Brassica juncea by introgression of the Brassica rapa genome

Abstract

Brassica juncea is an important vegetable and oil crop cultivated worldwide. To increase its genetic variation, we introgressed the A genome of Brassica rapa into B. juncea. We used three each of heading and semi-heading B. juncea accessions as recipient parents and a B. rapa line, B9008, as the donor parent. We obtained 101 BC₁S₁ lines in total with expanded phenotypic variations such as leafy head shapes. We developed 132 single-nucleotide polymorphism (SNP) markers that could distinguish the A genome of B. juncea from the B. rapa genome, and tracked the introgression of B. rapa segments in the new B. juncea germplasm. On average, 59.2% of the B. juncea A genome in the B. juncea introgression lines was covered by the donor segments. We also identified three markers whose donor genotype frequencies were significantly lower than the theoretical value, suggesting strong selection of the recipient genotype during the introgression process. We provide an effective strategy to evaluate the diversity of the new germplasm based on the combination of parental resequencing data and marker genotyping results. Further genetic analysis of 1642 SNPs showed that the genetic diversity of the new B. juncea germplasm with the introgressed B. rapa genome was significantly increased. This study illustrates the potential for expanding the genetic diversity of B. juncea through the introgression of the B. rapa genome.

References

[1] Nagaharu U . Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization . Jpn J Bot. 1935; 7: 389-452.

[2] Chen S, Wan Z, Nelson MN et al. Evidence from genome-wide simple sequence repeat markers for a polyphyletic origin and secondary centers of genetic diversity of Brassica juncea in China and India. J Hered. 2013; 104: 416-27.

[3] YAO QL, Chen FB, Fang P et al. Genetic diversity of Chinese vegetable mustard (Brassica juncea Coss) landraces based on SSR data . Biochem Syst Ecol. 2012; 45: 41-8.

[4] Chen FB, Liu HF, Yao QL et al. Evolution of mustard (Brassica juncea Coss) subspecies in China: evidence from the chalcone synthase gene . Genet Mol Res. 2016; 15: 2.

[5] Gao LW, Lyu SW, Tang J et al. Genome-wide analysis of auxin transport genes identifies the hormone responsive patterns associated with leafy head formation in Chinese cabbage. Sci Rep. 2017; 7: 42229.

[6] Yu J, Gao L, Liu W et al. Transcription coactivator ANGUSTIFOLIA3 (AN3) regulates leafy head formation in Chinese cabbage. Front Plant Sci. 2019; 10: 520.

[7] Cheng F, Sun R, Hou X et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea . Nat Genet. 2016; 48: 1218-24.

[8] Payal B, Shashi B, Banga SS et al. Heterosis as investigated in terms of polyploidy and genetic diversity using designed Brassica juncea amphiploid and its progenitor diploid species. PLoS One. 2012; 7: 2.

[9] Yadava SK, Arumugam N, Mukhopadhyay A et al. QTL mapping of yield-associated traits in Brassica juncea: meta-analysis and epistatic interactions using two different crosses between east European and Indian gene pool lines. Theor Appl Genet. 2012; 125: 1553-64.

[10] Li Q, Mei J, Zhang Y et al. A large-scale introgression of genomic components of Brassica rapa into B. napus by the bridge of hexaploid derived from hybridization between B. napus and B. oleracea . Theor Appl Genet. 2013; 126: 22073-80.

[11] Hasan MJ, Rahman H . Resynthesis of Brassica juncea for resistance to Plasmodiophora brassicae pathotype 3 . Breed Sci. 2018; 68: 385-91.

[12] Chen ZJ . Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu Rev Plant Biol. 2007; 58: 377-406.

[13] Chauhan JS, Singh KH, Singh VV et al. Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian J Agric Sci. 2011; 81: 1093-109.

[14] Kumar S, Hilario E, Deng CH et al. Turbocharging introgression breeding of perennial fruit crops: a case study on apple. Hortic Res. 2020; 7: 1-7.

[15] Bennett RA, Seguin-Swartz G, Rahman H . Broadening genetic diversity in canola using the C-genome species Brassica oleracea L. Crop Sci. 2012; 52: 2030-2039.

[16] Chen S, Nelson MN, Chèvre AM et al. Trigenomic bridges for Brassica improvement. Crit Rev Plant Sci. 2011; 30: 524-47.

[17] Mason AS, Batley J . Creating new interspecific hybrid and polyploid crops. Trends Biotechnol. 2015; 33: 436-41.

[18] Karim MM, Siddika A, Tonu N et al. Production of high yield short duration Brassica napus by interspecific hybridization between B. oleracea and B. rapa . Breed Sci. 2014; 63: 495-502.

[19] Chatterjee D, Banga S, Gupta M et al. Resynthesis of Brassica napus through hybridization between B. juncea and B. carinata . Theor Appl Genet. 2016; 129: 977-90.

[20] Narasimhulu SB, Kirti PB, Prakash S et al. Resynthesis of Brassica carinata by protoplast fusion and recovery of a novel cytoplasmic hybrid. Plant Cell Rep. 1992; 11: 428-32.

[21] Mei J, Shao C, Yang R et al. Introgression and pyramiding of genetic loci from wild Brassica oleracea into B. napus for improving Sclerotinia resistance of rapeseed . Theor Appl Genet. 2020; 133: 1313-9.

[22] Guo N, Cheng F, Wu J et al. Anthocyanin biosynthetic genes in Brassica rapa . BMC Genomics. 2011; 15: 135-42.

[23] Su TB, Wang W, Li P et al. Genomic variation map provides insights into the genetic basis of spring Chinese cabbage (Brassica rapa ssp. pekinensis) . Plant Commun. 2018; 11: 1360-76.

[24] Ramesh P, Mallikarjuna G, Sameena S et al. Advancements in molecular marker technologies and their applications in diversity studies. J Biosci. 2020; 45: 123.

[25] Navabi ZK, Stead KE, Pires JC et al. Analysis of B-genome chromosome introgression in interspecific hybrids of Brassica napus × B. carinata . Genetics. 2011; 187: 659-73.

[26] Badaeva ED, Amosova AV, Goncharov NP et al. A set of cytogenetic markers allows the precise identification of all A-genome chromosomes in diploid and polyploid wheat. Cytogenet Genome Res. 2015; 146: 71-9.

[27] Megyeri M, Mikó P, Farkas A et al. Cytomolecular discrimination of the A^m chromosomes of Triticum monococcum and the A chromosomes of Triticum aestivum using microsatellite DNA repeats . J Appl Genet. 2017; 58: 67-70.

[28] Li F, Kitashiba H, Inaba K et al. A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits. DNA Res. 2009; 16: 311-23.

[29] Wang X, Yu K, Li H et al. High-density SNP map construction and QTL identification for the apetalous character in Brassica napus L. Front Plant Sci. 2015; 6: 1164.

[30] Steele KA, Quinton-Tulloch MJ, Amgai RB et al. Accelerating public sector rice breeding with high-density KASP markers derived from whole genome sequencing of indica rice . Mol Breed. 2018; 38: 1-13.

[31] Yang JH, Liu D, Wang X et al. The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat Genet. 2016; 48: 1225-32.

[32] Cui YN, Zhuang M, Wu J et al. Segmental translocation contributed to the origin of the Brassica S-locus. Hortic Plant J. 2020; 6: 167-78.

[33] Gupta CM, Banga SS . Cytogenetic stability and genome size variations in newly developed derived Brassica juncea allopolyploid lines. J Oilseed Brassica. 2016; 5: 118-27.

[34] Zhang K, Wang XW, Cheng F . Plant polyploidy: origin, evolution, and its influence on crop domestication. Hortic Plant J. 2019; 5: 231-9.

[35] Sharma D, Nanjundan J, Singh L et al. Genetic diversity in leafy mustard (Brassica juncea var. rugosa) as revealed by agromorphological traits and SSR markers. Physiol Mol Biol Plants. 2020; 26: 2005-18.

[36] Cai X, Chang L, Zhang T et al. Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa . Genome Biol. 2021; 22: 166.

[37] Wen J, Tu J, Fu TD et al. Improving ovary and embryo culture techniques for efficient resynthesis of Brassica napus from reciprocal crosses between yellow-seeded diploids B. rapa and B. oleracea . Euphytica. 2008; 162: 81-9.

[38] Porebski S, Bailey LG, Baum BR . Modification of CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components. Plant Mol Biol Report. 1997; 15: 8-15.

[39] Chen SF, Zhou YQ, Chen Y et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018; 34: i884-90.

[40] Wang XW, Wang H, Wang J et al. The genome of the mesopolyploid crop species Brassica rapa . Nat Genet. 2011; 43: 1035-9.

[41] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[42] Li H, Durbin R . Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010; 26: 589-95.

[43] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[44] Liu J, Bo L, Feng C et al. A high density linkage map facilitates QTL mapping of flowering time in Brassica rapa . Hortic Plant J. 2016; 2: 217-23.

[45] Zhang L, Cai X, Wu J et al. Improved Brassica rapa reference genome by single-molecule sequencing and chromosome conformation capture technologies. Hortic Res. 2018; 5: 50.

[46] Letunic I, Bork P . Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 2019; 47: W256-9.

[47] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64.