Limitations of Climatic Data for Inferring Species Boundaries: Insights from Speckled Rattlesnakes

Phenotypes, DNA, and measures of ecological differences are widely used in species delimitation. Although rarely defined in such studies, ecological divergence is almost always approximated using multivariate climatic data associated with sets of specimens (i.e., the “climatic niche”); the justification for this approach is that species-specific climatic envelopes act as surrogates for physiological tolerances. Using identical statistical procedures, we evaluated the usefulness and validity of the climate-as-proxy assumption by comparing performance of genetic (nDNA SNPs and mitochondrial DNA), phenotypic, and climatic data for objective species delimitation in the speckled rattlesnake (Crotalus mitchellii) complex. Ordination and clustering patterns were largely congruent among intrinsic (heritable) traits (nDNA, mtDNA, phenotype), and discordance is explained by biological processes (e.g., ontogeny, hybridization). In contrast, climatic data did not produce biologically meaningful clusters that were congruent with any intrinsic dataset, but rather corresponded to regional differences in atmospheric circulation and climate, indicating an absence of inherent taxonomic signal in these data. Surrogating climate for physiological tolerances adds artificial weight to evidence of species boundaries, as these data are irrelevant for that purpose. Based on the evidence from congruent clustering of intrinsic datasets, we recommend that three subspecies of C. mitchellii be recognized as species: C. angelensis, C. mitchellii, and C. Pyrrhus.     

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.     

Ecological niche modelling and nDNA sequencing support a new, morphologically cryptic beetle species unveiled by DNA barcoding

Background

DNA sequencing techniques used to estimate biodiversity, such as DNA barcoding, may reveal cryptic species. However, disagreements between barcoding and morphological data have already led to controversy. Species delimitation should therefore not be based on mtDNA alone. Here, we explore the use of nDNA and bioclimatic modelling in a new species of aquatic beetle revealed by mtDNA sequence data.

Methodology/Principal Findings

The aquatic beetle fauna of Australia is characterised by high degrees of endemism, including local radiations such as the genus Antiporus. Antiporus femoralis was previously considered to exist in two disjunct, but morphologically indistinguishable populations in south-western and south-eastern Australia. We constructed a phylogeny of Antiporus and detected a deep split between these populations. Diagnostic characters from the highly variable nuclear protein encoding arginine kinase gene confirmed the presence of two isolated populations. We then used ecological niche modelling to examine the climatic niche characteristics of the two populations. All results support the status of the two populations as distinct species. We describe the south-western species as Antiporus occidentalis sp.n.

Conclusion/Significance

In addition to nDNA sequence data and extended use of mitochondrial sequences, ecological niche modelling has great potential for delineating morphologically cryptic species.     

Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus)

Genealogical data are an important source of evidence for delimiting species, yet few statistical methods are available for calculating the probabilities associated with different species delimitations. Bayesian species delimitation uses reversible-jump Markov chain Monte Carlo (rjMCMC) in conjunction with a user-specified guide tree to estimate the posterior distribution for species delimitation models containing different numbers of species. We apply Bayesian species delimitation to investigate the speciation history of forest geckos (Hemidactylus fasciatus) from tropical West Africa using five nuclear loci (and mtDNA) for 51 specimens representing 10 populations. We find that species diversity in H. fasciatus is currently underestimated, and describe three new species to reflect the most conservative estimate for the number of species in this complex. We examine the impact of the guide tree, and the prior distributions on ancestral population sizes (θ) and root age (τ₀), on the posterior probabilities for species delimitation. Mis-specification of the guide tree or the prior distribution for θ can result in strong support for models containing more species. We describe a new statistic for summarizing the posterior distribution of species delimitation models, called speciation probabilities, which summarize the posterior support for each speciation event on the starting guide tree.     

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

Background

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal “barcoding gene”. An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

Methods and Findings

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a “barcode gap” using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

Conclusions

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a “species”.     

Genetic characterization of Amur sturgeon <Emphasis Type="Italic">Acipenser schrenckii</Emphasis> and its hybrid caught around Hokkaido

Genetic characterization was performed in five individuals of wild Amur sturgeon Acipenser schrenckii, and/or its presumed hybrid caught around Hokkaido, using a mitochondrial DNA (mtDNA) marker and two markers of nuclear DNA (nDNA). Genetic analyses indicated that two of the five fish had the mtDNA haplotype of Kaluga, Huso dauricus, whereas the nDNA markers indicated signs for both A. schrenckii and H. dauricus genotypes, referring to a hybrid origin. The other three fish were plausibly pure A. schrenckii. The results indicated the importance of combined usage of mtDNA and nDNA markers for correct species identification in sturgeon.     

Speciation on the Rocks: Integrated Systematics of the Heteronotia spelea Species Complex (Gekkota; Reptilia) from Western and Central Australia

The isolated uplands of the Australian arid zone are known to provide mesic refuges in an otherwise xeric landscape, and divergent lineages of largely arid zone taxa have persisted in these regions following the onset of Miocene aridification. Geckos of the genus Heteronotia are one such group, and have been the subject of many genetic studies, including H. spelea, a strongly banded form that occurs in the uplands of the Pilbara and Central Ranges regions of the Australian arid zone. Here we assess the systematics of these geckos based on detailed examination of morphological and genetic variation. The H. spelea species complex is a monophyletic lineage to the exclusion of the H. binoei and H. planiceps species complexes. Within the H. spelea complex, our previous studies based on mtDNA and nine nDNA loci found populations from the Central Ranges to be genetically divergent from Pilbara populations. Here we supplement our published molecular data with additional data gathered from central Australian samples. In the spirit of integrative species delimitation, we combine multi-locus, coalescent-based lineage delimitation with extensive morphological analyses to test species boundaries, and we describe the central populations as a new species, H. fasciolatus sp. nov. In addition, within the Pilbara there is strong genetic evidence for three lineages corresponding to northeastern (type), southern, and a large-bodied melanic population isolated in the northwest. Due to its genetic distinctiveness and extreme morphological divergence from all other Heteronotia, we describe the melanic form as a new species, H. atra sp. nov. The northeastern and southern Pilbara populations are morphologically indistinguishable with the exception of a morpho-type in the southeast that has a banding pattern resembling H. planiceps from the northern monsoonal tropics. Pending more extensive analyses, we therefore treat Pilbara H. spelea as a single species with phylogenetic structure and morphological heterogeneity.     

Deep intra‐island divergence of a montane forest endemic: phylogeography of the Puerto Rican frog Eleutherodactylus portoricensis (Anura: Eleutherodactylidae)

Aim Hypotheses proposed for lineage diversification of tropical montane species have rarely been tested within oceanic islands. Our goal was to understand how basin barriers and Pleistocene climatic fluctuations shaped the distribution of diversity in Eleutherodactylus portoricensis (Eleutherodactylidae), a frog endemic to the montane rain forests of Puerto Rico. Location The north‐eastern (Luquillo) and south‐eastern (Cayey) mountains of Puerto Rico. Methods We generated mitochondrial DNA (mtDNA) control region sequences (c. 565 bp) from 144 individuals of E. portoricensis representing 16 localities, and sequenced 646 bp of cytochrome b and 596 bp of nuclear DNA (nDNA) rhodopsin exon and intron 1 from a subset of individuals. We conducted a phylogenetic analysis on the mtDNA sequence data and explored population substructure with maximum parsimony networks, a spatial analysis of molecular variance, and pairwise F_ST analysis. Coalescent simulations were performed to test alternative models of population divergence in response to late Pleistocene interglacial periods. Historical demography was assessed through coalescent analyses and Bayesian skyline plots. Results We found: (1) two highly divergent groups associated with the disjunct Luquillo and Cayey Mountains, respectively; (2) a shallow mtDNA genetic discontinuity across the La Plata Basin within the Cayey Mountains; (3) phylogeographic congruence between nDNA and mtDNA markers; (4) divergence dates for both mtDNA and nDNA pre‐dating the Holocene interglacial (c. 10 ka), and nDNA suggesting divergence in the penultimate interglacial (c. 245 ka); and (5) historical demographic stability in both lineages. Main conclusions The low‐elevation Caguas Basin is a long‐term barrier to gene flow between the two montane frog populations. Measures of genetic diversity for mtDNA were similar in both lineages, but lower nDNA diversity in the Luquillo Mountains lineage suggests infrequent dispersal between the two mountain ranges and colonization by a low‐diversity founder population. Population divergence began prior to the Holocene interglacial. Stable population sizes over time indicate a lack of demonstrable demographic response to climatic changes during the last glacial period. This study highlights the importance of topographic complexity in promoting within‐island vicariant speciation in the Greater Antilles, and indicates long‐term persistence and lineage diversification despite late Pleistocene climatic oscillations.     

Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

The Role of Integrative Taxonomy in the Conservation Management of Cryptic Species: The Taxonomic Status of Endangered Earless Dragons (Agamidae: Tympanocryptis) in the Grasslands of Queensland,Australia

Molecular phylogenetics is increasingly highlighting the prevalence of cryptic species, where morphologically similar organisms have long independent evolutionary histories. When such cryptic species are known to be declining in numbers and are at risk of extinction due to a range of threatening processes, the disjunction between molecular systematics research and conservation policy becomes a significant problem. We investigate the taxonomic status of Tympanocryptis populations in Queensland, which have previously been assigned to T. tetraporophora, using three species delimitation approaches. The taxonomic uncertainties in this species-group are of particular importance in the Darling Downs Earless Dragon (T. cf. tetraporophora), which is ranked as an endangered ‘species’ of high priority for conservation by the Queensland Department of Environment and Heritage Protection. We undertook a morphological study, integrated with a comprehensive genetic study and species delimitation analyses, to investigate the species status of populations in the region. Phylogenetic analyses of two gene regions (mtDNA: ND2; nuclear: RAG1) revealed high levels of genetic divergence between populations, indicating isolation over long evolutionary time frames, and strongly supporting two independent evolutionary lineages in southeastern Queensland, from the Darling Downs, and a third in the Gulf Region of northern Queensland. Of the three species delimitation protocols used, we found integrative taxonomy the most applicable to this cryptic species complex. Our study demonstrates the utility of integrative taxonomy as a species delimitation approach in cryptic complexes of species with conservation significance, where limited numbers of specimens are available.     

Delimiting Species without Nuclear Monophyly in Madagascar's Mouse Lemurs

Background

Speciation begins when populations become genetically separated through a substantial reduction in gene flow, and it is at this point that a genetically cohesive set of populations attain the sole property of species: the independent evolution of a population-level lineage. The comprehensive delimitation of species within biodiversity hotspots, regardless of their level of divergence, is important for understanding the factors that drive the diversification of biota and for identifying them as targets for conservation. However, delimiting recently diverged species is challenging due to insufficient time for the differential evolution of characters—including morphological differences, reproductive isolation, and gene tree monophyly—that are typically used as evidence for separately evolving lineages.

Methodology

In this study, we assembled multiple lines of evidence from the analysis of mtDNA and nDNA sequence data for the delimitation of a high diversity of cryptically diverged population-level mouse lemur lineages across the island of Madagascar. Our study uses a multi-faceted approach that applies phylogenetic, population genetic, and genealogical analysis for recognizing lineage diversity and presents the most thoroughly sampled species delimitation of mouse lemur ever performed.

Conclusions

The resolution of a large number of geographically defined clades in the mtDNA gene tree provides strong initial evidence for recognizing a high diversity of population-level lineages in mouse lemurs. We find additional support for lineage recognition in the striking concordance between mtDNA clades and patterns of nuclear population structure. Lineages identified using these two sources of evidence also exhibit patterns of population divergence according to genealogical exclusivity estimates. Mouse lemur lineage diversity is reflected in both a geographically fine-scaled pattern of population divergence within established and geographically widespread taxa, as well as newly resolved patterns of micro-endemism revealed through expanded field sampling into previously poorly and well-sampled regions.     

The African timber tree <Emphasis Type="Italic">Entandrophragma congoense</Emphasis> (Pierre ex De Wild.) A.Chev. is morphologically and genetically distinct from <Emphasis Type="Italic">Entandrophragma angolense</Emphasis> (Welw.) C.DC

Interpreting morphological variability in terms of species delimitation can be challenging. However, correcting species delineation can have strong implications for the sustainable management of exploited species. Up to now, species delimitation between two putative timber species from African forests, Entandrophragma congoense and E. angolense, remained unclear. To investigate their differences, we applied an integrated approach which combines morphological traits and genetic markers. We defined 13 morphological characters from 81 herbarium specimens and developed 15 new polymorphic microsatellite markers to genotype 305 samples (herbarium samples and specimens collected in the field across the species distribution ranges). Principal component analysis (PCA) of morphological data and the Bayesian clustering analyses of genetic data were used to assess differentiation between putative species. These analyses support two well-differentiated groups (F_ST?=?0.30) occurring locally in sympatry. Moreover, these two groups present distinct morphological characters at the level of the trunk, leaflets, and seeds. Our genetic markers identified few individuals (4%) that seem to be hybrids, though there is no evidence of genetic introgression from geographic patterns of genetic variation. Hence, our results provide clear support to recognize E. congoense as a species distinct from E. angolense, with a much lower genetic diversity than the latter, and that should be managed accordingly. This work highlights the power of microsatellite markers in resolving species boundaries.     

Congruence of Microsatellite and Mitochondrial DNA Variation in Acrobat Ants (Crematogaster Subgenus Decacrema,Formicidae: Myrmicinae) Inhabiting Macaranga (Euphorbiaceae) Myrmecophytes

A previously reported mitochondrial DNA (mtDNA) phylogeny of Crematogaster (subgenus Decacrema) ants inhabiting Macaranga myrmecophytes indicated that the partners diversified synchronously and their specific association has been maintained for 20 million years. However, the mtDNA clades did not exactly match morphological species, probably owing to introgressive hybridization among younger species. In this study, we determined the congruence between nuclear simple sequence repeat (SSR, also called microsatellite) genotyping and mtDNA phylogeny to confirm the suitability of the mtDNA phylogeny for inferring the evolutionary history of Decacrema ants. Analyses of ant samples from Lambir Hills National park, northeastern Borneo, showed overall congruence between the SSR and mtDNA groupings, indicating that mtDNA markers are useful for delimiting species, at least at the local level. We also found overall high host-plant specificity of the SSR genotypes of Decacrema ants, consistent with the specificity based on the mtDNA phylogeny. Further, we detected cryptic genetic assemblages exhibiting high specificity toward particular plant species within a single mtDNA clade. This finding, which may be evidence for rapid ecological and genetic differentiation following a host shift, is a new insight into the previously suggested long-term codiversification of Decacrema ants and Macaranga plants.     

Delimiting species using DNA and morphological variation and discordant species limits in spiny lizards (Sceloporus)

Haplotype phylogenies based on DNA sequence data are increasingly being used to test traditional species-level taxonomies based on morphology. However, few studies have critically compared species limits based on morphological and DNA data, and the methods used to delimit species using either type of data are only rarely explained. In this paper, we review three approaches for species delimitation (tree-based with DNA data and tree-based and character-based with morphological data) and propose explicit protocols for each. We then compare species limits inferred from these approaches, using morphological and mtDNA data for the Yarrow's spiny lizard (Sceloporus jarrovii), a traditionally polytypic species from the southwestern United States and Mexico. All three approaches support division of S. jarrovii into five species, but only two species are the same among the three approaches. We find the greatest support for the five species that are delimited based on mtDNA data, and we argue that mtDNA data may have important (and previously unappreciated) advantages for species delimitation. Because different data and approaches can disagree so extensively, our results demonstrate that the methodology of species delimitation is a critical issue in systematics.     

Evaluating multilocus Bayesian species delimitation for discovery of cryptic mycorrhizal diversity

The increasing availability of DNA sequence data enables exciting new opportunities for fungal ecology. However, it amplifies the challenge of how to objectively classify the diversity of fungal sequences into meaningful units, often in the absence of morphological characters. Here, we test the utility of modern multilocus Bayesian coalescent-based methods for delimiting cryptic fungal diversity in the orchid mycorrhiza morphospecies Serendipita vermifera. We obtained 147 fungal isolates from Caladenia, a speciose clade of Australian orchids known to associate with Serendipita fungi. DNA sequence data for 7 nuclear and mtDNA loci were used to erect competing species hypotheses by clustering isolates based on: (a) ITS sequence divergence, (b) Bayesian admixture analysis, and (c) mtDNA variation. We implemented two coalescent-based Bayesian methods to determine which species hypothesis best fitted our data. Both methods found strong support for eight species of Serendipita among our isolates, supporting species boundaries reflected in ITS divergence. Patterns of host plant association showed evidence for both generalist and specialist associations within the host genus Caladenia. Our findings demonstrate the utility of Bayesian species delimitation methods and suggest that wider application of these techniques will readily uncover new species in other cryptic fungal lineages.     

Mating type gene analyses in the genus Diplodia: From cryptic sex to cryptic species

Cryptic species are common in Diplodia, a genus that includes some well-known and economically important plant pathogens. Thus, species delimitation has been based on the phylogenetic species recognition approach using multigene genealogies. We assessed the potential of mating type (MAT) genes sequences as phylogenetic markers for species delimitation in the genus Diplodia. A PCR-based mating type diagnostic assay was developed that allowed amplification and sequencing of the MAT1-1-1 and MAT1-2-1 genes, and determination of the mating strategies used by different species. All species tested were shown to be heterothallic. Phylogenetic analyses were performed on both MAT genes and also, for comparative purposes, on concatenated sequences of the ribosomal internal transcribed spacer (ITS), translation elongation factor 1-alpha (tef1-α) and beta-tubulin (tub2). Individual phylogenies based on MAT genes clearly differentiated all species analysed and agree with the results obtained with the commonly used multilocus phylogenetic analysis approach. However, MAT genes genealogies were superior to multigene genealogies in resolving closely related cryptic species. The phylogenetic informativeness of each locus was evaluated revealing that MAT genes were the most informative loci followed by tef1-α. Hence, MAT genes can be successfully used to establish species boundaries in the genus Diplodia.     

Exploring Species Level Taxonomy and Species Delimitation Methods in the Facultatively Self-Fertilizing Land Snail Genus Rumina (Gastropoda: Pulmonata)

Delimiting species in facultatively selfing taxa is a challenging problem of which the terrestrial pulmonate snail genus Rumina is a good example. These snails have a mixed breeding system and show a high degree of shell and color variation. Three nominal species (R. decollata, R. saharica and R. paivae) and two color morphs within R. decollata (dark and light) are currently recognized. The present study aims at evaluating to what extent these entities reflect evolutionary diverging taxonomic units, rather than fixed polymorphisms due to sustained selfing. Therefore, a phylogenetic analysis of nuclear (ITS1, ITS2) and mitochondrial DNA (COI, CytB, 12S rDNA, 16S rDNA) sequences was performed. Putative species in Rumina, inferred from the mitochondrial DNA phylogeny, were compared with those proposed on the basis of the COI gene by (1) DNA barcoding gap analysis, (2) Automatic Barcode Gap Discovery, (3) the species delimitation plug-in of the Geneious software, (4) the Genealogical Sorting Index, and (5) the General Mixed Yule Coalescent model. It is shown that these methods produce a variety of different species hypotheses and as such one may wonder to what extent species delimitation methods are really useful. With respect to Rumina, the data suggest at least seven species, one corresponding to R. saharica and six that are currently grouped under the name R. decollata. The species-level status of R. paivae is rejected.     

Examination of evolutionary relationships in the Cape fossorial skink species complex (Acontinae: Acontias meleagris meleagris) reveals the presence of five cryptic lineages

Variable morphological characters have obscured genealogical relationships in the Cape fossorial skink Acontias meleagris meleagris species complex. Currently the species complex contains four dubious operational taxonomic units (A. meleagris meleagris, A. m. orientalis, A. percivali tasmani and the morph lineicauda) with poorly defined species boundaries. In the present study we examine the evolutionary relationships within the species complex by sampling 24 skink populations from the known geographical distribution in the Western and Eastern Cape provinces of South Africa, representing a total of 119 specimens. We used partial sequence data derived from two mitochondrial DNA genes, 16S rRNA and COI, and one nuclear DNA gene, intron β‐fibrinogen (β‐fibint 7), to examine evolutionary relationships. Phylogenetic relationships were determined using both Maximum Parsimony (MP) and Bayesian inference (BI) from the combined mtDNA, nDNA and the total evidence data. Additionally we employed Maximum likelihood (ML) analyses on the total evidence data that comprised ～1.5 kb. Topologies derived from the combined mtDNA analyses were congruent with the total evidence analyses (mtDNA + nDNA) and retrieved five major clades with strong statistical support inferred from bootstrapping and posterior probabilities. The five clades were genealogically and geographically exclusive, diagnostic at both the mtDNA and nDNA level and characterized by pronounced sequence divergence, with no shared haplotypes between clades. Collectively these results suggest the presence of five putative cryptic operational taxonomic units within the A. meleagris meleagris species complex. Constraining the traditionally recognized taxa always retrieved a statistically worse topology suggesting that considerable taxonomic revision is required. Our results indicate that traditional morphological characters need to be reassessed to define the five novel lineages in the A. meleagris meleagris species complex. The phylogeographic pattern for the fossorial skinks we retrieved was novel compared to phylogeographic studies for codistributed above ground living taxa. These results suggest that the abiotic and biotic factors that impact subterranean taxa may differ from supraterranean taxa.