mtDNA analysis reveals the ongoing speciation on Greek populations of Microtus (Terricola) thomasi (Arvicolidae,Rodentia)

The present article extends our previous work on the phylogenetic history of Microtus (Terricola) thomasi, analysing cytochrome b, 12S rRNA, 16S rRNA and the control region in 65 Greek populations. The analysis revealed three clades: one grouping the populations from Peloponnisos (Southern Greece); the second, the populations from Agios (Ag.) Stefanos and Evvoia island (Central East Greece); and the third, all the remaining populations with no geographical substructure. Genetic distances were low for most populations, with only the populations of Evvoia and Ag. Stefanos being substantially distant. Thus, although this species has a recent colonization history and probably descends from a highly polymorphic ancestor, a monophyletic and highly differentiated lineage is formed in Greece and is distributed in Ag. Stefanos and Evvoia. Molecular differentiation, distinct geographical distribution and restriction of gene flow between this lineage and the rest of the Greek populations provide evidence for its probable subspecific status, Microtus (Tericola) thomasi atticus. A possible mechanism leading the differentiation process of the proposed subspecies is suggested, based on the displacement of this species in central Greece by its congeneric, probably better‐fitted Microtus (Microtus) guentheri and the subsequent separation of Ag. Stefanos and Evvoia from the remaining Greek populations. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 117–130.     

Old Aegean palaeoevents driving the diversification of an endemic isopod species (Oniscidea,Trachelipodidae)

The Greek endemic isopod species Trachelipus aegaeus is distributed in Aegean islands and the adjacent coastal parts of the Greek mainland. Major palaeogeographic events of the Aegean archipelago, such as the formation of the mid‐Aegean trench and the Messinian Salinity Crisis, have been often employed as major causal factors of evolutionary events and phylogeographic patterns exhibited by several taxa. Herein, we infer phylogenetic relationships among T. aegaeus populations using partial cytochrome oxidase subunit I (COI) and 16S rRNA sequences. Due to the poor preservation of the specimens, we propose a modified DNA extraction protocol, which returned highly positive results in terms of the quality of the total extracted DNA. We implement a calibrated molecular clock and path sampling analysis, using alternative palaeogeographic events and rates of substitution, to evaluate the biogeographic history of the species and to estimate the chronology of diversification events among its populations. Our results are clearly in favour of the scenario of the MAT triggering vicariance among most T. aegaus populations. Moreover, the large intraspecific genetic divergence (0–19% for COI and 0–20.3% for the 16S rRNA) and the overall phylogeographic patterns depicted herein seem not to have been obscured by more recent palaeogeological events. A role of dispersal, probably human‐aided, is assumed for certain ‘deviant’ cases.     

Phylogeographic diversification of antelope squirrels (Ammospermophilus) across North American deserts

We investigated the biogeographic history of antelope squirrels, genus Ammospermophilus, which are widely distributed across the deserts and other arid lands of western North America. We combined range‐wide sampling of all currently recognized species of Ammospermophilus with a multilocus data set to infer phylogenetic relationships. We then estimated divergence times within identified clades of Ammospermophilus using fossil‐calibrated and rate‐calibrated molecular clocks. Lastly, we explored generalized distributional changes of Ammospermophilus since the last glacial maximum using species distribution models, and assessed responses to Quaternary climate change by generating demographic parameter estimates for the three wide‐ranging clades of A. leucurus. From our phylogenetic estimates we inferred strong phylogeographic structure within Ammospermophilus and the presence of three well‐supported major clades. Initial patterns of historical divergence were coincident with dynamic alterations in the landscape of western North America, and the formation of regional deserts during the Late Miocene and Pliocene. Species distribution models and demographic parameter estimates revealed patterns of recent population expansion in response to glacial retreat. When combined with evidence from co‐distributed taxa, the historical biogeography of Ammospermophilus provides additional insight into the mechanisms that impacted diversification of arid‐adapted taxa across the arid lands of western North America. We propose species recognition of populations of the southern Baja California peninsula to best represent our current understanding of evolutionary relationships among genetic units of Ammospermophilus. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 949–967.     

Extensive genetic divergence in the East Asian natricine snake,Rhabdophis tigrinus (Serpentes: Colubridae), with special reference to prominent geographical differentiation of the mitochondrial cytochrome b gene in Japanese populations

We investigated intraspecific phylogenetic relationships in the natricine snake, Rhabdophis tigrinus. A partial sequence of mitochondrial cytochrome b gene (990 bp) was sequenced for 220 individuals from 112 populations. The phylogeny indicated monophyly of the Japanese populations against the continental and Taiwanese populations, sister relationships of the Japanese and continental populations, and monophyly of the whole species. The results strongly suggested substantial genetic divergences among population assemblages from those three regions. We thus consider both lateralis from the continent, which is often synonymized to R. tigrinus, and formosanus from Taiwan, which is usually regarded as a subspecies of the latter, as distinct full species based on the evolutionary species concept. In the Japanese populations, haplotypes were classified to in two major clades (I and II) that were parapatric to each other. Clade I consisted of three distinct subclades (I‐A, I‐B, and I‐C), of which the former two were parapatric with each other, whereas the latter was sympatric with each of the former two subclades. The geographical haplotype structure exhibited by the Japanese populations is likely to have resulted from a series of allopatric differentiations with rapid range extensions of resultant lineages, leading to secondary contact or further admixture of mitochondrial haplotype clades and subclades. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 395–408.     

Molecular phylogeny of the Eremias persica complex of the Iranian plateau (Reptilia: Lacertidae), based on mtDNA sequences

The Persian racerunner Eremias persica Blanford, 1875 is confined to the Iranian plateau, and forms one of the most widespread but rarely studied species of the family Lacertidae. With many local populations inhabiting a variety of habitats, and exhibiting considerable morphological, genetic, and ecological variations, it represents a species complex. We analysed sequences of mitochondrial cytochrome b and 12S ribosomal RNA (rRNA) genes derived from 13 geographically distant populations belonging to the E. persica complex. Using our knowledge of palaeogeographical events, a molecular clock was calibrated to assess the major events in fragmentation, radiation, and intraspecific variation. The sequence data strongly support a basal separation of the highland populations of western Iran from those of the open steppes and deserts, occurring in the east. The subsequent radiation, fragmentation, and evolution of these major assemblages have led to several discernable geographical lineages across the wide area of the Iranian plateau. The results indicate a middle‐Miocene origin for the clade as a whole. The first split, isolating the western and eastern clades, appears to have occurred 11–10 Mya. Further fragmentations and divergence within the major clades began about 8 Mya, with an evolutionary rate of 1.6% sequence divergence per million years among the lineages in the genes studied (combined data set). Molecular and morphological data strongly support a taxonomic revision of this species complex. At least four of the discovered clades should be raised to species, and two to subspecies, rank. © 2009 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 158 , 641–660.     

Congruence between molecular phylogeny and cuticular design in Echiniscoidea (Tardigrada,Heterotardigrada)

Although morphological characters distinguishing echiniscid genera and species are well understood, the phylogenetic relationships of these taxa are not well established. We thus investigated the phylogeny of Echiniscidae, assessed the monophyly of Echiniscus, and explored the value of cuticular ornamentation as a phylogenetic character within Echiniscus. To do this, DNA was extracted from single individuals for multiple Echiniscus species, and 18S and 28S rRNA gene fragments were sequenced. Each specimen was photographed, and published in an open database prior to DNA extraction, to make morphological evidence available for future inquiries. An updated phylogeny of the class Heterotardigrada is provided, and conflict between the obtained molecular trees and the distribution of dorsal plates among echiniscid genera is highlighted. The monophyly of Echiniscus was corroborated by the data, with the recent genus Diploechiniscus inferred as its sister group, and Testechiniscus as the sister group of this assemblage. Three groups that closely correspond to specific types of cuticular design in Echiniscus have been found with a parsimony network constructed with 18S rRNA data. © 2013 The Linnean Society of London     

Cryptic speciation along a bathymetric gradient

The deep ocean supports a highly diverse and mostly endemic fauna, yet little is known about how or where new species form in this remote ecosystem. How speciation occurs is especially intriguing in the deep sea because few obvious barriers exist that would disrupt gene flow. Geographic and bathymetric patterns of genetic variation can provide key insights into how and where new species form. We quantified the population genetic structure of a protobranch bivalve, Neilonella salicensis, along a depth gradient (2200–3800 m) in the western North Atlantic using both nuclear (28S and calmodulin intron) and mitochondrial (cytochrome c oxidase subunit I) loci. A sharp genetic break occurred for each locus between populations above 2800 m and below 3200 m, defining two distinct clades with no nuclear or mitochondrial haplotypes shared between depth regimes. Bayesian phylogenetic analyses provided strong support for two clades, separated by depth, within N. salicensis. Although no morphological divergence was apparent, we suggest that the depth‐related population genetic and phylogenetic divergence is indicative of a cryptic species. The frequent occurrence of various stages of divergence associated with species formation along bathymetric gradients suggests that depth, and the environmental gradients that attend changes in depth, probably play a fundamental role in the diversification of marine organisms, especially in deep water. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 897–913.     

Rapid diversification of colouration among populations of a poison frog isolated on sky peninsulas in the central cordilleras of Peru

Aim Comparison of Epipedobates bassleri ( Myers, 1987 ), which occurs on high‐altitude mountain ridges (‘sky peninsulas’) in the Andean transition zone and demonstrates high levels of divergence in colouration among populations, and Epipedobates hahneli ( Schulte, 1999 ), which occurs throughout the lowland regions of the Amazon basin and is morphologically conserved, using phylogenetic analysis of mitochondrial sequence data and comparison of colour pattern. Location Central cordilleras of Peru (near Tarapoto, San Martin). Methods DNA was extracted from individuals of E. bassleri from the central cordilleras of Peru, and from individuals of E. hahneli from across Peru. The cytochrome b mitochondrial gene region was amplified and sequenced for individuals of each species, and phylogenetic analysis was carried out using Bayesian inference. Genetic distances among populations and geographic distances of each species were examined and compared using Mantel tests. Parametric bootstrapping was used to test the monophyly of E. bassleri. Results Epipedobates bassleri formed a well‐supported monophyletic group and showed higher levels of genetic divergence among populations than was shown among populations of E. hahneli from the same region. Distinct clades representing different geographic regions were recovered for E. hahneli. Levels of divergence among more geographically distant populations of E. hahneli were higher than levels of divergence among E. bassleri populations. We found a significant correlation between genetic divergence and geographic distance as measured along a 1000‐m contour line, but not as measured by direct routes (crossing putative biogeographical barriers). Main conclusions Levels of genetic divergence were higher among populations of morphologically conservative E. hahneli than among populations of morphologically variable E. bassleri, suggesting rapid divergence in colouration among populations of E. bassleri. These patterns support previous arguments concerning the role of the montane transition zone between the high mountains and lowlands in divergence and speciation. High levels of both genetic and phenotypic divergence among populations of E. bassleri indicate that ecological or behavioural factors may be responsible for the high levels of colour variation seen among E. bassleri, but not among E. hahnleli, populations.     

Are there two cryptic species of the lesser jerboa Jaculus jaculus (Rodentia: Dipodidae) in Tunisia? Evidence from molecular,morphometric, and cytogenetic data

Jerboas belonging to the genus Jaculus are widely distributed rodents inhabiting Palearctic desert and semi‐desert areas. Previous studies on the lesser Egyptian jerboa Jaculus jaculus showed the existence of various morphological forms of controversial taxonomic status. They were sometimes related to two different species, J. jaculus and Jaculus deserti, although this has not been recognized in recent taxonomic updates. To clarify the systematic status of J. jaculus in Tunisia, we performed molecular (phylogenetic analyses of cytochrome b sequences), morphological (multivariate analyses of 13 skull measurements) and karyotypic (standard preparations from bone marrow cells) analyses on a number of specimens collected from ten localities. Our analyses revealed two monophyletic, well differentiated clades, with a mean genetic divergence value (K2P = 10.9 ± 0.01%), which is within the range of distances generally observed between rodent species. Morphometric analyses clearly separated populations of the two genetic clades from each other. However, karyotypes of individuals from both clades appeared similar. Individuals from both molecular clades/morphometric groups were found in sympatry in most of the localities sampled. These results, as obtained from a restricted area of the total distribution, suggest that there are two separate species within the currently accepted J. jaculus in Tunisia. Alternative hypotheses such as the occurrence of a strong, ancient phylogeographic structure, or the presence of pseudogenes, are also considered to account for the results obtained. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 673–686.     

Phylogeography of the freshwater crayfish Cherax destructor Clark (Parastacidae) in inland Australia: historical fragmentation and recent range expansion

The yabby, Cherax destructor Clark, is the most widespread species in the most widespread genus of Australian freshwater crayfish. It has a distribution that spans several distinct drainage basins and biogeographical regions within semiarid and arid inland Australia. Here we report a study designed to investigate patterns of genetic variation within the species and hypotheses put forward to account for its extensive distribution using DNA sequences from the mitochondrial 16S rRNA gene region. Results of phylogenetic analyses contradicted previous allozyme data and revealed relatively deep phylogenetic structure in the form of three geographically correlated clades. The degree of genetic divergences between clades (8–15 bp) contrasted with the relatively limited haplotype diversity within clades (1–3 bp). Network-based analyses confirmed these results and revealed genetic structure on both larger and more restricted geographical scales. Nevertheless some haplotypes and 1-step clades had large distributions, some of which crossed boundaries between river basins and aquatic biogeographical regions. Thus both older and more recent historical processes, including fragmentation on a larger geographical scale and more recent range expansion on a local scale, appear to be responsible for the observed pattern of genetic variation within C. destructor . These results support elements of alternative hypotheses previously put forward to account for the evolutionary history of C. destructor and the origin of its large distribution.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 83 , 539–550.     

New data on the distribution and genetic structure of Greek moles of the genus Talpa (Mammalia,Talpidae)

In this work, we employed molecular markers and confirmed the occurrence of two mole species, Talpa stankovici and Talpa europaea in Greece. For the first species, all analyses revealed three major phylogroups, exhibiting great genetic divergence, possibly due to a vicariant event. This event was probably the Peloponnisos’ insulation during the Pliocene and the following sea level fluctuations. A scenario of stable, large‐sized populations for a long period rather than rapid growth from small‐sized populations or founder events is supported. Additionally, T. stankovici’s distribution area is extended to a southernmost limit. For T. europaea there is evidence of low genetic divergence between Greek and Central and North European populations.     

Phylogenetic relationships and rate of early diversification of Australian Sphenomorphus group scincids (Scincoidea, Squamata)

The Australian Sphenomorphus group is a morphologically and ecologically diverse clade of lygosomine scincids, collectively comprising more than one‐half of the Australian scincid fauna. A previous phylogenetic analysis of mitochondrial 12S and 16S rRNA, and ND4 and adjacent tRNA sequences for a series of Australian Sphenomorphus group scincids recovered several well‐supported, major clades, although these were generally separated by relatively short branches associated with low support values. Applying a recently described methodology for inferring lineage‐level polytomies, I employ ATP synthetase‐β subunit intron sequences and the existing mitochondrial (mt)DNA data set (with sequences for additional taxa) to assess the hypothesis that the poorly resolved basal relationships within the Australian Sphenomorphus group are a consequence of the major clades having originated essentially simultaneously. Phylogenetic analyses of the separate mtDNA and intron sequence data reveal a number of congruent clades, including Anomalopus, Calyptotis, Ctenotus, Lerista, the Eulamprus quoyii group, the Glaphyromorphus crassicaudis group (including Glaphyromorphus cracens, Glaphyromorphus darwiniensis, and Glaphyromorphus fuscicaudis), Glaphyromorphus gracilipes + Hemiergis, Coeranoscincus reticulatus + Ophioscincus truncatus + Saiphos, and Eulamprus amplus + Eulamprus tenuis + Gnypetoscincus + Nangura. The relationships among these clades indicated by the two data sets, however, are generally incongruent. Although this may be partially ascribed to error in estimating phylogenetic relationships due to insufficient data, some incongruence is evident when uncertainty in inferred relationships is allowed for. Moreover, the congruent clades are typically separated by very short branches, several having a length insignificantly different from zero. These results suggest that initial diversification of Australian Sphenomorphus group scincids was rapid relative to the substitution rates of the mtDNA and intron fragments considered, if not essentially simultaneous. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 92 , 347–366.     

Molecular phylogeny and phylogeography of the Greek populations of the genus Orthometopon (Isopoda, Oniscidea) based on mitochondrial DNA sequences

We infer phylogenetic relationships among isopod species of the genus Orthometopon distributed in the Greek area, comparing partial mitochondrial DNA sequences for cytochrome oxidase I (COI). All phylogenetic analyses produced topologically identical trees that revealed a well-resolved phylogeny. These trees support the monophyly of the genus Orthometopon , and suggest two clades that correspond to separate geographical regions (west and east of the mid-Aegean trench). However, the phylogenetic relationships among Greek populations of Orthometopon spp. are different from the presumed pattern on the basis of morphological evidence. The distinct geographical distribution of the major clades of the phylogenetic tree and its topology suggest a spatial and temporal sequence of phylogenetic separations, which coincide with some major palaeogeographical separations during the geological history of the Aegean Sea. The results stress the need for a reconsideration of the evolutionary history of Orthometopon species, which will help overcome difficulties encountered in classical taxonomy at the species level. © 2008 The Linnean Society of London, Zoological Journal of the Linnean Society , 2008, 152 , 707–715.     

Bioacoustic and genetic divergence in a frog with a wide geographical distribution

The study of intraspecific variation of acoustic signals and its relationship with genetic divergence is important for understanding the origin of divergence in communication systems. We studied geographical variation in the acoustic structure of advertisement calls from five populations of the four‐eyed frog, Pleurodema thaul, and its relationship with the genetic divergence among these populations. By analyzing temporal and spectral parameters of the advertisement calls, we report that the signals of northern, central, and southern populations have remarkable differences between them. A phylogeographical analysis from a mitochondrial DNA fragment demonstrated three phylogenetic groups coincident with those found with the bioacoustics analysis. Furthermore, bioacoustic and genetic distances show significant correlations after controlling for geographical distance. These results suggest that behavioural divergence among populations of P. thaul has a phylogenetic basis, supporting three evolutionary units within this species, as well as prompting the exploration of divergence processes in the sound communication system of this species. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 142–155.     

Phylogenetic analysis and phylogeography of the tetraploid rodent Tympanoctomys barrerae (Octodontidae): insights on its origin and the impact of Quaternary climate changes on population dynamics

The allotetraploid Tympanoctomys barrerae has a broad, patchy distribution around salt flats in western Argentina. To gain insights into its phylogenetic relationships, phylogeographical patterns, and origin, seven populations of T. barrerae and its allied taxa were studied through a 1075‐bp fragment of cytochrome b and cytochrome oxidase I. Matrilineal phylogenetic relationships were explored with maximum likelihood and Bayesian approaches. The intraspecific genealogy was inferred by median‐joining networks. The populational structure was assessed by molecular variance, and the demographic history through coalescence. The tree topology depicted sister‐group relationship between Octomys mimax to the pair T. barrerae–Pipanacoctomys aureus, suggesting P. aureus belongs to the maternal lineage that gave rise to T. barrerae. High degrees of intrapopulational variation and the several instances of interpopulational polyphyly suggest range shifts and secondary contact. Consistent with the phylogenetic results, the network analysis revealed two haplotypic lineages depicting genetic admixtures unrelated to the current geographical distribution, and a deep split with the southernmost lineage of Chubut. The origin of T barrerae was estimated at approximately 2.52 Mya, whereas the divergence estimate for Chubut coincides with the end of largest Patagonian glaciation, at 1.47 Mya. Apparently, this population remained isolated during a northward Pleistocenic range shift. The historical demographic patterns of the lineages of T. barrerae fit with a contraction–expansion model coincident with Quaternary cycling events. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 453–469.     

Hidden diversity of Euscorpius (Scorpiones: Euscorpiidae) in Greece revealed by multilocus species‐delimitation approaches

Phylogenetic analysis of the genus Euscorpius (Scorpiones: Euscorpiidae) across the Mediterranean region (86 specimens, 77 localities, four DNA markers: 16S rDNA, COI, COII, and ITS1), focusing on Greek fauna, revealed high variation, deep clade divergences, many cryptic lineages, paraphyly at subgenus level, and sympatry of several new and formerly known lineages. Numerous specimens from mainland and insular Greece, undoubtedly the least studied region of the genus' distribution, have been included. The reconstructed phylogeny covers representative taxa and populations across the entire genus of Euscorpius. The deepest clades detected within Euscorpius correspond (partially) to its current subgeneric division, outlining subgenera Tetratrichobothrius and Alpiscorpius. The rest of the genus falls into several clades, including subgenus Polytrichobothrius and a paraphyletic subgenus Euscorpius s.s. Several cryptic lineages are recovered, especially on the islands. The inadequacy of the morphological characters used in the taxonomy of the genus to delineate species is discussed. Finally, the time frame of differentiation of Euscorpius in the study region is estimated and the distributional patterns of the lineages are contrasted with those of other highly diversified invertebrate genera occurring in the study region. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 728–748.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Phylogeography of a pearl oyster (Pinctada maxima) across the Indo‐Australian Archipelago: evidence of strong regional structure and population expansions but no phylogenetic breaks

This study investigates the genetic structure and phylogeography of a broadcast spawning bivalve mollusc, Pinctada maxima, throughout the Indo‐West Pacific and northern Australia. DNA sequence variation of the mitochondrial cytochrome oxidase subunit I (COI) gene was analysed in 367 individuals sampled from nine populations across the Indo‐West Pacific. Hierarchical AMOVA indicated strong genetic structuring amongst populations (Φ_ST = 0.372, P < 0.001); however, sequence divergence between the 47 haplotypes detected was low (maximum 1.8% difference) and no deep phylogenetic divergence was observed. Results suggest the presence of genetic barriers isolating populations of the South China Sea and central Indonesian regions, which, in turn, show patterns of historical separation from northern Australian regions. In P. maxima, historical vicariance during Pleistocene low sea levels is likely to have restricted planktonic larval transport, causing genetic differentiation amongst populations. However, low genetic differentiation is observed where strong ocean currents are present and is most likely due to contemporary larval transport along these pathways. Geographical association with haplotype distributions may indicate signs of early lineage sorting arising from historical population separations, yet an absence of divergent phylogenetic clades related to geography could be the consequence of periodic pulses of high genetic exchange. We compare our results with previous microsatellite DNA analysis of these P. maxima populations, and discuss implications for future conservation management of this species. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 632–646.     

Molecular phylogeny and phylogeography of the Cuban cave-fishes of the genus Lucifuga: evidence for cryptic allopatric diversity

Underground environments are increasingly recognized as reservoirs of faunal diversity. Extreme environmental conditions and limited dispersal ability of underground organisms have been acknowledged as important factors promoting divergence between species and conspecific populations. However, in many instances, there is no correlation between genetic divergence and morphological differentiation. Lucifuga Poey is a stygobiotic fish genus that lives in Cuban and Bahamian caves. In Cuba, it offers a unique opportunity to study the influence of habitat fragmentation on the genetic divergence of stygobiotic species and populations. The genus includes four species and one morphological variant that have contrasting geographical distributions. In this study, we first performed a molecular phylogenetic analysis of the Lucifuga Cuban species using mitochondrial and nuclear markers. The mitochondrial phylogeny revealed three deeply divergent clades that were supported by nuclear and morphological characters. Within two of these main clades, we identified five lineages that are candidate cryptic species and a taxonomical synonymy between Lucifuga subterranea and Lucifuga teresinarum. Secondly, phylogeographic analysis using a fragment of the cytochrome b gene was performed for Lucifuga dentata, the most widely distributed species. We found strong geographical organization of the haplotype clades at different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. At a larger temporal scale, these processes could also explain the diversification and the distribution of the different species.     

Phylogeny of the spider genus Ixchela Huber, 2000 (Araneae: Pholcidae) based on morphological and molecular evidence (CO1 and 16S), with a hypothesized diversification in the Pleistocene

The genus Ixchela Huber is composed of 20 species distributed from north‐eastern Mexico to Central America, including the five new species described here from Mexico: I xchela azteca sp. nov. , I xchela jalisco sp. nov. , I xchela mendozai sp. nov. , I xchela purepecha sp. nov. and I xchela tlayuda sp. nov. We test the monophyly and investigate the phylogenetic relationships among species of the genus Ixchela using morphological and molecular data. Parsimony (PA) analysis of 24 taxa and 40 morphological characters with equal and implied weights supported the monophyly of Ixchela with eight morphological synapomorphies. The PA analyses with equal and implied weights, and separate Bayesian inference (BI) analyses for the CO1 gene (506 characters), concatenated gene fragments CO1 + 16S (885 characters), morphology + CO1 (546 characters) and the combined evidence data set (morphology + CO1 + 16S) (925 characters) support the monophyly of Ixchela. Our preferred topology shows two large clades; clade 1 has a natural distribution in the Mesoamerican biotic component, whereas clade 2 predominates in the Mexican Montane biotic component. The genus Ixchela diverged in the late Miocene, and the divergence between the internal clades in the genus occurred in the late Pliocene; by contrast, most of the speciation events seem to have occurred mainly during the Pleistocene, where climatic changes brought on by repeated glaciations played an important role in the diversification of the genus. © 2015 The Linnean Society of London