Chinese Root-type Mustard Provides Phylogenomic Insights into the Evolution of the Multi-use Diversified Allopolyploid Brassica Juncea.
Jinghua Yang,Changtong Zhang,Ning Zhao,Lili Zhang,Zhongyuan Hu,Sheng Chen,Mingfang Zhang
DOI: https://doi.org/10.1016/j.molp.2017.11.007
IF: 27.5
2018-01-01
Molecular Plant
Abstract:Allopolyploid mustard (Brassica juncea L.) is an important part of the U's triangle model of Brassica comprising vegetable, oilseed, and condiment crops. Vegetable B. juncea varieties are widely distributed in China and have been cultivated for a long period of 6000–7000 years (Yang et al., 2016Yang J.H. Liu D.Y. Wang X.W. Ji C.M. Cheng F. Liu B.N. Hu Z.Y. Chen S. Pental D. Ju Y.H. et al.The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection.Nat. Genet. 2016; 48: 1225-1232Crossref PubMed Scopus (324) Google Scholar). Tremendous phenotypic variation can be observed in the vegetable varieties of B. juncea, which include leaf-type, stem-type, root-type, and stalk-type varieties according to their edible organs (Qi et al., 2007Qi X.H. Zhang M.F. Yang J.H. Molecular phylogeny of Chinese vegetable mustard (Brassica juncea) based on the internal transcribed spacers (ITS) of nuclear ribosomal DNA.Genet. Resour. Crop Evol. 2007; 54: 1709-1716Crossref Scopus (17) Google Scholar). Oilseed B. juncea, as one of three major sources of edible oil in Brassica, is mainly distributed in the Indian subcontinent and Northwest China (Chen et al., 2013Chen S. Wan Z.J. Nelson M.N. Chauhan J.S. Redden R. Burton W.A. Lin P. Salisbury P.A. Fu T.D. Cowling W.A. 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-427Crossref PubMed Scopus (45) Google Scholar, Yang et al., 2016Yang J.H. Liu D.Y. Wang X.W. Ji C.M. Cheng F. Liu B.N. Hu Z.Y. Chen S. Pental D. Ju Y.H. et al.The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection.Nat. Genet. 2016; 48: 1225-1232Crossref PubMed Scopus (324) Google Scholar). It is also used for canola crops in Canada and Australia (Woods et al., 1991Woods D.L. Capcara J.J. Downey R.K. The potential of mustard (Brassica juncea (L) Coss) as an edible oil crop on the Canadian prairies.Can. J. Plant Sci. 1991; 71: 195-198Crossref Google Scholar, Burton et al., 2004Burton W.A. Ripley V.L. Potts D.A. Salisbury P.A. Assessment of genetic diversity in selected breeding lines and cultivars of canola quality Brassica juncea and their implications for canola breeding.Euphytica. 2004; 136: 181-192Crossref Scopus (47) Google Scholar), and condiment crops in Europe, China, and other regions (Vaughan and Hemingway, 1959Vaughan J.G. Hemingway J.S. The utilization of mustards.Econ. Bot. 1959; 13: 196-204Crossref Scopus (22) Google Scholar). About half a century ago, Vavilov proposed that Afghanistan and its adjoining regions were the primary origin center of B. juncea and that Asia Minor, central/western China, and eastern India were secondary centers of diversity (Vavilov, 1949Vavilov N.I. The origin, variation, immunity and breeding of cultivated plants.Chromica Botanica. 1949; 13: 1-364Google Scholar). Some subsequent studies proposed that B. juncea originated from several independent hybridizations based on cytogenetic and biochemical evidence (Olsson, 1960Olsson G. Species crosses within the genus Brassica. 1. Artificial Brassica juncea Coss.Hereditas. 1960; 46: 171-223Crossref Scopus (49) Google Scholar, Vaughan et al., 1963Vaughan J.G. Hemingway J.S. Schofield H.J. Contributions to a study of variation in Brassica juncea Coss. & Czern.Bot. J. Linn. Soc. 1963; 58: 435-447Crossref Scopus (37) Google Scholar, Prakash and Hinata, 1980Prakash S. Hinata K. Taxonomy, cytogenetics and origin of crop Brassicas, a review.Opera Botanica. 1980; 55: 1-57Google Scholar). Recently, simple sequence repeat (SSR) marker analysis of oilseed B. juncea suggested a polyphyletic origin and secondary centers of genetic diversity in China and India (Chen et al., 2013Chen S. Wan Z.J. Nelson M.N. Chauhan J.S. Redden R. Burton W.A. Lin P. Salisbury P.A. Fu T.D. Cowling W.A. 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-427Crossref PubMed Scopus (45) Google Scholar). However, there is also another view that China is the primary origin or center of B. juncea, which is supported by the discovery of wild B. rapa (Chinese variety) and B. nigra in Northwest China, the abundance of B. juncea varieties in China and archaeological evidence of a long history of cultivation (Yang et al., 2016Yang J.H. Liu D.Y. Wang X.W. Ji C.M. Cheng F. Liu B.N. Hu Z.Y. Chen S. Pental D. Ju Y.H. et al.The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection.Nat. Genet. 2016; 48: 1225-1232Crossref PubMed Scopus (324) Google Scholar). According to genomic analysis based on the re-sequencing of A-subgenomes of Brassica species, we suggested a monophyletic origin and evolution into vegetable and oilseed varieties for B. juncea (Yang et al., 2016Yang J.H. Liu D.Y. Wang X.W. Ji C.M. Cheng F. Liu B.N. Hu Z.Y. Chen S. Pental D. Ju Y.H. et al.The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection.Nat. Genet. 2016; 48: 1225-1232Crossref PubMed Scopus (324) Google Scholar). However, the phylogenetic evolutionary history of B. juncea crops has been debated for decades. A population of 109 B. juncea accessions collected from China, India, Australia, and European countries were sorted into vegetable and oilseed groups based on their agricultural usages (Supplemental Table 1). Using a specific-locus amplified fragment sequencing (SLAF-seq) strategy, all 109 accessions were sequenced, and a total of 891 289 single nucleotide polymorphisms (SNPs) were called from each accession (Supplemental Tables 2–4). Finally, we identified 170 181 SNPs after filtration and used them to investigate the genomic evolution of the population of 109 B. juncea accessions (Supplemental Figure 1 and Supplemental Table 4). A neighbor-joining phylogenetic tree was constructed, in which the 109 accessions were mainly clustered into four groups, including two vegetable groups and two oilseed groups (Figure 1A and Supplemental Figure 2). The accessions with varied leaves (28 accessions) or expanded stems (27 accessions) as edible organs from China formed one group, designated as Vegetable (Leaf/Stem) hereafter. A total of 15 root-type B. juncea accessions collected from China were clustered into an independent group designated as Vegetable (Root). The oilseed B. juncea accessions from Australia (9 accessions), China (7 accessions), and Europe (6 accessions) formed a group designated as Oil (A,C,E), whereas all the oilseed accessions from India (14 accessions) were distinctly clustered into a separate group named Oil (IN) (Figure 1A and Supplemental Figure 2). Accessions of the Oil (A,C,E) group were collected from Australia, Northwest China, and European countries and had brown or yellow seeds, whereas the Indian oilseeds were light or dark brown in color (Figure 1D). We used principal component analysis to assess the genetic variance. The first three components explained 30.14% of the total genetic variation, of which the first two components contributed 15.27% (PC1) and 9.21% (PC2) (Figure 1B and Supplemental Figure 3). Four main groups were identified according to the first three components: Vegetable (Leaf/Stem), Vegetable (Root), Oil (IN), and Oil (A,C,E) (Figure 1B). Population structure analysis was used to estimate the individual ancestry and admixture proportions, assuming the existence of certain populations. By analyzing the number of clusters (K value), we found clear evolutionary divergence among the four groups at K = 4 (Figure 1C and Supplemental Figures 4 and 5). When K = 4, we observed a few introgressions in Oil (IN) from Oil (A,C,E), in Vegetable (Root) from Vegetable (Leaf/Stem) and Oil (A,C,E), and in Vegetable (Leaf/Stem) from Oil (IN) and Vegetable (Root) (Figure 1C). We then analyzed the SNP variations in the B-subgenomes of vegetable and oilseed B. juncea varieties. The results showed that the number of polymorphic SNPs was greater than the number of fixed SNPs in the B-subgenomes (Supplemental Table 5 and Supplemental Figure 6), of which this pattern was similar to that in the A-subgenomes, supporting a monophyletic origin for B. juncea (Yang et al., 2016Yang J.H. Liu D.Y. Wang X.W. Ji C.M. Cheng F. Liu B.N. Hu Z.Y. Chen S. Pental D. Ju Y.H. et al.The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection.Nat. Genet. 2016; 48: 1225-1232Crossref PubMed Scopus (324) Google Scholar). Considering the earlier diversification of the A-subgenomes of Vegetable (Root) varieties identified in genomic analysis of the A-subgenomes of 17 vegetable and oilseed varieties and their ancestor B. rapa, we propose that the Vegetable (Root) group diversified first, the Oil (A,C,E) group diversified second, and the Oil (IN) group diversified last from Vegetable (Leaf/Stem) B. juncea (Figure 1D). We also calculated the divergence times for these four groups of B. juncea. The Vegetable (Root) group diverged 0.067 Mya, followed by the Oil (A,C,E) group 0.065 Mya, and the Oil (IN) and Vegetable (Leaf/Stem) groups 0.032 Mya (Figure 1D and Supplemental Figure 7). In this study, we found that Vegetable (Root) accessions diversified at an earlier stage before the advent of oilseed and other vegetable mustard crops, implying that China was the primary origin and diversity center of B. juncea. The much earlier diversification and close relationship between the Vegetable (Root) and Oil (A,C,E) groups also supports the hypothesis of China as the primary diversity center. We found that the oilseed B. juncea accessions were distinctly clustered into two groups, Oil (A,C,E) and Oil (IN), which was consistent with the results of SSR marker analysis (Chen et al., 2013Chen S. Wan Z.J. Nelson M.N. Chauhan J.S. Redden R. Burton W.A. Lin P. Salisbury P.A. Fu T.D. Cowling W.A. 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-427Crossref PubMed Scopus (45) Google Scholar). The Vegetable (Leaf/Stem) group diverged from Oil (IN) group, which are widely cultivated in China and India, respectively, but much later than the Vegetable (Root) and Oil (A,C,E) groups. These results suggest that the Vegetable (Leaf/Stem) and Oil (IN) groups were newly formed during domestication, indicating the secondary diversification of these two groups in China and India. In conclusion, we propose that China was the primary origin and diversity center of B. juncea, and that more than one migration event occurred from the primary to the secondary centers of domestication in both China and India. These findings provide a novel perspective on the evolutionary trajectory of this multiple agricultural use of B. juncea and will shed new light on functional gene mining by repositioning the phylogenetic role of root-type landrace in B. juncea. Recently, we reported a new class of Turnip mosaic virus (TuMV) resistance gene, eukaryotic translation initiation factor 2B-beta (eIF2Bβ), which is isolated from root-type B. juncea (Shopan et al., 2017Shopan J. Mou H.P. Zhang L.L. Zhang C.T. Ma W.W. Walsh J.A. Hu Z.Y. Yang J.H. Zhang M.F. Eukaryotic translation initiation factor 2B-beta (eIF2Bβ), a new class of plant virus resistance gene.Plant J. 2017; 90: 929-940Crossref PubMed Scopus (44) Google Scholar). This work is supported by grants from the Ear-marked Special Fund from the Ministry of Agriculture of China (Grant # 091721301354051001).