Mitochondrial phylogeography of the long-eared bats (Plecotus) in the Mediterranean Palaearctic and Atlantic Islands

Long-eared bats of the genus Plecotus are widespread and common over most of the western Palaearctic. Based on recent molecular evidence, they proved to represent a complex of several cryptic species, with three new species being described from Europe in 2002. Evolutionary relationships among the different lineages are still fragmentary because of the limited geographic coverage of previous studies. Here we analyze Plecotus mitochondrial DNA sequences from the entire Mediterranean region and Atlantic Islands. Phylogenetic reconstructions group these western Palaearctic Plecotus into two major clades which split at least 5 Myr ago and that are each subdivided into further subgroups. An 'auritus group' includes the traditional P. auritus species and its sister taxon P. macrobullaris (=P. alpinus) plus related specimens from the Middle East. P. auritus and P. macrobullaris have broadly overlapping distributions in Europe, although the latter is apparently more restricted to mountain ranges. The other major clade, the 'austriacus group,' includes the European species P. austriacus and at least two other related taxa from North Africa (including P. teneriffae from the Canary Islands), the Balkans and Anatolia (P. kolombatovici). The sister species of this 'austriacus group' is P. balensis, an Ethiopian endemic. Phylogenetic reconstructions further suggest that P. austriacus reached Madeira during its relatively recent westward expansion through Europe, while the Canary Islands were colonized by a North African ancestor. Although colonization of the two groups of Atlantic Islands by Plecotus bats followed very distinct routes, neither involved lineages from the 'auritus group.' Furthermore, the Strait of Gibraltar perfectly segregates the distinct lineages, which confirms its key role as a geographic barrier. This study also stresses the biogeographical importance of the Mediterranean region, and particularly of North Africa, in understanding the evolution of the western Palaearctic biotas.     

Intraspecific evolution of Canary Island Plecotine bats,based on mtDNA sequences

Island differentiation and relationships with congenerics were investigated in the endemic Canary Island bat Plecotus teneriffae, based on approximately 1 kb of mtDNA from the 16S rRNA and cytochrome b genes. P. teneriffae had closer affinities with P. austriacus than with P. auritus. Levels of differentiation between Canary Islands were quite high relative to Pipistrelle-like bats, consistent with philopatric behaviour in the Plecotus genus. Cladogenesis within P. teneriffae appears to have occurred after the emergence of the islands of El Hierro and La Palma during the Pleistocene. An intraspecific network shows that haplotypes from the younger islands of La Palma and El Hierro are connected to the Tenerife haplotype by a similarly large number of mutational steps. This suggests that they were both colonised at a similar time from the much older island of Tenerife. The other Plecotine bat species, Barbastellus barbastellus shows close affinities with B. barbastellus from mainland Spain, with levels of mtDNA divergence being comparable with intraspecific variation within other mammal species.     

Conflicting molecular phylogenies of European long-eared bats (Plecotus) can be explained by cryptic diversity

Conflicting phylogenetic signals of two data sets that analyse different portions of the same molecule are unexpected and require an explanation. In the present paper we test whether (i) differential evolution of two mitochondrial genes or (ii) cryptic diversity can better explain conflicting results of two recently published molecular phylogenies on the same set of species of long-eared bats (genus Plecotus). We sequenced 1714bp of three mitochondrial regions (16S, ND1, and D-loop) of 35 Plecotus populations from 10 European countries. A likelihood ratio test revealed congruent phylogenetic signals of the three data partitions. Our phylogenetic analyses demonstrated that the existence of a previously undetected Plecotus lineage caused the incongruities of previous studies. This lineage is differentiated on the species level and lives in sympatry with its sister lineage, Plecotus auritus, in Switzerland and Northern Italy. A molecular clock indicates that all European Plecotus species are of mid or late Pliocene origin. Plecotus indet. was previously described as an intergrade between P. auritus and Plecotus austriacus since it shares morphological characters with both. It is currently known from elevations above 800 m a.s.l. in the Alps, the Dinarian Alps and the Pindos mountains in Greece. Since we could demonstrate that incongruities of two molecular analyses simply arose from the mis-identification of one lineage, we conclude that molecular phylogenetic analyses do not free systematists from a thorough inclusion of morphological and ecological data.     

Identification and distribution of bats of the genus Plecotus in England

The chief aim of this work was to find characters which would enable positive identifications to be made of living Plecotus auritus , Linn., 1758, and Plecotus austriacus Fischer, 1829. Previous work has related to museum material and usually did not give measurements which could be obtained from living bats. Sixteen indices per bat were at first measured, but only 12 of these were found to be significant and necessary for identification. The results were reasonably reproducible which allows their comparison with European data. Other diagnostic characters are described together with the known distribution in the British Isles. A summary of the known distribution and ecological data from Europe is given.     

Deep sympatric mitochondrial divergence without reproductive isolation in the common redstart Phoenicurus phoenicurus

Mitochondrial DNA usually shows low sequence variation within and high sequence divergence among species, which makes it a useful marker for phylogenetic inference and DNA barcoding. A previous study on the common redstart (Phoenicurus phoenicurus) revealed two very different mtDNA haplogroups (5% K2P distance). This divergence is comparable to that among many sister species; however, both haplogroups coexist and interbreed in Europe today. Herein, we describe the phylogeographic pattern of these lineages and test hypotheses for how such high diversity in mtDNA has evolved. We found no evidence for mitochondrial pseudogenes confirming that both haplotypes are of mitochondrial origin. When testing for possible reproductive barriers, we found no evidence for lineage‐specific assortative mating and no difference in sperm morphology, indicating that they are not examples of cryptic species, nor likely to reflect the early stages of speciation. A gene tree based on a short fragment of cytochrome c oxidase subunit 1 from the common redstart and 10 other Phoenicurus species, showed no introgression from any of the extant congenerics. However, introgression from an extinct congeneric cannot be excluded. Sequences from two nuclear introns did not show a similar differentiation into two distinct groups. Mismatch distributions indicated that the lineages have undergone similar demographic changes. Taken together, these results confirm that deeply divergent mitochondrial lineages can coexist in biological species. Sympatric mtDNA divergences are relatively rare in birds, but the fact that they occur argues against the use of threshold mtDNA divergences in species delineation.     

Homogeneity at nuclear microsatellite loci masks mitochondrial haplotype diversity in the endangered fanshell pearlymussel (Cyprogenia stegaria)

We report on multiple patterns of differentiation and connectivity in the fanshell pearlymussel (Cyprogenia stegaria), based on different markers. Knowledge of genetic variation and genetic connectivity among remaining populations of this federally endangered species is needed to initiate implementation of the species recovery plan. We collected tissue samples from 96 specimens from the Green, Rolling Fork, and Licking Rivers, tributaries to the Ohio River, and the Clinch River, a tributary to the Tennessee River, providing broad coverage of the current distributional range of the species. Results from 7 nuclear DNA microsatellite markers suggested minimal population-level differentiation, whereas a mitochondrial DNA (mtDNA) marker (ND1) exhibited significant differentiation between C. stegaria in the Clinch River and the Ohio River populations. The ND1 data also confirm the existence of 2 distinct mtDNA lineages in the genus that transcends species boundaries. Further analyses suggest that the disproportionally strong signal from 2 very divergent ND1 lineages possibly masks finer-grained structure in the Ohio River population, based on one of the mtDNA lineages only. We recommend further sampling to confirm the absence of one lineage from the upper Clinch River drainage and suggest that provisional management guidelines should limit reciprocal exchanges among C. stegaria populations from the Clinch River and those in the Ohio River system.     

Discordant patterns of morphological and molecular change in broadtail madtoms (genus Noturus)

Two morphologically distinct forms of an undescribed madtom catfish (Noturus sp.) occur in the rivers and lakes of southeastern USA. 'Lake' broadtail madtoms are endemic to Lake Waccamaw and are probably related to nearby 'river' broadtail populations. To investigate phylogenetic relationships, we surveyed mitochondrial DNA (mtDNA) sequence variation in 'lake' and 'river' broadtails and other members of the genus Noturus. Mitochondrial rDNA data suggest a sister group relationship between broadtail madtoms and N. insignis, not N. leptacanthus as posited previously. Population-level analyses using additional mtDNA characters (NADH dehydrogenase subunit 2 (ND2) and cytochrome b (Cytb)) identified two highly divergent genetic lineages within broadtail madtoms that do not correspond to the morphological designations 'river' and 'lake'.     

Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

Georgian and kurd mtDNA sequence analysis shows a lack of correlation between languages and female genetic lineages

Mitochondrial DNA sequences from Georgians and Kurds were analyzed in order to test the possible correlation between female lineages and languages in these two neighboring West Eurasian groups. Mitochondrial sequence pools in both populations are very similar despite their different linguistic and prehistoric backgrounds. Both populations present mtDNA lineages that clearly belong to the European gene pool, as shown by 1) similar nucleotide and sequence diversities; 2) a large number of sequences shared with the rest of European samples; 3) nonsignificant genetic distances; and 4) classification of the present lineages into the major European mtDNA haplogroups already described. The outlier position of the populations from the Caucasus according to classical genetic markers is not recognized in the present Georgian mtDNA sequence pool. This result suggests that the differentiation of mtDNA sequences in West Eurasia and the outlier features of Caucasian populations should be attributed to different processes. Moreover, the putative linguistic relationship between Caucasian groups and the Basques, another outlier population within Europe for classical genetic markers, is not detected by the analysis of mtDNA sequences.     

The genetic relationship between the Finns and the Finnish Saami (Lapps): analysis of nuclear DNA and mtDNA.

The genetic relationships between two Finno-Ugric-speaking populations, the Finns and the Finnish Saami (Lapps), were studied by using PCR for six nuclear-DNA marker loci, mitochondrial restriction-site polymorphism, and sequence variation of a 360-bp segment of the mitochondrial control region. The allele frequencies of each of the nuclear-DNA marker loci and the frequencies of mtDNA restriction haplotypes were significantly different between the populations. The Saami showed exceptionally low variation in their mtDNA restriction sites. The 9-bp deletion common in East Asian populations was not observed, nor did the haplotype data fit into the haplogroup categorization of Torroni et al. The average number of nucleotide substitutions from the mtDNA haplotype data indicated that the Finnish Saami may be closer to the Finns than to the other reference populations, whereas nuclear DNA suggested that the Finns are more closely related to the European reference populations than to the Finnish Saami. The similarity of the Finns to the other Europeans was even more pronounced according to the sequence data. We were unable to distinguish between the Finns and either the Swiss or Sardinian reference populations, whereas the Finnish Saami clearly stood apart. The Finnish Saami are distinct from other Circumarctic populations, although two of the lineages found among the Saami showed closer relationship to the Circumarctic than to the European lineages. The sequence data indicated an exceptionally high divergence for the Saami mtDNA control lineages. The distribution of the pairwise nucleotide differences in the Saami suggested that this population has not experienced an expansion similar to what was indicated for the Finns and the reference populations.     

Incongruence in the pattern and timing of intra-specific diversification in bronze frogs and bullfrogs (Ranidae)

We compare patterns of lineage divergence in mitochondrial DNA (mtDNA) sequences of two protein-encoding mitochondrial genes (cyt b and ND2) in two ecologically similar, co-distributed, and closely related ranid frogs (Rana clamitans and Rana catesbeiana), that are geographically widespread, and frequently syntopic. We identified three lineages in R. clamitans, separated by 0.5% to 2.1% net corrected sequence divergence, comparable to two R. catesbeiana lineages separated by 0.6%. The geographic pattern of lineage distribution differed notably between the two species. In R. clamitans, we found a Coastal Plain-Appalachian (CPA) lineage restricted to south and east of the Appalachian Mountains and a widespread lineage that encompassing nearly all the sampled range. A third distinct and divergent lineage was detected in one location in the southwest portion of the range (Louisiana). This pattern contrasts with the east-west pattern in R. catesbeiana, and reflects possible differences in refugial dynamics and patterns of range expansion. Although both species have undergone range expansion and population growth, coalescent reconstruction of N(e) reflects larger lineages but more recent divergence in R. clamitans relative to R. catesbeiana, reflecting significant differences in population history or divergent patterns of molecular evolution at mtDNA.     

Molecular systematics of Xenocyprinae (teleostei: cyprinidae): taxonomy, biogeography, and coevolution of a special group restricted in East Asia

We surveyed mitochondrial DNA (mtDNA) sequence variation in the subfamily Xenocyprinae from China and used these data to estimate intraspecific, interspecific, and intergeneric phylogeny and assess biogeographic scenarios underlying the geographic structure of lineages. We sequenced 1140 bp of cytochrome b from 30 individuals of Xenocyprinae and one putative outgroup (Myxocypris asiaticus) and also sequenced 297 bp of ND4L, 1380 bp of ND4, 68 bp of tRNA(His), and 69 bp of tRNA(Ser) from 17 individuals of Xenocyprinae and the outgroup (M. asiaticus). We detected high levels of nucleotide variation among populations, species, and genera. The phylogenetic analysis suggested that Distoechodon hupeinensis might be transferred to the genus Xenocypris, the taxonomic status of the genus Plagiognathops might be preserved, and species of Xenocypris and Plagiognathops form a monophyletic group that is sister to the genus Distoechodon and Pseudobrama. The introgressive hybridization might occur among the populations of X. argentea and X. davidi, causing the two species to not be separated by mtDNA patterns according to their species identification, and the process and direction of hybridization are discussed. The spatial distributions of mtDNA lineages among populations of Xenocypris were compatible with the major geographic region, which indicated that the relationship between Hubei + Hunan and Fujian is closer than that between Hubei + Hunan and Sichuan. From a perspective of parasite investigation, our data suggested that the fauna of Hexamita in Xenocyprinae could be used to infer the phylogeny of their hosts.     

So close and so different: comparative phylogeography of two small mammal species, the yellow-necked fieldmouse (Apodemus flavicollis) and the woodmouse (Apodemus sylvaticus) in the Western Palearctic region

In Europe, concordant geographical distribution among genetic lineages within different species is rare, which suggests distinct reactions to Quaternary ice ages. This study aims to determine whether such a discrepancy also affects a pair of sympatric species, which are morphologically and taxonomically closely related but which have slight differences in their ecological habits. The phylogeographic structures of two European rodents, the Yellow-necked fieldmouse (A. flavicollis) and the woodmouse (Apodemus sylvaticus) were, therefore, compared on the basis of mitochondrial DNA cytochrome b (mtDNA cyt b) sequences (965 base pairs) from 196 specimens collected from 59 European localities spread throughout the species distributions. The results indicate that the two species survived in different ways through the Quaternary glaciations. A. sylvaticus survived in the Iberian Peninsula from where it recolonized almost all Europe at the end of the last glaciation. Conversely, the refuge from which A. flavicollis recolonized Europe, including northern Spain, during the Holocene corresponds to the Italo-Balkan area, where A. sylvaticus suffered a serious genetic bottleneck. This study confirms that even closely related species may have highly different phylogeographic histories and shows the importance of ecological plasticity of the species for their survival through climate change. Finally, it suggests that phylogeographic distinctiveness may be a general feature of European species.     

Congruent Avian Phylogenies Inferred from Mitochondrial and Nuclear DNA Sequences

Recent molecular studies addressing the phylogenetic relationships of avian orders have had conflicting results. While studies using nuclear DNA sequences tend to support traditional taxonomic views, also supported by morphological data [(paleognaths (galloanseres (all other birds)))], with songbirds forming a clade within Neoaves (all other birds), analyses with complete mtDNA genomes have resulted in topologies that place songbirds as one of the earliest-diverging avian lineages. Considering that over half of the extant bird species are songbirds, these different results have very different implications for our understanding of avian evolution. We analyzed data sets comprising nearly 4 kb of mitochondrial DNA (mtDNA) (complete 12S, ND1, ND2, and cytochrome b) plus 600 bp of the nuclear gene c-mos for 15 birds that were chosen to represent all major avian clades and to minimize potential long-branch attraction problems; we used a partition-specific maximum likelihood approach. Our results show congruence with respect to the ingroup among phylogenies obtained with mtDNA and the nuclear gene c-mos, separately or combined. The data sets support a traditional avian taxonomy, with paleognaths (ratites and tinamous) occupying a basal position and with songbirds more derived and forming a monophyletic group. We also show that, for mtDNA studies, turtles may be a better outgroup for birds than crocodilians because of their slower rate of sequence evolution.     

High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species

Sympatric cryptic species, characterized by low morphological differentiation, pose a challenge to understanding the role of interspecific competition in structuring ecological communities. We used traditional (morphological) and novel molecular methods of diet analysis to study the diet of two cryptic bat species that are sympatric in southern England (Plecotus austriacus and P. auritus) (Fig. 1). Using Roche FLX 454 (Roche, Basel, CH) high-throughput sequencing (HTS) and uniquely tagged generic arthropod primers, we identified 142 prey Molecular Operational Taxonomic Units (MOTUs) in the diet of the cryptic bats, 60% of which were assigned to a likely species or genus. The findings from the molecular study supported the results of microscopic analyses in showing that the diets of both species were dominated by lepidopterans. However, HTS provided a sufficiently high resolution of prey identification to determine fine-scale differences in resource use. Although both bat species appeared to have a generalist diet, eared-moths from the family Noctuidae were the main prey consumed. Interspecific niche overlap was greater than expected by chance (O(jk) = 0.72, P < 0.001) due to overlap in the consumption of the more common prey species. Yet, habitat associations of nongeneralist prey species found in the diets corresponded to those of their respective bat predator (grasslands for P. austriacus, and woodland for P. auritus). Overlap in common dietary resource use combined with differential specialist prey habitat associations suggests that habitat partitioning is the primary mechanism of coexistence. The performance of HTS is discussed in relation to previous methods of molecular and morphological diet analysis. By enabling species-level identification of dietary components, the application of DNA sequencing to diet analysis allows a more comprehensive comparison of the diet of sympatric cryptic species, and therefore can be an important tool for determining fine-scale mechanisms of coexistence.     

Heteroplasmy and Ancient Translocation of Mitochondrial DNA to the Nucleus in the Chinese Horseshoe Bat (Rhinolophus sinicus) Complex

The utility and reliability of mitochondrial DNA sequences in phylogenetic and phylogeographic studies may be compromised by widespread and undetected nuclear mitochondrial copies (numts) as well as heteroplasmy within individuals. Both numts and heteroplasmy are likely to be common across diverse taxa yet few studies have characterised their frequencies and variation at the intra-specific level. Here we report the presence of both numts and heteroplasmy in the mitochondrial control region of the Chinese horseshoe bat Rhinolophus sinicus. In total we generated 123 sequences from 18 bats, which contained two different numt clades (i.e. Numt-1 and Numt-2) and one mtDNA clade. The sequence divergence between Numt-1 and Numt-2 was 16.8% and each numt type was found in all four R. sinicus taxa, suggesting either two ancient translocations of mitochondrial DNA into the nucleus from the same source taxon, or a single translocation from different source taxa that occurred before the split of R. sinicus into different lineages. Within the mtDNA clade, phylogenetic relationships among the four taxa of R. sinicus were similar to those seen in previous results. Based on PCR comparisons, heteroplasmy was inferred between almost all individuals of R. sinicus with respect to sequence variation. Consistent with introgression of mtDNA between Central sinicus and septentrionalis, individuals from these two taxa exhibited similar signatures of repeated sequences in the control region. Our study highlights the importance of testing for the presence of numts and heteroplasmy when applying mtDNA markers to phylogenetic studies.     

Transfer of mitochondrial DNA from the northern red-backed vole (Clethrionomys rutilus) to the bank vole (C. glareolus)

Summary Using a silver staining method to detect DNA fragments produced by restriction enzymes, it was possible to compare mitochondrial DNAs (mtDNAs) from 85 individuals of the bank vole (Clethrionomys glareolus) trapped at 25 localities in Fennoscandia. There are two distinctly different mtDNA lineages, one occurring in southern and central Fennoscandia and the other in the northern parts. A fragment comparison method shows about 12.7% nucleotide sequence divergence between these two lineages. This major difference between animals of the same species could theoretically be explained by intraspecific lineage survivorship independent of species hybridization, or by introduction of an atypical mtDNA via hybridization with a closely related species. Analysis of mtDNAs from the two otherClethrionomys species present in Fennoscandia (C. rutilus andC. rufocanus) shows that the mtDNA of northernC. glareolus is very similar to that ofC. rutilus and that the mtDNA lineages of these two species cluster together in a phenogram, with small genetic distances among them. By contrast, electrophoresis of proteins encoded by 17 nuclear loci reveals fixed allelic differences between these two species at 8 loci. Hence the presence of two distinctly different mtDNA lineages withinC. glareolus may be a consequence of a limited episode of hybridization betweenC. glareolus andC. rutilus, probably during the postglacial recolonization of Fennoscandia 8000–13,000 years ago.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Genetic diversity in the North American captive African elephant collection

One hundred and eight mitochondrial DNA (mtDNA) sequences were used in combination with multilocus genotypes generated from 19 nuclear microsatellite loci to assess the genetic diversity in the North American captive African elephant collection. The analysis of 4488 bp of mitochondrial nicotinamide adenine dinucleotide hydrogenenase subunit 4-cytochrome b (ND4-cyt b ) sequence revealed the existence of three distinct mtDNA lineages (Haplogroups I–III). The genetic distance between Haplogroups I/II and Haplogroup III was determined to be 0.046. The nucleotide diversity for the elephant mtDNA ND4-cyt b was <0.004. Analysis of multilocus genotypes indicated an average observed heterozygosity of 0.53 compared with the expected heterozygosity of 0.64. No subdivision in the data was detected including data partitioned by mtDNA haplogroup. Relatedness among the 108 individuals did not differ from that expected in random mating populations.     

The hidden diversity of Coleodactylus amazonicus (Sphaerodactylinae,Gekkota) revealed by molecular data

Coleodactylus amazonicus, a small leaf-litter diurnal gecko widely distributed in Amazon Basin, has been considered a single species with no significant morphological differences between populations along its range. A recent molecular study, however, detected large genetic differences between populations of central Amazonia and those in the easternmost part of the Amazon Basin, suggesting the presence of taxonomically unrecognised diversity. In this study, DNA sequences of three mitochondrial (16S, cytb, and ND4) and two nuclear genes (RAG-1, c-mos) were used to investigate whether the species currently identified as C.amazonicus contains morphologically cryptic species lineages. The present phylogenetic analysis reveals further genetic subdivision including at least five potential species lineages, restricted to northeastern (lineage A), southeastern (lineage B), central-northern (lineage E) and central-southern (lineages C and D) parts of Amazon Basin. All clades are characterized by exclusive groups of alleles for both nuclear genes and highly divergent mitochondrial haplotype clades, with corrected pairwise net sequence divergence between sister lineages ranging from 9.1% to 20.7% for the entire mtDNA dataset. Results of this study suggest that the real diversity of “C. amazonicus” has been underestimated due to its apparent cryptic diversification.