Genetic Diversity of Cultivated and Wild Radish and Phylogenetic Relationships Among Raphanus and Brassica Species Revealed by the Analysis of Trnk/Matk Sequence
Total Page:16
File Type:pdf, Size:1020Kb
Breeding Science 58 : 15–22 (2008) Genetic diversity of cultivated and wild radish and phylogenetic relationships among Raphanus and Brassica species revealed by the analysis of trnK/matK sequence Na Lü*1), Kyoko Yamane2) and Ohmi Ohnishi1) 1) Laboratory of Crop Evolution, Graduate School of Agriculture, Kyoto University, 1 Nakajo, Mozume, Mukoh, Kyoto 617-0001, Japan 2) Plant Resource Laboratory, Graduate School of Life and Environmental Sciences, Osaka Prefecture University, 1-1 Gakuen, Naka, Sakai, Osaka 599-8531, Japan Nucleotide sequence variations of the trnK/matK region were investigated in cultivated and wild radish spe- cies and related Brassica species. Two insertions/deletions and 9 substitutions were detected among the Raphanus accessions. The value of the nucleotide diversity (π) was found to be higher in cultivated radish (0.00184) than in the wild species (0.00134). Based on the nucleotide diversity, the phylogenetic relation- ships of Raphanus and its related species were inferred by constructing a Neighbor-Joining tree. Raphanus species and Brassica barrelieri formed a sister clade located between the Rapa/Oleracea group and the Nigra group of Brassica. These results were in complete agreement with those obtained by Warwick and Black. The placement of Raphanus species at this position showed the existence of a paraphyletic relationship among the Brassica species. Each of the three varieties of cultivated radish, R. sativus var. sativus (European small radish), var. hortensis (East Asian big radish) and var. niger (black radish), belonged to a different clus- ter of the phylogenetic tree, suggesting the existence of independent multiple origins of these varieties. Based on the phylogenetic tree, problems related to the identification of the wild ancestral species of cultivated rad- ish and original birthplaces of cultivated radish varieties were discussed. Key Words: Raphanus, Brassica, trnK/matK, nucleotide diversity, phylogenetic relationships, origin of culti- vated radish. Introduction wild Raphanus species, is distributed in the Far East, mainly in the coastal areas of Japan, Korea and China (Makino Raphanus (2n = 18), the radish genus includes one cul- 1961). However, most of the radish scientists were not con- tivated species and several wild species. Cultivated radish vinced that the East Asian wild species was an endemic true has been classified into many varieties based on the mor- wild species. Since this wild species is morphologically phology of its edible root and different uses (Kitamura 1958) similar to cultivated radish in East Asia, it was assumed to be as follows: R. sativus var. sativus L. (syn. var. radicula an escaped form of cultivated radish, and it was classified as Pers.) (European small radish), var. hortensis Becker (East R. sativus var. hortensis f. raphanistroides Kitamura Asian big long radish), var. niger Kerner (black Spanish (Kitamura 1958). However, analysis of mitochondrial DNA radish), var. chinensis Gallizioli (Chinese oil radish) and (Yamagishi and Terachi 2001, Yamagishi 2004) did not con- var. caudatus Hooker & Anderson (tail-podded radish or rat- firm this assumption. tail radish). Among the wild species, R. raphanistrum L., The phylogenetic relationships among cultivated radish R. landra DC. and R. martimus Smith are distributed in the and its wild relatives have been studied using molecular coastal areas of the Mediterranean Sea. Both R. raphanistrum markers. For example, Lewis-Jones et al. (1982) who stud- and R. landra were once considered to be candidates for the ied allozyme variations failed to detected any critical dif- wild ancestor of cultivated radish. They cannot be discrimi- ferentiation among the wild species and concluded that nated from each other morphologically at the species level, R. raphanistrum and R. landra could not be identified as dis- and some authors considered that R. landra was a subspecies tinct species. However, differentiation between wild radish of R. raphanistrum (Munoz and Bermejo 1978). On the species from the Mediterranean area and the Far East was other hand, R. sativus var. raphanistroides Makino, another difficult because no sample of East Asian wild radish was included in their study. A RAPD survey on Raphanus Communicated by T. Sato (Yamagishi et al. 1998) suggested that East Asian wild rad- Received April 9, 2007. Accepted October 23, 2007. ish was more closely related to cultivated radish in East Asia *Corresponding author (e-mail: [email protected]) than to the wild radish in Europe, R. raphanistrum. As for the 16 Lü, Yamane and Ohnishi wild ancestor of cultivated radish, Lewis-Jonas et al. (1982) among wild radish species, including the Far Eastern wild suggested that a variant of the raphanistrum-landra complex species R. sativus var. raphanistroides, which has been from Italy or the East Mediterranean area might be a candi- scarcely investigated simultaneously with European wild date for the wild ancestor. Recently, Yamagishi and Terachi species; 2) to identify the wild ancestral species of cultivated (2003) have analyzed mtDNA variation in Raphanus. They radish. To achieve these objectives, we analyzed 9 acces- showed that haplotypes of cultivated radish were derived sions of three wild species and 9 accessions of different cul- from multiple sources of wild radish and concluded that cul- tivated varieties of R. sativus. Since it had been postulated tivated radish displayed multiple origins. A similar conclu- that a wild Brassica species was involved in the origin of sion was drawn from the studies on mtDNA conducted by cultivated radish, we also included several Brassica species Yamagishi (2004) and on chloroplast DNA by Yamane et al. in the present investigation. This is the first study in which (2005). However, it remains to be determined which wild wild and cultivated radish and related Brassica species was species provided haplotype genomes to cultivated radish. analyzed to elucidate the origin of cultivated radish. Brassica is the genetically closest genus to Raphanus, and cultivated forms of this genus are important vegetables, Materials and Methods such as cabbage, mustard and Chinese cabbage. The basic diploid species include B. oleracea L. (n = 9), B. nigra (L.) Plant materials Koch (n = 8) and B. rapa L. (syn. B. campestris (L.) In the present study, 18 accessions of Raphanus were Clapman, n = 10), and the allotetraploid species, B. napus L. used. They consisted of 9 cultivated radish accessions (n = 19), B. juncea (L.) Czern (n = 18) and B. carinata Braun (R. sativus), 2 accessions of East Asian wild radish, 3 of (n = 17) were derived by interspecific hybridization among R. raphanistrum, and 4 of R. landra. Seven accessions of these basic groups (U 1935). The interspecific relationships Brassica (B. barrelieri, B. rapa (syn. B. campestris), of Brassica species including Raphanus have been studied B. oleracea, B. nigra, B. juncea, B. napus and B. carinata) based on chloroplast DNA diversity (Warwick and Black and one accession of Sinapis alba were also used to analyze 1991, 1993, 1997 and Warwick et al. 1992). The results sug- the phylogenetic relationships of Raphanus and Brassica gested that Brassica could be divided into two lineages (Table 1). The Isatis tinctoria L. species was used as an out- (Rapa/Oleracea lineage and Nigra lineage), and that the group in our study, because I. tinctoria formed a sister group Raphanus species were located in the Rapa/Oleracea lineage. with the group that included Brassica in the phylogenetic On the other hand, the results from nuclear DNA analysis in- tree of Beilstein et al. (2006). All the R. landra, Brassica, dicated that Raphanus was more closely related to the Nigra S. alba and I. tinctoria samples were lines stored at the gene lineage (RFLP: Song et al. 1990, 5S rRNA sequence: Yang bank of the Institute of Plant Genetics and Crop Plant et al. 1998). When using mtDNA, two contradictory conclu- Research (IPK), Gatersleben, Germany (Table 1). The other sions have been obtained. Palmer and Herbon (1988) sug- accessions were collected by one of us (O. Ohnishi) during gested that Raphanus was located in the Rapa/Oleracea lin- expeditions conducted from 1990 to 2005 (Table 1). The eage, whereas Pradhan et al. (1992) indicated that Raphanus samples of wild radish populations, consisting of one was closer to B. nigra. As early as in 1990, Song et al. hundred to one thousand siliques, were originally collected (1990) suggested that R. sativus might have been derived from natural populations. For the cultivated species, mass from species that belonged to two different lineages by in- seed samples were bought at local markets. trogression or hybridization. However, they could not iden- tity the ancestral wild species because no wild species was Polymerase chain reaction and sequencing utilized in their study. The plants were grown in a greenhouse and total DNA Chloroplast DNA, which is characterized by a maternal was isolated from young leaf tissues of an arbitrarily chosen in- inheritance, the absence of recombination and a slow evolu- dividual for each accession, according to the plant DNAzol tionary rate of base substitutions, has been used to analyze Reagent protocol (Gibco BRL, Grand Island, NY, USA). all the types of plant species. One of the chloroplast genes, Using the two universal primers trnK/3914F and trnK/2R matK exhibits a higher nucleotide substitution rate than oth- (Johnson and Soltis 1994), the trnK/matK region was ampli- er chloroplast genes (Olmstead and Palmer 1994), and its fied for all the plants described above. PCR amplifications substitution rate is 3 times higher than that of the rbcL gene were performed in a final reaction solution of 50 µl contain- (Johnson and Soltis 1995). It is considered to be a useful ing 50 ng DNA template, 10 mM Tris-HCl (pH 8.3), 20 mM gene to resolve phylogenetic and evolutionary problems in a KCl, 1.5 mM MgCl2, 100 mM each of dATP, dCTP, dGTP, wide range of taxonomic levels, especially at the closely re- dTTP, 0.2 mM of each primer and 2.5 units of Taq DNA lated taxonomic levels, such as inter- and intra- species polymerase (Takara Shuzo Co.