Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

Genetic divergence between cave and surface populations of Astyanax in Mexico (Characidae,Teleostei)

A study of genetic diversity at microsatellite loci and the mitochondrial DNA (mtDNA) cytochrome b gene was carried out to assess genetic relationships among four Mexican cave (Pachon, Sabinos, Tinaja, Chica) and four surface populations of Astyanax fasciatus (Characidae) from northeast Mexico and the Yucatan. With the exception of Chica, the cave populations were all characterized by extremely low microsatellite variability, which most likely resulted from bottleneck events. Population analyses of the microsatellite data indicated no measurable levels of gene flow between all cave and surface populations (F(ST) > 0.0707). Phylogenetic analyses of mtDNA data showed that only two cave populations - Sabinos and Tinaja - group together to the exclusion of surface populations. From the microsatellite data these cave populations cluster with the Pachon cave fish population. The mtDNA thus appears to have been replaced in Pachon because of introgressive hybridization. It is likely that these three cave populations have descended from a surface ancestor in common with current surface populations, rather than evolving recently from one of the extant surface populations. Like Pachon, the Chica population clustered with the surface populations according to mtDNA data, but was not clearly associated with either the surface or the other cave populations according to the microsatellite data. Our data indicate that the Chica population evolved recently from a surface population, and subsequently hybridized with a phylogenetically older cave population. In conclusion, both the microsatellite and mtDNA data suggest multiple origins of cave populations and the Chica and Sabinos/Tinaja/Pachon were founded after at least two independent invasions from surface populations.     

Phylogeography of the pine processionary moth Thaumetopoea wilkinsoni in the Near East

Phylogeographic structure of the eastern pine processionary moth Thaumetopoea wilkinsoni was explored in this study by means of nested clade phylogeographic analyses of COI and COII sequences of mitochondrial DNA and Bayesian estimates of divergence times. Intraspecific relationships were inferred and hypotheses tested to understand historical spread patterns and spatial distribution of genetic variation. Analyses revealed that all T. wilkinsoni sequences were structured in three clades, which were associated with two major biogeographic events, the colonization of the island of Cyprus and the separation of southwestern and southeastern Anatolia during the Pleistocene. Genetic variation in populations of T. wilkinsoni was also investigated using amplified fragment length polymorphisms and four microsatellite loci. Contrasting nuclear with mitochondrial data revealed recurrent gene flow between Cyprus and the mainland, related to the long-distance male dispersal. In addition, a reduction in genetic variability was observed at both mitochondrial and nuclear markers at the expanding boundary of the range, consistent with a recent origin of these populations, founded by few individuals expanding from nearby localities. In contrast, several populations fixed for one single mitochondrial haplotype showed no reduction in nuclear variability, a pattern that can be explained by recurrent male gene flow or selective sweeps at the mitochondrial level. The use of both mitochondrial and nuclear markers was essential in understanding the spread patterns and the population genetic structure of T. wilkinsoni, and is recommended to study colonizing species characterized by sex-biased dispersal.     

MATRIARCHAL POPULATION GENETIC STRUCTURE IN AN AVIAN SPECIES WITH FEMALE NATAL PHILOPATRY

We employ mitochondrial (mt) DNA markers to examine the matrilineal component of population genetic structure in the snow goose Chen caerulescens. From banding returns, it is known that females typically nest at their natal or prior nest site, whereas males pair with females on mixed wintering grounds and mediate considerable nuclear gene flow between geographically separate breeding colonies. Despite site philopatry documented for females, mtDNA markers show no clear distinctions between nesting populations across the species' range from Wrangel Island, USSR to Baffin Island in the eastern Canadian Arctic. Two major mtDNA clades (as well as rare haplotypes) are distributed widely and provide one of the few available examples of a phylogeographic pattern in which phylogenetic discontinuity in a gene tree exists without obvious geographic localization within a species' range. The major mtDNA clades may have differentiated in Pleistocene refugia, and colonized current nesting sites through recent range expansion via pulsed or continual low-level dispersal by females. The contrast between results of banding returns and mtDNA distributions in the snow goose raises general issues regarding population structure: direct contemporary observations on dispersal and gene flow can in some cases convey a misleading impression of phylogeographic population structure, because they fail to access the evolutionary component of population connectedness; conversely, geographic distributions of genetic markers can provide a misleading impression of contemporary dispersal and gene flow because they retain a record of evolutionary events and past demographic parameters that may differ from those of the present. An understanding of population structure requires integration of both evolutionary (genetic) and contemporary (direct observational) perspectives.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Parallel speciation in Astyanax cave fish (Teleostei) in Northern Mexico

We investigated differentiation processes in the Neotropical fish Astyanax that represents a model system for examining adaptation to caves, including regressive evolution. In particular, we analyzed microsatellite and mitochondrial data of seven cave and seven surface populations from Mexico to test whether the evolution of the cave fish represents a case of parallel evolution. Our data revealed that Astyanax invaded northern Mexico across the Trans-Mexican Volcanic Belt at least three times and that populations of all three invasions adapted to subterranean habitats. Significant differentiation was found between the cave and surface populations. We did not observe gene flow between the strongly eye and pigment reduced old cave populations (Sabinos, Tinaja, Pachon) and the surface fish, even when syntopically occurring like in Yerbaniz cave. Little gene flow, if any, was found between cave populations, which are variable in eye and pigmentation (Micos, Chica, Caballo Moro caves), and surface fish. This suggests that the variability is due to their more recent origin rather than to hybridization. Finally, admixture of the young Chica cave fish population with nuclear markers from older cave fish demonstrates that gene flow between populations that independently colonized caves occurs. Thus, all criteria of parallel speciation are fulfilled. Moreover, the microsatellite data provide evidence that two co-occurring groups with small sunken eyes and externally visible eyes, respectively, differentiated within the partly lightened Caballo Moro karst window cave and might represent an example for incipient sympatric speciation.     

Polymorphic microsatellites isolated from the cave fish Astyanax fasciatus

Six polymorphic microsatellite loci from a library of the tetra Astyanax fasciatus from Mexico were isolated. Amplification and heterozygosity were tested in four cave and four surface populations. These loci were developed for population genetic study to detect gene flow between cave and surface populations.     

Conservation genetics of eastern hellbenders <Emphasis Type="Italic">Cryptobranchus alleganiensis alleganiensis</Emphasis> in the Tennessee Valley

The eastern hellbender has declined across much of its range and is a candidate for listing under the Endangered Species Act. Some of the most viable remaining populations exist in the Southern Appalachian Region in the Tennessee Valley watershed; however these populations are highly isolated and fragmented, and occupy several physiographic provinces, suggesting they may exhibit significant levels of genetic differentiation. We investigated genetic and phylogeographic relationships among eastern hellbender populations across the Tennessee Valley, using nuclear microsatellite markers and mitochondrial sequence data. Our population genetic analyses of microsatellite data revealed a strong pattern of isolation by stream distance, and 4 genetically distinct populations. These four populations were mainly associated with major watersheds, although middle Tennessee samples were difficult to assign to any particular population. Our phylogeographic analysis of mtDNA resulted in a strongly supported monophyletic ingroup containing nine largely allopatric clades, which also largely corresponded to major watersheds. Our findings suggest that hellbenders from different watersheds in the Tennessee Valley should be recognized as genetically distinct populations, and care should be taken to balance the needs of rescuing declining populations with translocation or headstart programs, while also preserving genetic diversity across the region.     

The Effect of Pleistocene Climate Fluctuations on Distribution of European Abalone (<Emphasis Type="Italic">Haliotis tuberculata</Emphasis>), Revealed by Combined Mitochondrial and Nuclear Marker Analyses

The genetic differentiation among the populations of the European abalone Haliotis tuberculata was investigated using different markers to better understand the evolutionary history and exchanges between populations. Three markers were used: mitochondrial cytochrome oxidase I (COI), the sperm lysin nuclear gene, and eight nuclear microsatellites. These markers present different characteristics concerning mutation rate and inheritance, which provided complementary information about abalone history and gene diversity. Genetic diversity and relationships among subspecies were calculated from a sample of approximately 500 individuals, collected from 17 different locations in the north-eastern Atlantic Ocean, Macaronesia, and Mediterranean Sea. COI marker was used to explore the phylogeny of the species with a network analysis and two phylogenetic methods. The analysis revealed 18 major haplotypes grouped into two distinct clades with a pairwise sequence divergence up to 3.5 %. These clades do not correspond to subspecies but revealed many contacts along Atlantic coast during the Pleistocene interglaciations. The sperm lysin gene analysis separated two different subtaxa: one associated to Macaronesian islands, and the other to all other populations. Moreover, a small population of the northern subtaxon was isolated in the Adriatic Sea—probably before the separation of the two lineages—and evolved independently. Microsatellites were analyzed by different genetics methods, including the Bayesian clustering method and migration patterns analysis. It revealed genetically distinct microsatellite patterns among populations from Mediterranean Sea, Brittany and Normandy, Morocco, and Canary and Balearic islands. Gene flow is asymmetric among the regions; the Azores and the Canary Islands are particularly isolated and have low effective population sizes. Our results support the hypothesis that climate changes since the Pleistocene glaciations have played a major role in the geographic distribution of the European abalone. Traces of these events related to maternal inheritance were shown on COI marker.     

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.     

Speciation with gene flow and the genetics of habitat transitions

Whether speciation can advance to completion in the face of initially high levels of gene flow is a very controversial topic in evolutionary biology. Extensive gene exchange is generally considered to homogenize populations and counteract divergence. Moreover, the role of introgressive hybridization in evolution remains largely unexplored in animals, particularly in freshwater zooplankton in which allopatric speciation is considered to be the norm. Our work investigates the genetic structure of two young ecological species: the pond species, Daphnia pulex and the lake species, Daphnia pulicaria. Phylogenetic and population genetics analyses were conducted on mitochondrial NADH dehydrogenase 5 (ND5) gene, the nuclear Lactate dehydrogenase (Ldh) gene and 21 nuclear microsatellite markers in 416 individuals from habitats with various degrees of permanence. The strong and consistent phylogenetic discordance between nuclear and mitochondrial markers suggests a complex evolutionary history of multiple independent habitat transition events that involved hybridization and introgression between lake and pond Daphnia. On the other hand, the low level of contemporary gene flow between adjacent populations indicates the presence of effective habitat isolating barriers. The Daphnia system provides strong evidence for a divergence-with-gene flow speciation model that involves multiple habitat transition events.     

Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.     

Deciphering the global phylogeography of a coastal shrub (Scaevola taccada) reveals the influence of multiple forces on contemporary population structure

The phylogeography of coastal plant species is heavily influenced by past sealevel fluctuations, dispersal barriers, and life-history traits, such as long-distance dispersal ability of the propagules. Unlike the widely studied mangroves, phylogeographic patterns have remained mostly obscure for other coastal plant species. In this study, we sampled 42 populations of Scaevola taccada (Gaertn.) Roxb., a coastal shrub of the family Goodeniaceae, from 17 countries across its distribution range. We used five chloroplast DNA (cpDNA) and 14 nuclear microsatellite (simple sequence repeat [SSR]) markers to assess the influence of abiotic factors and population genetic processes on the phylogeographic pattern of the species. Geographical distribution of cpDNA haplotypes suggests that the species originated in Australia, followed by historical dispersal and expansion of its geographic range. Multiple abiotic factors, including the sealevel changes during the Pleistocene, the presence of landmasses like the Malay Peninsula, and contemporary oceanic circulation patterns, restricted gene flow between geographically distinct populations, thereby creating low haplotype diversity and a strong population structure. Population genetic processes acted on these isolated populations, leading to high nuclear genetic diversity and population differentiation, as revealed from analyzing the polymorphic SSR loci. Although genetic divergence was mostly concordant between cpDNA and SSR data, asymmetrical gene flow and ancestral polymorphism could explain the discordance in the detailed genetic structure. Overall, our findings indicate that abiotic factors and population genetic processes interactively influenced the evolutionary history and current phylogeographic pattern of S. taccada across its distribution range.     

Phylogeographic breaks without geographic barriers to gene flow

Abstract.— The spatial distribution of genetic markers can be useful both in estimating patterns of gene flow and in reconstructing biogeographic history, particularly when gene genealogies can be estimated. Genealogies based on nonrecombining genetic units such as mitochondrial and chloroplast DNA often consist of geographically separated clades that come into contact in narrow regions. Such phylogeographic breaks are usually assumed to be the result of long-term barriers to gene flow. Here I show that deep phylogeographic breaks can form within a continuously distributed species even when there are no barriers to gene flow. The likelihood of observing phylogeographic breaks increases as the average individual dispersal distance and population size decrease. Those molecular markers that are most likely to show evidence of real geographic barriers are also most likely to show phylogeographic breaks that formed without any barrier to gene flow. These results might provide an explanation as to why some species, such as the greenish warblers ( Phylloscopus trochiloides ), have phylogeographic breaks in mitochondrial or chloroplast DNA that do not coincide with sudden changes in other traits.     

Genetic diversity and phylogeography of Daphnia similoides sinensis located in the middle and lower reaches of the Yangtze River

Geographical patterns, climate, and environmental change have important influences on the distribution and spread of aquatic organisms. However, the relationships between the geographical pattern and phylogenetics of Daphnia as well as environmental change are not well known. The genetic diversity and phylogeography of seven D. similoides sinensis populations located in the middle and lower reaches of the Yangtze River were investigated based on the combination of mitochondrial (COI gene) and nuclear (14 microsatellite primers) markers. Based on the mitochondrial gene markers, D. similoides sinensis from the middle and lower reaches of the Yangtze River had one ancestral haplotype and two evolutionary clades. In addition, D. similoides sinensis population deviated from neutral evolution, showing signs of a bottleneck effect followed by population expansion. Based on the microsatellite markers, the seven D. similoides sinensis populations formed three main groups. The dendrogram (NJ/ME) showed that D. similoides sinensis based on the mitochondrial genes marker were obviously clustered two main clades, whereas there were three clades based on the microsatellite markers. Our results suggested that the habitat fragmentation due to the barrier of the dams and sluices promoted the genetic differentiation and phylogeography of D. similoides sinensis populations in the middle and lower reaches of the Yangtze River.     

Assessing genetic structure with multiple classes of molecular markers: a case study involving the introduced fire ant Solenopsis invicta.

We used 30 genetic markers of 6 different classes to describe hierarchical genetic structure in introduced populations of the fire ant Solenopsis invicta. These included four classes of presumably neutral nuclear loci (allozymes, codominant random amplified polymorphic DNAs (RAPDs), microsatellites, and dominant RAPDs), a class comprising two linked protein-coding nuclear loci under selection, and a marker of the mitochondrial DNA (mtDNA). Patterns of structure revealed by F statistics and exact tests of differentiation were highly concordant among the four classes of neutral nuclear markers, although the microsatellites were the most effective markers for detecting structure. The results from the mtDNA complemented those from the neutral nuclear markers by revealing that strong limitations to female-mediated gene flow were the cause of the local structure registered by the nuclear markers. The pattern of structure inferred from the selected nuclear loci was markedly different from the patterns derived from the other sets of markers but was predictable on the basis of the presumed mode of selection acting on these loci. In general, the results for all six classes of markers can be explained by known features of the social and reproductive biology of fire ants. Thus, the results from these diverse sets of markers, combined with detailed natural history data, provide an unusually complete picture of how the fundamental evolutionary forces of gene flow, drift, and selection govern the distribution of genetic variation within and between fire ant populations.     

Identifying main evolutionary mechanisms shaping genetic variation of Leuciscus cephalus L. 1758 (Cyprinidae) in western Greece: discordance between methods.

Numerous methods can be used in intraspecific phylogeographic studies to infer the evolutionary mechanisms that shaped observed genetic variation in populations. However, these methods are scarcely used jointly, and the evolutionary outcomes they could propose are not fully compared. In this study, using a chub (Leuciscus cephalus; Cyprinidae) mitochondrial DNA data set (13 populations in Western Greece, 14 haplotypes), we compare three distinct 'historical' methods that could possibly infer relative importance of basic evolutionary mechanisms (isolation vs migration) shaping genetic variation: the nested clade analysis, the psi-test and the 'mismatch distributions'. Taking together, interpretations of these analyses allow to draw a picture of the evolutionary history of chub in Western Greece based on isolation and genetic drift for higher clades. However, results issued each method can differ for low differentiated clades. We discuss such differences and suggest that methods should be used jointly in phylogeographic studies for a better evaluation of the evolutionary mechanisms that shaped genetic variation.     

Patterns of molecular genetic variation among African elephant populations

The highly threatened African elephants have recently been subdivided into two species, Loxodonta africana (savannah or bush elephant) and L. cyclotis (forest elephant) based on morphological and molecular studies. A molecular genetic assessment of 16 microsatellite loci across 20 populations (189 individuals) affirms species level genetic differentiation and provides robust genotypic assessment of species affiliation. Savannah elephant populations show modest levels of phylogeographic subdivision based on composite microsatellite genotype, an indication of recent population isolation and restricted gene flow between locales. The savannah elephants show significantly lower genetic diversity than forest elephants, probably reflecting a founder effect in the recent history of the savannah species.     

Riverine habitat specificity constrains dispersion in a Neotropical fish (Characidae) along Southern Brazilian drainages

Freshwater fishes often display a marked phylogeographic structure strongly associated with historical and ecological changes in the aquatic environment. Different ecological conditions in the same river drainage may act as permeable barriers to dispersion and gene flow. Previous studies recognized two discrete spatial components for the ichthyofauna in the freshwater coastal drainages of southern Brazil: the lowland fish fauna in the lagoons and the fish fauna of the rivers flowing in the valleys. In order to test if the coastal lagoons may limit the dispersion of a riverine species, we describe the phylogeographic structure among populations of Cyanocharax itaimbe, a species endemic to this region. We analysed 55 specimens characterized for two mitochondrial and one nuclear genes. Sequences were analysed using gene trees and species tree approaches, together with standard population genetics methods. Molecular analyses indicated three evolutionary groups which diverged from each other between an estimated 1,600,000 and 450,000 years before the present. However, two currently isolated river systems share the same evolutionary clade, whereas a single drainage contains two different lineages. Our results indicate strong genetic structure among populations along with generally conserved morphology. The strong genetic structure among populations living in the same drainage system may be explained by ecological differences between lagoons and rivers (or palaeochannels) that act as barriers to dispersion.     

Phylogeography of surface and cave Astyanax (Teleostei) from Central and North America based on cytochrome b sequence data

Astyanax fasciatus has become a model organism for the study of regressive and adaptive evolution in cave animals. To fully understand these processes, it is important to have background information on the systematics and phylogeography of surface and cave populations of this species. Here we investigate the phylogeography of A. fasciatus in North and Central America and also the historical biogeography of this region. Phylogenetic analysis of part of the mtDNA cytochrome b gene from 26 surface and nine cave A. fasciatus populations revealed seven major clades, which, in principle, represent geographical patterns of distribution. However, the four strongly eye and pigment reduced cave populations, Piedras, Sabinos, Tinaja, and Curva, form a separate cluster, which is not sister group to the surface populations from the same locality. Similarly the Belizean populations do not cluster with their geographic neighbors from the Yucatan. The analyses indicate that there have been recurrent invasions of surface Astyanax from the south, that were most likely influenced by major climate changes during the Pleistocene. During this period, ancestors of the strongly eye and pigment reduced cave populations were able to survive underground as thermophilic relics when the surface populations became extinct. The high level of genetic divergence among the different clades shows that differing haplotype lineages must have reinvaded the surface waters from the south and/or back-colonized them from residual habitats and also penetrated into the caves. Nested clade analyses show that recurrent gene flow as well as historic processes like past fragmentation and range expansion have influenced current populations of A. fasciatus in Central and North America. Different haplotype clades of the phylogeny are not compatible with the present taxonomy of Astyanax and, therefore, we propose the application of a single systematic unit, called A. fasciatus.