Bangladesh J. Plant Taxon. 27(2): 261-271, 2020 (December) © 2020 Bangladesh Association of Plant Taxonomists THE IMPORTANCE IN DNA BARCODING OF THE REGIONS WHICH IS COVERING rRNA GENES AND ITS SEQUENCES IN THE GENUS QUERCUS L. YILMAZ AYKUT Faculty of Science and Arts, Department of Molecular Biology and Genetics, Uşak University, Uşak, Turkey Keywords: DNA Barcoding; IGS; rRNA gene; ITS; Quercus; Turkey. Abstract Turkey with 18 oak (Quercus) species is one of the richest country according to species number and diversity. The most important reason of the species diversity in Turkey is its location and geomorphological structure which increase climatic effects and seperate Turkey into different phytogeographic regions. Furthermore, hybridization behaviours which frequently observed between oak species, genetic drift, gene flow and ecological factors cause morphological variations in the plants species. All of these factors make it difficult to define the species concept for plant groups like oaks. Therefore, the region covering 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene and secondly intergenic spacer (IGS)/ 5S rRNA gene for barcoding were obtained from genbank and used as a useful tool for the determination and solution of the phylogenetic relations of taxonomically problematic species, also these barcoding regions were compared with each other according to species recognition ability for oak species. As a result, it can be stated that both barcoding regions have high variable sites based on sequence information to identify the species and evaluate relationships of species studied. Introduction DNA barcoding is very important molecular approach for definition biodiversity, evolutionary studies and especially for identifying species with taxonomically problems. DNA sequences prefered in DNA barcoding must have sufficient variability in grouping of species according to common characteristics and separation of taxonomically closely related species. Therefore in past years, different sequences regions belonging to genomic and plastid DNA are experienced and tried to find the best regions for DNA barcoding. Universal barcoding system would be avaluable resource for recognition of unambiguous species and in terms of speed, low cost, reliability (Piredda et al., 2011). Short DNA sequence that contain sufficient sequence variation to distinguish species is used for DNA barcode as molecular marker (Kress and Erickson, 2007). Especially internal transcribed spacer (ITS) regions of rDNA genes in genomic DNA are the sequences prefered the most commonly for plant molecular systematic studies (Baldwin et al., 1995; Alvarez and Wendel, 2003; Bailey, 2003; Sramko, 2008; Sramko et al., 2014). Also the external transcribed spacers (ETS) and the intergenic spacer (IGS) are widely utilized in phylogenetics in addition to ITS region. The cytochrome c oxidase-1 (CO1) gene from the mitochondria has enough nucleotide differentiation rates to identify many groups of animals and is routinely used to identify new species as an universal barcode (Hebert et al., 2003, 2004; Greenstone et al., 2005; Ward et al., 2005; Smith et al., 2006; Piredda et al., 2011; Hürkan, 2017) but this rate is relatively low and *Corresponding author, E-mail: aykut.yilmaz@usak.edu.tr mailto:aykut.yilmaz@usak.edu.tr 262 AYKUT unsuitable in plants (Chase et al., 2005; Kress et al., 2005; Fazekas et al., 2008; Hollingsworth et al., 2009). Therefore alternative barcode regions should be screened that can be universally successful in all species, however such a barcode region has not been found yet (Chase and Fay, 2009; Hollingsworth et al., 2009). Many regions of chloroplast genome for plant species as effective strategy for barcoding are recently used to resolve problems and the relationships in species level. Nevertheless there is still much debate related to the most suitable regions to be used in chloroplast genome. There is no barcoding region available to be used for all plant groups. Barcoding regions used together or whole chloroplast genome could provide enormous data and specificity for universal barcoding in plants. The regions and the region combinations belonging to chloroplast genome like rbcL, matK, trnK, trnH-psbA, atpB-rbcL, trnT-trnF are commonly and effectively used for plant phylogenetic analysis. As a result, barcoding is used as a useful tool for the determination and solution of the phylogenetic relations of taxonomically problematic species. The genus Quercus represented by over 500 species in the northern hemisphere show high phenotypic variation with natural hybrids (Manos et al., 2001; Borazan and Babaç, 2003; Yılmaz, 2018a). Turkey with 18 oak species belonging to three subgeneric sections (Quercus, Cerris and Ilex) is among the richest country with species number and diversity (Yaltırık, 1984). Section Quercus L. is characterized by the widest distribution and the greatest number of species among the sections which is presented in Turkey: Q. frainetto Ten., Q. petraea (Mattuschka) Lieb., Q. pontica C. Koch., Q. robur L., Q. infectoria Oliver, Q. hartwissiana Steven., Q. vulcanica (Boiss. Heldr. ex) Kotschy, Q. macranthera subsp. syspirensis (C. Koch.) Menitsky, Q. pubescens Willd and Q. virgiliana Ten. (Yaltirik, 1984). Section Cerris Loudon. is the second largest section with five species: Q. libani Olivier, Q. trojana Webb, Q. cerris L., Q. brantii Lindl. and Q. ithaburensis subsp. macrolepis (Kotschy) Hedge et Yalt. (Yaltirik, 1984). Section Ilex Loudon is represented by three species: Q. ilex L., Q. coccifera L. and Q. aucheri Jaub. et Spach. (Yaltirik, 1984). The most important reason of the high species diversity in Turkey is its location and geomorphological structure which increase climatic effects and seperate Turkey into different phytogeographic regions (Uslu and Bakış, 2012; Yılmaz, 2018b). Turkey is between the Asian and European continents that is used an important migration route for many plants and animals. Another factor on species diversity and number in Turkey is the Anatolian Diagonal which divides Anatolia as eastern and western parts (Davis, 1971; Çıplak et al., 1993; Borazan and Babaç, 2003; Yılmaz, 2018 a,b). Furthermore oak species can spread across wide geographic regions via wind and grow in mixed populations that increase the hybridization between species belonging to same or different sections (Hokanson et al., 1993; Kremer and Petit, 1993; Bacilieri et al., 1996). In addition to all factors, insufficient diagnostic morphological characters that it is sometimes not possible to identify oak species due to high morphological variation (Denk and Grimm, 2010; Simeone et al., 2013) and the lack of investigations such as ecological, historical and genetic descriptors make problematic the genus Quercus in Turkey and similarly in the world. Hybridization behaviours, gene flow, genetic drift, ecological factors and epigenetic mechanisms cause morphological variation in the plants species. The classical taxonomic system that is based on the morphological similarity of individuals makes it difficult to define the concept of biological species, especially for plant groups like oaks. Therefore molecular markers instead of morphological characters are frequently prefered to identify the oak species and understand the THE IMPORTANCE IN DNA BARCODING OF THE REGIONS 263 oak evolution (Oh and Manos, 2008; Denk and Grimm, 2010; Simeone et al., 2013; Yılmaz et al., 2013; Yılmaz, 2016). Especially DNA barcoding has been used as the most useful tool in solving these problems. The objective of this study is to evaluate phylogenetic relationships of Quercus species by using the 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene and intergenic spacer (IGS)/ 5S rRNA gene from genbank and compare these barcoding regions according to species recognition ability. Materials and Methods Study materials and DNA regions Sequences analysis for Quercus taxa was seperately done for two regions of rDNA containing 18S rRNA gene(partial)/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene (partial) and secondly intergenic spacer (IGS)/ 5S rRNA gene (partial). Informations related to studied taxa were obtained from National Centre of Biotechnology Information (NCBI). Studied taxa and genbank codes for the rRNA gene regions analysed in this study are presented in Table 1 and 5. Sixteen taxa for first region containing ITS1 and ITS2 together with related genes of rRNA and 15 taxa for the other containing IGS and 5S rRNA gene were prefered and analysed for phylogenetic relations. Almost all of taxa selected for this study belong to Turkey except a few species. While the locations of 15 studied taxa for first region analysed (18S rRNA gene ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA) belong to completely Turkey except three species, all taxa prefered for second region belong to Turkey. Table 1. Studied taxa, sections of the species studied and their Gen Bank accession numbers for 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene region. Genus Section Species GenBank Acc. No. Quercus Quercus Q. hartwissiana Q. vulcanica Q.infectoria subsp. boisseri Q. infectoria subsp. infectoria Q. macranthera Q. frainetto Q. petraea subsp. petraea Q. pubescens Q. pontica FM244036 FM244264 FM243942 FM244072 FM244101 FM244015 FM244134 FM244253 FM244159 Cerris Q. ithaburensis subsp. macrolepis Q. trojana Q.cerris Q. brantii FM243873 FM243920 FM243851 FM243826 Ilex Q. ilex Q. coccifera Q. aucheri FM244455 FM244318 FM244282 Sequence alignment and Phylogenetic analysis Multiple sequence alignments for both regions were seperately performed by using Molecular Evolutionary Genetics Analysis (MEGA). The probabilities of substitution from one base to 264 AYKUT another base, transition/transversion ratios for purines-pyrimidines and overall, nucleotide frequencies were computed by using alignment sequences that were edited (Tables 2-4 & 6-8). Neighbour-joining dendrograms that bootstrap values are reported above branches for two regions such as 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene and secondly IGS/ 5S rRNA gene were obtained with MEGA X program (Figs 1- 2). All positions containing gaps and missing data were eliminated (complete deletion option). Consequently, evolutionary analyses were conducted by using a total of 283 positions in the final dataset for IGS/ 5S rRNA gene and a total of 690 positions in the final dataset for 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene. Fig. 1. Neighbor-Joining dendrogram given by the 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene for 16 Quercus taxa. Bootstrap values are reported in the branches. All positions containing gaps and missing data were eliminated (complete deletion option). There were a total of 690 positions in the final dataset. Fig. 2. Neighbor-Joining dendrogram given by the IGS/5S rRNA gene for 15 Quercus taxa. Bootstrap values are reported in the branches. All positions containing gaps and missing data were eliminated (complete deletion option). THE IMPORTANCE IN DNA BARCODING OF THE REGIONS 265 Results and Discussion Eighteen oak species belonging to three subgeneric sections (Quercus, Cerris and Ilex) currently occur in Turkey that is one of the richest country with species diversity and number (Yaltirik, 1984). Two rDNA regions that sequence information is provided from taxonomy database of NCBI were used for phylogenetic analysis. The foundamental aims of the study is firstly to evaluate the taxa from Turkey belonging to genus Quercus according to phylogenetic relations and to contribute the solution of taxonomic problems, secondly to compare two rDNA regions frequently used in DNA barcoding and evaluate the region that gives the best results for barcoding. Analysis results for region covering 18S rRNA gene/ ITS1/5.8S rRNA gene/ ITS2/25S rRNA gene: The valuable informations about the taxonomy of studied taxa were provided from analysis of the first region. This genomic DNA region has the quite wide sequence data covering three rRNA gene (18S rRNA gene, 5.8S rRNA gene and 25S rRNA gene) and two spacer regions (ITS1 and ITS2) giving the information useful for plant systematics in species and generic level. This DNA region has alignment length of 697 bp for taxa studied and showed 94 variable sites. Studied taxa and accession numbers obtained from NCBI are given in Table 1. Table 2 shows the probability of substitution (r) from one base to another base. For simplicity, the sum of r values is made equal to 100. Rates of different transitional substitutions are shown in bold and those of transversional substitutions are shown as italics in Table 2. This analysis involved 16 nucleotide sequences. All positions containing gaps and missing data were eliminated (complete deletion option). There were a total of 690 positions in the final data set. Evolutionary analyses were conducted in MEGA X. Table 2. The probability of substitution (r) from one base (row) to another base (column). A T C G A - 1.67 3.02 12.19 T 1.85 - 39.33 2.84 C 1.85 21.8 - 2.84 G 7.91 1.67 3.02 - Transitional substitutions with the rate of 81,23 % are much higher than transversional substitutions according to total base substitutions showing in Table 2. Moreover, transitional substitutions of the pyrimidines are higher than purines (Table 2). In the comparison of purines (k1) and pyrimidines (k2) according to transition/transversion ratio, pyrimidines with 13,02 show the higher value from purines (Table 3). Overall transition/transversion ratio (R) is 4,01 (R = [A*G*k1 + T*C*k2]/[(A+G)*(T+C)]) in the evaulation of all positions in the final dataset (Table 3). The nucleotide frequencies are 19.66% (A), 17.84% (T/U), 32.19% (C), and 30.30% (G) (Table 4). It can be stated that the percentage of G and C bases for all studied Quercus taxa for the DNA region containing 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene is higher than the percentage of A and T/U bases (Table 4). Neighbor-Joining (NJ) dendrogram was drawn to show the phylogenetic relations of 16 Quercus taxa (Fig. 1). The evolutionary history was inferred using the Neighbor-Joining method (Saitou and Nei, 1987). The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method (Tamura et al., 2004) and are in the 266 AYKUT units of the number of base substitutions per site. The differences in the composition bias among sequences were considered in evolutionary comparisons (Tamura and Kumar, 2002). It can be stated as a result of the examination of NJ tree that the DNA region of interest has sufficient information for species separation and sectional grouping of species. Sequence datas provided from studied taxa seperates the species to three sections as Quercus, Ilex and Cerris (Fig. 1). All samples are clearly differentiated from each other. Furthermore, NJ tree showed that this barcoding region for the Quercus taxa has the enough sequence information with 94 variable sites to compare and evaulate especially taxonomically closely related species. Table 3. The transition/transversion rates for purines and pyrimidines. Transition/transversion ratio Purines (k1) 4.29 Pyrimidines (k2) 13.02 Overall (R) 4.01 Table 4. Nucleotide frequencies for each base. Transition/transversion ratio Nucleotide Frequence (%) A 19.66 T/U 17.84 C 32.19 G 30.30 Analysis results for region covering IGS and 5SrRNA gene IGS/5S rRNA gene region has the alignment length of 398 bp and 115 variable sites for taxa studied. Studied taxa and accession numbers obtained from NCBI are given in Table 5. Variable region of the IGS/5S rRNA gene sequences is widest than the variable sites of 18S rRNA gene/ ITS1/ 5.8S rRNA gene/ ITS2/ 25S rRNA gene sequences. In other words, it can be stated that although IGS/5S rRNA gene has shorter DNA sequences, it has more distinctive information for Quercus taxa. The probability of substitution (r) from one base to another base was shown in Table 6. This analysis involved nucleotide sequences belonging to 15 Quercus taxa. All positions containing gaps and missing data were eliminated and evolutionary analyses were conducted in MEGA X. The rate of transitional substitutions with 63.78 are higher than transversional substitutions according to total base substitutions (Table 6). Furthermore, 65% of transitional substitutions is caused by base substitutions of pyrimidines with each other. Transition/transversion ratios of purines and pyrimidines are 2.87 and 4.00, respectively (Table 7). In other words, transitional substitutions of both base group are higher than transversional substitutions. Overall transition/transversion ratio is 1.86 in the evaulation of all positions in the final data set (Table 7). The nucleotide frequencies for the DNA region containing IGS/5S rRNA gene are 17.74% (A), 28.93% (T/U), 28.32% (C), and 25.02% (G) (Table 8). Finally, phylogenetic relations of 15 Quercus taxa was showed with Neighbor-Joining (NJ) dendrogram (Fig. 2). The evolutionary distances were computed using the Maximum Composite Likelihood method. It can be stated that NJ tree of interested DNA region seperated the studied taxa to three group as sectional and besides had sufficient information for species separation and THE IMPORTANCE IN DNA BARCODING OF THE REGIONS 267 phylogenetic relations of species. Furthermore, NJ tree showed that this barcoding region like other studied region has the enough sequence information to evaulate taxonomically problematic species like members of the genus Quercus. Table 5. Studied taxa , sections of the species studied and their Gen Bank accession numbers for IGS/5S rRNA gene region. Genus Section Species GenBank Acc. No. Quercus Quercus Q. vulcanica Q.infectoria subsp. boisseri Q. infectoria subsp. infectoria Q. frainetto Q. petraea subsp. petraea Q. pubescens Q. pontica FM243389 FM243123 FM243193 FM243158 FM243232 FM243345 FM243258 Cerris Q. ithaburensis subsp. macrolepis Q. trojana Q.cerris Q. brantii Q. libani FM242972 FM243104 FM242924 FM242906 FM242961 Ilex Q. ilex Q. coccifera Q. aucheri FM243603 FM243492 FM243431 Table 6. The probability of substitution (r) from one base (row) to another base (column) for IGS/5S rRNA gene region. A T C G A - 5.24 5.13 13.02 T 3.21 - 20.54 4.53 C 3.21 20.99 - 4.53 G 9.23 5.24 5.13 - Table 7. The transition/transversion rates belonging to IGS/5S rRNA gene region for purines and pyrimidines. The transition/transversion ratio Purines 2.87 Pyrimidines 4.00 Overall 1.86 Table 8. Nucleotide frequencies provided from IGS/5S rRNA gene region for each base. Nucleotide Frequence (%) A 17.74 T/U 28.93 C 28.32 G 25.02 268 AYKUT In Turkey, oaks which are represented by 18 species have wide geographical ditribution and dominated the most of forests. Oaks having such a wide geographical spread and variety of species has been used many purposes because of economically importance, such as foods, furniture and especially fuel wood. This situations increases the taxonomic problems in the genus and make it difficult the species definition. Additionally, location of Turkey between the Asian and European continents serve as a migration route for many plants such as oaks. Weak reproductive barriers and mixed populations in many regions are observed between oak species. All of these factors may be reason of the extensive hybridization behaviours, morphological variation in the species level and also taxonomic problems. The determination of succesful barcoding regions for the genus Quercus would have a considerable effect in improving available taxonomic problems and in the species level identification. For this reason, two genomic DNA region containing 18S rRNA gene/ITS1/5.8S rRNA gene/ITS2/25S rRNA gene and intergenic spacer (IGS)/5S rRNA gene proposed by the Consortium for the Barcode of Life (CBOL) were used as molecular markers and compared with each other. As a result, it can be stated that both barcoding regions have high variable sites based on sequence information to identify the species and evaluate relationships of species studied. Furthermore, it was observed that neighbour-joining dendrograms containing the full oak data for both barcoding regions seperated the species to three group as sectional and besides studied taxa from each other. When it is evaluated the phylogenetic relationships of oaks which are completely similarly grouped as sectional by both NJ dendrograms; it can be stated that the evolutionary distances among Q. ithaburensis subsp. macrolepis, Q. brantii, Q. trojana and Q. cerris belonging to section Cerris showed similarity for each barcoding region. Also it is observed that section Ilex and section Cerris is phylogenetically more close than section Quercus for both barcoding region. The comparisons of three species (Q. coccifera, Q. ilex and Q. aucheri) belonging to section Ilex show to us that Q. ilex and Q. aucheri are closer two taxa than Q. coccifera. Similarly, Yılmaz et al., (2013) stated in previous report on DNA comparison of related three species from section Ilex that Q. ilex and Q. aucheri were observed as close two separate groups and populations of Q. coccifera showed more differences than populations of Q. ilex and Q. aucheri. Besides that, other study on the based the all chromosomal parameters such as length range, haploid complement, A1 and A2 values of these three taxa show similarity the results provided from barcoding regions and supports the study results (Yılmaz, 2018b). In a previous study; Denk and Grimm (2010) used the ITS and 5S-IGS data to recognize the major infrageneric groups and the phylogenetic relationships among the species of Quercus from western Eurasia. However sequence regions encoding for the 18S, 5.8S and 25S rRNA were excluded from the analyses by Denk and Grimm (2010) on the contrary of this study. While the individuals of Q. pontica formed a distinct group in the study of Denk and Grimm (2010), in this study NJ dendrograms showed that Q. pontica evaluated within the section Quercus is the outmost species in the comparison to other species belonging to section Quercus. The comparisons of alignment lengths and variable sites of the barcoding regions studied show to us that although IGS-5S rRNA gene region with the 398 bp alignment length is smaller than other barcoding region, it exhibit more sequence variation with the range of 28.29%. However, when the sites with missing/ambiguous data and gaps were excluded for effective analyses, IGS-5S rRNA gene for studied taxa show 16.83% variation range. In other words, it has high missing data in the comparison to other barcoding region containing the sequence variation of 13.48%. THE IMPORTANCE IN DNA BARCODING OF THE REGIONS 269 Analyses especially for the species whose sequence lengths differ due to regions containing wide deletions exhibit missing/ambiguous data and gaps in sequence alignment. All species of the section Quercus analyzed using the IGS-5S rRNA gene sequence information have regions of such deletion in comparison to other species belonging to section Cerris and Ilex. Denk and Grimm (2010) states that "A number of newly assembled and gene bank sequences include missing data due to the fact that a guanine-rich region within the 5′ ITS1 region can be difficult to sequence". For this reason, Denk and Grimm whose added the sequence information to NCBI GenBank used by us re-run the sequencing to guarantee at least one completely sequenced ITS clone per individual. It can be said that both barcoding regions have important sequence information for species identification and evaluation of evolutionary relations in oaks, also these are recommended for further studies. In Turkey, another important reason that makes it difficult to understant the oaks besides hybridization is the lack of adequate conservation programs. Turkey is a very valuable country with 11000 taxon and 35% endemism rate in terms of plant diversity (Vural, 2003). Therefore, the results of the study are important for the determination of plant diversity and the conservation of genetic resources. Especially, Q. aucheri, Q. vulcanica and Q. macranthera subsp. syspirensis which are endemic taxa are valuable resources due to restricted distribution area. While distribution area of Q. aucheri is restricted to south-west Anatolia in Turkey, Q. macranthera subsp. syspirensis which is distributed in north and north-east regions of Anatolia has shown wider distribution than Q. aucheri. Other endemic species, Q. vulcanica distributed from 1200 to 2000 m altitude has more restricted area and isolated habitats such as Isparta-Eğirdir (Yukari Gokdere village), Konya-Sultan Mountains and Kutahya-Turkmen Mountains when compared with other species. However, Q. vulcanica has been faced with the threat of extinction. Furthermore, there is not enough protection program for conservation of oak biodiversity. This study contributes to understant the biodiversity and genetic resources of oaks besides the understanding the phylogenetic relationships of the oaks. Acknowledgments The author would like to thank to NCBI for sequence informations used in this study and the authors sharing the sequence informations in NCBI. References Alvarez, I. and Wendel, J.F. 2003. Ribosomal ITS sequences and plant phylogenetic inference. Mol. Phylogen. Evol. 29: 417 ̶434. Bacilieri, R., Ducousso, A., Petit, R.J. and Kremer, A. 1996. Mating system and asymmetric hybridization in a mixed stand of european oaks. Evolution. 50: 900 ̶908. Bailey, C.D. 2003. Characterization of angiosperm nrDNA polymorphism, paralogy, and pseudogenes. Mol. Phylogen. Evol. 29: 435 ̶455. Baldwin, B.G., Sanderson, M.J., Porter, J.M., Wojciechowski, M.F., Campbell, C.S. and Donoghue, M.J. 1995. The its Region of Nuclear Ribosomal DNA: A Valuable Source of Evidence on Angiosperm Phylogeny. Ann. Missouri Bot. Gard. 82(2): 247 ̶ 277. Borazan, A. and Babaç, M.T. 2003. Morphometric leaf variation in oaks (Quercus) of Bolu, Turkey. Ann Bot Fennici. 40: 233 ̶ 242. Chase, M.W. and Fay, M.F. 2009. Barcoding of plants and fungi. Science. 325(5941): 682 ̶ 683. 270 AYKUT Chase, M.W., Salamin, N. and Wilkinson, M. 2005. Land plants and DNA barcodes: short-term and long- term goals. Philosophical transactions of the Royal Society of London, Series B, Biological sciences. 360: 1889 ̶ 1895. Çıplak, B., Demirsoy, A. and Bozcuk, N. 1993. Distribution in Orthoptera in relation to the Anatolian Diagonal in Turkey. Articulata. 8: 1 ̶20. Davis, P.H. 1971. Distribution patterns in Anatolia with particular reference to endemism. In: Davis PH, Harper PC, Hedge IC, editors. Plant life of South West Asia. Edinburgh: Botanical Society of Edinburgh p. 15–27. Denk, T. and Grimm, G.W. 2010. The oaks of western Eurasia: Traditional classifications and evidence from two nuclear markers. Taxon. 59(2): 351 ̶366. Fazekas, A.J., Burgess, K.S., Kesanakurti, P.R., Graham, S.W., Newmaster, S.G., Husband, B.C., Percy, D.M., Hajibabaei, M. and Barrett, S.C.H. 2008. Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well. PLoS ONE. 3: e282. Greenstone, M.H., Rowley, D.L., Heimbach, U., Lundgren, J.G., Pfannenstiel, R.S. and Rehner, S.A. 2005. Barcoding generalist predators by polymerase chain reaction: carabids and spiders. Mol. Ecol. 14: 2347 ̶ 3266. Hebert, P.D.N., Ratnasingham, S. and deWaard, J.R. 2003. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society of London, series B. 270: 96-99. Hebert, P.D.N., Penton, E.H., Burns, J.M., Janzen, D.H. and Hallwachs, W. 2004. Ten species in one: DNA barcoding reveals cryptic species in the Neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the USA. 101: 14812 ̶ 14817. Hokanson, S.C., Isebrands, J.G., Jensen, R.J. and Hancock, J.F. 1993. Isozyme variation in oaks of the Apostle Islands in Wisconsin: Genetic structure and levels of inbreeding in Quercus rubra and Quercus ellipsoidalis (Fagaceae). Amer. J. Bot. 80: 1349 ̶1357. Hollingsworth, M.L., Clark, A.A., Forrest, L.L., Richardson, J., Pennington, R.T., Long, D.G., Cowan, R., Chase, M.W., Gaudeul, M. and Hollingsworth, P.M. 2009. Selecting barcoding loci for plants: evaluation of seven candidate loci with species-level sampling in three divergent groups of land plants. Mol. Ecol. Res. 9: 439 ̶457. Hürkan, K. 2017. Karasal Bitkilerde DNA Barkodlama: Bazı DNA Barkod Bölgelerinin İncelenmesi. The International Journal of Innovative Approaches in Science Research. 1(1): 57 ̶67. Kremer, A. and Petit, R.J. 1993. Gene diversity in natural populations of oak species. Ann. For Sci. 50: 186 ̶ 202. Kress, W.J. and Erickson, D.L. 2007. A two-locus global DNA barcode for land plants: the coding rbcL gene complements the non-coding trnH-psbA spacer region. PLoS ONE. 2: e508. Kress, W.J., Wurdack, K.J., Zimmer, E.A., Weight, L.A. and Janzen, D.H. 2005. Use of DNA barcodes to identify flowering plants. Proceedings of the National Academy of Sciences of the USA. 102: 8369 ̶ 8374. Manos, P.S., Zhou, Z. and Cannon, C.H. 2001. Systematics of Fagaceae: phylogenetic tests of reproductive trait evolution. Int. J. Plant Sci. 162: 1361 ̶ 1379. NCBI, National Centre of Biotechnology Information, https://www.ncbi.nlm.nih.gov/genbank. Oh, S.H. and Manos, P.S. 2008. Molecular phylogenetics and cupule evolution in Fagaceae as inferred from nuclear CRABS CLAW sequences. Taxon. 57: 434 ̶451. Piredda, R., Simeone, M.C., Attimonelli, M., Bellarosa, R. and Schirone, B. 2011. Prospects of barcoding the Italian wild dendroflora: oaks reveal severe limitations to tracking species identity. Mol. Ecol. Res. 11: 72 ̶83. Saitou, N. and Nei, M. 1987. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4: 406 ̶425. https://www.ncbi.nlm.nih.gov/genbank. THE IMPORTANCE IN DNA BARCODING OF THE REGIONS 271 Simeone, M.C., Piredda, R., Papini, A., Vessella, F. and Schirone, B. 2013. Application of plastid and nuclear markers to DNA barcoding of Euro-Mediterranean oaks (Quercus, Fagaceae): Problems, prospects and phylogenetic implications. Bot. J. Linn. Soc. 172: 478 ̶499. Smith, M.A., Woodley, N.E., Janzen, D.H., Hallwachs, W. and Hebert, P.D.N. 2006. DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: Tachinidae). Proceedings of the National Academy of Sciences of the USA. 103: 3657 ̶3662. Sramko, G. 2008. Sequence variability of the nrITS in the Ophrys fuciflora species-complex of the Mediterranean bee-orchid (Ophrys L.) genus. Department of Botany. University of Debrecen, Debrecen. Sramko, G., Molnar, A.V., Hawkins, J.A. and Bateman, R.M. 2014. Molecular phylogeny and evolutionary history of the Eurasiatic orchid genus Himantoglossum s.l. (Orchidaceae). Ann. Bot. 114: 1609 ̶ 1626. Tamura, K., Nei, M. and Kumar, S. 2004. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proceedings of the National Academy of Sciences (USA). 101: 11030 ̶ 11035. Tamura, K. and Kumar, S. 2002. Evolutionary distance estimation under heterogeneous substitution pattern among lineages. Mol. Biol. Evol. 19: 1727 ̶1736. Uslu, E. and Bakış, Y. 2012. Geographic distribution of Turkish Oaks. Dendrobiology. 67: 41 ̶48. Vural, M. 2003. Türkiye’nin tehlike altındaki bitkileri. FAO/BM Tematik Grubu, Türkiye’de Biyolojik Çeşitlilik ve Organik Tarım Çalıştay Raporu, 15-16 Nisan 2003. D168-183. Ward, R.D., Zemlak, T.S., Innes, B.H., Last, P.R. and Hebert, P.D.N. 2005. DNA barcoding Australia’s fish species. Philosophical transactions of the Royal Society of London, Series B, Biological sciences. 360: 1847 ̶ 1857. Yaltirik, F. 1984. Türkiye meşeleri teşhis kılavuzu. İstanbul: Yenilik Basımevi. Yılmaz, A. 2016. Phylogenetic relationships of the genus Quercus L. (Fageceae) from three different sections. African J. Biotech. 15(40): 2265 ̶ 2271. Yılmaz, A. 2018a. Karyomorphology of some Quercus L. species from section Quercus and Cerris in Turkey. Caryologia: International Journal of Cytology, Cytosystem. Cytogenet. 71(3): 210 ̶216. Yılmaz, A. 2018b. Cytogenetic Relationships of Turkish Oaks. Cytogenetics - Past, Present and Further Perspectives, Chapter 2. Yılmaz, A., Uslu, E. and Babaç, M.T. 2013. Molecular diversity among Turkish oaks (QUERCUS) using random amplified polymorphic DNA (RAPD) analysis. African J. Biotech. 12(45): 6358 ̶6365. (Manuscript received on 11 May 2019; revised on 10 July 2020)