Ecological speciation

Allopatric and Sympatric Drivers of Speciation in Alviniconcha Hydrothermal Vent Snails

Despite significant advances in our understanding of speciation in the marine environment, the mechanisms underlying evolutionary diversification in deep-sea habitats remain poorly investigated. Here, we used multigene molecular clocks and population genetic inferences to examine processes that led to the emergence of the six extant lineages of Alviniconcha snails, a key taxon inhabiting deep-sea hydrothermal vents in the Indo-Pacific Ocean. We show that both allopatric divergence through historical vicariance and ecological isolation due to niche segregation contributed to speciation in this genus. The split between the two major Alviniconcha clades (separating A. boucheti and A. marisindica from A. kojimai, A. hessleri, and A. strummeri) probably resulted from tectonic processes leading to geographic separation, whereas the splits between co-occurring species might have been influenced by ecological factors, such as the availability of specific chemosynthetic symbionts. Phylogenetic origin of the sixth species, Alviniconcha adamantis, remains uncertain, although its sister position to other extant Alviniconcha lineages indicates a possible ancestral relationship. This study lays a foundation for future genomic studies aimed at deciphering the roles of local adaptation, reproductive biology, and host–symbiont compatibility in speciation of these vent-restricted snails.

Niche specificity influences gene flow across fine-scale habitat mosaics in Succulent Karoo plants

While the tempo of diversification in biodiversity hotspots has received much attention, the spatial scale of diversification has often been overlooked. Addressing this deficiency requires understanding the drivers of population divergence and the spatial scales at which they operate in species-rich clades and ecosystems. South Africa's Succulent Karoo (SK) hotspot provides an excellent system for such research, being both compact (ca. 110,000 km² ) and home to spectacular in-situ radiations, such as the ruschioid Aizoaceae. Here we use GBS to document genetic structure in two co-occurring ruschioid species, at both coarse (>10 km) and fine (<500 m) spatial scales. Where Ruschia burtoniae shows strong between-population genetic differentiation and no gene flow, Conophytum calculus shows weak differentiation, with high levels of admixture suggesting recent or ongoing gene flow. Community analysis and transplant experiments reveal that R. burtoniae occupies a narrow, low-pH edaphic niche, and at scales of a few hundred metres, areas of elevated genetic turnover correspond to patches of edaphically unsuitable habitat. In contrast, C. calculus occupies a broader niche and exhibits isolation-by-distance without a habitat effect. We suggest that edaphic specialisation, coupled with highly restricted seed and pollen dispersal in heterogeneous landscapes, has played a major role in driving rapid diversification at small spatial scales in this system. However, the contrasting patterns in our study species show that these factors do not influence all organisms uniformly, being strongly modulated by lineage-specific traits that influence both the spatial scale of gene flow and habitat specificity.

Environmental heterogeneity determines the ecological processes that govern bacterial metacommunity assembly in a floodplain river system

How diversity is structured has been a central goal of microbial ecology. In freshwater ecosystems, selection has been found to be the main driver shaping bacterial communities. However, its relative importance compared with other processes (dispersal, drift, diversification) may depend on spatial heterogeneity and the dispersal rates within a metacommunity. Still, a decrease in the role of selection is expected with increasing dispersal homogenization. Here, we investigate the main ecological processes modulating bacterial assembly in contrasting scenarios of environmental heterogeneity. We carried out a spatiotemporal survey in the floodplain system of the Paraná River. The bacterioplankton metacommunity was studied using both statistical inferences based on phylogenetic and taxa turnover as well as co-occurrence networks. We found that selection was the main process determining community assembly even at both extremes of environmental heterogeneity and homogeneity, challenging the general view that the strength of selection is weakened due to dispersal homogenization. The ecological processes acting on the community also determined the connectedness of bacterial networks associations. Heterogeneous selection promoted more interconnected networks increasing β-diversity. Finally, spatiotemporal heterogeneity was an important factor determining the number and identity of the most highly connected taxa in the system. Integrating all these empirical evidences, we propose a new conceptual model that elucidates how the environmental heterogeneity determines the action of the ecological processes shaping the bacterial metacommunity.

Ecological and morphological determinants of evolutionary diversification in Darwin's finches and their relatives

Darwin's finches are a classic example of adaptive radiation, a process by which multiple ecologically distinct species rapidly evolve from a single ancestor. Such evolutionary diversification is typically explained by adaptation to new ecological opportunities. However, the ecological diversification of Darwin's finches following their dispersal to Galápagos was not matched on the same archipelago by other lineages of colonizing land birds, which diversified very little in terms of both species number and morphology. To better understand the causes underlying the extraordinary variation in Darwin's finches, we analyze the evolutionary dynamics of speciation and trait diversification in Thraupidae, including Coerebinae (Darwin's finches and relatives) and, their closely related clade, Sporophilinae. For all traits, we observe an early pulse of speciation and morphological diversification followed by prolonged periods of slower steady-state rates of change. The primary exception is the apparent recent increase in diversification rate in Darwin's finches coupled with highly variable beak morphology, a potential key factor explaining this adaptive radiation. Our observations illustrate how the exploitation of ecological opportunity by contrasting means can produce clades with similarly high diversification rate yet strikingly different degrees of ecological and morphological differentiation.

Moniliophthora perniciosa, the mushroom causing witches' broom disease of cacao: Insights into its taxonomy, ecology and host range in Brazil

Witches' broom caused by Moniliophthora perniciosa is the main disease of cacao (Theobroma cacao) in Brazil. The fungus is known to occur on other host families and these populations have been addressed in the literature as biotypes: C (Malvaceae); H (Malpighiaceae); L (Bignoniaceae) and S (Solanaceae). No complete elucidation of the phylogenetic relationships of isolates obtained from this disparate host range appears in the literature. One member of H (ex Heteropterys acutifolia) has been described as a distinct species. But should other biotypes be also recognized as distinct taxa? In the present study, a survey yielding 24 isolates of M. perniciosa from ten hosts and covering a wide range of geographic regions in Brazil was undertaken. These isolates were compared with those from T. cacao using three DNA regions for the phylogenetic analyses: ITS, LSU and RPB1. Morphology was also examined. All isolates in this study were found to belong to M. perniciosa, including the population from H. acutifolia, formerly treated as Moniliophthora brasiliensis but reduced here to a synonym of M. perniciosa. This species ranged from pathogenic to a previously unreported occurrence as a non-pathogenic endophyte in the Atlantic rainforest tree Allophylus edulis (Sapindaceae). M. perniciosa was recorded on a range of solanaceous hosts (16 species) over a wide variety of ecosystems. The ecological and evolutionary significance of these novel findings are discussed.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Genomic landscape of reproductive isolation in Lucania killifish: The role of sex loci and salinity

Adaptation to different environments can directly and indirectly generate reproductive isolation between species. Bluefin killifish (Lucania goodei) and rainwater killifish (L. parva) are sister species that have diverged across a salinity gradient and are reproductively isolated by habitat, behavioural, extrinsic and intrinsic post-zygotic isolation. We asked if salinity adaptation contributes indirectly to other forms of reproductive isolation via linked selection and hypothesized that low recombination regions, such as sex chromosomes or chromosomal rearrangements, might facilitate this process. We conducted QTL mapping in backcrosses between L. parva and L. goodei to explore the genetic architecture of salinity tolerance, behavioural isolation and intrinsic isolation. We mapped traits relative to a chromosome that has undergone a centric fusion in L. parva (relative to L. goodei). We found that the sex locus appears to be male determining (XX-XY), was located on the fused chromosome and was implicated in intrinsic isolation. QTL associated with salinity tolerance were spread across the genome and did not overly co-localize with regions associated with behavioural or intrinsic isolation. This preliminary analysis of the genetic architecture of reproductive isolation between Lucania species does not support the hypothesis that divergent natural selection for salinity tolerance led to behavioural and intrinsic isolation as a by-product. Combined with previous studies in this system, our work suggests that adaptation as a function of salinity contributes to habitat isolation and that reinforcement may have contributed to the evolution of behavioural isolation instead, possibly facilitated by linkage between behavioural isolation and intrinsic isolation loci on the fused chromosome.

Divergent Selection and Primary Gene Flow Shape Incipient Speciation of a Riparian Tree on Hawaii Island

A long-standing goal of evolutionary biology is to understand the mechanisms underlying the formation of species. Of particular interest is whether or not speciation can occur in the presence of gene flow and without a period of physical isolation. Here, we investigated this process within Hawaiian Metrosideros, a hypervariable and highly dispersible woody species complex that dominates the Hawaiian Islands in continuous stands. Specifically, we investigated the origin of Metrosideros polymorpha var. newellii (newellii), a riparian ecotype endemic to Hawaii Island that is purportedly derived from the archipelago-wide M. polymorpha var. glaberrima (glaberrima). Disruptive selection across a sharp forest-riparian ecotone contributes to the isolation of these varieties and is a likely driver of newellii's origin. We examined genome-wide variation of 42 trees from Hawaii Island and older islands. Results revealed a split between glaberrima and newellii within the past 0.3-1.2 My. Admixture was extensive between lineages within Hawaii Island and between islands, but introgression from populations on older islands (i.e., secondary gene flow) did not appear to contribute to the emergence of newellii. In contrast, recurrent gene flow (i.e., primary gene flow) between glaberrima and newellii contributed to the formation of genomic islands of elevated absolute and relative divergence. These regions were enriched for genes with regulatory functions as well as for signals of positive selection, especially in newellii, consistent with divergent selection underlying their formation. In sum, our results support riparian newellii as a rare case of incipient ecological speciation with primary gene flow in trees.

The signature of competition in ecomorphological traits across the avian radiation

Competition for shared resources represents a fundamental driver of biological diversity. However, the tempo and mode of phenotypic evolution in deep-time has been predominantly investigated using trait evolutionary models which assume that lineages evolve independently from each other. Consequently, the role of species interactions in driving macroevolutionary dynamics remains poorly understood. Here, we quantify the prevalence for signatures of competition between related species in the evolution of ecomorphological traits across the bird radiation. We find that mechanistic trait models accounting for the effect of species interactions on phenotypic divergence provide the best fit for the data on at least one trait axis in 27 out of 59 clades ranging between 21 and 195 species. Where it occurs, the signature of competition generally coincides with positive species diversity-dependence, driven by the accumulation of lineages with similar ecologies, and we find scarce evidence for trait-dependent or negative diversity-dependent phenotypic evolution. Overall, our results suggest that the footprint of interspecific competition is often eroded in long-term patterns of phenotypic diversification, and that other selection pressures may predominantly shape ecomorphological diversity among extant species at macroevolutionary scales.

Model-Based Inference of Punctuated Molecular Evolution

In standard models of molecular evolution, DNA sequences evolve through asynchronous substitutions according to Poisson processes with a constant rate (called the molecular clock) or a rate that can vary (relaxed clock). However, DNA sequences can also undergo episodes of fast divergence that will appear as synchronous substitutions affecting several sites simultaneously at the macroevolutionary timescale. Here, we develop a model, which we call the Relaxed Clock with Spikes model, combining basal, clock-like molecular substitutions with episodes of fast divergence called spikes arising at speciation events. Given a multiple sequence alignment and its time-calibrated species phylogeny, our model is able to detect speciation events (including hidden ones) cooccurring with spike events and to estimate the probability and amplitude of these spikes on the phylogeny. We identify the conditions under which spikes can be distinguished from the natural variance of the clock-like component of molecular substitutions and from variations of the clock. We apply the method to genes underlying snake venom proteins and identify several spikes at gene-specific locations in the phylogeny. This work should pave the way for analyses relying on whole genomes to inform on modes of species diversification.

Appendicular skeletal morphology of North American Martes reflect independent modes of evolution in conjunction with Pleistocene glacial cycles

Pleistocene glacial cycles are thought to have driven ecological niche shifts, including novel niche formation. North American pine martens, Martes americana and M. caurina, are exemplar taxa thought to have diverged molecularly and morphologically during Pleistocene glaciation. Previous research found correlations between Martes limb morphology with biome and climate, suggesting that appendicular evolution may have occurred via adaptation to selective pressures imposed by novel and shifting habitats. Such variation can also be achieved through non-adaptive means such as genetic drift. Here, we evaluate whether regional genetic differences reflect limb morphology differences among populations of M. americana and M. caurina by analyzing evolutionary tempo and mode of six limb elements. Our comparative phylogenetic models indicate that genetic structure predicts limb shape better than size. Marten limb size has low phylogenetic signal, and the best supported model of evolution is punctuational (kappa), with morphological and genetic divergence occurring simultaneously. Disparity through time analysis suggests that the tempo of limb evolution in Martes tracks Pleistocene glacial cycles, such that limb size may be responding to shifting climates rather than population genetic structure. Contrarily, we find that limb shape is strongly tied to genetic relationships, with high phylogenetic signal and a lambda mode of evolution. Overall, this pattern of limb size and shape variation may be the result of geographic isolation during Pleistocene glacial advance, while declines in disparity suggest hybridization during interglacial periods. Future inclusion of extinct populations of Martes, which were more morphologically and ecologically diverse, may further clarify Martes phenotypic evolution.

Host-symbiont coevolution, cryptic structure, and bleaching susceptibility, in a coral species complex (Scleractinia; Poritidae)

The 'species' is a key concept for conservation and evolutionary biology, yet the lines between population and species-level variation are often blurred, especially for corals. The 'Porites lobata species complex' consists of branching and mounding corals that form reefs across the Pacific. We used reduced representation meta-genomic sequencing to examine genetic relationships within this species complex and to identify candidate loci associated with colony morphology, cryptic genetic structure, and apparent bleaching susceptibility. We compared existing Porites data with bleached and unbleached colonies of the branching coral P. compressa collected in Kāne'ohe Bay Hawai'i during the 2015 coral bleaching event. Loci that mapped to coral, symbiont, and microbial references revealed genetic structure consistent with recent host-symbiont co-evolution. Cryptic genetic clades were resolved that previous work has associated with distance from shore, but no genetic structure was associated with bleaching. We identified many candidate loci associated with morphospecies, including candidate host and symbiont loci with fixed differences between branching and mounding corals. We also found many loci associated with cryptic genetic structure, yet relatively few loci associated with bleaching. Recent host-symbiont co-evolution and rapid diversification suggests that variation and therefore the capacity of these corals to adapt may be underappreciated.

A non-exhaustive survey revealed possible genetic similarity in mitochondrial adaptive evolution of marine fish species in the northwestern Pacific

Mitochondrial coding genes involved in the oxidative phosphorylation pathway play vitally important roles in energy production and thermal adaptation. Investigating the underlying molecular mechanism of mitochondrial adaptive evolution is crucial for understanding biodiversity and ecological radiation. In this study, we collated population genetic studies of marine fish species in the northwestern Pacific based on mitochondrial cytochrome b gene sequences, to investigate whether similar patterns could be detected in mitochondrial adaptive evolution. After filtering, nine studies containing eight marine fish species (Ammodytespersonatus, Boleophthalmuspectinirostris, Larimichthyspolyactis, Mugilcephalus, Pampusargenteus, Platycephalus sp.1, Sebastiscusmarmoratus, and Trachidermusfasciatus) belonging to eight different families were retained. Multiple codon-based approaches were used to identify potential sites under selection in each species. By comparison, our results showed that the posterior part of the mitochondrial cytochrome b gene (particularly codon 372 and its neighboring sites) seemed to be involved in the adaptive evolution process, suggesting potential genetic similarity among distantly related species. We also summarized four types of adaptive patterns in the reviewed species, and suggest that the level of genetic differentiation and mitochondrial adaptive evolution might be correlated. Further studies are needed to confirm such relationship by detecting RNA-level evidence and investigating more species and samples.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

Dynamic shifts in social network structure and composition within a breeding hybrid population

Mating behaviour and the timing of reproduction can inhibit genetic exchange between closely related species; however, these reproductive barriers are challenging to measure within natural populations. Social network analysis provides promising tools for studying the social context of hybridization, and the exchange of genetic variation, more generally. We test how social networks within a hybrid population of California Callipepla californica and Gambel's quail Callipepla gambelii change over discrete periods of a breeding season. We assess patterns of phenotypic and genotypic assortment, and ask whether altered associations between individuals (association rewiring), or changes to the composition of the population (individual turnover) drive network dynamics. We use genetic data to test whether social associations and relatedness between individuals correlate with patterns of parentage within the hybrid population. To achieve these aims, we combine RFID association data, phenotypic data and genomic measures with social network analyses. We adopt methods from the ecological network literature to quantify shifts in network structure and to partition changes into those due to individual turnover and association rewiring. We integrate genomic data into networks as node-level attributes (ancestry) and edges (relatedness, parentage) to test links between social and parentage networks. We show that rewiring of associations between individuals that persist across network periods, rather than individual turnover, drives the majority of the changes in network structure throughout the breeding season, and that the traits involved in phenotypic/genotypic assortment were highly dynamic over time. Social networks were randomly assorted based on genetic ancestry, suggesting weak behavioural reproductive isolation within this hybrid population. Finally, we show that the strength of associations within the social network, but not levels of genetic relatedness, predicts patterns of parentage. Social networks play an important role in population processes such as the transmission of disease and information, yet there has been less focus on how networks influence the exchange of genetic variation. By integrating analyses of social structure, phenotypic assortment and reproductive outcomes within a hybrid zone, we demonstrate the utility of social networks for analysing links between social context and gene flow within wild populations.

Maladaptive migration behaviour in hybrids links to predator-mediated ecological selection

Different migratory species have evolved distinct migratory characteristics that improve fitness in their particular ecological niches. However, when such species hybridize, migratory traits from parental species can combine maladaptively and cause hybrids to fall between parental fitness peaks, with potential consequences for hybrid viability and species integrity.

Here, we take advantage of a natural cross‐breeding incident to study migratory behaviour in naturally occurring hybrids as well as in their parental species and explore links between migratory traits and predation risk.

To achieve this, we used electronic tags and passive telemetry to record detailed individual migration patterns (timing and number of migratory trips) in two common freshwater fish species, roach Rutilus rutilus, common bream Abramis brama as well as their hybrids. Next, we scanned for tags regurgitated by a key avian predator (great cormorant Phalacrocorax carbo) at nearby roosting sites, allowing us to directly link migratory behaviour to predation risk in the wild.

We found that hybrid individuals showed a higher number of short, multi‐trip movements between lake and stream habitats as compared to both parental species. The mean date of first lake departure differed between bream and roach by more than 10 days, while hybrids departed in two distinct peaks that overlapped with the parental species' averages. Moreover, the probability of cormorant predation increased with multi‐trip movement frequency across species and was higher for hybrids.

Our data provide novel insights into hybrid viability, with links to predator‐mediated ecological selection. Increased exposure to predators via maladaptive migratory behaviour reduces hybrid survival and can thereby reinforce species integrity.

Speciation patterns in the Forficula auricularia species complex: cryptic and not so cryptic taxa across the western Palaearctic region

Forficula auricularia (the European earwig) is possibly a complex of cryptic species. To test this hypothesis, we performed: (1) a phylogeographic study based on fragments of the mitochondrial COI and the nuclear ITS2 markers on a wide geographic sampling, (2) morphometric analyses of lineages present in Spain and (3) niche overlap analyses. We recovered five reciprocally monophyletic ancient phylogroups with unique historical patterns of distribution, climatic niches and diversification. External morphology was conserved and not correlated with speciation events, except in one case. Phylogenetic placement of the morphologically distinct taxon renders F. auricularia paraphyletic. Based on the congruence of the phylogenetic units defined by mtDNA and nuclear sequence data, we conclude that phylogroups have their own historical and future evolutionary trajectory and represent independent taxonomic units. Forficula auricularia is a complex of at least four species: the morphologically diagnosable Forficula aeolica González-Miguéns & García-París sp. nov., and the cryptic taxa: Forficula mediterranea González-Miguéns & García-París sp. nov., Forficula dentataFabricius, 1775stat. nov. and Forficula auriculariaLinnaeus, 1758s.s. We also provide new synonymy for F. dentata.

A generalized framework for AMOVA with multiple hierarchies and ploidies

The analysis of molecular variance (AMOVA) is a widely used statistical method in population genetics and molecular ecology. The classic framework of AMOVA only supports haploid and diploid data, in which the number of hierarchies ranges from two to four. In practice, natural populations can be classified into more hierarchies, and polyploidy is frequently observed in extant species. The ploidy level may even vary within the same species, and/or within the same individual. We generalized the framework of AMOVA such that it can be used for any number of hierarchies and any level of ploidy. Based on this framework, we present four methods to account for data that are multilocus genotypic and allelic phenotypic (with unknown allele dosage). We use simulated datasets and an empirical dataset to evaluate the performance of our framework. We make freely available our methods in a new software package, polygene, which is freely available at https://github.com/huangkang1987/polygene.

The effectiveness of pseudomagic traits in promoting divergence and enhancing local adaptation

"Magic traits," in which the same trait is both under divergent ecological selection and forms the basis of assortative mating, have been sought after due to their supposed unique ability to promote divergence with gene flow. Here, we ask how unique magic traits are, by exploring whether a tightly linked complex of a locus under divergent selection and a locus that acts as a mating cue can mimic a magic trait in its divergence. We find that these "pseudomagic traits" can be very effective in promoting divergence; with tight linkage they are essentially as effective as a magic trait and with loose linkage, and even no linkage, divergence can still be enhanced. Distinguishing between magic and pseudomagic traits in empirical cases may thus not be important when inferring their role in divergence. The ability of divergence in the mating trait to drive divergence in the ecological trait by lowering the effective migration rate, which occurs somewhat even without linkage, is particularly striking; magic traits are typically considered to have the other direction of causality. Our results thus suggest that divergence in a mating trait can at least modestly increase local adaption by allowing more ecological divergence, particularly with tighter linkage.

Morphological and trophic divergence of lake and stream minnows (Phoxinus phoxinus)

Phenotypic divergence in response to divergent natural selection between environments is a common phenomenon in species of freshwater fishes. Intraspecific differentiation is often pronounced between individuals inhabiting lakes versus stream habitats. The different hydrodynamic regimes in the contrasting habitats may promote a variation of body shape, but this could be intertwined with morphological adaptations to a specific foraging mode. Herein, I studied the divergence pattern of the European minnow (Phoxinus phoxinus), a common freshwater fish that has received little attention despite its large distribution. In many Scandinavian mountain lakes, European minnows are considered as being invasive and were found to pose threats to the native fish populations due to resource competition. Minnows were recently found to show phenotypic adaptations in lake versus stream habitats, but the question remained if this divergence pattern is related to differences in resource use. I therefore studied the patterns of minnow divergence in morphology (i.e., using geometric morphometrics) and trophic niches (i.e., using stomach content analyses) in the lake Ånnsjön and its tributaries to link the changes in body morphology to the feeding on specific resources. Lake minnows showed a strong reliance on benthic Cladocera and a more streamlined body shape with a more upward facing snout, whereas stream minnows fed on macroinvertebrates (larvae and adults) to a higher degree and had a deeper body with a snout that was pointed down. Correlations showed a significant relationship of the proportion of macroinvertebrates in the gut and morphological features present in the stream minnows. The results of this study highlight the habitat-specific divergence pattern in morphology and resource use in this ubiquitous freshwater fish. Consequently, interspecific interactions of invasive minnows and the native fish population could differ in the respective food webs and resource competition could target different native fish species in the contrasting habitats.

A Genomewide Scan for Genetic Structure and Demographic History of Two Closely Related Species, Rhododendron dauricum and R. mucronulatum (Rhododendron, Ericaceae)

Understanding the processes of divergence and speciation is an important task for evolutionary research, and climate oscillations play a pivotal role. We estimated the genetic structure and demographic history of two closely related species of Rhododendron, R. dauricum, and R. mucronulatum, distributed in northeastern China using 664,406 single nucleotide polymorphic loci of specific-locus amplified fragment sequencing (SLAF-seq) and 4 chloroplast DNA (cpDNA) fragments, sampling 376 individuals from 39 populations of these two species across their geographic distributions. The geographical distribution of cpDNA haplotypes revealed that R. dauricum and R. mucronulatum have different spatial genetic structures and haplotype diversity. Analysis of molecular variance (AMOVA) results showed that these two species have significant genetic differentiation and that the phylogeny demonstrates that these two species clustered a monophyletic group based on SLAF data, respectively, but not in cpDNA data. The evidence of significant gene flow was also detected from R. mucronulatum to R. dauricum. A deep divergence between the two species was observed and occurred during the early Oligocene. The niche models showed that the two species have different demographic histories. Thus, our results imply that geography and climate changes played important roles in the evolutionary process of R. dauricum and R. mucronulatum, and although there was an interspecific gene flow, the divergence was maintained by natural selection.

The genomic view of diversification

The process of species diversification is traditionally summarized by a single tree, the species tree, whose reconstruction from molecular data is hindered by frequent conflicts between gene genealogies. Here, we argue that instead of seeing these conflicts as nuisances, we can exploit them to inform the diversification process itself. We adopt a gene-based view of diversification to model the ubiquitous presence of gene flow between diverging lineages, one of the most important processes explaining disagreements among gene trees. We propose a new framework for modelling the joint evolution of gene and species lineages relaxing the hierarchy between the species tree and gene trees inherent to the standard view, as embodied in a popular model known as the multispecies coalescent (MSC). We implement this framework in two alternative models called the gene-based diversification models (GBD): (a) GBD-forward following all evolving genomes through time and (b) GBD-backward based on coalescent theory. They feature four parameters tuning colonization, gene flow, genetic drift and genetic differentiation. We propose an inference method based on differences between gene trees. Applied to two empirical data sets prone to gene flow, we find better support for the GBD-backward model than for the MSC model. Along with the increasing awareness of the extent of gene flow, this work shows the importance of considering the richer signal contained in genomic histories, rather than in the mere species tree, to better apprehend the complex evolutionary history of species.

Repeated horizontal gene transfers triggered parallel evolution of magnetotaxis in two evolutionary divergent lineages of magnetotactic bacteria

Under the same selection pressures, two genetically divergent populations may evolve in parallel toward the same adaptive solutions. Here, we hypothesized that magnetotaxis (i.e., magnetically guided chemotaxis) represents a key adaptation to micro-oxic habitats in aquatic sediments and that its parallel evolution homogenized the phenotypes of two evolutionary divergent clusters of freshwater spirilla. All magnetotactic bacteria affiliated to the Magnetospirillum genus (Alphaproteobacteria class) biomineralize the same magnetic particle chains and share highly similar physiological and ultrastructural features. We looked for the processes that could have contributed at shaping such an evolutionary pattern by reconciling species and gene trees using newly sequenced genomes of Magnetospirillum related bacteria. We showed that repeated horizontal gene transfers and homologous recombination of entire operons contributed to the parallel evolution of magnetotaxis. We propose that such processes could represent a more parsimonious and rapid solution for adaptation compared with independent and repeated de novo mutations, especially in the case of traits as complex as magnetotaxis involving tens of interacting proteins. Besides strengthening the idea about the importance of such a function in micro-oxic habitats, these results reinforce previous observations in experimental evolution suggesting that gene flow could alleviate clonal interference and speed up adaptation under some circumstances.

Warning signals promote morphological diversification in fossorial uropeltid snakes (Squamata: Uropeltidae)

Many species possess warning colourations that signal unprofitability to predators. Warning colourations are also thought to provide prey with a ‘predator-free space’ and promote niche expansion. However, how such strategies release a species from environmental constraints and facilitate niche expansion is not clearly understood. Fossoriality in reptiles imposes several morphological limits on head and body size to facilitate burrowing underground, but many fossorial snakes live close to the surface and occasionally move above ground, exposing them to predators. In such cases, evolving antipredator defences that reduce predation on the surface could potentially relax the morphological constraints associated with fossoriality and promote morphological diversification. Fossorial uropeltid snakes possess varying degrees of conspicuous warning colourations that reduce avian predation when active above ground. We predicted that species with more conspicuous colourations will exhibit more robust body forms and show faster rates of morphological evolution because constraints imposed by fossoriality are relaxed. Using a comparative phylogenetic approach on the genus Uropeltis, we show that more conspicuous species tend to have more robust morphologies and have faster rates of head-shape evolution. Overall, we find that the evolution of warning colourations in Uropeltis can facilitate niche expansion by influencing rates of morphological diversification.

Genetic consequences of being a dwarf: do evolutionary changes in life-history traits influence gene flow patterns in populations of the world's smallest goldenrod?

BACKGROUND AND AIMS

Contrasting life-history traits can evolve through generations of dwarf plant ecotypes, yet such phenotypic changes often involve decreased plant size and reproductive allocation, which can configure seed dispersal patterns and, subsequently, population demography. Therefore, evolutionary transitions to dwarfism can represent good study systems to test the roles of life-history traits in population demography by comparing genetic structure between related but phenotypically divergent ecotypes.

METHODS

In this study, we examined an ecotypic taxon pair of the world's smallest goldenrod (stem height 2.6 cm) in alpine habitats and its closely related lowland taxon (30-40 cm) found on Yakushima Island, Japan. Genetic variation in chloroplast DNA sequences, nuclear microsatellites and genome-wide single-nucleotide polymorphisms were used to investigate 197 samples from 16 populations, to infer the population genetic demography and compare local genetic structure of the ecotypes.

KEY RESULTS

We found a pronounced level of genetic differentiation among alpine dwarf populations, which were much less geographically isolated than their lowland counterparts. In particular, several neighbouring dwarf populations (located ~500 m apart) harboured completely different sets of chloroplast haplotypes and nuclear genetic clusters. Demographic modelling revealed that the dwarf populations have not exchanged genes at significant levels after population divergence.

CONCLUSIONS

These lines of evidence suggest that substantial effects of genetic drift have operated on these dwarf populations. The low-growing stature and reduced fecundity (only 3.1 heads per plant) of the dwarf plants may have reduced gene flow and rare long-distance seed dispersal among habitat patches, although the effects of life-history traits require further evaluation using ecological approaches.

Effects of temperature changes on life parameters, oxidative stress, and antioxidant defense system in the monogonont marine rotifer Brachionus plicatilis

Global warming is a big concern for all organisms and many efforts have been made to reveal the potential effects of temperature elevation on aquatic organisms. However, limited studies on molecular mechanistic approaches on physiological effects due to temperature changes are available. Here, we investigated the effects of temperature changes on life parameters (e.g., population growth [total number of rotifers], and lifespan), oxidative stress levels and antioxidant activities (e.g., glutathione S-transferase [GST], catalase [CAT], superoxide dismutase [SOD]) with expression levels in the monogonont marine rotifer Brachionus plicatilis. The changes in temperatures led to significant reduction (P < 0.05) in lifespan, possibly due to significant decrease (P < 0.05) in antioxidant activities, reducing the potential to cope with significant elevation in the temperature-induced oxidative stress in B. plicatilis. To further assess the actual induction and clearance of reactive oxygen species (ROS), N-acetyl-L-cysteine was used to examine whether the temperature-induced oxidative stress could be successfully scavenged. Furthermore, expression patterns of the antioxidant-related genes (GSTs, SODs, and CATs) were down- or upregulated (P < 0.05) in response to different temperatures in B. plicatilis. Overall, these findings indicate that ROS-mediated oxidative stress led to cellular damage and antioxidant defense system, resulting in deleterious effects on life parameters in rotifer.

Discrepancies between genetic and ecological divergence patterns suggest a complex biogeographic history in a Neotropical genus

Phylogenetic patterns and the underlying speciation processes can be deduced from morphological, functional, and ecological patterns of species similarity and divergence. In some cases, though, species retain multiple similarities and remain almost indistinguishable; in other cases, evolutionary convergence can make such patterns misleading; very often in such cases, the "true" picture only emerges from carefully built molecular phylogenies, which may come with major surprises. In addition, closely related species may experience gene flow after divergence, thus potentially blurring species delimitation. By means of advanced inferential methods, we studied molecular divergence between species of the Virola genus (Myristicaceae): widespread Virola michelii and recently described, endemic V. kwatae, using widespread V. surinamensis as a more distantly related outgroup with different ecology and morphology-although with overlapping range. Contrary to expectations, we found that the latter, and not V. michelii, was sister to V. kwatae. Therefore, V. kwatae probably diverged from V. surinamensis through a recent morphological and ecological shift, which brought it close to distantly related V. michelii. Through the modeling of the divergence process, we inferred that gene flow between V. surinamensis and V. kwatae stopped soon after their divergence and resumed later, in a classical secondary contact event which did not erase their ecological and morphological differences. While we cannot exclude that initial divergence occurred in allopatry, current species distribution and the absence of geographical barriers make complete isolation during speciation unlikely. We tentatively conclude that (a) it is possible that divergence occurred in allopatry/parapatry and (b) secondary contact did not suppress divergence.

Lotka-Volterra approximations for evolutionary trait-substitution processes

A set of axioms is formulated characterizing ecologically plausible community dynamics. Using these axioms, it is proved that the transients following an invasion into a sufficiently stable equilibrium community by a mutant phenotype similar to one of the community's finitely many resident phenotypes can always be approximated by means of an appropriately chosen Lotka–Volterra model. To this end, the assumption is made that similar phenotypes in the community form clusters that are well-separated from each other, as is expected to be generally the case when evolution proceeds through small mutational steps. Each phenotypic cluster is represented by a single phenotype, which we call an approximate phenotype and assign the cluster’s total population density. We present our results in three steps. First, for a set of approximate phenotypes with arbitrary equilibrium population densities before the invasion, the Lotka–Volterra approximation is proved to apply if the changes of the population densities of these phenotypes are sufficiently small during the transient following the invasion. Second, quantitative conditions for such small changes of population densities are derived as a relationship between within-cluster differences and the leading eigenvalue of the community’s Jacobian matrix evaluated at the equilibrium population densities before the invasion. Third, to demonstrate the utility of our results, the ‘invasion implies substitution’ result for monomorphic populations is extended to arbitrarily polymorphic populations consisting of well-recognizable and -separated clusters.

Phylogenetic niche conservatism and plant diversification in South American subtropical grasslands along multiple climatic dimensions

Phylogenetic niche conservatism can be investigated at multiple scales on an explicit geographical context. Haplotype-based comparative analyses of lineages occupying the same region, and thus subjected to similar environmental factors, allow decoupling shared evolutionary and ecological patterns, as well as multiple dimensions of adaptive diversification. Here we aimed to assess the role of environmental drivers on diversification of subtropical grassland, based on haplotypic diversity of two plant genera. We sampled two closely related and co-distributed grassland plant genera, Petunia and Calibrachoa, across their entire distribution area. Eigenvectors extracted from pairwise distances based on chloroplast DNA haplotypes were used to fit Phylogenetic Signal-Representation (PSR) curves to estimate evolutionary patterns in 19 bioclimatic variables and altitude. The PSR curves showed that altitude, precipitation, and temperature variables changed at different rates with haplotype differentiation. Altitude and temperature traits evolved under conditions closer to a neutral dynamics, whereas precipitation traits differentiated following more complex models. Our results indicated that the diversification in the two genera was more limited by precipitation conditions. Based on these novel findings, we suggest that future studies should test the possible impact of precipitation variables on the process of ecological differentiation in these genera.

Morphological variation associated with trophic niche expansion within a lake population of a benthic fish

Ecological theory suggests that generalist species should have traits with multiple adaptive peaks. Consequently, in heterogeneous environments such adaptive landscapes may lead to phenotypic divergence that becomes fixed in populations via reproductive isolation, thus driving speciation. However, contrary to this expectation, the process of ecological diversification in wild populations is not always associated with obvious trait divergence and reproductive isolation due to some ecological and geographic constrains. To examine the ecological conditions that promote (or inhibit) divergence is quite important to improve our understanding of the underlying mechanisms. Here we examine how the patterns of trait variation (divergence/non-divergence) are determined in relation to ecological niche expansion and gene flow using a benthic fish, Pseudogobio esocinus, in the Lake Biwa system, Japan. The fish exhibited various patterns of morphological variation in mouth parts among populations. Lake fish tended to have a smaller mouth compared with river fish and also showed remarkable individual variations within some local samples. Lake fish utilized chironomid larvae as the primary prey, as in riverine fish. But, fish with smaller and narrower mouths utilized significantly higher proportions of amphipods (a novel prey unique to the lake) as their secondary prey. Microsatellite analysis detected no genetic structuring in the Lake Biwa catchment, suggesting no reproductive separation among eco-morphologically divergent individuals. Our results exemplify population niche expansion associated with continuous eco-morphological variation without divergence, and provide insights into the role of non-discrete diversification for thriving in heterogeneous environments.

Whole Genome Scan Reveals Molecular Signatures of Divergence and Selection Related to Important Traits in Durum Wheat Germplasm

The first breeding program in the world for durum wheat was conceived in Italy in the early 1900s. Over the decades, pressure exerted by natural and artificial selection could have progressively reduced the genetic diversity of the durum wheat germplasm. In the present study, a large panel of Italian durum wheat accessions that includes landraces, old and modern cultivars was subjected to genotyping using the Illumina iSelect 15K wheat SNP array. The aim was to assess the impact that selection has in shaping Italian durum wheat genetic diversity and to exploit the patterns of genetic diversity between populations to identify molecular signatures of divergence and selection. Relatively small differences in genetic diversity have been observed among accessions, which have been selected and cultivated in Italy over the past 150 years. Indeed, directional selection combined with that operated by farmers/breeders resulted in the increase of linkage disequilibrium (LD) and in changes of the allelic frequencies in DNA regions that control important agronomic traits. Results from this study also show that major well-known genes and/or QTLs affecting plant height (RHT), earliness (VRN, PPD) and grain quality (GLU, PSY, PSD, LYC, PPO, LOX3) co-localized with outlier SNP loci. Interestingly, many of these SNPs fall in genomic regions where genes involved in nitrogen metabolism are. This finding highlights the key role these genes have played in the transition from landraces to modern cultivars. Finally, our study remarks on the need to fully exploit the genetic diversity of Italian landraces by intense pre-breeding activities aimed at introducing a new source of adaptability and resistance in the genetic background of modern cultivars, to contrast the effect of climate change. The list of divergent loci and loci under selection associated with useful agronomic traits represents an invaluable resource to detect new allelic variants for target genes and for guiding new genomic selection programs in durum wheat.

Temporally consistent species differences in parasite infection but no evidence for rapid parasite-mediated speciation in Lake Victoria cichlid fish

Parasites may have strong eco‐evolutionary interactions with their hosts. Consequently, they may contribute to host diversification. The radiation of cichlid fish in Lake Victoria provides a good model to study the role of parasites in the early stages of speciation. We investigated patterns of macroparasite infection in a community of 17 sympatric cichlids from a recent radiation and 2 older species from 2 nonradiating lineages, to explore the opportunity for parasite‐mediated speciation. Host species had different parasite infection profiles, which were only partially explained by ecological factors (diet, water depth). This may indicate that differences in infection are not simply the result of differences in exposure, but that hosts evolved species‐specific resistance, consistent with parasite‐mediated divergent selection. Infection was similar between sampling years, indicating that the direction of parasite‐mediated selection is stable through time. We morphologically identified 6 Cichlidogyrus species, a gill parasite that is considered a good candidate for driving parasite‐mediated speciation, because it is host species‐specific and has radiated elsewhere in Africa. Species composition of Cichlidogyrus infection was similar among the most closely related host species (members of the Lake Victoria radiation), but two more distantly related species (belonging to nonradiating sister lineages) showed distinct infection profiles. This is inconsistent with a role for Cichlidogyrus in the early stages of divergence. To conclude, we find significant interspecific variation in parasite infection profiles, which is temporally consistent. We found no evidence that Cichlidogyrus‐mediated selection contributes to the early stages of speciation. Instead, our findings indicate that species differences in infection accumulate after speciation.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

Rapid adaptation of the Irish potato famine pathogen Phytophthora infestans to changing temperature

Temperature plays a multidimensional role in host-pathogen interactions. As an important element of climate change, elevated world temperature resulting from global warming presents new challenges to sustainable disease management. Knowledge of pathogen adaptation to global warming is needed to predict future disease epidemiology and formulate mitigating strategies. In this study, 21 Phytophthora infestans isolates originating from seven thermal environments were acclimated for 200 days under stepwise increase or decrease of experimental temperatures and evolutionary responses of the isolates to the thermal changes were evaluated. We found temperature acclimation significantly increased the fitness and genetic adaptation of P. infestans isolates at both low and high temperatures. Low-temperature acclimation enforced the countergradient adaptation of the pathogen to its past selection and enhanced the positive association between the pathogen's intrinsic growth rate and aggressiveness. At high temperatures, we found that pathogen growth collapsed near the maximum temperature for growth, suggesting a thermal niche boundary may exist in the evolutionary adaptation of P. infestans. These results indicate that pathogens can quickly adapt to temperature shifts in global warming. If this is associated with environmental conditions favoring pathogen spread, it will threaten future food security and human health and require the establishment of mitigating actions.

Genetic structure is stronger across human-impacted habitats than among islands in the coral Porites lobata

We examined genetic structure in the lobe coral Porites lobata among pairs of highly variable and high-stress nearshore sites and adjacent less variable and less impacted offshore sites on the islands of Oahu and Maui, Hawaii. Using an analysis of molecular variance framework, we tested whether populations were more structured by geographic distance or environmental extremes. The genetic patterns we observed followed isolation by environment, where nearshore and adjacent offshore populations showed significant genetic structure at both locations (AMOVA F ST = 0.04∼0.19, P < 0.001), but no significant isolation by distance between islands. Strikingly, corals from the two nearshore sites with higher levels of environmental stressors on different islands over 100 km apart with similar environmentally stressful conditions were genetically closer (FST = 0.0, P = 0.73) than those within a single location less than 2 km apart (FST = 0.04∼0.08, P < 0.01). In contrast, a third site with a less impacted nearshore site (i.e., less pronounced environmental gradient) showed no significant structure from the offshore comparison. Our results show much stronger support for environment than distance separating these populations. Our finding suggests that ecological boundaries from human impacts may play a role in forming genetic structure in the coastal environment, and that genetic divergence in the absence of geographical barriers to gene flow might be explained by selective pressure across contrasting habitats.

Evolutionary history of two cryptic species of northern African jerboas

BACKGROUND

Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence.

RESULTS

The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization.

CONCLUSIONS

Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

Biogeographic barriers, Pleistocene refugia, and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species-specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.