Recent divergence with gene flow in Tennessee cave salamanders (Plethodontidae: Gyrinophilus) inferred from gene genealogies

Selection Maintains the Phenotypic Divergence of Cave and Surface Fish

AbstractGenetic divergence in the presence of gene flow has been well documented, but there is little information on the specific factors maintaining divergence. The present study investigates this in the Mexican tetra (Astyanax mexicanus), an excellent model for studying this question because surface and cave populations differ markedly in phenotype and genotype but are interfertile. Previous population studies documented significant gene flow among cave and surface populations, but they focused on analyses of neutral markers whose evolutionary dynamics likely differ from those of genes involved in cave adaptation. The present study advances our understanding of this question by focusing specifically on the genetics responsible for eye and pigmentation reduction, signature traits of cave populations. Direct observations of two cave populations over the course of 63 years verify that surface fish frequently move into the caves and even hybridize with the cave fish. Importantly, however, historical records show that surface alleles for pigmentation and eye size do not persist but are rapidly eliminated from the cave gene pool. It has been argued that regression of eyes and pigmentation was driven by drift, but the results of this study suggest that strong selection actively eliminates surface alleles from the cave populations.

Sensory evolution in a cavefish radiation: patterns of neuromast distribution and associated behaviour in Sinocyclocheilus (Cypriniformes: Cyprinidae)

The genus Sinocyclocheilus, comprising a large radiation of freshwater cavefishes, are well known for their presence of regressive features (e.g. variable eye reduction). Fewer constructive features are known, such as the expansion of the lateral line system (LLS), which is involved in detecting water movements. The precise relationship between LLS expansion and cave adaptation is not well understood. Here, we examine morphology and LLS-mediated behaviour in Sinocyclocheilus species characterized by broad variation in eye size, habitat and geographical distribution. Using live-staining techniques and automated behavioural analyses, we examined 26 Sinocyclocheilus species and quantified neuromast organ number, density and asymmetry within a phylogenetic context. We then examined how these morphological features may relate to wall-following, an established cave-associated behaviour mediated by the lateral line. We show that most species demonstrated laterality (i.e. asymmetry) in neuromast organs on the head, often biased to the right. We also found that wall-following behaviour was distinctive, particularly among eyeless species. Patterns of variation in LLS appear to correlate with the degree of eye loss, as well as geographical distribution. This work reveals that constructive LLS evolution is convergent across distant cavefish taxa and may mediate asymmetric behavioural features that enable survival in stark subterranean microenvironments.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species-all plethodontid and proteid salamanders-are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Evidence for ecological processes driving speciation among endemic lizards of Madagascar

Although genetic patterns produced by population isolation during speciation are well documented, the biogeographic and ecological processes that trigger speciation remain poorly understood. Alternative hypotheses for the biogeography and ecology of speciation include geographic isolation combined with niche conservation (soft allopatry) or parapatric distribution on an environmental gradient with niche divergence (ecological speciation). Here, we use species' distributions, environmental data, and two null models (the Random Translation and Rotation and the Background Similarity Test) to test these alternative hypotheses among 28 sister pairs of microendemic lizards in Madagascar. Our results demonstrate strong bimodal peaks along a niche divergence-conservation spectrum, with at least 25 out of 28 sister pairs exhibiting either niche conservation or divergence, and the remaining pairs showing weak ecological signals. Yet despite these significant results, we do not find strong associations of niche conservation with allopatric distributions or niche divergence with parapatric distributions. Our findings thus provide strong evidence of a role for ecological processes driving speciation, rather than the classic expectation of speciation through geographic isolation, but demonstrate that the link between ecological speciation and parapatry is complex and requires further analysis of a broader taxonomic sample to fully resolve.

Repeated colonization of caves leads to phenotypic convergence in catfishes (Siluriformes: Trichomycterus) at a small geographical scale

Across various animal groups, adaptation to the challenging conditions of cave environments has resulted in convergent evolution. We document a Neotropical cavefish system with ample potential to study questions related to convergent adaptation to cave environments at the population level. In the karstic region of the Andes of Santander, Colombia, cave-dwelling catfish in the genus Trichomycterus exhibit variable levels of reduction of eyes and body pigmentation relative to surface congeners. We tested whether cave-dwelling, eye-reduced, depigmented Trichomycterus from separate caves in Santander were the result of a single event of cave colonization and subsequent dispersal, or of multiple colonizations to caves by surface ancestors followed by phenotypic convergence. Using mitochondrial DNA sequences to reconstruct phylogenetic relationships, we found that caves in this region have been colonized independently by two separate clades. Additional events of cave colonization – and possibly recolonization of surface streams – may have occurred in one of the clades, where surface and cave-dwelling populations exhibit shallow differentiation, suggesting recent divergence or divergence with gene flow. We also identify potentially undescribed species and likely problems with the circumscription of named taxa. The system appears promising for studies on a wide range of ecological and evolutionary questions.

Speciation in Daphnia

The microcrustacean Daphnia is arguably one of the most studied zooplankton species, having a well understood ecology, life history, and a relatively well studied evolutionary history. Despite this wealth of knowledge, species boundaries within closely related species in this genus often remain elusive and the major evolutionary forces driving the diversity of daphniids remain controversial. This genus contains more than 80 species with multiple cryptic species complexes, with many closely related species able to hybridize. Here, we review speciation research in Daphnia within the framework of current speciation theory. We evaluate the role of geography, ecology, and biology in restricting gene flow and promoting diversification. Of the 253 speciation studies on Daphnia, the majority of studies examine geographic barriers (55%). While evidence shows that geographic barriers play a role in species divergence, ecological barriers are also probably prominent in Daphnia speciation. We assess the contribution of ecological and nonecological reproductive isolating barriers between closely related species of Daphnia and found that none of the reproductive isolating barriers are restricting gene flow completely. Research on reproductive isolating barriers has disproportionally focused on two species complexes, the Daphnia pulex and Daphnia longispina species complexes. Finally, we identify areas of research that remain relatively unexplored and discuss future research directions that build our understanding of speciation in daphniids.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

Living in isolation for almost 40 years: molecular divergence of the 28S rDNA and COI sequences between French and Polish populations of the cave beetle Speonomus normandi hydrophilus (Jeannel, 1907)

The paper gives the results of the first studies on the molecular divergence between native and non-native populations of Speonomus normandi hydrophilus (Jeannel, 1907). This species is endemic to Massif Arize in the Central Pyrenees (France), and represents highly specialised organisms that live underground. In 1982, one hundred specimens of S. normandi hydrophilus had been experimentally introduced into the Dzwonnica Cave (Poland). Since then, a numerous population has developed in the Towarna-Dzwonnica cave system, and the neighbouring Cabanowa Cave. After almost 40 years of isolation between native and non-native populations, the genetic variations were examined using the COI and 28S rDNA genes. Analyses of the haplotypes of 28S showed one common haplogroup, which confirms the origin of the Polish group. The differentiation of haplotypes for the COI marker was high for both the French and Polish populations. Altogether 18 haplotypes of this marker have been detected, 12 in the French population and 9 in the Polish. However, only a portion of the haplotypes is shared between the native and introduced population.

Higher Gene Flow in Sex-Related Chromosomes than in Autosomes during Fungal Divergence

Nonrecombining sex chromosomes are widely found to be more differentiated than autosomes among closely related species, due to smaller effective population size and/or to a disproportionally large-X effect in reproductive isolation. Although fungal mating-type chromosomes can also display large nonrecombining regions, their levels of differentiation compared with autosomes have been little studied. Anther-smut fungi from the Microbotryum genus are castrating pathogens of Caryophyllaceae plants with largely nonrecombining mating-type chromosomes. Using whole genome sequences of 40 fungal strains, we quantified genetic differentiation among strains isolated from the geographically overlapping North American species and subspecies of Silene virginica and S. caroliniana. We inferred that gene flow likely occurred at the early stages of divergence and then completely stopped. We identified large autosomal genomic regions with chromosomal inversions, with higher genetic divergence than the rest of the genomes and highly enriched in selective sweeps, supporting a role of rearrangements in preventing gene flow in genomic regions involved in ecological divergence. Unexpectedly, the nonrecombining mating-type chromosomes showed lower divergence than autosomes due to higher gene flow, which may be promoted by adaptive introgressions of less degenerated mating-type chromosomes. The fact that both mating-type chromosomes are always heterozygous and nonrecombining may explain such patterns that oppose to those found for XY or ZW sex chromosomes. The specific features of mating-type chromosomes may also apply to the UV sex chromosomes determining sexes at the haploid stage in algae and bryophytes and may help test general hypotheses on the evolutionary specificities of sex-related chromosomes.

Phenotypic diversity assessment within a major ex situ collection of wild endemic coffees in Madagascar

BACKGROUND AND AIMS

Like other clades, the Coffea genus is highly diversified on the island of Madagascar. The 66 endemic species have colonized various environments and consequently exhibit a wide diversity of morphological, functional and phenological features and reproductive strategies. The trends of interspecific trait variation, which stems from interactions between genetically defined species and their environment, still needed to be addressed for Malagasy coffee trees.

METHODS

Data acquisition was done in the most comprehensive ex situ collection of Madagascan wild Coffea. The structure of endemic wild coffees maintained in an ex situ collection was explored in terms of morphological, phenological and functional traits. The environmental (natural habitat) effect was assessed on traits in species from distinct natural habitats. Phylogenetic signal (Pagel's λ, Blomberg's K) was used to quantify trait proximities among species according to their phylogenetic relatedness.

KEY RESULTS

Despite the lack of environmental difference in the ex situ collection, widely diverging phenotypes were observed. Phylogenetic signal was found to vary greatly across and even within trait categories. The highest values were exhibited by the ratio of internode mass to leaf mass, the length of the maturation phase and leaf dry matter content (ratio of dry leaf mass to fresh leaf mass). By contrast, traits weakly linked to phylogeny were either constrained by the original natural environment (leaf size) or under selective pressures (phenological traits).

CONCLUSIONS

This study gives insight into complex patterns of trait variability found in an ex situ collection, and underlines the opportunities offered by living ex situ collections for research characterizing phenotypic variation.

Phenotypic plasticity as a mechanism of cave colonization and adaptation

A widely accepted model for the evolution of cave animals posits colonization by surface ancestors followed by the acquisition of adaptations over many generations. However, the speed of cave adaptation in some species suggests mechanisms operating over shorter timescales. To address these mechanisms, we used Astyanax mexicanus, a teleost with ancestral surface morphs (surface fish, SF) and derived cave morphs (cavefish, CF). We exposed SF to completely dark conditions and identified numerous altered traits at both the gene expression and phenotypic levels. Remarkably, most of these alterations mimicked CF phenotypes. Our results indicate that many cave-related traits can appear within a single generation by phenotypic plasticity. In the next generation, plasticity can be further refined. The initial plastic responses are random in adaptive outcome but may determine the subsequent course of evolution. Our study suggests that phenotypic plasticity contributes to the rapid evolution of cave-related traits in A. mexicanus.

The Mexican tetra is a fish that has two forms: a surface-dwelling form, which has eyes and silvery grey appearance, and a cave-dwelling form, which is blind and has lost its pigmentation. Recent studies have shown that the cave-dwelling form evolved rapidly within the last 200,000 years from an ancestor that lived at the surface. The recent evolution of the cave-dwelling form of the tetra poses an interesting evolutionary question: how did the surface-dwelling ancestor of the tetra quickly adapt to the new and challenging environment found in the caves?

‘Phenotypic plasticity’ is a phenomenon through which a single set of genes can produce different observable traits depending on the environment. An example of phenotypic plasticity occurs in response to diet: in animals, poor diets can lead to an increase in the size of the digestive organs and to the animals eating more. To see if surface-dwelling tetras can quickly adapt to cave environments through phenotypic plasticity, Bilandžija et al. have exposed these fish to complete darkness (the major feature of the cave environment) for two years. After spending up to two years in the dark, these fish were compared to normal surface-dwelling and cave-dwelling tetras.

Results revealed that surface-dwelling tetras raised in the dark exhibited traits associated with cave-dwelling tetras. These traits included changes in the activity of many genes involved in diverse processes, resistance to starvation, metabolism, and levels of hormones and molecules involved in neural signaling, which could lead to changes in behavior. However, the fish also exhibited traits, including an increase in the cells responsible for pigmentation, that would have no obvious benefit in the darkness. Even though the changes observed require no genetic mutations, they can help or hinder the fish’s survival once they occur, possibly determining subsequent evolution. Thus, a trait beneficial for surviving in the dark that appears simply through phenotypic plasticity may eventually be selected for and genetic mutations that encode it more reliably may appear too.

These results shed light on how species may quickly adapt to new environments without accumulating genetic mutations, which can take hundreds of thousands of years. They also may help to explain how colonizer species succeed in challenging environments. The principles described by Bilandžija et al. can be applied to different organisms adapting to new environments, and may help understand the role of phenotypic plasticity in evolution.

Geographically structured genetic diversity in the cave beetle Darlingtonea kentuckensis Valentine, 1952 (Coleoptera, Carabidae, Trechini, Trechina)

Cave beetles of the eastern USA are one of many poorly studied groups of insects and nearly all previous work delimiting species is based solely on morphology. This study assesses genetic diversity in the monotypic cave carabid beetle genus DarlingtoneaValentine 1952, to test the relationship between putative geographical barriers to subterranean dispersal and the boundaries of genetically distinct groups. Approximately 400bp of the mitochondrial cytochrome oxidase I (COI) gene was sequenced from up to four individuals from each of 27 populations, sampled from caves along the escarpments of the Mississippian and Cumberland plateaus in eastern Kentucky, USA. The 81 individuals sequenced yielded 28 unique haplotypes. Hierarchical analyses of molecular variance (AMOVA) within and among geographically defined groups tested two a priori hypotheses of structure based on major and minor river drainages, as well as genetic distance clusters defined a posteriori from an unrooted analysis. High genetic differentiation (FST) between populations was found across analyses. The influence of isolation by distance could potentially account for much but not all of the variation found among geographically defined groups at both levels. High variability among the three northernmost genetic clusters (FCT), low variability among populations within clusters (FSC), and low within-cluster Mantel correlations indicate the importance of unidentified likely intra-karst barriers to gene flow separating closely grouped cave populations. Overall phylogeographic patterns are consistent with previous evidence of population isolation among cave systems in the region, revealing geographically structured cryptic diversity in Darlingtonea over its distribution. The landscape features considered a priori in this study were not predictive of the genetic breaks among the three northern clusters, which are genetically distinct despite their close geographic proximity.

Process-based species delimitation leads to identification of more biologically relevant species

Most approaches to species delimitation to date have considered divergence-only models. Although these models are appropriate for allopatric speciation, their failure to incorporate many of the population-level processes that drive speciation, such as gene flow (e.g., in sympatric speciation), places an unnecessary limit on our collective understanding of the processes that produce biodiversity. To consider these processes while inferring species boundaries, we introduce the R-package delimitR and apply it to identify species boundaries in the reticulate taildropper slug (Prophysaon andersoni). Results suggest that secondary contact is an important mechanism driving speciation in this system. By considering process, we both avoid erroneous inferences that can be made when population-level processes such as secondary contact drive speciation but only divergence is considered, and gain insight into the process of speciation in terrestrial slugs. Further, we apply delimitR to three published empirical datasets and find results corroborating previous findings. Finally, we evaluate the performance of delimitR using simulation studies, and find that error rates are near zero when comparing models that include lineage divergence and gene flow for three populations with a modest number of Single Nucleotide Polymorphisms (SNPs; 1500) and moderate divergence times (<100,000 generations). When we apply delimitR to a complex model set (i.e., including divergence, gene flow, and population size changes), error rates are moderate (∼0.15; 10,000 SNPs), and, when present, misclassifications occur among highly similar models.

Polyploid plants have faster rates of multivariate niche differentiation than their diploid relatives

Polyploid speciation entails substantial and rapid postzygotic reproductive isolation of nascent species that are initially sympatric with one or both parents. Despite strong postzygotic isolation, ecological niche differentiation has long been thought to be important for polyploid success. Using biogeographic data from across vascular plants, we tested whether the climatic niches of polyploid species are more differentiated than their diploid relatives and if the climatic niches of polyploid species differentiated faster than those of related diploids. We found that polyploids are often more climatically differentiated from their diploid parents than the diploids are from each other. Consistent with this pattern, we estimated that polyploid species generally have higher rates of multivariate niche differentiation than their diploid relatives. In contrast to recent analyses, our results confirm that ecological niche differentiation is an important component of polyploid speciation and that niche differentiation is often significantly faster in polyploids.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Morphometrics and phylogeography of the cave-obligate land snail Helicodiscus barri (Gastropoda, Stylommatophora, Helicodiscidae)

Molecular studies have recently led to the detection of many cryptic species complexes within morphologically ambiguous species formerly undescribed by the scientific community. Organisms such as land snails are at a particularly higher risk of species misidentification and misinterpretation, in that gastropod systematics are based almost entirely on external shell morphology. Subterranean ecosystems are associated with especially high degrees of cryptic speciation, largely owing to the abiotic similarities of these systems. In this study, we attempt to diagnose the potential cryptic diversity in the troglobitic land snail Helicodiscusbarri. Land snails are generally associated with having low vagility, and as such this species’ broad, mosaic distribution indicates the misdiagnosis of this organism as a single species. We analyze both mitochondrial (16S, CO1) and nuclear (28S, H3) genetic data for 23 populations. Phylogeny for H.barri was reconstructed using both maximum-likelihood and Bayesian approaches to assess relationships among populations, and two species delimitation methods (mPTP and ABGD) were used to detect the presence of unique molecular operational taxonomic units (MOTUs). Species delimitation results revealed seven and sixteen MOTUs respectively, suggesting the presence of several cryptic lineages within H.barri. To assess how external shell morphology corresponds with patterns of genetic and environmental variation, two morphometric approaches were used incorporating 115 shells from 31 populations. Both morphometric approaches reveal a significant environmental influence on shell morphology, and one approach showed the significance of MOTU groups. We discuss the delimitation and morphometric results and additionally provide discussion on the taxonomic and conservation implications of this study.

Population Genetic Analyses Using 10 New Polymorphic Microsatellite Loci Confirms Genetic Subdivision within the Olm, Proteus anguinus

We provide a comparative population genetic study of the elusive amphibian, Proteus anguinus, by comparing the genetic diversity and divergence among 4 cave populations (96 individuals) sampled in the Dinaric Karst of Croatia. We developed 10 variable microsatellite markers using pyrosequencing and applied them to the 4 selected populations belonging to 4 different cave systems. The results showed strong genetic differentiation between the 4 caves corroborating with previous findings suggesting that Proteus might comprise several unrecognized taxa. Our results confirmed that gene flow should be high within the caves, whereas it is low between hydrographic systems since geological periods. Finally, we conclude that the high genetic subdivision suggests the necessity of treating the 4 studied Proteus populations as evolutionary significant units.

Mitochondrial sequence data indicate “Vicariance by Erosion” as a mechanism of species diversification in North American Ptomaphagus (Coleoptera, Leiodidae, Cholevinae) cave beetles

Small carrion beetles (Coleoptera: Leiodidae: Cholevinae) are members of cave communities around the world and important models for understanding the colonization of caves, adaptation to cave life, and the diversification of cave-adapted lineages. We developed a molecular phylogeny to examine the diversification of the hirtus-group of the small carrion beetle genus Ptomaphagus. The hirtus-group has no surface-dwelling members; it consists of 19 short-range endemic cave- and soil-dwelling species in the central and southeastern United States of America. Taxonomic, phylogenetic and biogeographic data were previously interpreted to suggest the hirtus-group diversified within the past 350,000 years through a series of cave colonization and speciation events related to Pleistocene climate fluctuations. However, our time-calibrated molecular phylogeny resulting from the analysis of 2,300 nucleotides from five genes across three mitochondrial regions (cox1, cytb, rrnL-trnL-nad1) for all members of the clade paints a different picture. We identify three stages of diversification in the hirtus-group: (1) ~10 million years ago (mya), the lineage that develops into P.shapardi, a soil-dwelling species from the Ozarks, diverged from the lineage that gives rise to the 18 cave-obligate members of the group; (2) between 8.5 mya and 6 mya, seven geographically distinct lineages diverged across Kentucky, Tennessee, Alabama and Georgia; six of these lineages represent a single species today, whereas (3) the ‘South Cumberlands’ lineage in Tennessee and Alabama diversified into 12 species over the past ~6 my. While the events triggering diversification during the first two stages remain to be determined, the distributions, phylogenetic relationships and divergence times in the South Cumberlands lineage are consistent with populations being isolated by vicariant events as the southern Cumberland Plateau eroded and fragmented over millions of years.

At the confluence of vicariance and dispersal: Phylogeography of cavernicolous springtails (Collembola: Arrhopalitidae, Tomoceridae) codistributed across a geologically complex karst landscape in Illinois and Missouri

The processes of vicariance and dispersal are central to our understanding of diversification, yet determining the factors that influence these processes remains a significant challenge in evolutionary biology. Caves offer ideal systems for examining the mechanisms underlying isolation, divergence, and speciation. Intrinsic ecological differences among cavernicolous organisms, such as the degree of cave dependence, are thought to be major factors influencing patterns of genetic isolation in caves. Using a comparative phylogeographic approach, we employed mitochondrial and nuclear markers to assess the evolutionary history of two ecologically distinct groups of terrestrial cave-dwelling springtails (Collembola) in the genera Pygmarrhopalites (Arrhopalitidae) and Pogonognathellus (Tomoceridae) that are codistributed in caves throughout the Salem Plateau-a once continuous karst region, now bisected by the Mississippi River Valley in Illinois and Missouri. Contrasting phylogeographic patterns recovered for troglobiotic Pygmarrhopalites sp. and eutroglophilic Pogonognathellus sp. suggests that obligate associations with cave habitats can restrict dispersal across major geographic barriers such as rivers and valleys, but may also facilitate subterranean dispersal between neighboring cave systems. Pygmarrhopalites sp. populations spanning the Mississippi River Valley were estimated to have diverged 2.9-4.8 Ma, which we attribute to vicariance resulting from climatic and geological processes involved in Mississippi River Valley formation beginning during the late Pliocene/early Pleistocene. Lastly, we conclude that the detection of many deeply divergent, morphologically cryptic, and microendemic lineages highlights our poor understanding of microarthropod diversity in caves and exposes potential conservation concerns.

From top to bottom: Do Lake Trout diversify along a depth gradient in Great Bear Lake, NT, Canada?

Depth is usually considered the main driver of Lake Trout intraspecific diversity across lakes in North America. Given that Great Bear Lake is one of the largest and deepest freshwater systems in North America, we predicted that Lake Trout intraspecific diversity to be organized along a depth axis within this system. Thus, we investigated whether a deep-water morph of Lake Trout co-existed with four shallow-water morphs previously described in Great Bear Lake. Morphology, neutral genetic variation, isotopic niches, and life-history traits of Lake Trout across depths (0-150 m) were compared among morphs. Due to the propensity of Lake Trout with high levels of morphological diversity to occupy multiple habitat niches, a novel multivariate grouping method using a suite of composite variables was applied in addition to two other commonly used grouping methods to classify individuals. Depth alone did not explain Lake Trout diversity in Great Bear Lake; a distinct fifth deep-water morph was not found. Rather, Lake Trout diversity followed an ecological continuum, with some evidence for adaptation to local conditions in deep-water habitat. Overall, trout caught from deep-water showed low levels of genetic and phenotypic differentiation from shallow-water trout, and displayed higher lipid content (C:N ratio) and occupied a higher trophic level that suggested an potential increase of piscivory (including cannibalism) than the previously described four morphs. Why phenotypic divergence between shallow- and deep-water Lake Trout was low is unknown, especially when the potential for phenotypic variation should be high in deep and large Great Bear Lake. Given that variation in complexity of freshwater environments has dramatic consequences for divergence, variation in the complexity in Great Bear Lake (i.e., shallow being more complex than deep), may explain the observed dichotomy in the expression of intraspecific phenotypic diversity between shallow- vs. deep-water habitats. The ambiguity surrounding mechanisms driving divergence of Lake Trout in Great Bear Lake should be seen as reflective of the highly variable nature of ecological opportunity and divergent natural selection itself.

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open- and closed-canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade-off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open- and closed-canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.

Social interactions predict genetic diversification: an experimental manipulation in shorebirds

Mating strategy and social behavior influence gene flow and hence affect levels of genetic differentiation and potentially speciation. Previous genetic analyses of closely related plovers Charadrius spp. found strikingly different population genetic structure in Madagascar: Kittlitz’s plovers are spatially homogenous whereas white-fronted plovers have well segregated and geographically distinct populations. Here, we test the hypotheses that Kittlitz’s plovers are spatially interconnected and have extensive social interactions that facilitate gene flow, whereas white-fronted plovers are spatially discrete and have limited social interactions. By experimentally removing mates from breeding pairs and observing the movements of mate-searching plovers in both species, we compare the spatial behavior of Kittlitz’s and white-fronted plovers within a breeding season. The behavior of experimental birds was largely consistent with expectations: Kittlitz’s plovers travelled further, sought new mates in larger areas, and interacted with more individuals than white-fronted plovers, however there was no difference in breeding dispersal. These results suggest that mating strategies, through spatial behavior and social interactions, are predictors of gene flow and thus genetic differentiation and speciation. Our study highlights the importance of using social behavior to understand gene flow. However, further work is needed to investigate the relative importance of social structure, as well as intra- and inter-season dispersal, in influencing the genetic structures of populations.

Ecological speciation in darkness? Spatial niche partitioning in sibling subterranean spiders (Araneae : Linyphiidae : Troglohyphantes)

Speciation in subterranean habitats is commonly explained as the result of divergent selection in geographically isolated populations; conversely, the contribution of niche partitioning in driving subterranean species diversification has been rarely quantified. The present study integrated molecular and morphological data with a hypervolume analysis based on functional traits to investigate a potential case of parapatric speciation by means of niche differentiation in two sibling spiders inhabiting contiguous subterranean habitats within a small alpine hypogean site. Troglohyphantes giachinoi, sp. nov. and T. bornensis are diagnosed by small details of the genitalia, which are likely to be involved in a reproductive barrier. Molecular analysis recovered the two species as sister, and revealed a deep genetic divergence that may trace back to the Messinian (~6 million years ago). The hypervolume analysis highlighted a marginal overlap in their ecological niches, coupled with morphological character displacement. Specifically, T. giachinoi, sp. nov. exhibits morphological traits suitable for thriving in the smaller pores of the superficial network of underground fissures (Milieu Souterrain Superficiel, MSS), whereas T. bornensis shows a greater adaptation to the deep subterranean habitat. Our results suggest that different selective regimes within the subterranean environment, i.e. deep caves v. MSS, may either drive local speciation or facilitate contiguous distributions of independently subterranean adapted species.

Hybridization patterns in two contact zones of grass snakes reveal a new Central European snake species

Recent studies found major conflicts between traditional taxonomy and genetic differentiation of grass snakes and identified previously unknown secondary contact zones. Until now, little is known about gene flow across these contact zones. Using two mitochondrial markers and 13 microsatellite loci, we examined two contact zones. One, largely corresponding to the Rhine region, involves the western subspecies Natrix natrix helvetica and the eastern subspecies N. n. natrix, whereas in the other, more easterly, contact zone two lineages meet that are currently identified with N. n. natrix and N. n. persa. This second contact zone runs across Central Europe to the southern Balkans. Our analyses reveal that the western contact zone is narrow, with parapatrically distributed mitochondrial lineages and limited, largely unidirectional nuclear gene flow. In contrast, the eastern contact zone is very wide, with massive nuclear admixture and broadly overlapping mitochondrial lineages. In combination with additional lines of evidence (morphology, phylogeny, divergence times), we conclude that these differences reflect different stages in the speciation process and that Natrix helvetica should be regarded as a distinct species. We suggest a nomenclatural framework for presently recognized grass snake taxa and highlight the need for reconciling the conflicts between genetics and taxonomy.

Divergence and evolution of assortative mating in a polygenic trait model of speciation with gene flow

Assortative mating is an important driver of speciation in populations with gene flow and is predicted to evolve under certain conditions in few-locus models. However, the evolution of assortment is less understood for mating based on quantitative traits, which are often characterized by high genetic variability and extensive linkage disequilibrium between trait loci. We explore this scenario for a two-deme model with migration, by considering a single polygenic trait subject to divergent viability selection across demes, as well as assortative mating and sexual selection within demes, and investigate how trait divergence is shaped by various evolutionary forces. Our analysis reveals the existence of sharp thresholds of assortment strength, at which divergence increases dramatically. We also study the evolution of assortment via invasion of modifiers of mate discrimination and show that the ES assortment strength has an intermediate value under a range of migration-selection parameters, even in diverged populations, due to subtle effects which depend sensitively on the extent of phenotypic variation within these populations. The evolutionary dynamics of the polygenic trait is studied using the hypergeometric and infinitesimal models. We further investigate the sensitivity of our results to the assumptions of the hypergeometric model, using individual-based simulations.

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow - a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co-occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment-length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host-associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

Post-glacial habitat release and incipient speciation in the genus Delphinus

The role of ecological and changing environmental factors in the radiation of species diversity is a fundamental question in evolutionary biology. Of particular interest is the potential for these factors to determine the boundary between what we would consider differentiation among populations and incipient speciation. Dolphins in the genus Delphinus provide a useful test case, exhibiting morphological variation in beak length, coloration and body size across their wide geographic distribution, and in particular among coastal and more pelagic habitats. Two species have been proposed, D. delphis and D. capensis, but morphologically similar allopatric populations are not monophyletic, indicating that the mostly coastal 'long-beaked' D. capensis form is not a single globally distributed species. However, the sympatric populations in the Eastern North Pacific currently designated as these two species are both morphologically and genetically differentiated. Here we use microsatellite DNA and mitochondrial DNA markers to investigate the evolutionary mechanisms that led to this incipient speciation event. We used coalescent and assignment methods to investigate the timing and extent of reproductive isolation. Our data indicate that although there is some level of on-going gene flow, the putative species found in the Eastern North Pacific are reciprocally monophyletic. The timing of isolation appears to be associated with regional changes in paleoceanographic conditions within the Holocene timeframe.

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site-associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.

Predicting the Occurrence of Cave-Inhabiting Fauna Based on Features of the Earth Surface Environment

One of the most challenging fauna to study in situ is the obligate cave fauna because of the difficulty of sampling. Cave-limited species display patchy and restricted distributions, but it is often unclear whether the observed distribution is a sampling artifact or a true restriction in range. Further, the drivers of the distribution could be local environmental conditions, such as cave humidity, or they could be associated with surface features that are surrogates for cave conditions. If surface features can be used to predict the distribution of important cave taxa, then conservation management is more easily obtained. We examined the hypothesis that the presence of major faunal groups of cave obligate species could be predicted based on features of the earth surface. Georeferenced records of cave obligate amphipods, crayfish, fish, isopods, beetles, millipedes, pseudoscorpions, spiders, and springtails within the area of Appalachian Landscape Conservation Cooperative in the eastern United States (Illinois to Virginia and New York to Alabama) were assigned to 20 x 20 km grid cells. Habitat suitability for these faunal groups was modeled using logistic regression with twenty predictor variables within each grid cell, such as percent karst, soil features, temperature, precipitation, and elevation. Models successfully predicted the presence of a group greater than 65% of the time (mean = 88%) for the presence of single grid cell endemics, and for all faunal groups except pseudoscorpions. The most common predictor variables were latitude, percent karst, and the standard deviation of the Topographic Position Index (TPI), a measure of landscape rugosity within each grid cell. The overall success of these models points to a number of important connections between the surface and cave environments, and some of these, especially soil features and topographic variability, suggest new research directions. These models should prove to be useful tools in predicting the presence of species in understudied areas.

Divergence history of the Carpathian and smooth newts modelled in space and time

Information about demographic history is essential for the understanding of the processes of divergence and speciation. Patterns of genetic variation within and between closely related species provide insights into the history of their interactions. Here, we investigated historical demography and genetic exchange between the Carpathian (Lissotriton montandoni, Lm) and smooth (L. vulgaris, Lv) newts. We combine an extensive geographical sampling and multilocus nuclear sequence data with the approximate Bayesian computation framework to test alternative scenarios of divergence and reconstruct the temporal and spatial pattern of gene flow between species. A model of recent (last glacial period) interspecific gene flow was favoured over alternative models. Thus, despite the relatively old divergence (4-6 mya) and presumably long periods of isolation, the species have retained the ability to exchange genes. Nevertheless, the low migration rates (ca. 10(-6) per gene copy per generation) are consistent with strong reproductive isolation between the species. Models allowing demographic changes were favoured, suggesting that the effective population sizes of both species at least doubled as divergence reaching the current ca. 0.2 million in Lm and 1 million in Lv. We found asymmetry in rates of interspecific gene flow between Lm and one evolutionary lineage of Lv. We suggest that intraspecific polymorphism for hybrid incompatibilities segregating within Lv could explain this pattern and propose further tests to distinguish between alternative explanations. Our study highlights the importance of incorporating intraspecific genetic structure into the models investigating the history of divergence.

Incipient speciation with gene flow on a continental island: Species delimitation of the Hainan Hwamei (Leucodioptron canorum owstoni, Passeriformes, Aves)

Because of their isolation, continental islands (e.g., Madagascar) are often thought of as ideal systems to study allopatric speciation. However, many such islands have been connected intermittently to their neighboring continent during recent periods of glaciation, which may cause frequent contact between the diverging populations on the island and continent. As a result, the speciation processes on continental islands may not meet the prerequisites for strictly allopatric speciation. We used multiple lines of evidence to re-evaluate the taxonomic status of the Hainan Hwamei (Leucodioptron canorum owstoni), which is endemic to Hainan, the largest continental island in the South China Sea. Our analysis of mitochondrial DNA and twelve nuclear loci suggests that the Hainan Hwamei can be regarded as an independent species (L. owstoni); the morphological traits of the Hainan Hwamei also showed significant divergence from those of their mainland sister taxon, the Chinese Hwamei (L. canorum). We also inferred the divergence history of the Hainan and Chinese Hwamei to see whether their divergence was consistent with a strictly allopatric model. Our results suggest that the two Hwameis split only 0.2 million years ago with limited asymmetrical post-divergence gene flow. This implies that the Hainan Hwamei is an incipient species and that speciation occurred through ecologically divergent selection and/or assortative mating rather than a strictly allopatric process.

Phylogenetic species delimitation for crayfishes of the genus Pacifastacus

Molecular genetic approaches are playing an increasing role in conservation science by identifying biodiversity that may not be evident by morphology-based taxonomy and systematics. So-called cryptic species are particularly prevalent in freshwater environments, where isolation of dispersal-limited species, such as crayfishes, within dendritic river networks often gives rise to high intra- and inter-specific genetic divergence. We apply here a multi-gene molecular approach to investigate relationships among extant species of the crayfish genus Pacifastacus, representing the first comprehensive phylogenetic study of this taxonomic group. Importantly, Pacifastacus includes both the widely invasive signal crayfish Pacifastacus leniusculus, as well as several species of conservation concern like the Shasta crayfish Pacifastacus fortis. Our analysis used 83 individuals sampled across the four extant Pacifastacus species (omitting the extinct Pacifastacus nigrescens), representing the known taxonomic diversity and geographic distributions within this genus as comprehensively as possible. We reconstructed phylogenetic trees from mitochondrial (16S, COI) and nuclear genes (GAPDH), both separately and using a combined or concatenated dataset, and performed several species delimitation analyses (PTP, ABGD, GMYC) on the COI phylogeny to propose Primary Species Hypotheses (PSHs) within the genus. All phylogenies recovered the genus Pacifastacus as monophyletic, within which we identified a range of six to 21 PSHs; more abundant PSHs delimitations from GMYC and ABGD were always nested within PSHs delimited by the more conservative PTP method. Pacifastacus leniusculus included the majority of PSHs and was not monophyletic relative to the other Pacifastacus species considered. Several of these highly distinct P. leniusculus PSHs likely require urgent conservation attention. Our results identify research needs and conservation priorities for Pacifastacus crayfishes in western North America, and may inform better understanding and management of P. leniusculus in regions where it is invasive, such as Europe and Japan.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.