Phylogeographic analysis of mitochondrial gene flow and introgression in the salamander, Plethodon shermani

Plethodon shermani comprises a series of geographically disjunct populations occupying high-elevation mountain isolates. These populations hybridize at their borders with salamanders of the Plethodon glutinosus species complex, and past range expansions inferred from Pleistocene climatic cycles may have increased the possible genetic interactions between P. shermani and species of the P. glutinosus complex. Because mitochondrial DNA haplotypes often show introgression across species borders, we survey mtDNA variation for evidence of past and ongoing genetic interactions between P. shermani, its close relative Plethodon cheoah, and species of the P. glutinosus complex. Ongoing hybridization with the P. glutinosus-complex species Plethodon teyahalee is accompanied by extensive mitochondrial introgression in some Unicoi populations of P. shermani, but it has little genetic impact on P. shermani populations outside hybrid zones at three other isolates (Tusquitee, Wayah Bald, Standing Indian). Some Unicoi populations of P. shermani exhibit mtDNA evidence of past hybridization with diverse lineages from P. aureolus and P. glutinosus. The Tusquitee isolate of P. shermani is also characterized by mtDNA haplotypes most closely related to Plethodon aureolus and P. glutinosus, presumably introduced by past genetic contact with these species or with introgressed populations of Unicoi P. shermani. The mtDNA variation in sampled populations of the Wayah Bald and Standing Indian isolates of P. shermani appears largely unaffected by ongoing hybridization. Principal components analyses of allozymic data indicate that P. shermani isolates collectively form a genetically homogeneous unit clearly demarcated from species with which they have had current or past genetic interactions. Rapid mtDNA introgression associated with transient contacts between P. shermani and other species permits a fine-level resolution of evolutionary lineages not evident from allozymic data.     

Phylogeny of South American Bufo (Anura: Bufonidae) inferred from combined evidence

Despite the considerable research that has focused on the evolutionary relationships and biogeography of the genus Bufo, an evolutionary synthesis of the entire group has not yet emerged. In the present study, almost 4 kb of DNA sequence data from mitochondrial (12S, tRNA^Val, and 16S) and nuclear (POMC; Rag-1) genes, and 83 characters from morphology were analysed to infer a phylogeny of South American toads. Phylogenies were reconstructed with parsimony and maximum likelihood and Bayesian model-based methods. The results of the analysis of morphological data support the hypothesis that within Bufo , some skull characters (e.g. frontoparietal width), correlated with the amount of cranial ossification, are prone to homoplasy. Unique and unreversed morphological synapomorphies are presented that can be used to diagnose recognized species groups of South American toads. The results of all phylogenetic analyses support the monophyly of most species groups of South American Bufo . In most DNA-only and combined analyses, the South American (minus the B. guttatus and part of the ' B. spinulosus ' groups), North American, Central American, and African lineages form generally well-supported clades: ((((((((South America) (North America + Central America)) Eurasia) Africa) Eurasia) South America) West Indies) South America). This result confirms and extends prior studies recovering South American Bufo as polyphyletic. The biogeographical results indicate that: (1) The origin of Bufo predates the fragmentation of Gondwana; (2) Central and North American species compose the sister group to a large, 'derived' clade of South American Bufo ; and (3) Eurasian species form the sister group to the New World clade.  © 2006 The Linnean Society of London, Zoological Journal of the Linnean Society , 2006, 146 , 407–452.     

Phylogeography and conservation genetics of endangered European Margaritiferidae (Bivalvia: Unionoidea)

Margaritifera margaritifera and M. auricularia are among the most endangered freshwater mussels in the world, and the only species of the genus found in Europe. Our genetic study explores allozymic variability (27 loci) and differentiation at the mitochondrial sequence level (partial COI and 16S rRNA gene sequences). The Spanish M. auricularia population showed genetic parameters of variation that were of the same order as those of other freshwater molluscs (though at the lower end of the range), probably permitting its potential recovery. The difference between this species and M. margaritifera was clearly established (ten diagnostic allozymic loci, Nei = 0.462, and mean nucleotide divergence around 9.4%). The M. margaritifera populations analysed showed a certain degree of population genetic structure (according to allozyme data) that was not, however, related to a geographical cline. Nevertheless, two mitochondrial lineages (albeit very closely related) were identified: a northern lineage extending from Ireland to the Kola Peninsula including the western Atlantic coast, and a second cluster distributed from Ireland to the Iberian Peninsula. The phylogenetic relationships between these two species and other related taxa were established. The putative M. m. durrovensis could be considered an 'ecophenotype'. Palaeobiogeographical scenarios are presented and indicate unexpected 'recent' gene flow between M. margaritifera populations that were theoretically isolated in the early Tertiary.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society 2003, 78, 235–252.     

Is the last glaciation the only relevant event for the present genetic population structure of the meadow brown butterfly Maniola jurtina (Lepidoptera: Nymphalidae)?

Phylogeographical studies are available for a considerable number of European species, but few analyses exist for temperate species with very large and fairly continuous populations that are also absent from Northern Europe. Therefore, we studied the butterfly Maniola jurtina as a model for this group. The species has two major genetic lineages (mean genetic distance between lineages: 0.033; F _CT: 0.052), most probably evolving in glacial differentiation centres in the western and eastern Mediterranean. The onset of this differentiation might have been the beginning of the last glacial stage maximum some 40 kyr bp . A hybrid zone between these two lineages exists in western Central Europe. No genetic substructures have been found within the two lineages ( F _SC: 0.017) and average genetic distances are very small. Therefore, it is highly probable that postglacial expansion was of the phalanx type. There is, at most, very limited differentiation at regional and local scales. However, the genetic diversity within populations is high (means: A : 2.68; H _E: 17.2%; P : 78%), as would be predicted for such a common species. Comparison of these results with a published allozyme analysis revealed a similar phylogeographical pattern, but lower genetic diversity in the latter. Morphological patterns of wings and genitalia show similar geographical patterns as allozyme data.   © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 419–431.     

Phylogeography of Xiphorhynchus fuscus (Passeriformes, Dendrocolaptidae): vicariance and recent demographic expansion in southern Atlantic forest

Knowledge of the evolutionary processes that shaped a biota is important for both academic and conservation purposes. The objective of the present study is to analyse the mitochondrial genetic variation of Xiphorhynchus fuscus (Aves: Dendrocolaptidae) from the southern Atlantic forest in Brazil and Argentina, and to discuss whether the results support different hypotheses regarding the local intraspecific diversification of this species . We sequenced 575 bp of the control region of 114 specimens collected in the Brazilian states of Bahia, Minas Gerais, Rio de Janeiro, São Paulo, Paraná, and Santa Catarina, and in the province of Misiones in Argentina. We studied the population genetic structure with analysis of molecular variance and the demographic history with multiple regression analysis, coalescence simulations, and demographic tests. Xiphorhynchus fuscus presented a significant population genetic structure (Φ_st = 0.57). Three mitochondrial lineages were described, one associated with Xiphorhynchus fuscus tenuirostris and the others with Xiphorhynchus fuscus fuscus. The data did not support the primary influence of geographical barriers or rivers in the intraspecific diversification of X. fuscus in the southern Atlantic forest. Instead, the data supported the influence of isolation by geographical distance, recent vicariance events, and demographic expansions apparently related to Pleistocene and Holocene forest dynamics.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 73–84.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Cryptic speciation in liverworts – a case study in the Aneura pinguis complex

Bryophytes are amongst the most ancestral terrestrial plants and often have large distribution ranges across continents. Recent biochemical and molecular studies have suggested that many worldwide morphological species of bryophyte may represent genetically divergent and reproductively isolated cryptic species. We tested the cryptic species hypothesis in the thalloid liverwort Aneura pinguis complex. We applied analyses of chloroplast DNA (cpDNA) sequence variation and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods to discriminate between specimens of A. pinguis derived from various habitats in five distant geographical regions in Poland. Of the 19 specimens sequence characterized for the cpDNA tRNA^Leu region, seven haplotypes were identified divided into three nonmonophyletic clusters. The application of developed PCR-RFLP markers confirmed the existence of three tRNA^Leu types of A. pinguis (A–C) within the specimens derived from 21 populations. Sympatric populations of different tRNA^Leu types were found in lowland and mountain regions. No clear correlation between stand type and the presence of two tRNA^Leu types (A, B) was observed, as both were growing on soil, humus, and rocks. The tRNA^Leu type C was found only on humus and its distribution was restricted to low-lying northern populations. The above results indicate that the A. pinguis complex is highly differentiated at the molecular level and may represent three cryptic species.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 273–282.     

Mitochondrial differentiation, introgression and phylogeny of species in the Tegenaria atrica group (Araneae: Agelenidae)

The relationships between the three members of the Tegenaria atrica group ( T. atrica , T. saeva and T. gigantea ) were examined with DNA sequence data from mitochondrial CO1, 16S rRNA, tRNA^leu(CUN) and ND1 genes. Members of this group of large house spiders have overlapping distributions in western Europe and hybridize with each other to a variable degree. The close relatedness of all three species was supported by all analyses. T. saeva and T. gigantea are more closely affiliated than either is to T. atrica . Haplotypes clearly assignable to T. gigantea were also present in many specimens of T. saeva , suggesting asymmetrical introgression of mtDNA from T. gigantea into T. saeva . Molecular clock calibrations (CO1) suggest that deeper divisions within the genus Tegenaria may be in excess of 10 million years old, and that the evolutionary history of the T. atrica group has been moulded by Quaternary glacial–interglacial cycles.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 81 , 79–89.     

Falling apart and merging: diversification of slender salamanders (Plethodontidae: Batrachoseps) in the American West

The plethodontid genus Batrachoseps , the slender salamanders, is the most diverse clade of salamanders in western North America, but it has posed taxonomic difficulties because it contains many morphologically cryptic species. A segment of the mitochondrial DNA gene cytochrome b was studied for 278 individuals densely sampled from throughout the range of all 18 described species and several undescribed species. Phylogenetic analyses of the mtDNA data identify six major clades, one corresponding to the subgenus Plethopsis and five within a monophyletic subgenus Batrachoseps. All major clades and most species within these clades display strong phylogeographic structuring. Comparisons of mtDNA and allozyme data show that several allozymically cohesive groups are not monophyletic with respect to mtDNA. We suggest that this phenomenon results from fragmentation of populations, divergence in allopatry, and then recontact and gradual merging of units caused predominantly by male-mediated gene flow. The mtDNA offers evidence that populations were once more isolated than they are now, while the patterns of allozyme variation reflect recent and current interactions among populations. The complex patterns of morphological, allozymic and mtDNA variation associated with the constantly changing geological landscape give insight into the nature of processes responsible for species formation in Batrachoseps .  © 2002 The Linnean Society of London. Biological Journal of the Linnean Society , 2002, 76 , 361–391.     

Systematics and evolutionary relationships of the mountain lizard Liolaemus monticola (Liolaemini): how morphological and molecular evidence contributes to reveal hidden species diversity

The delimitation of species is a major issue in systematic biology and has been a re-emerging discipline in the last decade. A number of studies have shown that the use of multiple data sets is critical for the identification of cryptic species, particularly in groups with complex evolutionary histories. Liolaemus monticola is a montane lizard species distributed in central Chile (32°–42°S), with four described subspecies in a latitudinal gradient from north to south: L. m. monticola , L. m. chillanensis , L. monticola ssp. and L. m. villaricensis . In order to test the systematic status and phylogenetic relationships of the taxa included in the L. monticola group, we analysed morphological (morphometric and meristic) and molecular (allozyme and mitochondrial DNA) data sets. The results of the morphological analyses showed that meristic variables correctly assigned individuals with higher accuracy than did morphometric characters. The results of the analyses of allozyme data revealed eight diagnostic loci that are evidence for significant differences among the four L. monticola subspecies. Phylogenetic analyses with mitochondrial DNA data, including additional species, showed that the L. monticola group is polyphyletic. We postulate that the four current subspecies represent independent evolutionary lineages and must be raised to the specific level as L. monticola , L. chillanensis and L. villaricensis . The taxonomic status of the unnamed L. monticola ssp. remains unresolved, although we provide a preliminary proposal.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 635–650.     

Hidden Mediterranean biodiversity: molecular evidence for a cryptic species complex within the reef building vermetid gastropod Dendropoma petraeum (Mollusca: Caenogastropoda)

The reef building vermetid gastropod Dendropoma petraeum inhabits the warmest waters of the Mediterranean Sea and is considered a threatened marine species. The aim of this study was to characterize its genetic structure throughout its whole distribution range using mitochondrial and nuclear sequence data. Because of its sessile adult lifestyle and lack of a pelagic larval stage, we expected a markedly subdivided population structure with limited levels of gene flow. Fragments of the mitochondrial genes cytochrome c oxidase subunit I (COI) and 16S rRNA (16S), were sequenced, along with the nuclear ribosomal cluster (internal transcribed spacer; ITS) in specimens from 18 localities. Our analyses identified four highly distinct phylogroups separated by a mean divergence of > 14% according to the COI sequence data or > 9% according to 16S, but differing only slightly in morphology. The nuclear data (ITS) indicated a lower substitution rate (divergence among groups of around 1%). These large genetic distances among the four lineages clearly point to the existence of a cryptic species complex within D. petraeum comprising at least four species. Differences in the characteristics of intracapsular larval development and protoconch were also detected among these lineages. The allopatric distribution of these cryptic species supports a predominantly vicariant-based cladogenetic pattern for the genus Dendropoma in the Mediterranean.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 898–912.     

Population genetics of Dioon edule Lindl. (Zamiaceae, Cycadales): biogeographical and evolutionary implications

Dioon edule Lindl. (Zamiaceae) is a cycad endemic to Mexico, that occurs as one species D. edule and the geographical variety D. edule var. angustifolium (Miq.) Miq. Dioon edule has a north to south distribution in eastern Mexico. In this study, we analysed 14 allozymic loci in eight populations of D. edule from its total distribution range by sampling all known populations. Patterns of diversity and genetic variability, within and among populations, were obtained. The mean number of alleles per locus ( A ) was 1.44 and the percentage of polymorphic loci was relatively high ( P  = 54.78). The mean observed ( H _O ) and expected heterozygosity ( H _E ) were 0.27 and 0.24, respectively. F -statistics revealed an excess of heterozygous genotypes, locally and globally ( F  = −0.17 and f  = −0.27, respectively). The genetic variation explained by differences among populations was only 7.5%. We also detected a negative relationship between genetic diversity and latitude. On average, the gene flow between population pairs was relatively high ( Nm  = 2.98); furthermore, gene flow between population pairs was significantly correlated with geographical distances ( r  = −0.38, P  = 0.025). Therefore, patterns of genetic diversity in D. edule appear to be associated with the post-Pleistocene spread of the species, from its southerly (origin) to its northerly range (derived populations, including its central distribution). The biogeographical and evolutionary aspects of the results of this study are discussed. We recognize Dioon angustifolium Miq. for the northernmost disjunct populations.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 457–467     

Historical diversification of floodplain forest specialist species in the Amazon: a case study with two species of the avian genus Xiphorhynchus (Aves: Dendrocolaptidae)

Phylogeographical and population genetics methods are used to reconstruct the diversification history of two species of the genus Xiphorhynchus (Aves: Dendrocolaptidae) associated with seasonally flooded forest types in Amazonia. Sequences of the mitochondrial gene cytochrome b were assessed for 21 and 30 individuals, belonging to eight and ten populations, of Xiphorhynchus kienerii and Xiphorhynchus obsoletus , respectively. Uncorrected genetic distances among unique haplotypes recovered ranged only from 0.01% to 0.4% for both species. Over 90% of the genetic variation detected in both species was partitioned within populations, and therefore was not structured geographically. Mismatch distributions and values of Tajima's D -tests indicate that both X. kienerii and X. obsoletus have had small evolutionary effective population sizes, but experienced a recent demographic expansion. These demographic expansions are tentatively dated as occurring over the last 18 000 years BP, a time frame which coincides with the establishment of the early and mid-Holocene age floodplain forest in most of central and eastern Amazonia, following a period of increased river stages throughout the basin. Based on phylogenetic, phylogeographical, and populations genetics data obtained for X. kienerii and X. obsoletus , an evolutionary scenario is proposed to account for the historical diversification of floodplain specialist species in Amazonia.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 383–395.     

Divergent and reticulate processes in evolution of Ethiopian Lophuromys flavopunctatus species complex: evidence from mitochondrial and nuclear DNA differentiation patterns

Molecular study of mitochondrial and nuclear genes and cytogenetic analysis were performed to examine possible patterns of speciation in the diverse Lophuromys flavopunctatus species complex of Ethiopia. Phylogenetic analysis of mtDNA data resulted in an unresolved bush of ten deeply diverged haplotype groups corresponding to potential species either well supported by various types of character or 'cryptic'. The cytogenetic analysis showed representatives of five of these mtDNA lineages to share an identical karyotype (2 n  = 70, NFa = 84), that has not been found previously in Ethiopia. One of them, L.  cf.  sikapusi , being a member of the L. flavopunctatus species complex, demonstrates remarkable morphological similarity to representatives of another species complex, L. sikapusi s.l ., which might be considered as a result of convergent evolution in analogous environments. Analysis of RAPD data suggests that at least two mtDNA types might have been subject to interspecific transfer due to hybridization. In the case of two sympatric haplotypes of L. brunneus we may assume that the contemporary pattern of variation between them can be explained by relatively recent hybridization with another distinct species, L. flavopunctatus . The formation of two groups belonging to distinct mitochondrial lineages within northern populations could be associated with more complex processes including ancient hybridization.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 83 , 301–316.     

The geography of diversification in the mormoopids (Chiroptera: Mormoopidae)

The traditional explanation of the distribution of the Mormoopidae is that this family originated in southern Central America or northern South America, later expanding its range north to Mexico and the West Indies, and differentiating into eight species. An alternative fossil-based hypothesis argues that the family originated in the northern Neotropics, reached the Caribbean early in its history, and dispersed to South America after the completion of the Isthmus of Panama. The present study analyses new and previously published sequence data from the mitochondrial 12S, tRNA^val, 16S, and cytochrome b , and the nuclear Rag 2, to evaluate species boundaries and infer relationships among extant taxa. Fixed differences in cytochrome b often coincide with published morphological characters and show that the family contains at least 13 species. Two additional, morphologically indistinct, lineages are restricted to Suriname and French Guiana. Phylogeny-based inferences of ancestral area are equivocal on the geographical origin of mormoopids, in part because several internal nodes are not resolved with the available data. Divergences between Middle American and Antillean populations are greater than those between Mexico/Central America and South America. This suggests that mormoopids diversified in northern Neotropics before entering South America. A northern neotropical origin for mormoopids is congruent with both the Tertiary fossil record and recent phylogenetic hypotheses for the sister family to the Mormoopidae, the Phyllostomidae.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 101–118.     

Clinal genetic variation and isolation by distance in the European eel Anguilla anguilla (L.)

The genetic variability and structure of the European eel ( Anguilla anguilla L.) in populations throughout Europe was reassessed using 15 allozymic loci, seven of which were polymorphic. Seven sites were sampled on a latitudinal gradient across the natural continental range, extending from southern France to southern Norway. Heterozygosity ( H _e = 0.05) and level of polymorphism (P = 0.43) were comparable to other marine fish. Populations were poorly differentiated ( G _ST = 0.014, F _ST = 0.002), which is not surprising considering the high dispersal capability of the European eel. However, a significant geographical cline was detected at two alleles ( IDH-1 * 100 and GPI-1 * 110 ), and genetic distances ( D _CE) were concordant with geographical coastal distances. Mantel tests, pairwise F _ST's and multidimensional scaling analyses identify three distinct groups: Northern Europe, Western Europe and the Mediterranean Sea. We propose that the clinal genetic structure in the European eel may be due to (1) isolation by distance (as recently detected with microsatellites), (2) temporal reproductive separation, (3) post-larval selective forces, (4) contact between formerly separated groups or (5) some combination thereof.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 77 , 509–521     

The evolutionary pathway of light emission in myodocopid Ostracoda

The evolutionary history of bioluminescence and iridescence in myodocopid ostracods was estimated by phylogenetic analysis of mitochondrial 16S ribosomal RNA sequences. The inferred phylogeny of the myodocopids suggests that the common ancestor of Myodocopida evaluated in this study exhibits iridescence. This type of light emission was once lost and recaptured independently in the descendant lineages. Bioluminescent species also evolved from non-luminous ancestral species. In the suborder Myodocopina, all the bioluminescent species form a monophyletic group, suggesting that bioluminescence evolved only once. Structural differences between two bioluminescent groups in the order Myodocopida suggests independent origins for bioluminescence.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 87 , 449–455.     

Rapid diversification of mating systems in ciliates

Ciliates are a diverse group of microbial eukaryotes that exhibit tremendous variety in several aspects of their mating systems. To understand the evolutionary forces driving mating system diversification in ciliates, we use a comparative approach synthesizing data from many ciliate species in light of recent phylogenetic analyses. Specifically, we investigate the evolution of number of mating types, mode of mating type inheritance, and the molecular determinants of mating types across the taxonomic diversity of ciliates, with an emphasis on three well-studied genera: Tetrahymena , Paramecium , and Euplotes . We find that there have been many transitions in the number of mating types, and that the requirement of nuclear reorganization may be a more important factor than genetic exchange in determining the optimum number of mating types in a species. We also find that the molecular determinants of mating types and mode of inheritance are evolving under different constraints in different lineages of ciliates. Our results emphasize the need for further detailed examination of mating systems in understudied ciliate lineages.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 187–197.     

Phylogeography and species limits in the Gymnodactylus darwinii complex (Gekkonidae, Squamata): genetic structure coincides with river systems in the Brazilian Atlantic Forest

Phylogenetic analyses based on mtDNA cytochrome  b were performed in 42 lizards of the Gymnodactylus darwinii complex from three regions within Brazil's Atlantic Forest. Mainland regions and continental shelf islands in the south-eastern range and mainland areas from the north-east were sampled. The criteria of maximum parsimony (MP), maximum likelihood (ML) and Bayesian methods were explored, with the robustness for nodes assessed by bootstrapping (MP and ML) and posterior probabilities (Bayesian searches). By all methods, three distinctive phylogroups were recovered: a south-eastern clade (SE) and two clades from northern regions (NE₁ and NE₂). The pattern of genetic structure of the major clades coincided with the presence of river systems in the Atlantic Forest, and based on corrected genetic distances between those clades, divergence times were tentatively estimated using mtDNA rates calibrated for squamate reptiles. The putative role of Atlantic Forest rivers in generating differentiation is discussed. We present a hypothesis of species limits for G. darwinii , based on concordant lines of evidence including cytogenetic and mtDNA analyses. Two chromosome races (cytotype A, 2n = 38; and cytotype B, 2n = 40) had distributions concordant with clades SE and NE₁ + NE₂, respectively. These races are interpreted to be full species on the basis of a number of empirical criteria.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 13–26.