Scientifica / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 952395 |

Bharti Sarin, Juan Pedro Martín, Babeeta Chrungu Kaula, Aparajita Mohanty, "Chloroplast DNA Variations in Wild Brassicas and Their Implication in Breeding and Population Genetics Studies", Scientifica, vol. 2015, Article ID 952395, 7 pages, 2015.

Chloroplast DNA Variations in Wild Brassicas and Their Implication in Breeding and Population Genetics Studies

Academic Editor: Karl-Josef Dietz
Received25 Jun 2015
Accepted22 Jul 2015
Published10 Aug 2015


Evaluation of chloroplast DNA (cpDNA) diversity in wild relatives of crop brassicas is important for characterization of cytoplasm and also for population genetics/phylogeographic analyses. The former is useful for breeding programs involving wide hybridization and synthesis of alloplasmic lines, while the latter is important for formulating conservation strategies. Therefore, PCR-RFLP (Polymerase Chain Reaction-Restriction Fragment Length Polymorphism) technique was applied to study cpDNA diversity in 14 wild brassicas (including 31 accessions) which revealed a total of 219 polymorphic fragments. The combination of polymorphisms obtained by using only two primer pair-restriction enzyme combinations was sufficient to distinguish all 14 wild brassicas. Moreover, 11 primer pairs-restriction enzyme combinations revealed intraspecific polymorphisms in eight wild brassicas (including endemic and endangered species, B. cretica and B. insularis, resp.). Thus, even within a small number of accessions that were screened, intraspecific polymorphisms were observed, which is important for population genetics analyses in wild brassicas and consequently for conservation studies.

1. Introduction

The wild relatives of crop brassicas are repositories of genes conferring resistance to several biotic and abiotic stresses [1] and also source of male sterility-inducing cytoplasm in cultivars [2]. As maternal inheritance of chloroplast and mitochondrial genomes has been observed in Brassica species [3], evaluation of chloroplast genome diversity in wild brassicas can demonstrate the maternal lineage of related species [4]. This is important for breeding programs, because the type of cytoplasm/maternal lineage in brassicas can influence the direction of cross and extent of success achieved in wide hybridization [5, 6]. Also, analysis of chloroplast DNA (cpDNA) variations can reveal genetic relatedness within and between wild and cultivated species [7, 8]. Studies on cpDNA diversity are also important for population genetics and phylogeographic analyses of rare, endemic, and endangered species. Many of the wild relatives (e.g., Brassica insularis and B. cretica) are endemic and/or are endangered species [9, 10] and therefore the population genetics studies of such species are essential for formulating conservation strategies. Therefore, for carrying out such genetic and conservation studies, the first step is to assess chloroplast genome of wild brassicas for intergeneric/interspecific and intraspecific polymorphisms.

PCR-RFLP (Polymerase Chain Reaction-Restriction Fragment Length Polymorphism) is a simple, rapid, and reproducible technique [11] that uses universal primers [12] to amplify chloroplast genome regions followed by digestion with restriction enzymes to reveal fragment length polymorphisms [13]. There are very few studies in brassicas, where PCR-RFLP of chloroplast genome has been analyzed [4, 14, 15]. Cunha et al. [15] used PCR-RFLP of cpDNA to discriminate three diploid cultivars of brassicas. However, this technique did not reveal any intraspecific or interspecific polymorphisms in wild and cultivated B. oleracea members [14]. Yamane et al. [4] used PCR-RFLP technique to detect interspecific polymorphisms in Raphanus sp. which facilitated the understanding of maternal lineage of cultivated radish. Use of simple sequence repeat (SSR) markers of cpDNA and sequencing of short noncoding regions of cpDNA and dCAPS (derived cleaved amplified polymorphic sequences) markers have detected polymorphisms that have been used in phylogenetic and genetic diversity analyses in brassicas [3, 8, 10, 16, 17]. These techniques need either Polyacrylamide Gel Electrophoresis (PAGE) with silver staining or sequencing facility. On the other hand, PCR-RFLP also known as CAPS is a simpler (only agarose gels required), reliable, and fast technique and can encompass a large region (as many universal primers are available) of chloroplast genome for analyses.

Keeping this in view, our objective was to assess suitability of PCR-RFLP technique to study cpDNA variations in some wild brassicas belonging to the same cytodeme, which can facilitate (i) population genetics and phylogeographic studies for conservation purposes and (ii) analyses of maternal lineage and genetic relatedness, which is essential for breeding of brassicas. To the best of our knowledge, the present investigation for the first time reports interspecific and intraspecific variations in cpDNA regions of 14 wild brassicas (including endemic and endangered species, B. cretica and B. insularis, resp.), using PCR-RFLP technique.

2. Materials and Methods

2.1. Plant Material

Seeds of 31 accessions belonging to 14 wild relatives of brassicas of tribe Brassiceae were obtained from National Bureau of Plant Genetic Resources (NBPGR), New Delhi, India (Table 1). The seeds of each accession were sown in two replicates in pots and the plants thus obtained were maintained in Botanical Garden of Zakir Husain College, University of Delhi. Although the identification of germplasm is accurately maintained at NBPGR (a national level germplasm bank of India), the plant species were further confirmed by morphology based classification using Flora Europaea [18]. Fresh leaves of individuals (3–5) of each accession/species were collected, frozen, and stored at −80°C till DNA extraction. In addition, leaf material of Cardamine flexuosa (tribe Cardamineae) was collected from naturally growing population. This species was used as outgroup, since it belongs to tribe Cardamineae and is expected to be genetically distant from the rest of the wild species (which belong to tribe Brassiceae).

SpeciesAccessionCodeNumber of individuals

Brassica barrelieri EC662175Bba15

B. cretica EC662183Bcr15

B. elongata EC662178Bel14

B. gravinae EC662191Bgr15

B. insularis EC675815Bin15

B. maurorum EC662185Bma15

B. souliei EC662188Bso14

B. tournefortii IC560703Bto15

B. villosa EC675808Bvi15

Diplotaxis assurgens EC662192Das13

D. catholica EC662194Dca15

D. erucoides EC662197Der14

D. virgata EC662198Dvi13

Hirschfeldia incana EC675810Hin14

Cardamine flexuosa Cfl5


2.2. DNA Extraction, Amplification, and Digestion

Total genomic DNA was extracted in replicates of two from each individual, following protocol by Torres et al. [19]; subsequently quantified and working dilutions of 5 ng/μL were made.

For PCR amplification, six pairs of universal cpDNA primers (CD, DT, HK, K1K2, TF, and VL) described in Dumolin-Lapegue et al. [12] were used. Three replicates of PCR were carried out for each primer pair. The amplification was carried out in 30 μL of reaction mixture containing 0.2 μM of each primer, 200 μM of each of the four dNTPs, 2 mM MgCl2, 1 U of Taq DNA Polymerase in 1x buffer, provided by the manufacturer (Merck) of the enzyme, and 15 ng of genomic DNA. The PCR was set with an initial cycle of 4 min at 94°C, followed by 30 cycles of 45 s at 94°C, 45 s at 50–54°C, 2 min–4 min 30 s at 72°C, and finally 10 min extension at 72°C (Table 2). Agarose gel (1.2%) was used to run PCR products in 1X TBE buffer, along with 1 kilobase (kb) ladder as molecular size marker. Two restriction enzymes, HinfI and TaqI, chosen based on report of Cunha et al. [15] were used to digest the amplified products. Following digestion, the fragments were separated on 2.4% agarose gels, run at 3 V/cm for 3 h with 100-base pair (bp) and 50 bp ladders, as molecular size markers. All restriction digestions and electrophoresis were repeated thrice. Negative controls for PCR amplifications and restriction digestions were also set and run on gels along with the samples. The gels were stained with ethidium bromide, photographed, and documented using Gel Doc XR+ (BioRad) with Image Lab TM software.

Abbreviation of cpDNA primersaPCR conditionsSize of amplified fragment (bp)
Annealing temperatureExtension time

CD52°C4 min 30 sec2400
DT52°C2 min1200
HK51°C2 min 30 sec1750
K1K253°C3 min 2600
TF50°C2 min 30 sec1700
VL54°C4 min 30 sec3800

Abbreviation and primer sequences as in Dumolin-Lapegue et al. (1997) [12].
2.3. Data Analysis

All clearly resolved polymorphic restriction fragments were scored as 1 (present) or 0 (absent). A matrix of similarities between every pair of samples was created using Jaccard’s similarity coefficient [20], ; and are the total number of fragments analyzed in individuals and , respectively, and is the number of fragments shared by the two individuals. The similarity matrix was employed to construct a UPGMA dendrogram, using the SAHN-clustering and TREE programs from NTSYS-pc, version 2.2 [21]. A cophenetic matrix was produced from the tree matrix to test the goodness of fit of the cluster analysis to the similarity matrix on which it was based, by comparing the two matrices using the Mantel matrix correspondence test [22] in the MXCOMP program of the NTSYS-pc package.

3. Results

Six pairs of universal cpDNA primers were used to amplify 13.5 kb (approx.) region of chloroplast genome from wild relatives of brassicas. The size of the amplified fragments with each primer pair was same in all the species (Table 2). Of the 12 combinations (i.e., 6 primer pairs × 2 restriction enzymes), fragments obtained with HK-TaqI could not be clearly resolved and therefore not included in the analyses. A total of 219 restriction fragments (between 1 kb and 100 bp) were scored for analysis of polymorphisms. The interspecific/intergeneric variations between the 14 wild brassicas were scored by the combination of presence or absence of polymorphic fragments that were obtained from the PCR-RFLP patterns of each of the wild species. Two (K1K2-TaqI and DT-TaqI) out of 11 primer pair-restriction enzyme combinations were sufficient to distinguish all wild brassicas. The PCR-RFLP patterns obtained with DT-TaqI that distinguished some wild brassicas are shown in Figure 1.

For intraspecific variations, the size differences of the polymorphic fragments (with all the primer pair-restriction enzyme combinations) were assessed (Table 3). All individuals belonging to the same accession within a species were monomorphic. Of the 14 wild brassicas, 11 were represented by two or three accessions. Eight (B. barrelieri, B. cretica, B. elongata, B. insularis, B. maurorum, D. assurgens, D. catholica, and D. virgata) of the 11 species revealed intraspecific polymorphisms (Table 3) whereas three wild brassicas (B. tournefortii, B. villosa, and H. incana) showed no intraspecific polymorphism. B. gravinae, B. souliei, and D. erucoides were represented by one accession each (Table 1) and no intraspecific polymorphism was observed in these species. The combinations, DT-HinfI, DT-TaqI, K1K2-HinfI, and K1K2-TaqI, showed intraspecific polymorphisms in six wild species (B. cretica, B. elongata, B. insularis, B. maurorum, D. assurgens, and D. virgata; see Table 3). In B. barrelieri and D. assurgens only one primer pair-restriction enzyme combination could detect intraspecific polymorphisms, that is, HK-HinfI and K1K2-HinfI, respectively. Intraspecific polymorphisms in B. cretica and B. insularis were revealed with five and nine primer pair-restriction enzyme combinations, respectively (Table 3).

SpeciesAccessionPolymorphic fragmentsa (bp) in various primer pair-restriction enzyme combinations

B. elongata EC6621783300940
EC662180310150580 + 210 + 150
B. insularis EC675815110
D. catholica EC662194330550


B. cretica EC662183500310
B. elongata EC6621780220
B. insularis EC675815470140230140
B. maurorum EC662185170230


B. barrelieri EC6621750
D. catholica EC662194190


B. cretica EC662183370200330
B. elongata EC662178970520410340
EC662179600 + 370510410230
EC662180600 + 370510400340
EC662182600 + 370520410230
B. insularis EC675815200320 + 300510290
EC675816200320 + 300510290
D. assurgens EC662192220
D. virgata EC662198230210


B. cretica EC662183510360
B. insularis EC675815610350200


B. insularis EC675815230420

The numbers (1 to 4) shown besides each primer pair-restriction enzyme combination indicate the order of polymorphic fragments (highest to lowest size).

The genetic relatedness between the wild brassicas is represented in the dendrogram (Figure 2). The Mantel test revealed a high and significant cophenetic correlation (; ), thus, showing a very good fit to Jaccard’s similarity matrix. Two major clusters were observed in the dendrogram with C. flexuosa as the outgroup. The accessions of each species grouped together. The dendrogram separated the wild brassicas into two groups: group I, consisting of five species (B. barrelieri, B. cretica, B. insularis, B. villosa, and D. erucoides), and group II, including nine species (B. gravinae, B. elongata, B. maurorum, B. souliei, B. tournefortii, D. assurgens, D. catholica, D. virgata, and H. incana). All species in group I belonged to oleracea lineage. Group II was represented by members included in nigra lineage, except B. elongata and B. gravinae (which belong to oleracea group).

4. Discussion

In the present investigation, all the wild Brassica species were grown in field up to flowering and/or fruit set. Morphological characteristics were studied in the field grown plants and identification of all species was confirmed with the help of Flora Europaea [18]. All 31 accessions from 14 wild brassicas were subjected to PCR-RFLP of six cpDNA regions (with 11 primer pair-restriction enzyme combinations), which revealed intergeneric/interspecific and intraspecific polymorphisms. In an earlier study, Panda et al. [14] used PCR-RFLP technique (with 36 primer pair-restriction enzyme combinations) and reported absence of interspecific and intraspecific polymorphisms in members of oleracea group, that is, wild B. oleracea, B. alboglabra, B. bourgeaui, and B. montana. Here, it may be noted that the species used in both studies are different, although they belong to the same diploid cytodeme described under Brassica coenospecies. Cunha et al. [15] found 38 of 110 combinations (10 primer pairs × 11 restriction enzymes) useful in distinguishing only four diploid cultivated crucifers (B. nigra, B. oleracea, B. rapa, and Raphanus sativus). In the present study, of the 11 primer pair-restriction enzyme combinations, the PCR-RFLP patterns of two combinations (DT-TaqI and K1K2-TaqI; see Figure 1 for distinct PCR-RFLP patterns of some wild brassicas with DT-TaqI) were sufficient to distinguish the 14 wild species. Here, it may be suggested that, since DT-TaqI and K1K2-TaqI can distinguish the cytoplasm (maternal lineage) of wild brassicas (used in the present study), they may be used to assess natural hybridization processes within the Brassica coenospecies. Earlier, PCR-RFLP of cpDNA regions has been used to analyze maternal lineage of only cultivated radish [4]. Also, they can be used to characterize and/or confirm the alloplasmic or cytoplasmic male sterile lines of crop brassicas developed through various breeding programs.

Intraspecific polymorphisms were also observed in eight wild species with different primer pair-restriction enzyme combinations (Table 3). This information can facilitate population genetics and phylogeographic studies which is crucial for formulating conservation strategies. Five and nine primer pair-restriction enzyme combinations could reveal intraspecific variations in B. cretica and B. insularis, respectively. Although, earlier, cpDNA SSR markers have been used for understanding population genetic structure of B. cretica [10], the present set of polymorphisms obtained using PCR-RFLP technique provide additional set of cpDNA markers for similar studies. The number of accessions analyzed per species ranged between one and three. It is worthwhile to note that, even within this small number of accessions that were screened, intraspecific polymorphisms were observed. This result encourages the extension of PCR-RFLP technique for chloroplast genome analyses in larger number of wild species and their accessions.

The dendrogram showing the genetic relatedness of the wild brassicas was mostly in agreement with previous studies [7, 23, 24] except for B. elongata and B. gravinae which belong to oleracea lineage [24, 25]. Here, it may be suggested that PCR-RFLP of larger number of noncoding regions of cpDNA (which can detect more number of relevant interspecific polymorphisms) can further help in understanding genetic relationships. The psbD-trnT sequence which corresponds to amplicon of DT primer pair (used in the present investigation) and trnT-trnF sequence (corresponds to amplicon of primer pair TF) have been used in earlier studies, along with additional noncoding sequences [7, 23], and have revealed reliable genetic relationships amongst Brassicaceae members.

5. Conclusion

PCR-RFLP of cpDNA can reveal interspecific and intraspecific polymorphisms in wild brassicas. The primer-restriction enzyme combinations which have revealed intraspecific polymorphisms (as detailed in Table 3) in the wild brassicas including B. cretica and B. insularis (endemic and endangered species) can be useful for assessment of their population genetic structure and phylogeographic studies, which is important to formulate conservation strategies. Appropriate combinations of PCR-RFLP which can reveal that numerous interspecific polymorphisms (e.g., DT-TaqI and K1K2-TaqI) may be used for characterizing or confirming maternal lineage of natural hybrids and alloplasmic lines developed by cross breeding wild and crop brassicas. Thus, PCR-RFLP of cpDNA can be used for marker assisted Brassica breeding programs as well.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was supported under Star College Scheme funded by Department of Biotechnology, Ministry of Science and Technology, Government of India. The sanction number is SAN/no. 102/IFD/DBT/SAN/1911/2008-2009.


  1. S. I. Warwick, A. Francis, and R. K. Gugel, Guide to Wild Germplasm of Brassica and Allied Crops (tribe Brassiceae, Brassicaceae), part 4, 3rd edition, 2009.
  2. S. Bang, K. Tsutsui, S. Shim, and Y. Kaneko, “Production and characterization of the novel CMS line of radish (Raphanus sativus) carrying Brassica maurorum cytoplasm,” Plant Breeding, vol. 130, no. 3, pp. 410–412, 2011. View at: Publisher Site | Google Scholar
  3. Y. Zhang, Z. Fang, Q. Wang et al., “Chloroplast subspecies-specific SNP detection and its maternal inheritance in Brassica oleracea L. by using a dCAPS marker,” Journal of Heredity, vol. 103, no. 4, pp. 606–611, 2012. View at: Publisher Site | Google Scholar
  4. K. Yamane, N. Lü, and O. Ohnishi, “Chloroplast DNA variations of cultivated radish and its wild relatives,” Plant Science, vol. 168, no. 3, pp. 627–634, 2005. View at: Publisher Site | Google Scholar
  5. R. Ayotte, P. M. Harney, and V. S. Machado, “The transfer of triazine resistance from Brassica napus L. to B. oleracea L. I. Production of F1 hybrids through embryo rescue,” Euphytica, vol. 36, no. 2, pp. 615–624, 1987. View at: Publisher Site | Google Scholar
  6. L. A. Chapman and D. R. Goring, “Pollen-pistil interactions regulating successful fertilization in the Brassicaceae,” Journal of Experimental Botany, vol. 61, no. 7, pp. 1987–1999, 2010. View at: Publisher Site | Google Scholar
  7. Y.-W. Yang, P.-Y. Tai, Y. Chen, and W.-H. Li, “A study of the phylogeny of Brassica rapa, B. nigra, Raphanus sativus, and their related genera using noncoding regions of chloroplast DNA,” Molecular Phylogenetics and Evolution, vol. 23, no. 2, pp. 268–275, 2002. View at: Publisher Site | Google Scholar
  8. S. I. Warwick and C. A. Sauder, “Phylogeny of tribe Brassiceae (Brassicaceae) based on chloroplast restriction site polymorphisms and nuclear ribosomal internal transcribed spacer and chloroplast trnL intron sequences,” Canadian Journal of Botany, vol. 83, no. 5, pp. 467–483, 2005. View at: Publisher Site | Google Scholar
  9. S. Glémin, L. Vimond, J. Ronfort, T. Bataillon, and A. Mignot, “Marker-based investigation of inbreeding depression in the endangered species Brassica insularis,” Heredity, vol. 97, no. 4, pp. 304–311, 2006. View at: Publisher Site | Google Scholar
  10. K. Edh, B. Widén, and A. L. F. Ceplitis, “Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae),” Molecular Ecology, vol. 16, no. 23, pp. 4972–4983, 2007. View at: Publisher Site | Google Scholar
  11. L. Mondini, A. Noorani, and M. A. Pagnotta, “Assessing plant genetic diversity by molecular tools,” Diversity, vol. 1, no. 1, pp. 19–35, 2009. View at: Publisher Site | Google Scholar
  12. S. Dumolin-Lapegue, M.-H. Pemonge, and R. J. Petit, “An enlarged set of consensus primers for the study of organelle DNA in plants,” Molecular Ecology, vol. 6, no. 4, pp. 393–397, 1997. View at: Publisher Site | Google Scholar
  13. B. Sarin, J. P. M. Clemente, and A. Mohanty, “PCR–RFLP to distinguish three Phyllanthus sp., commonly used in herbal medicines,” South African Journal of Botany, vol. 88, pp. 455–458, 2013. View at: Publisher Site | Google Scholar
  14. S. Panda, J. P. Martín, and I. Aguinagalde, “Chloroplast and nuclear DNA studies in a few members of the Brassica oleracea L. group using PCR-RFLP and ISSR-PCR markers: a population genetic analysis,” Theoretical and Applied Genetics, vol. 106, no. 6, pp. 1122–1128, 2003. View at: Google Scholar
  15. C. Cunha, M. Tonguç, and P. D. Griffiths, “Discrimination of diploid crucifer species using PCR-RFLP of chloroplast DNA,” HortScience, vol. 39, no. 7, pp. 1575–1577, 2004. View at: Google Scholar
  16. C. J. Allender, J. Allainguillaume, J. Lynn, and G. J. King, “Simple sequence repeats reveal uneven distribution of genetic diversity in chloroplast genomes of Brassica oleracea L. and (n = 9) wild relatives,” Theoretical and Applied Genetics, vol. 114, no. 4, pp. 609–618, 2007. View at: Publisher Site | Google Scholar
  17. L. Liu, B. Zhao, D. Tan, and J. Wang, “Phylogenetic relationships of Brassicaceae in China: insights from a non-coding chloroplast, mitochondrial, and nuclear DNA data set,” Biochemical Systematics and Ecology, vol. 39, no. 4–6, pp. 600–608, 2011. View at: Publisher Site | Google Scholar
  18. “Psilotaceae to platanaceae,” in Flora Europaea, T. G. Tutin, N. A. Burges, A. O. Chater et al., Eds., vol. 1, Cambridge University Press, Cambridge, UK, 2nd edition, 2010. View at: Google Scholar
  19. A. M. Torres, N. F. Weeden, and A. Martín, “Linkage among isozyme, RFLP and RAPD markers in Vicia faba,” Theoretical and Applied Genetics, vol. 85, no. 8, pp. 937–945, 1993. View at: Publisher Site | Google Scholar
  20. P. Jaccard, “Nouvelles recherches sur la distribution florale,” Bulletin de la Société Vaudoise des Sciences Naturelles, vol. 44, pp. 223–270, 1908. View at: Google Scholar
  21. F. J. Rohlf, NTSYS-pc: Numerical Taxonomy and Multivariate Analysis System, version 2.20, Exeter Software, Setauket, NY, USA, 2005.
  22. N. Mantel, “The detection of disease clustering and a generalized regression approach,” Cancer Research, vol. 27, no. 2, pp. 209–220, 1967. View at: Google Scholar
  23. T. Arias and J. C. Pires, “A fully resolved chloroplast phylogeny of the brassica crops and wild relatives (Brassicaceae: Brassiceae): novel clades and potential taxonomic implications,” Taxon, vol. 61, no. 5, pp. 980–988, 2012. View at: Google Scholar
  24. S. I. Warwick and J. C. Hall, “Phylogeny of Brassica and wild relatives,” in Biology and Breeding of Crucifers, S. K. Gupta, Ed., pp. 19–36, CRC Press, Taylor & Francis, 2009. View at: Google Scholar
  25. A. K. Pradhan, S. Prakash, A. Mukhopadhyay, and D. Pental, “Phytogeny of Brassica and allied genera based on variation in chloroplast and mitochondrial DNA patterns: molecular and taxonomic classifications are incongruous,” Theoretical and Applied Genetics, vol. 85, no. 2-3, pp. 331–340, 1992. View at: Publisher Site | Google Scholar

Copyright © 2015 Bharti Sarin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.