Phylogeography of anterastes serbicus species group (orthoptera, Tettigoniidae): phylogroups correlate with mountain belts, but not with the morphospecies.
Discovering and describing species diversity is a fundamental part of evolution and biodiversity studies, and this task is far from complete, especially among poorly studied groups, i.e., almost all invertebrates. This is certainly the case for the bush cricket genus Anterastes Brunner von Wattenwyl, 1882 (Orthoptera, Tettigoniidae), particularly for the A. serbicus group distributed in Anatolia and the Balkans.
The complex includes three described species, Anterastes serbicus Brunner von Wattenwyl, 1882, A. burri Karaba?, 1951 and A. antitauricus Ciplak, 2004a. Despite the small number of species, and an effort to discover diagnostic morphological characters distinguishing members of this complex from each other and from others in the genus (Ciplak 2004a), the A. serbicus species complex appears to represent evolutionarily distinct lineages. Additional work is clearly needed to delimit the several members of these lineages, on morphological criteria alone, to test the morphospecies approach, to reveal the number, phylogenetic relationships and geographic distribution of evolutionary lineages, to propose historical biogeography and demography scenarios that accommodate the genetic variation observed among and within lineages and to detect, if any exists, hidden cryptic species diversity.
More accurate phylogenies have allowed authors to test previously suggested taxa (Elmer et al. 2007, Stephenes & Wiens 2009, Cicconardi etal. 2010) or to reconstruct their speciation in correlation with their historical ranges (Avis 2009). There is an overwhelming literature on phylogeography of West Eurasian lineages which suggests that present European taxa evolved mainly from the ancestral stocks inhabiting glacial refugia such as Iberia, Italy and the Balkans (Schmitt 2007). Although Anatolia (Asia Minor or Asian Turkey) has been considered to be a fourth glacial refugium in some studies (Hewitt 1996, 2004; Ciplak 2004a, 2004b, 2008), only a few regional or peninsula-wide phylogenetic studies showing different phylogeographic patterns in different taxa have been conducted (Bohlen et al. 2006; Dubey et al. 2007; Gunduz et al. 2005, 2007; Veith et al. 2003, 2008).
Any further phylogenetic study that may provide information on the phylogeographic history of the Anatolian fauna and flora is of interest. This is because 1) Anatolia has a very active tectonic history through the Neogene in the Mediterranean (Bozkurt 2001) which possibly has repeatedly led to dispersal/vicariance of the inhabiting populations (Kosswig 1955, Ciplak 2004b, Hberk 2002, Hberk et al. 2003), creating major genetic breaks within species. 2) It is an important crossroad at present (possibly also in the past) for the faunal/floral exchange between Europe, Asia and even Africa, via the Arabian Peninsula. 3) It was not covered by an ice sheet (the 52nd latitude is considered the border of permafrost; Hewitt 2000) during glacial periods; being located just south of the border of the permafrost, it may have played a refugial role. 4) It has a very heterogeneous topography/climate that leads to habitat fragmentation, geographic barriers and a biodiversity that is not uniformly distributed throughout the peninsula.
Large scale geographic definitions may indicate the importance of Anatolia as a general refugium, but these do not allow us to elucidate the importance of this region as a system of refugia and speciation areas. Firstly, Anatolia cannot be considered as a homogenous refugium, because most taxa could not spread into the whole Anatolian landmass; rather each taxon could extend its distribution only to particular districts, using suitable corridors during range shifts from the northern territories or vice versa (Hewitt 1996, Tarkhnishvili et al. 2000, Veith et al. 2003, Seddon et al. 2001, Ciplak 2004a). Secondly, large scale geographical definitions focus mainly on latitudinal range changes and these are trivially related to vertical changes. Anatolian topography might have played an important role in Pleistocene range changes both by constituting distribution corridors for latitudinal dispersals, and by providing suitable habitats during altitudinal shifts (Ciplak 2004a). Thirdly, the southern geographic location of the Anatolian refugium requires paying special attention to its role in the matter of rear-edge concept. According to the rear-edge concept (Hampe & Petit 2005), the cold-preferring lineages or populations should be fragmented from the south of their range. The expectation is that such forms disappeared in lowlands, and were confined to some highlands, as isolated populations, the fate of which is erosion in genetic diversity (Hewitt 2004) or drift (Knowles & Richards 2005). Thus, designation of a meaningful phylogeographic study or testing validity of present hypotheses requires considering this dual refugial role of Anatolia in combination with the influence of its extremely complex topography.
There are a few studies focusing on the horizontal distribution patterns around Anatolia (Rokas et al. 2003, Gunduz et al. 2007, Challis et al. 2007, Stone et al. 2007); however, the generality of these patterns needs to be tested using different lineages. Anatolia harbors a highly endemic biodiversity at different taxonomic levels (for Orthoptera see, Ciplak & Demirsoy 1995, 1996; Ciplak 2003) and radiation of the endemic mountainous forms is suggested to occur mainly due to climatic shifts of the Quaternary (see Ciplak 2003, 2004a, 2008). Possibly the vertical range changes were more important for the species or lineages endemic to Anatolia.
The present study, ofa four-species lineage, mainly aims to address the following questions: 1) What is the portion of Anatolian forms in establishment of present populations of nonrefugial northern areas and how far could these forms have extended their range during the Holocene? 2) Are there lineages, reflecting past vertical range changes (in altitudinal heterogeneities) or fragmentations where this kind of range change could be defined by using the phylogeography of such lineages? Anatolia is one of the southernmost refugia and the Mediterranean Sea constitutes a strong barrier in the south for terrestrial species/populations. Therefore, Anatolia is a perfect area to test the rear-edge concept. 3) Is the rear-edge concept applicable to Anatolian lineages? The starting point for conservation activities is defining evolutionarily independent units and their genetic diversity. Speaking for Anatolian taxa, this is especially important for cold-preferring taxa isolated on its southern highlands or sky islands (Ciplak 2008). 4) Is each highland population isolated and is the isolation more prominent for those in the south?
The Anatolian-Balkanian genus Anterastes has recently been revised by Ciplak (2004a) who presented a morphology-based phylogeny and a biogeographic account. The genus includes nine species that are endemic to Anatolia, and one (A. serbicus) that also extends its range to the Balkanian highlands. As presently known, their populations prefer cold climates and are present on alpine/ subalpine regions of highlands around/above 1500 m altitude. Some of the species are only known from a single upland, while others are known from more than one highland, but each also occurs as an isolated population. For example, A. serbicus that is present both in North Anatolia and the Balkans has not been recorded from lowlands of the Balkans and Anatolia (Ciplak 2004a, Ciplak unpublished observations). Since this is the distribution pattern/habitat preference observed in all species, the genus is a model candidate lineage to answer the questions just listed.
On the other hand, as mentioned, species of the genus differ in minor diagnostic characters (Ciplak 2004a). One of the species group in the genus is the A. serbicus group which consists of three species: A. serbicus (Anatolia+Balkans), A. burri (Northwest Anatolia) and A. antitauricus (South Anatolia). Each of these species is suitable for testing the above statements, since they are presently known from three or more different altitudes. On the other hand, A. tolunayi, a species not in the A. serbicus group, but very close through shared similarities (Ciplak 2004a), known from a single altitude, was also included in the study in order to locate the species of the A. serbicus complex within a phylogenetic frame. Using 16S rDNA sequences obtained from different populations of these four species we produced a phylogeny with the aim of testing the above addressed questions.
Materials and Methods
Sampling and DNA extraction.--Twelve populations of four ingroup species (A. serbicus, A. burri, A. antitauricus and A. tolunayi) and one population of an outgroup (A. uludaghensis), were sampled from Anatolia and the Balkans (Tables 1, 2; Fig. 1). Morphological species identification of the nominal species was based on differences given by Ciplak (2004). Total DNA was extracted from 153 adult specimens preserved in 96% ethanol, kept at +4[degrees]C (25 belonging to A. serbicus, 56 to A. burri, 43 to A. antitauricus, 20 to A. tolunayi and 8 to outgroup A. uludaghensis) with the salt-extraction method (Aljanabi & Martinez 1997). The 16S rDNA gene was amplified using the universal primers LR-J-12887 and LR-N-13398 described by Simon et al. (1994). Amplifications were performed in a 50-ul volume containing 0.3 [micro]l of each primer (100 [micro]M), 1 [micro]l dNTP mix (10 mM), 2 [micro]l 50 mM MgCl2, 5 [micro]l 10X Platinum PCR buffer (containing 200 mM Tris-HCl [pH 8.4], 500 mM KCl), 1.25 U Platinium Taq DNA polymerase (Invitrogen), and 0.5-1 [micro]l of 50-70 ng template DNA. Temperature cycling was carried out in an Eppendorf Mastercycler Personal. Cycling conditions were 1 min denaturing at 95[degrees]C; (45 s at 95[degrees]C, 40 s at 49[degrees]C, and 50 s at 72[degrees]C) x 35. PCR products were purified and sequenced in both directions. The sequencing products were loaded onto an ABI 3730 XL automated sequencer.
Phylogenetic analyses.--The sequences were aligned manually in BioEdit [version 5.09] (Hall 1999). DnaSP [version 5] (Librado & Rozas 2009) was used to determine unique haplotypes and to calculate their frequency. Sequences were deposited in GeneBank. The phylogenetic signal in the data partitions was estimated by maximum likelihood mapping method (Strimmer & von Haeseler 1997) using the program TREE-PUZZLE [version 5.2] (Schmidt et al. 2002). Initially, the evolutionary relationships between haplotypes were examined with the software Network version 184.108.40.206 (http: //fluxus-engineering.com) (Bandelt et al. 1999). Later, the aligned sequences were analyzed by maximum parsimony (MP) and maximum likelihood (ML), as implemented in PAUP [version 4.0b 10] (Swofford 2000) as well as by using Bayesian inference of phylogeny as implemented in MrBayes [version 3.1.2] (Ronquist & Huelsenbeck 2003). Bayesian search was carried out using four simulations of Markov chains, five million generations and sampling every 100th generation. The software TRACER [version 1.5] (Rambout & Drummond 2003) was used to examine the parameters and determine the number of trees needed to reach stationarity (burn-in). Bayesian posterior branch probabilities were obtained by taking the majority rule consensus of the sampled trees, excluding the first 5000 trees as burn-in. Prior to calculations all ambiguously aligned positions and gap-loaded positions were excluded from the dataset. The MP analyses were carried out 10 times with the heuristic search approach using the TBR algorithm. The confidence of branching was assessed using 100 nonparametric bootstrap resamplings (Felsenstein 1985). The parameters and the best-fit model were estimated using MODELTEST [version 3.06] (Posada & Crandall 1998). The selected model was implemented in the ML analysis performed with PAUP. The ML tree search was conducted using the heuristic search approach, the 'as is' addition replicate. Nonparametric bootstrapping (Felsenstein 1985) was used to evaluate the support of nodes based on 1000 pseudoreplicates.
[FIGURE 1 OMITTED]
Estimation of divergence times.--To estimate divergence times and to correlate them with the historical events, a generalized clock was applied. Divergence dates among the main mitochondrial clades were estimated using a Bayesian Markov Chain Monte Carlo (MCMC) approach, as implemented in BEAST [version 1.5.2] (Drummond & Rambaut 2007). Prior to running BEAST all ambiguously aligned positions and gap-loaded positions were excluded from the dataset. Several related analyses were run for the same alignment with the model selected by MODELTEST. The haplotype network in Fig. 3 and the phylogenetic tree in Fig. 4 were considered to constrain the analysis with regard to the topologies sampled. Dates of the divergence were inferred using a relaxed molecular clock, following the uncorrelated relaxed lognormal clock as implemented. The dataset was calibrated by general substitution rate (1.4%, per lineage, per 700 kyr) as suggested by Brower (1994) for 16S rDNA in Laupala (Orthoptera). All analyses were performed assuming a Yule process of diversification. BEAST was run for 20 million generations, sampling every 1000 and the convergence to stationarity and the effective sample size (ESS) of model parameters were checked by TRACER. The maximum clade-credibility tree built with TREEANNOTATOR, discarding the initial 10% of samples as burn-in [FIG-TREE version 1.3.1] (Rambaut 2008) was used to visualize the results, including the confidence intervals.
Analyses of genetic diversity.--The number of haplotypes, either without gaps (k) or with gaps (K), haplotype diversity (h), the number of polymorphic sites (S) and nucleotide diversity (n) were calculated for geographical populations and species separately. With unequal numbers of individuals, varying from 4 to 28 per population in the study, Pearson's correlation tests were conducted to test value independence of the number of haplotypes (k), haplotype diversity (h) and nucleotide diversity (n), against the number of individuals (N). To evaluate the significance of differentiation between populations based on pairwise FST an exact test (Reymond & Rousset 1995) was performed with a Markov Chain length of 10,000 and 1,000 burn-in steps. Analysis of molecular variance (AMOVA) (Weir & Cockerham 1984) was conducted to measure the extent to which genetic variance is assigned to the three hierarchical levels of organization: 1) between morphospecies, 2) between geographic populations within morphospecies and 3) within geographic populations. The statistical significance of AMOVA was evaluated using 10,000 permutations. All above analyses were performed using Arlequin 3.01 software (Excoffier et al. 2005) under the TrN+r model, since it is suggested that it is the more general case of the HKY (Hasewaga et al. 1985).
[FIGURE 2 OMITTED]
The outgroup A. uludaghensis is sampled from a single locality (it is so far known only from Bursa, Uludag Mt). Of the ingroup species, A. tolunayi was found from one locality (Izmir), A. serbicus (Sivas, Giresun and Kastamonu provinces of Turkey, and Bulgaria) and A. burri from four (Usak, Bursa, Bolu and Cankiri provinces) and A. antitauricus from three (Karaman, Konya and Afyon provinces) (Tables 1, 2; Fig. 1). From the sequence analyses, 41 haplotypes were identified among the 153 samples belonging to one outgroup and four ingroup Anterastes species: single haplotype for A. uludaghensis (outgroup, n=8) and A. tolunayi (n= 20), 12 for A. serbicus (n= 25), 18 for A. burri (n= 56) and 9 for A. antitauricus (n= 43) (Table 1). Two of these 41 haplotypes (one from A. antitauricus and the other from A. burri) are definable by the gaps in three sites. The phylogenetic and time-estimation analyses were made using 523 bp of 39 haplotypes (excluding the gaps), but a dataset of all was used in comparative genetic analyses. Of the 523 sites, 443 are constant, 80 are variable (50 of which are parsimony informative). There is no haplotype shared either by two or more different morphospecies or by two or more geographic populations.
Phylogenetic analyses.--Results of likelihood mapping are presented in Fig. 2. Poor dichotomic phylogenetic signal was detected in the dataset.
The percentage of the quartets suggesting a star-or network-like phylogeny is 17.7%, an amount close to threshold, indicating that data arenotreliableforadichotomicphylogenetic analysis (Lemeyet al. 2009). Thus, we initially performed ahaplotype network analysis. This analysis suggested 1IZ and 2 LIS to be distinctly diverged by 20 and 17 mutations from respectively most similar haplotypes. Remaining haplotypes constituted four clusters which are definable by geography: 1) Bulgarian cluster (A. serbicus), 2) Bursa-TIludag cluster (topotypical A. burri), 3) Black Sea Region cluster excluding Bolu (Cankiri-A. burri; Sivas, Giresun and Kastamonu-A. serbicus) and 4) Mediterranean + Aegean Anatolia plus Bolu province (from west part of Turkish Black Sea area) cluster (Bolu-A. burri, Karaman, Konya and Afyon- A. antitauricus). The haplotype 15KA occurs as ancestor of others in the last cluster.
MODELTEST, applied to the dataset, selected the HKY model (Hasegewa et al. 1985), according to hLRT with proportion of invariable characters I, of 0.6537, gamma correction [GAMMA], of 0.6851 and ti/tv ratio of 3.9628. This model was implemented in BI and ML analyses. Each of BI, MP and ML analyses resulted in trees with different topologies; however, there are clades supported in all trees (Fig. 4).
The single haplotype from the lzmirpopulation (1IZ,A. tolunayi) occurs as either an independent (BI, Fig. 4B) or abasal branch leading to all others (MP, ML; Fig. 4A, C). The latter clade accounts for all samples of A. serbicus, A. burri and A. antitauricus. The single haplotype (2US) from the Usak-Murat Mt (Aegean Anatolia) population (A. burri) branches off basally from the remaining in all analyses (Fig. 4).
[FIGURE 3 OMITTED]
For the next step, MP and BI analyses suggested a polytomy of several branches, however, there are three internal phylogroups well supported in both trees (Fig. 4A, B). One of them accounts for the seven haplotypes from the Bulgarian population (A. serbicus), a second for the eight haplotypes from the Bursa-Uludag population (topotypical A. burri) and the last for the nine haplotypes from Sivas, Giresun, Kastamonu (A. serbicus) and Cankiri (A. burri). Network analysis to the contrary (Fig. 3), the remaining 12 haplotypes (A. burri + A. antitauricus) occurred as independent branches of single or dual haplotypes in BI and MP analyses.
ML analysis resulted in a different topology. After two successive branches (1IZ and 2US) there is a dichotomy of two haplotype clades: 1) seven haplotypes of A. serbicus from Bulgaria and 2) the remaining. On the next step the nine haplotypes from Sivas, Giresun, Kastamonu (A. serbicus) and Cankiri (A. burri) constitute a phylogroup, while the others make a second. The crown clade includes two phylogroups: 1) eight haplotypes from the BursaUludag population (topotypical A. burri) and 2) 12 haplotypes from Mediterranean + Aegean Anatolia plus Bolu. All of the gene trees supported monophyly of their species clades (A. serbicus, A. burri and A. antitauricus), but not that of each morphospecies.
[FIGURE 4 OMITTED]
Estimation of divergence times.--Considering both haplotype network (Fig. 3) and phylogenetic analyses (Fig. 4), BEAST analysis was constrained with the following topologies: 1) outgroup and others, 2) the outgroup+1IZ and others, 3) the outgroup+HZ+2TJS and others, and 4) three terminal groups, the first of which to include seven haplotypes from Bulgaria, the second seven haplotypes from Bursa-Uluda? and the last 19 haplotypes originatingfrom 15 KA (see Fig. 3). The BEAST analysis, under a HKY+I+ T model, revealed that this segment of DNA showed a clock-like evolution (ucld.mean= 0.52, ucld.stdev = 0.27; corresponding to a rate of 2% per million years; Fig. 5).
From BEAST, we estimated the time to most recent common ancestor (TMRCA) forthe entire ingroup and internal clades as follows: 1) to ingroup origin 2.27 Myr (95% HPD interval = 1.534-3.059), 2) to three ingroup species clade (including A. serbicus, A. burri and A. antitauricus) 2.01 (95% HPD interval = 1.385-2.659) and 3) to the clade including Bulgarian plus its sister Anatolian phylogroup 1.56 (95% HPD interval = 1.077-2.094). The first corresponds to Late Pliocene/Early Pleistocene and the other two to Early Pleistocene. The diversification of terminal phylogroups (those corresponding to haplotype clusters suggested by haplotype network analysis) correspond to a time prior to Gunz glacial (for TMRCA of each clade see Fig. 5).
Analyses of genetic diversity.--Strong geographicstructuring of genetic variation was observed at the 16S rDNA marker. All of the haplotypes were unique to a single geographic population. Diversity indices in A. tolunayi were zero, since 20 samples, represent just a single haplotype. For the three ingroup species, the haplotype diversity [h) varies between 0.8239-0.8933 and the nucleotide diversity ([xi]) between 0.007-0.028 (Tablel). Diversity indices for the geographic populations were calculated separately, since monophyly of morphospecies was not supported and a strong geographic structuring was found. These indices are zero for the Izmir (A. tolunayi, 20 samples), Usak (A. burri, 16 samples), Sivas (A. serbicus; 4 samples) and Giresun (A. serbicus, 5 samples) populations, since each have a single haplotype. This is the case for the Karaman population (A. antitauricus) also, since only a single individual was available for study (from our observation, this population seems to have gone extinct very recently when its restricted habitat was destroyed).
In the other geographic populations haplotype diversity [h) varies between 0.0053-0.9048 and nucleotide diversity ([xi]) between 0.0019-0.0136 (Table 2). The highest haplotype diversities were found in the Bursa, Kastamonu and Bulgarian populations, while nucleotide diversity was highest in the Bulgarian population. There is no evidence that heterogeneous levels of genetic variability were produced by the unequal sample sizes for each population (Pearson's correlation test between diversity indices and sample size: r= 0.196 and P= 0.564 for number of haplotypes, r= 0.108 and P= 0.752 for haplotype diversity, r= - 0.007 and P= 0.983 for nucleotide diversity).
A significantly high differentiation was observed for both inter-population, within species, and interhaplotypes, within population comparisons. Based on an AMOVA of the entire 16S rDNA dataset, 22.03% of genetic variation was found to be partitioned between morphospecies, 67.28% among sampling sites inside each species, and 10.68% inside each sampling site, the last two being statistically significant (Table 3). Possibly this is because each geographic population has its unique haplotypes (or there are no haplotypes shared by two or more populations). Mean interspecies pairwise FST values were in the range of 0.1-0.52. A. tolunayi is significantly divergent from all others. The interpopulations differentiation is even higher. Mean interpopulation pairwise FST range between 0.301.0. FST values are below 0.5 in a few comparisons, indicating highly diverged geographic populations (Table 4).
The first and possibly most striking aspect of the data is that the gene tree does not support morphospecies (see Ciplak 2004a for morphological data). Of the four species, A. tolunayi seems to constitute a genetic clade, while the other three morphospecies, A. serbicus, A. burri and A. antitauricus, are para(poly)phyletic. Consistent with phylogenetic relationships, AMOVA analysis suggests that genetic divergence between species is not significant, while that between populations is obviously high. Furthermore, several haplotypes diverge from each other by several substitutions and nearly all are joined to each other by hypothetical ancestral haplotypes (see Fig. 3).
There are two possible explanations of this situation. First, the ingroup species, other than A. tolunayi, may be considered as a single species. The small morphological differences between the three species in the A. serbicus group (see Ciplak 2004a) may support this assumption. Second, this three-species clade may include several undescribed, morphologically cryptic, species. Present data support the second assumption, since 1) there are no clues in the data indicating gene flow between morphospecies, since there is no haplotype shared by two or more of them; 2) there are strongly diverged geographic populations (pairwise FST values are mostly higher than 0.5 and run up to 1.0 in most of the comparisons), each harboring only its unique haplotypes; 3) there are some internal phylogroups, such as Bulgarian (A. serbicus), Usak-Murat Mt (A. burri), Bursa-Uludag Mt (topotypical A. burri), which seem to have been reproductively isolated since the Mid-Pleistocene (Figs 3, 5).
It is suggested that these animals are present only in highland habitats with restricted ecological conditions (Ciplak 2004a). If this habitat preference is the case in their history, it can be assumed that living in similar conditions may lead to a similar morphology because of similar selection pressures (Elmer et al. 2007, Stephenes & Wiens 2009, Cicconardi et al. 2010).
Although the Usak-Murat Mt population is genetically the most divergent population of these three species, it is a typical A. burri in morphology. This is also the case for several other populations (e.g., Bulgarian, A. serbicus), and all these are consistent with the second assumption. However, it should be kept in mind that gene trees are not species trees (Edwards 2009), thus a nomenclatural decision at this stage will be premature,
A second aspect of the results to be evaluated is that analyses of the data by different procedures did not result in a robust phylogeny. For example, haplotypes of A. antitauricus from a single population (Konya) occur as independent branches in all phylogenetic analyses. Contrary to the phylogenetic analyses, these haplotypes differ from each other in a few base pairs and they constitute a single cluster originated from the haplotype found in the Karaman population in the network analysis (Fig. 3). Similarly, the four haplotypes from Cankiri population either constitute a phylogroup with those from Kastamonu, Giresun and Sivas (BI and ML trees) or independent branches (MP tree).
One reason underlying this inefficiency relates partially to poor phylogenetic signs in mtDNA sequences due to A-T bias (Simon et al. 1994). Results of the likelihood mapping analysis and the percentage of A-T (69.58%) and that of G-C (30.42%) support this assumption. However, there are signs in the data indicating another reason: if the local populations were independent evolutionary units with no or limited gene flow through their divergence from a polymorphic ancestral population, this could have led to loss of ancestral haplotyes. Such mosaic evolution may cause a limited number of synapomorphic sites (Simon et al. 1994).
Since such a process is suggested to have been more severe in recent radiations, the unresolved upper parts of the tree, particularly the polytomy of haplotypes found in A. antitauricus, are consistent with this statement. Additionally, the geographic populations exhibit strong genetic isolation, as suggested from the compelling finding that this phylogeny reflects relatively well, i.e., corresponds to, geography (i.e., haplotypes from the same or nearby localities are always genealogically exclusive and very high pairwise FST values exist among populations). This is as might be expected for species with limited dispersal capabilities (e.g., Bond et al. 2001, Hedin & Wood 2002). All these issues require discussing the phylogeography of the group in detail.
Phylogeography and divergence times.--We now address a few questions to be answered in a biogeographic perspective in the light of the data obtained. Anatolia is suggested as a glacial refugium, playing an important role in the establishment of fauna of nonrefugial northern areas during the Holocene (Taberlet et al. 1998, Hewitt 2004).
The first question is: 'what is the allotted portion of cold-preferring Anatolian forms extending to nonrefugial northern areas during the Holocene, or during the last four glacial periods?' The Balkanian population of A. serbicus is the only representative of this lineage extending beyond Anatolia. La Greca (1999) suggested that this dispersal occurred recently during the Holocene. Present data contradict this statement. Bulgarian haplotypes constitute a distinct and highly divergent phylogroup which share the last common ancestor with a sister Anatolian clade 1.56 My ago (see Fig. 5). This corresponds to a time much earlier than Gunz glacial.
[FIGURE 5 OMITTED]
Species of Anterastes prefer cold mountainous habitats, hence a northward range expansion during warm interglacial periods, as suggested for such forms (Hewitt 1996; Ciplak 2004a, 2008), is an expectation. However, the range changes during glacial/interglacial periods are possibly related to dispersal ability. Anterastes are flightless insects and have strict habitat requirements (Ciplak 2004a): both may indicate reduced dispersal ability.
Although most of the hypotheses explaining the high biodiversity in the region refer to Anatolia as a glacial refugium (Taberlet et al. 1998; Hewitt 2004; Ciplak 2004a, b; Ciplak 2008), evidence from this study suggests that the Anatolian Peninsula is also home to old endemic biota. This situation is attested to by the high rate of specific endemism in several groups, such as Tettigoniinae with 83% endemism (Ciplak 2003), Pamphagidae with 75% endemism (Ciplak & Demirsoy 1996), and plants with 33% endemism (Davis 1965-1985, Medail & Quezel 1997, Medail & Diadema 2009). Evidence from phylogenetic studies and the presence of deep genetic divergence in populations (Cooper et al. 1995, Seddon et al. 2001, Rokas et al. 2003, Gunduz et al. 2007, Murienne et al. 2010) also suggest long biogeographic isolation. In this context, the long and complex palaeogeographic history, high habitat heterogeneity, and great topographic and climatic variations of the Anatolia are of fundamental importance in explaining its biodiversity.
The pattern discussed above considers the horizontal range changes in a latitudinal gradient. Although present data do not indicate such a pattern during the last four glacial periods, climatic shifts might have led to vertical range changes of cold-preferring Anatolian forms in an altitudinal gradient (Ciplak 2004a, 2008). Thus, our second question was 'are there any clues in the data reflecting vertical range changes during warming/cooling cycles of the Pleistocene?' Since the present populations were sampled during the Holocene warm period, an expectation for these mountainous forms is their restriction to uplands and disappearance from lowlands. Up to the present they have not been found in lowlands and each population is geographically isolated in its respective highland. So the distribution patterns fit the expectation.
The further consistent signs in the data are as follows. First, phylogroups or genetic clusters do not support monophyly of morphospecies, but they correlate well with altitudinal chains (see Fig. 3). Second, the populations from highlands within the same altitudinal chains, fall in the same phylogroups. Third, from our observations each local population is geographically isolated, as represented by unique haplotypes. Fourth, they differ from each other in very high values of pairwise FST; however, the lowest FST values were observed between populations of the same highland chains.
Another expectation related to restriction to uplands is reduction in population size and so loss of genetic diversity (Knowles & Richards 2005). Low haplotype diversity and high nucleotide diversity are suggested as clues to reduction in genetic diversity (Zhao et al. 2008). There are several populations exhibiting this pattern (Table 2). In the populations from Izmir, Usak, Giresun and Sivas, diversity indices are zero. Also, in Afyon, Konya and Cankiri populations, when compared to others, the nucleotide diversity is higher and the haplotype diversity is lower. Thus, both distribution pattern and population genetic parameters indicate a recent reduction in population size consistent with the expectation.
Anatolia is one of the most southern of refugia, and the Mediterranean Sea constitutes a strong barrier in Anatolia's south, limiting ranges of terrestrial lineages, whose ranges were possibly affected by climatic changes of the Pleistocene. Thus, this peninsula constitutes an ideal area in which to test the rear-edge concept (Hampe & Petit 2005). Although the rear-edge concept applies to fragments of a single species, it may be applicable to lineages radiated in a similar way due to past fragmentations. Thus, our third question: 'is the rear-edge concept applicable to Anatolian lineages?'
From the perspective of the concept (Hampe & Petit 2005), the fragmented populations on the southern edge are small, with reduced genetic diversity. In this case, the diversity indices should be lower for southern populations than for northern ones. From the latitudinal locations and the climate/vegetation characteristic of Anatolia (Ciplak 2003), of the 12 populations studied, the Konya, Karaman, Afyon, Izmir and Usak populations may be considered as southern and the others as northern. A comparison of the diversity indices (k, K, h and n) supports the prediction. The diversity indices are zero in Izmir and Usak and are lower in Konya, Afyon and Karaman, than in the north, although the highest sample size is from the Konya population. The highest diversity indices are higher for those populations from the north (Table 2).
Conservation activities require definition of the conservation units and the documenting of their genetic diversities to assess the plans for their future. Definition of conservation units can either be done on a taxonomic or evolutionarily independent unit basis (Frankham et al. 2006). For Anterastes, a taxon-based definition does not seem to work, since each species, whether valid or not when our results are considered, consists of populations isolated both geographically (on their sky island) and genetically (each has haplotypes not shared by any other). Thus, the conclusion is that each has its own evolutionary trend and needs to be considered separately in conservation plans.
Another issue relating to conservation is the genetic diversity per conservation unit, since the genetically homogenized populations are assumed to be more prone to extinction (Frankham et al. 2006). Haplotype and nucleotide diversities are zero in 5 of 12 populations, indicating a highly reduced genetic diversity. Combining the diversity values with the range width, suggests that some of these populations, especially that of Izmir (A. tolunayi) and Usak (A. burri), are threatened. Ciplak (2008) suggested that cold-preferring taxa/populations of Anatolia are in general, isolated on highlands, particularly in its southern part. Thus, conservation plans for southern and northern populations require different considerations: an area conservation perspective for those in the north, but a conservation genetic perspective for those in the south. Present data are consistent with this suggestion and validate the statements. Moreover, this study provides a very good example of species-level "paraphyly" and adds to a growing body of literature suggesting that traditional methods for delineating species boundaries in different groups may underestimate actual species-level diversity.
We thank Dragan Chobanov (Bulgaria) for providing specimens from Bulgaria. Dr Deniz Sirin (Turkey) and M. Sait Taylan (Antalya) joined in one of the field studies. The young orthopterist Zehra Boztepe helped in preparation of the manuscript. Our research was supported by the Scientific and Technical Research Council ofTurkey (Project No: 107T462). Studies were carried out in laboratories in the Department of Biology, Akdeniz University, and the paper was supported by the Akdeniz University Research Fund.
Submitted April 9, accepted May 31, 2010
Aljanabi S.M., Martinez I. 1997. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Research 25: 692-4693.
Avis J.C. 2009. Phylogeography: retrospect and prospect. Journal of Biogeography 36: 3-15.
Bandelt H.J., Forster F., Rohl A. 1999. Median-Joining Networks for Inferring Intraspecific Phylogenies. Molecular Biology and Evolution 16: 37-48.
Bohlen J., Perdices A., Doadrio I., Economidis P.S., 2006. Vicariance, colonisation, and fast local speciation in Asia Minor and the Balkans as revealed from the phylogeny ofspined loaches (Osteichthyes; Cobitidae). Molecular Phylogenetics and Evolution, 39: 552-561.
Bond J.E., Hedin M.C, Ramirez M.G., Opell B.D., 2001. Deep molecular divergence in the absence of morphological and ecological change in the Californian coastal dune endemic trapdoor spider Aptostichus simus. Molecular Ecology 10: 899-910.
Bozkurt E. 2001. Neotectonics of Turkey--a synthesis. Geodinamica Acta 14: 3-30.
Brower A.V.Z. 1994. Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from pattern of mitochondrial DNA evolution. Proceedings National Academy of Sciences United States of America 91: 6491-6495.
Challis R., Mutun S., Nieves-Aldrey J.L., Preuss S, Rokas A., Aebi A., Sadeghi E, Tavakoli M, Stone GN. 2007. Longitudinal range expansion and cryptic eastern species in the western Palaearctic oak gallwasp, Andricus coriarius. Molecular Ecology 16: 2103-2114.
Cicconardi F., Nardi F., Emerson B.C., Frati F., Fanciulli P.P. 2010. Deep phylogeograpic and long-term persistence of forest invertebrates (Hexopoda: Collembola) in the northwestern Mediterranenan basin. Molecular Ecology 19: 386-400.
Ciplak B., Demirsoy A. 1995. Turkiye'de Ensifera (Orthoptera, Insecta) alttak?m?n?n endemizm ac?s?ndan de?erlendirilmesi. Turkish Journal of Zoology 19: 213-220.
Ciplak B., Demirsoy A. 1996. Caelifera (Orthoptera, Insecta) alt tak?m?n?n Turkiye'de endemizm durumu. Turkish Journal of Zoology 20: 241246.
Ciplak B. 2003. Distribution of Tettigoniinae (Orthoptera, Tettigoniidae) bush-crickets in Turkey: the importance of the Anatolian Taurus Mountains in biodiversity and implications for conservation. Biodiversity and Conservation 12: 47-64.
Ciplak B. 2004c. Biogeography of Anatolia: the marker group Orthoptera. Memorie della Societa Entomologica Italiana 82: 357-372.
Ciplak B. 2004a. Systematics, phylogeny and biogeography of Anterastes (Orthoptera, Tettigoniidae, Tettigoniinae): evolution within a refugium. Zoologica Scripta 33: 19-44.
Ciplak B. 2008. The analogy between glacial cycles and global warming for the glacial relicts in a refugium: a biogeographic perspective for conservation of Anatolian Orthoptera, pp. 135-163. In: Fattorini, S.: Insect Ecology and Conservation (Chapter 6). Research Signpost.
Cooper S.J., Ibrahim K.M., Hewitt G.M. 1995. Postglacial expansion and genome subdivision in the European grasshopper Chorthippus parallelus. Molecular Ecology 10: 2187-2198.
Davis P.H. 1965-1985. Flora of Turkey and the East Aegean Islands 1-9. 1965 (Vol. 1); 1967 (2); 1970 (3); 1972 (4); 1975 (5); 1978 (6); 1982 (7); 1984 (8); 1985 (9). Edinburgh.
Drummond A.J., Rambaut A. 2007. BEAST: Bayesian evolutionary analyses by sampling trees. BMC Evolutionary Biology 7: 214-221.
Dubey S., Cosson J.-F., Vohralik V., Krystufek B., Diker E., Vogel P. 2007. Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae). Journal of Evolutionary Biology 20: 1799-1808.
Elmer K.R., Davila A.J., Lougheed S.C. 2007. Cryptic diversity and deep divergence in an upper Amazonian leaf litter frog, Eleutherodactylus ockendeni. BMC Evolutionary Biology 7: 247-260.
ExcoffIer L, Laval G, Schneider S. 2005. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1: 47-50.
Felsenstein J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783-791.
Frankham R., Ballou J.D., Briscoe D.A. 2006. Introduction to conservation genetics. Cambridge University Press, London.
Garland T., Bennett A.F., Rezende E.L. 2005. Phylogenetic approaches in comparative physiology. Journal of Experimental Biology 208: 30153035.
Gunduz I., Rambau R.V., Tez C., Searle J.B., 2005. Mitochondrial DNA variation in the western house mouse (A4us musculus domesticus) close to its site of origin: studies in Turkey. Biological Journal of Linnean Society 84: 473-485.
Gunduz I., Jaarola M., Tez C., Yeniyurt C., Polly P.D., Searle J.B. 2007. Multigenic and morphometric differentiation of ground squirrels (Spermophilus, Sciuridae, Rodentia) in Turkey. Molecular Phylogenetics and Evolution 43: 916-935.
Hall T. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT. Nucleic Acids Symposium Series 41: 95-98.
Hampe A., Petit R. 2005. Conserving biodiversity under climate change: the rear edge matters. Ecology Letters 8: 461-467.
Hasegewa M., Kishino H., Yano T., 1985. Dating of the human-ape splitting by a molecular cloack of mitochondrial DNA. Journal of Molecular Biology 22: 160-174.
Hedin M.C., Wood D.A., 2002. Genealogical exclusivity in geographically proximate populations of Hypochilus thorelli Marx (Araneae, Hypochilidae) on the Cumberland Plateau of North America. Molecular Ecology 11: 1975-1988.
Hewitt G.M. 1996. Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of Linnean Society 58: 247-276.
Hewitt G.M. 1999. Post-glacial re-colonization of European biota. Biological Journal of Linnean Society 68: 87-112.
Hewitt G.M. 2000. The genetic legacy of the Quaternary ice ages. Nature 405: 907-913.
Hewitt G.M. 2004. Genetic consequences of climatic oscillation in the Quaternary. Philosophical Transactions Royal Society London B 359: 183-195.
Hrbek T., Meyer A. 2002. Closing of the Tethys Sea and the phylogeny of Eurasian killifishes (Cyprinodontiformes : Cyprinodontidae). Journal of Evolutionary Biology 16: 17-36.
Hrbek T., Kucuk F., Frickey T., Stolting K.N., Wildekamp R.H., Meyer A. 2003. Molecular Phylogeny and historical biogeography of the Aphanius (Pisces, Cyprinodontiformes) species complex of central Anatolia, Turkey.
Molecular Phylogeneticst an Evolution. 25: 125-137. Knowles L.L., Richards C.L. 2005. Importance of genetic drift during Pleistocene divergence as revealed by analyses of genomic variation. Molecular Ecology 14: 4023-4032.
Kosswig C. 1955. Zoogeography of the Near East. Systematic Zoology, 4: 50-74.
La Greca M. 1999. Il contributo degli Ortotteri (Insecta) alla conoscenza della biogeografia dell' Anatolia: la componente gondwaniana. Biogeographia 20: 179-200.
Lemey P., Salemi M., Vandamme A.M. (Eds) 2009. The Phylogenetic Handbook: a Practical Approach to Phylogenetic Analysis and Hypothesis Testing, 2nd Edition, pp. 381-404, Cambridge University Press, Cambridge.
Librado P., Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452
Medail F., Quezel P. 1997. Hot-spots analysis for conservation of plant biodiversity in the Mediterranean basin. Annals of Missouri Botanical Garden 84: 112-127.
Medail F., Diadema K. 2009. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. Journal of Biogeography 36: 1333-1345.
Murienne J., Karaman I., Giribet G. 2010. Explosive evolution of an ancient group of Cyphophthalmi (Arachnida: Opiliones) in the Balkan Peninsula.
Journal of Biogeography 37: 90-102. Posada D., Krandal K.A. 1998. MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817-818.
Rambout A., Drummond A.J. 2003. Tracer v1.3. Available form: http: //evolve.zoo.ox.ac.uk/.
Rambaut A. 2008. FigTree v1.2. Available from: <http://tree.bio.ed.ac.uk/ software/ FigTree/>.
Raymond M., Rousset F. 1995. An exact test for population differentiation. Evolution 49: 1280-1283.
Rokas A., Atkinson R.J., Webster, M.I., Csoka, G., Stone G.N. 2003. Out of Anatolia: longitudinal gradients in genetic divesity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae. Molecular Ecology 12: 2153-2174.
Ronquist F., Huelsenbeck J.P. 2003. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574.
Schmidt H.A., Strimmer K., Vingron M., von Haeseler A. 2002. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics 18: 502-504.
Schmitt T. 2007. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Frontiers in Zoology 4: 11-23.
Seddon J.M., Santucci F., Reeve N.J., Hewitt G.M. 2001. DNA footprints of European hedgehogs, Erinaceus europaeus and E. concolor: Pleistocene refugia, postglacial expansion and colonization routes. Molecular Ecology 10: 2187-2198.
Simon C., Frati F., Beckenbach A., Crespi B., Liu H., Flook P. 1994. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Annals Entomological Society of America 87: 651-702.
Stephens P. R., Wiens J.J. 2009. Bridging the gap between community ecology and historical biogeography: niche conservatism and community structure in emydid turtles. Molecular Ecology 18: 4664-4679.
Stone G.N., Challis R.J., Atkinson R. J., Csoka G., Hayward A., Melika G., Mutun S., Preuss S., Rokas A., Sadeghi E., Schorogge K. 2007. The phylogeographical clade trade: tracing the impact of human-mediated dispersal on the colonization of northern Europe by the oak gallwasp Andricus kollari. Molecular Ecology 16: 2768-2781.
Strimmer K., von Haeseler A. 1997. Likelihood mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proceedings National Academy Science, USA 94: 6815-6819.
Swofford D.L. 2002. PAUP* Phylogenetic Analysis Using Parsimony (* and other methods). v. 4.0 beta. M.A. Sinauer Associates, Sunderland.
Sword G.A., Senior L.B., Gaskin J., Joern A. 2007. Double trouble for grasshopper molecular systematics: intra-individual heterogeneity of both mitochondrial 12S-valine-16S and nuclear internal transcribed spacer ribosomal DNA sequences in Hesperotettix viridis (Orthoptera: Acrididae). Systematic Entomology 32: 420-428.
Taberlet P., Fumagalli L., Wust-Saucy A.G., Cosson J.F. 1998. Comparative phylogeography and postglacial colonization routes in Europe. Molecular Ecology 7: 453-464.
Tarkhnishvili D. Thorpe R.S., Arntzen J.W. 2000. Pre-pleistocene refugia and differentiation between populations of the Caucasian salamander (Mertensiella caucasica). Molecular Phylogenetics and Evolution 14: 414-422.
Weir B.S., Cockerham C.C. 1984. Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370.
Veith M., Schmidtler J.F., Kosuch J., Baran I., Seitz A. 2003. Palaeoclimatic changes explain Anatolian mountain frog evolution: a test for alternating vicariance and dispersal events. Molecular Ecology 12: 185-199.
Veith M., Lipscher E., Oz M., Kiefer A., Baran I., Polymeni R.M., Steinfartz S., 2008. Cracking the nut: geographical adjacency of sister taxa supports vicariance in a polytomic salamander clade in the absence of node support. Molecular Phylogenetics and Evolution, 47: 916-931.
Zhao L., Zhang J., Liu Z.J., Funk S.M., Wei F.W., Xu M.Q., Li M. 2008,Complex population genetic and demographic history of the salangid, Neosalanx taihuensis, based on cytochrome b sequences. BMC Evolutionary Biology 8: 201-228.
Battal Ciplak, Sarp Kaya and ?slam Gunduz
(BC)Department of Biology, Faculty of Art & Science, Akdeniz University 07058 Antalya, Turkey. Email: firstname.lastname@example.org (SK)Department of Biology, Faculty of Art & Science, Akdeniz University Antalya, Turkey (IG)Department of Biology, Faculty of Art & Science, Ondokuz May?s University Samsun Turkey
Table 1. Sampling localities and genetic diversity indices of four species. Shown from left to right are: the sample sizes (n), the number of segregating sites (s), the number of haplotypes without gaps (k), the number of haplotypes including gaps (K), haplotype diversity (h), nucleotide diversity ([pi]); NA, not included in population genetics analyses. Altitude Location Coordinates (m) Konya (KO) 41[degrees]03'24"N 33[degrees]43'05"E 1560 Afyon (AF) 38[degrees]28'08"N 30[degrees]22'69"E 1802 Karaman (KA) 37[degrees]01'30"N 33[degrees]16'45"E 1386 Bursa (BU) 40[degrees]07'18"N 29[degrees]08'57"E 2014 Bolu (BO) 40[degrees]35'48"N 31[degrees]48'06"E 2026 Usak (US) 38[degrees]56'64"N 29[degrees]37'89"E 1849 Cankiri (CA) 40[degrees]41'21"N 32[degrees]44'06"E 1595 Kastamonu (KS) 41[degrees]03'24"N 33[degrees]43'05"E 2006 Giresun (GS) 40[degrees]27'96"N 38[degrees]42'48"E 2294 Sivas (SI) 40[degrees]09'95"N 37[degrees]49'73"E 1999 Bulgaria (BUL) Rila Mt, Ossagovska Planina Mts 2100-2850 Izmir (IZ) 38 21'15"N 28[degrees]06'14"E 1400-1550 Location Species n s k K Konya (KO) A. antitauricus 28 7 4 4 Afyon (AF) A. antitauricus 14 5 3 4 Karaman (KA) A. antitauricus 1 - - - Bursa (BU) A. burri 15 10 8 8 Bolu (BO) A. burri 12 7 4 5 Usak (US) A. burri 16 0 1 1 Cankiri (CA) A. burri 13 6 4 4 Kastamonu (KS) A. serbicus 4 8 3 3 Giresun (GS) A. serbicus 5 0 1 1 Sivas (SI) A. serbicus 4 0 1 1 Bulgaria (BUL) A. serbicus 12 23 7 7 Izmir (IZ) A. tolunayi 20 0 1 1 Location h [pi] Konya (KO) 0.6825 +/- 0.0446 0.005386 +/- 0.003259 Afyon (AF) 0.0053 +/- 0.0032 0.001938 +/- 0.001549 Karaman (KA) NA NA Bursa (BU) 0.9048 +/- 0.0456 0.004890 +/- 0.003115 Bolu (BO) 0.6667 +/- 0.1409 0.002537 +/- 0.001903 Usak (US) 0 0 Cankiri (CA) 0.6795 +/- 0.1116 0.003506 +/- 0.002417 Kastamonu (KS) 0.8333 +/- 0.2224 0.008733 +/- 0.006459 Giresun (GS) 0 0 Sivas (SI) 0 0 Bulgaria (BUL) 0.7727 +/- 0.1276 0.013601 +/- 0.007729 Izmir (IZ) 0 0 Table 2. Geographic information and indices of genetic diversity. Shown from left to right are; the sample sizes (n), the number of segregating sites (s), the number of haplotypes without gaps (k), the number of haplotypes including gaps (K), haplotype diversity (h), nucleotide diversity (n). Number of species localities n s k K A. antitauricus 3 43 13 8 9 A. burri 4 56 40 17 18 A. serbicus 4 25 43 12 12 A. tolunayi 1 20 0 1 1 species h [pi] A. antitauricus 0.8239 +/- 0.0276 0.007016 +/- 0.004018 A. burri 0.8851 +/- 0.0290 0.026057 +/- 0.013177 A. serbicus 0.8933 +/- 0.0378 0.028000 +/- 0.014452 A. tolunayi 0 0 Table 3. AMOVA of four species of Anterastes genus (* P<0.0001; df, degree of freedom). Source of variation df Sum of Variance squares Components Among groups 3 501.980 1.99267 Among populations 8 488.218 6.08465 Within groups Within populations 132 127.524 0.96609 Total 143 1117.722 9.04341 Source of variation Percentage of F variance (%) Among groups 22.03 0.22035 Among populations 67.28 0.86298 * Within groups Within populations 10.68 0.89317 * Total Table 4. Pairwise [F.sub.st] values among localities of four species (only significant values are represented). A. tolunayi A. antitauricus A. burri IZ KO AF BU BO IZ KO 0.92897 AF 0.98194 0.51899 BU 0.96189 0.77195 0.84740 BO 0.97802 0.43700 0.50236 0.82696 US 1.00000 0.91775 0.97874 0.94306 0.97367 CA 0.96829 0.70382 0.82865 0.86166 0.79885 KS 0.97247 0.72153 0.83814 0.83376 0.80214 GS 1.00000 0.78142 0.92672 0.89376 0.90374 SI 1.00000 0.72396 0.90326 0.87376 0.87361 BUL 0.90657 0.79851 0.81546 0.82014 0.79189 A. burri A. serbicus US CA KS GS SI BUL IZ KO AF BU BO US CA 0.96310 KS 0.96194 0.59287 GS 1.00000 0.77733 0.61710 SI 1.00000 0.67363 1.00000 BUL 0.89259 0.78499 0.72483 0.78960 0.75497
|Printer friendly Cite/link Email Feedback|
|Author:||Ciplak, Battal; Kaya, Sarp; Gunduz, Islam|
|Publication:||Journal of Orthoptera Research|
|Date:||Jan 1, 2010|
|Previous Article:||Ultrasound avoidance behaviors in two species of Neoconocephalus (Orthoptera, Tettigoniidae).|
|Next Article:||Studies on neotropical phasmatodea x: redescriptions of aplopocranidium Zompro, 2004 and Jeremia Redtenbacher, 1908, with a survey of the tribe...|