Research Article |
Corresponding author: Bekir Dogan ( doganbekir2000@yahoo.com ) Academic editor: Peter de Lange
© 2015 Bekir Dogan, Ahmet Duran, Özlem Çetin, Meryem Öztürk, Esra Martin.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Dogan B, Duran A, Çetin Ö, Öztürk M, Martin E (2015) Study of phylogenetic relationship of Turkish species of Klasea (Asteraceae) based on ISSR amplification. PhytoKeys 56: 29-40. https://doi.org/10.3897/phytokeys.56.5608
|
Klasea is a taxonomically complex genus in which there are many problems, mostly with Klasea kotschyi and K. hakkiarica. It is challenging to differentiate the genera based on morphological characters alone. Revision studies performed on the basis of molecular data obtained from studies conducted in recent years have made the phylogenetic relationships and systematic positions of the taxa more apparent and reliable. In this study, Klasea, Serratula, Jurinea and Centaurea species native to Turkey, were collected from different localities of Anatolia and DNA was isolated from the collected samples. The data were analyzed ordination analyses including UPGMA and PCA using NTSYSpc 2.1. The infrageneric and intergeneric phylogenetic relationships between Klasea and other related genera were also characterized. The Klasea species were grouped into three clusters. It was determined that taxa Klasea kotschyi and K. hakkiarica are separate but closely related. Moreover, it was observed that the Klasea lasiocephala a separate group within the genera. Clearly the genera Klasea, Serratula, Jurinea and Centaurea are phylogenetically differentiated on the dendogram.
Asteraceae , ISSR, Klasea , Serratula , Molecular systematics, Phylogeny
The tribe Cardueae (Asteraceae) is generally accepted to be classified into five subtribes named Echinopinae, Carlininae, Carduinae, Centaureinae and Cardopatiinae (
16 species were reported for the genus Serratula in Turkey (
Currently, morphological revisions of various plant taxa are often supported by molecular data (
The RAPD (Randomly Amplified Polymorphic DNA) fingerprinting method is widely used and has a wide range of applications (
Simple sequence repeats (SSRs), also known as microsatellites, are tandemly repeated di-, tri, tetra- or penta-nuclotide sequences (mainly within the range of 10–80 repeats of the core unit) that are abundant within eukaryotic genomes. A high level of genomic variation is generated by the more or less evenly distributed microsatellite sequences present within the plant and animal genomes. The high levels of genomic variation are widely used for genetic variation analysis of both wild plants (
In this study, Klasea species, which are difficult to delineate using morphological traits, were collected from their natural habitats in Turkey. DNA was isolated and fingerprinting was performed using a highly reliable and reproducible technique that mimics the application ease of RAPDs. The method employed to assess the genetic diversity and to resolve the genetic relationships among the species is a technique derived from SSR characterization based on PCR amplification of ISSR regions primed by a single oligonucleotide corresponding to the targeted repeat motif. The SSR-containing primers are usually 16-25 base pair long oligonucleotides anchored at the 3’- or 5’-end by two to four arbitrary, and often degenerate, nucleotides (
The aim of this study was to determine the genetic relationships among selected Anatolian-originated Klasea, Serratula, Jurinea and Centaurea species collected from diverse regions of Turkey and to use a DNA-based molecular marker system to resolve the unclear and controversial status of these species based on conventional morphological characters.
Silica gel dried plant leaf samples belonging to 15 Klasea taxa and Serratula tinctoria, and 2 out-group taxa (Jurinea and Centaurea) were collected from the natural flora of Turkey. The species and provinces of their localities are as follows: Klasea quinquefolia (Artvin), K. oligocephala (Kahramanmaraş), K. kotschyi (Bitlis), K. serratuloides (Van), K. erucifolia (Erzurum), K. lasiocephala (Antalya), K. cerinthifolia (Kahramanmaraş), K. grandifolia (Antalya), K. radiata subsp. radiata (Kars), K. hakkiarica (Hakkari), K. haussknechtii (Muş), K. coriaceae (Kars), K. radiata subsp. biebersteiniana (Kars), K. kurdica (Osmaniye), K. bornmuelleri (Malatya), Serratula tinctoria (Bolu), Centaurea ptosimopappoides (Adana), C. straminicephala (Erzurum), and Jurinea cataonica. For details see Table
List of sampled taxa. Including location data, collectors, and herbarium in which the voucher specimens are accessioned.
Species | Voucher |
---|---|
Klasea serratuloides | Turkey, Van: Van to Gurpinar, 2125 m, 38°24.434’N, 043°23.079’E, 19.07.2009, B.Doğan 2117 & A.Duran ( |
K. lasiocephala | Turkey, Antalya, Gazipasa, Çayıryaka mountain pasture, 1730 m, 36°.30.027'N, 032°32.181'E, 30.06.2009, B.Doğan 2105 & A.Duran ( |
K. bornmuelleri | Turkey, Malatya, Darende, near Akçatoprak, 1010 m, 38°30.064'N, 037°33.907'E, 17.07.2009, B.Doğan 2110 & A.Duran ( |
K. kurdica | Turkey, Osmaniye, Yarpuz, 1465 m, 37°00.774'N, 036°26.683'E, 15.07.2009, B.Doğan 2106 & A.Duran ( |
K. coriaceae | Turkey, Kars, Tuzluca to Kağızman, 1055 m, 40°06.399'N, 043°29.567'E, 20.07.2009, B.Doğan 2122 & A.Duran ( |
K. cerinthifolia | Turkey, Kahramanmaraş, Ahir mountain, 990 m, 37°36.470'N, 036°52.917'E, 16. 07.2009, B.Doğan 2107 & A.Duran ( |
K. grandifolia | Turkey, Antalya, Akseki, Süleymanlı village, 1425 m, 37°17.980'N, 031°46.520'E, 31.07.2009, B.Doğan 2130 & A.Duran ( |
K. haussknechtii | Turkey, Muş, Malazgirt, Karıncalı village, 1840 m, 39°21.219'N, 042°20.010'E, 18.07.2009, B.Doğan 2113 ( |
K. radiata subsp. biebersteiniana | Turkey, Kars, Kağızman, Akçay to Cumaçay, 1830 m, 20.07.2009, B.Doğan 2124 & A.Duran ( |
K. radiata subsp. radiata | Turkey, Kars, Arpaçay, Kardeşköy to Dağköy, 2190 m, 06.08.2010, 40°55.087'N, 043°11.209'E, B.Doğan 2283 & A.Duran |
K. hakkiarica | Turkey, Hakkari, Cilo mountain, Kırıkdağ, near dez stream, 2210 m, 37°32.974'N, 043°57.615'E, 07.08.2009, B.Doğan 2132 & A.Duran ( |
K. kotschyi | Turkey, Bitlis, Tatvan, Sapur village, 1965 m, 38°26.154'N, 042°24.413'E, 06.08.2009, B.Doğan 2131 & A.Duran ( |
K. quinquefolia | Turkey, Artvin, Ardanuç, Boyalı village, 1210 m, 41°06.967'N, 042°07.283'E, 11.08.2009, B.Doğan 2139 & A.Duran ( |
K. erucifolia | Turkey, Erzurum, Köprüköy, Eğirmez village, 1635 m, 39°57.056'N, 041°51.530'E, 09.08.2009, B.Doğan 2137 & A.Duran ( |
K. oligocephala | Turkey, Kahramanmaraş, Ahir mountain, 995 m, 37°36.475'N, 036°52.947'E, 16. 07.2009, B.Doğan 2108 & A.Duran ( |
Serratula tinctoria | Turkey, Bolu, Gerede to Bolu, 28. km, 1105 m, 09.08.2010, 40°45.340'N, 031°54.888'E, B.Doğan 2290 & A.Duran ( |
Jurinea cataonica | Turkey, Erzincan, Old Çayırlı road, 10. km, 1750 m, 39°47.954'N, 039°30.343'E, 07.08.2005, B.Doğan 1029 ( |
Centaurea ptosimopappoides | Turkey, Adana,Aladağ to Kızıldag, 890 m, 19. 06. 2010, A.Duran 9042 & M.Öztürk ( |
C. straminicephala | Turkey, Erzurum,Uzundere to Artvin, 1100 m, 26.07.2002, A.Duran 6048 & M.Sağıroğlu ( |
Nuclear DNA of silica gel dried leaf samples were extracted according to the instructions of the Nucleon phytopure plant DNA extraction kit (RPN 8510, Amersham Life Science, England). For each sample, DNA was extracted from 100 mg of leaf. After concentrations were determined using an Eppendorf BioPhotometer, DNA samples were diluted to the working concentration of 25 ng/µL. To better quantify the DNA and to assess the quality of the DNA, samples were run on an agarose gel (0.9%), stained with ethidium bromide, against a DNA standard with known concentrations. Stock DNA was kept at -86 °C.
Of the 20 primers investigated during our initial screening, the primers that gave the most informative patterns (in terms of repeatability, scorability, and the ability to distinguish between varieties) were selected for fingerprinting. The characteristics of the primers used are given in Table
Primer | Primer sequence | Tm (°C) | Size (bp) | GC (%) | Number of polymorphic bands |
---|---|---|---|---|---|
ISSR F1 | GAG CAA CAA CAA CAA CAA | 49.1 | 18 | 38.9 | 13 |
ISSR F2 | CTC GTG TGT GTG TGT GTG T | 56.7 | 19 | 52.6 | 11 |
ISSR F3 | AGA GAG AGA GAG AGA GCG | 56 | 18 | 55.6 | 14 |
ISSR F4 | AGA GAG AGA GAG AGA GTG | 53.7 | 18 | 50 | 12 |
ISSR F5 | AGA GAG AGA GAG AGA G | 49.2 | 16 | 50 | 10 |
ISSR F6 | CCA CCA CCA CCA CCA | 53.3 | 15 | 66.7 | 13 |
ISSR F7 | ACA CAC ACA CAC ACA C | 49.2 | 16 | 50 | 12 |
Each reaction contained 2.5 mM MgCl2, 10 mM Tris-HCl (pH 8.8), 50 mM KCl; 0.8% Nonidet P40, 200 mM of each of dNTP, 0.5 mM primer, 25 ng DNA template and 0.4 units of Taq DNA Polymerase (Bioron, Germany) in a final reaction volume of 25 µl. After a pre-denaturation step of 3 minutes at 94 °C, amplification reactions were cycled 40 times at 94 °C for 1 minute, at annealing temperature (Table
Amplifications were repeated at least twice at different time periods for each primer using the same reagents and procedures.
Amplified fragments were visualized under a UV transiluminator and photographed using a gel documentation system (Vilbert Lourmat, Infinity model). All of the amplified fragments were treated as dominant genetic markers. Each DNA band generated was visually scored as an independent character or locus (1 for presence and 0 for absence). Qualitative differences in band intensities were not considered. Every gel was scored in triplicate (independent scorings) and only the fragments consistently scored were considered for analysis. A rectangular binary data matrix was prepared and all the data analysis was performed using the Numerical Taxonomy System, NTSYS-pc version 2.1 (Applied Biostatistic, Exeter Software, Setauket, New York, USA).
In cluster analysis of the samples, the unweighted pair-group method with the arithmetic mean (UPGMA) procedure was followed (
Silica gel dried plants collected from 19 different natural habitats were taken to the laboratory. The total number of species collected and used in the phylogenetic analysis was 19. DNA extractions were first attempted using a standard 2X CTAB method. Due to the poor DNA quality produced by the CTAB procedure, a commercial kit (Nucleon phytopure) was used in all isolations and repeated extractions were conducted whenever necessary.
From an initial screening of 20 ISSR primers, seven primers revealed high levels of polymorphisms. These primers generated 85 highly polymorphic fragments that were consistently amplified in repeated experiments conducted on separate dates. The GC percentages of the selected primers were within the range of 38.8–66.7%. The characteristics as well as the sequences of the primers revealing a polymorphism are shown in Table
Representative agarose gels where PCR products were amplified with the primers ISSR F5 (highest number of polymorphic bands, top) and ISSR F5 (lowest level of polymorphic bands, down). 1 Serratula tinctoria 2 Klasea quinquefolia 3 K. oligocephala 4 K. kotschyi 5 K. serratuloides 6 K. erucifolia 7 K. lasiocephala 8 K. cerinthifolia 9 K. grandifolia 10K. radiata subsp. radiata11 K. hakkiarica 12 K. haussknechtii 13 K. coriaceae 14K. radiata subsp. radiata15 K. kurdica 16 K. bornmuelleri 17 Centaurea ptosimopappoides 18 C. straminicephala 19 Jurinea cataonica, M: marker.
A total of 15 Klasea, 1 Serratula, 1 Jurinea and 2 Centaurea taxa were used in the scoring analysis. The Jurinea and Centaurea taxa, which were used as the out-group, formed a cluster that was distinct from the Klasea and Serratula cluster in the constructed dendogram. Furthermore the Klasea and the Serratula taxa form clearly separate clusters among themselves (Figure
Dendrogram showing genetic relationship of Klasea, Serratula, Centaurea and Jurinea species as shown using inter simple sequence repeats. (Serratula tinctoria, Klasea quinquefolia, K. oligocephala, K. kotschyi, K. serratuloides, K. erucifolia, K. lasiocephala, K. cerinthifolia, K. grandifolia, K. radiata subsp. radiata, K. hakkiarica, K. haussknechtii, K. coriaceae, K. radiata subsp. radiata, K. kurdica, K. bornmuelleri, Centaurea ptosimopappoides, C. straminicephala, Jurinea cataonica)
The Klasea radiata subsp. radiata and Klasea radiata subsp. biebersteiniana taxa were observed to be very closely positioned in the dendogram. The K. kotschyi, K. hakkiarica and K. haussknechtii taxa, which have similar leaf characteristics, were also correlated in terms of their molecular data and were located in the same sub-cluster. K. coriaceae is a taxon, which spreads over distinct areas in the Eastern Anatolia and is distinguished owing to its distinctive height. The Klasea oligocephala and K. kurdica clustered together. The two taxa are very similar also in terms of morphological properties such as leaf, capitulum and pollen characteristics. K. serratuloides taxon has the largest capitulum among the genus. It has a similar profile to the K. radiata subspecies with respect to the leaf characteristics and the similarity of overspreading areas. Klasea serratuloides and K. radiata were also positioned close to one another on the dendogram.
K. lasiocephala is distinguished within the genus by its very short stems or the absent stems.. K. lasiocephala differs morphologically from other Klasea taxa in having absent or reduced stems and that it is also somewhat genetically distinct from other Klasea taxa, as the sole taxon in the cluster in which it is placed. The K. erucifolia and K. grandifolia taxa have similar leaf characteristics and were also located in the same sub-cluster owing to their molecular characteristics. K. bornmuelleri taxon does not have a morphologically close relative in the genus. Its position on the dendogram confirmed this classification. K. cerinthifolia is distinguished by its yellow flowers and semiamplexicaul leaf structure and was also molecularly identified to be distinct. All these findings were consistent with the morphological classifications made in the Flora of Turkey (
The inspection of the dendogram indicated that molecularly similar taxa were also morphologically similar. This separation was also shown in the PCA plot (Figure
Principal co-ordinate analysis (PCA) of Klasea, Serratula, Centaurea and Jurinea species. 1 Serratula tinctoria 2 Klasea quinquefolia 3 K. oligocephala 4 K. kotschyi 5 K. serratuloides 6 K. erucifolia 7 K. lasiocephala 8 K. cerinthifolia 9 K. grandifolia 10K. radiata subsp. radiata11 K. hakkiarica 12 K. haussknechtii 13 K. coriaceae 14K. radiata subsp. radiata15 K. kurdica 16 K. bornmuelleri 17 Centaurea ptosimopappoides 18 C. straminicephala 19 Jurinea cataonica
Our study has demonstrated that ISSR is a powerful tool in resolving the genetic relationships within problematic taxonomical entities. In conclusion, the morphologically close taxa were, in the molecular aspect, also located in the same clade. The genera used as out-groups (Serratula, Jurinea, and Centaurea) were clearly separate from the genus Klasea. According to our knowledge, this is the first report on the use of ISSR in Klasea.
This work was financially supported by The Scientific and Technical Research Council of Turkey (TUBITAK Project No. TBAG-109T243) and The Scientific Research Coordination Center (BAP) of Selcuk University (Project No. 11401048).