Genomic divergence in a ring species complex

Patterns of genetic, phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.

Genotyping-by-Sequencing in a Species Complex of Australian Hummock Grasses (Triodia): Methodological Insights and Phylogenetic Resolution

Next-generation sequencing is becoming increasingly accessible to researchers asking biosystematic questions, but current best practice in both choosing a specific approach and effectively analysing the resulting data set is still being explored. We present a case study for the use of genotyping-by-sequencing (GBS) to resolve relationships in a species complex of Australian arid and semi-arid grasses (Triodia R.Br.), highlighting our solutions to methodological challenges in the use of GBS data. We merged overlapping paired-end reads then optimised locus assembly in the program PyRAD to generate GBS data sets for phylogenetic and distance-based analyses. In addition to traditional concatenation analyses in RAxML, we also demonstrate the novel use of summary species tree analyses (taking gene trees as input) with GBS loci. We found that while species tree analyses were relatively robust to variation in PyRAD assembly parameters, our RAxML analyses resulted in well-supported but conflicting topologies under different assembly settings. Despite this conflict, multiple clades in the complex were consistently supported as distinct across analyses. Our GBS data assembly and analyses improve the resolution of taxa and phylogenetic relationships in the Triodia basedowii complex compared to our previous study based on Sanger sequencing of nuclear (ITS/ETS) and chloroplast (rps16-trnK spacer) markers. The genomic results also partly support previous evidence for hybridization between species in the complex. Our methodological insights for analysing GBS data will assist researchers using similar data to resolve phylogenetic relationships within species complexes.

Barriers to gene flow and ring species formation

Ring species are groups of organisms that dispersed along a ring-shaped region in such a way that the two ends of the population that meet after many generations are reproductively isolated. They provide a rare opportunity to understand the role of spatial structuring in speciation. Here, we simulate the evolution of ring species assuming that individuals become sexually isolated if the genetic distance between them is above a certain threshold. The model incorporates two forms of dispersal limitation: exogenous geographic barriers that limit the population range and endogenous barriers that result in genetic structuring within the population range. As expected, species' properties that reduce gene flow within the population range facilitate the evolution of reproductive isolation and ring species formation. However, if populations are confined to narrow ranges by geographic barriers, ring species formation increases when local mating is less spatially restricted. Ring species are most likely to form if a population expands while confined to a quasi-unidimensional range but preserving high mobility in the direction of the range expansion. These conditions are unlikely to be met or persist in real populations and may explain why ring species are rare.

Genomic evidence of demographic fluctuations and lack of genetic structure across flyways in a long distance migrant, the European turtle dove

BACKGROUND

Understanding how past climatic oscillations have affected organismic evolution will help predict the impact that current climate change has on living organisms. The European turtle dove, Streptopelia turtur, is a warm-temperature adapted species and a long distance migrant that uses multiple flyways to move between Europe and Africa. Despite being abundant, it is categorized as vulnerable because of a long-term demographic decline. We studied the demographic history and population genetic structure of the European turtle dove using genomic data and mitochondrial DNA sequences from individuals sampled across Europe, and performing paleoclimatic niche modelling simulations.

RESULTS

Overall our data suggest that this species is panmictic across Europe, and is not genetically structured across flyways. We found the genetic signatures of demographic fluctuations, inferring an effective population size (Ne) expansion that occurred between the late Pleistocene and early Holocene, followed by a decrease in the Ne that started between the mid Holocene and the present. Our niche modelling analyses suggest that the variations in the Ne are coincident with recent changes in the availability of suitable habitat.

CONCLUSIONS

We argue that the European turtle dove is prone to undergo demographic fluctuations, a trait that makes it sensitive to anthropogenic impacts, especially when its numbers are decreasing. Also, considering the lack of genetic structure, we suggest all populations across Europe are equally relevant for conservation.

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within-group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within-population selection-induced reductions in variation. An important implication of this 'sweep-before-differentiation' model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Multispecies Outcomes of Sympatric Speciation after Admixture with the Source Population in Two Radiations of Nicaraguan Crater Lake Cichlids

The formation of species in the absence of geographic barriers (i.e. sympatric speciation) remains one of the most controversial topics in evolutionary biology. While theoretical models have shown that this most extreme case of primary divergence-with-gene-flow is possible, only a handful of accepted empirical examples exist. And even for the most convincing examples uncertainties remain; complex histories of isolation and secondary contact can make species falsely appear to have originated by sympatric speciation. This alternative scenario is notoriously difficult to rule out. Midas cichlids inhabiting small and remote crater lakes in Nicaragua are traditionally considered to be one of the best examples of sympatric speciation and lend themselves to test the different evolutionary scenarios that could lead to apparent sympatric speciation since the system is relatively small and the source populations known. Here we reconstruct the evolutionary history of two small-scale radiations of Midas cichlids inhabiting crater lakes Apoyo and Xiloá through a comprehensive genomic data set. We find no signs of differential admixture of any of the sympatric species in the respective radiations. Together with coalescent simulations of different demographic models our results support a scenario of speciation that was initiated in sympatry and does not result from secondary contact of already partly diverged populations. Furthermore, several species seem to have diverged simultaneously, making Midas cichlids an empirical example of multispecies outcomes of sympatric speciation. Importantly, however, the demographic models strongly support an admixture event from the source population into both crater lakes shortly before the onset of the radiations within the lakes. This opens the possibility that the formation of reproductive barriers involved in sympatric speciation was facilitated by genetic variants that evolved in a period of isolation between the initial founding population and the secondary migrants that came from the same source population. Thus, the exact mechanisms by which these species arose might be different from what had been thought before.

A genomic perspective on hybridization and speciation

Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.

Mitochondrial phylogeny of the Eurasian/African reed warbler complex (Acrocephalus, Aves). Disagreement between morphological and molecular evidence and cryptic divergence: A case for resurrecting Calamoherpe ambigua Brehm 1857

A tree based on the mitochondrial cyt b gene for 278 samples from throughout the range of the Eurasian Reed Warbler Acrocephalus scirpaceus - African Reed Warbler A. baeticatus complex shows well-supported geographically structured divergence for eight distinct lineages. The phylogenetic structuring together with the clarification of priority, provided by sequence data from seven type specimens, suggests that both taxonomy and distribution boundaries are in need of revision. The Iberian and Moroccan populations form a well-supported clade, and we propose that these are treated as taxonomically distinct, under the name ambiguus (Brehm, 1857). We propose that the names scirpaceus, fuscus, avicenniae, ambiguus, minor, cinnamomeus, hallae and baeticatus are used for the well supported clades in the complex, which we recommend to treat as one polytypic species, A. scirpaceus, pending studies of gene flow and assortative mating in the contact zones.

Ring species as demonstrations of the continuum of species formation

In the mid-20th century, Ernst Mayr (1942) and Theodosius Dobzhansky (1958) championed the significance of 'circular overlaps' or 'ring species' as the perfect demonstration of the gradual nature of species formation. As an ancestral species expands its range, wrapping around a geographic barrier, derived taxa within the ring display interactions typical of populations, such as genetic and morphological intergradation, while overlapping taxa at the terminus of the ring behave largely as sympatric, reproductively isolated species. Are ring species extremely rare or are they just difficult to detect? What conditions favour their formation? Modelling studies have attempted to address these knowledge gaps by estimating the biological parameters that result in stable ring species (Martins et al. 2013), and determining the necessary topographic parameters of the barriers encircled (Monahan et al. 2012). However, any generalization is undermined by a major limitation: only a handful of ring species are known to exist in nature. In addition, many of them have been broken into multiple species presumed to be evolving independently, usually obscuring the evolutionary dynamics that generate diversity. A paper in this issue of Molecular Ecology by Fuchs et al. (2015), focused on the entire genealogy of a bulbul (Alophoixus) species complex, offers key insights into the evolutionary processes underlying diversification of this Indo-Malayan bird. Their findings fulfil most of the criteria that can be expected for ring species (Fig. ): an ancestor has colonized the mainland from Sundaland, expanded along the forested habitat wrapping around Thailand's lowlands, adjacent taxa intergrade around the ring distribution, and terminal taxa overlap at the ring closure. Although it remains unclear whether ring divergence has resulted in restrictive gene flow relative to that observed around the ring, their results suggest that circular overlaps might be more common in nature than currently recognized in the literature. Most importantly, this work shows that the continuum of species formation that Mayr and Dobzhansky praised in circular overlaps is found in biological systems currently described as 'rings of species', in addition to the idealized 'ring species'.

A Spatial Framework for Understanding Population Structure and Admixture

Geographic patterns of genetic variation within modern populations, produced by complex histories of migration, can be difficult to infer and visually summarize. A general consequence of geographically limited dispersal is that samples from nearby locations tend to be more closely related than samples from distant locations, and so genetic covariance often recapitulates geographic proximity. We use genome-wide polymorphism data to build "geogenetic maps," which, when applied to stationary populations, produces a map of the geographic positions of the populations, but with distances distorted to reflect historical rates of gene flow. In the underlying model, allele frequency covariance is a decreasing function of geogenetic distance, and nonlocal gene flow such as admixture can be identified as anomalously strong covariance over long distances. This admixture is explicitly co-estimated and depicted as arrows, from the source of admixture to the recipient, on the geogenetic map. We demonstrate the utility of this method on a circum-Tibetan sampling of the greenish warbler (Phylloscopus trochiloides), in which we find evidence for gene flow between the adjacent, terminal populations of the ring species. We also analyze a global sampling of human populations, for which we largely recover the geography of the sampling, with support for significant histories of admixture in many samples. This new tool for understanding and visualizing patterns of population structure is implemented in a Bayesian framework in the program SpaceMix.

Speciation Progress: A Case Study on the Bushcricket Poecilimon veluchianus

Different mechanisms such as selection or genetic drift permitted e.g. by geographical isolation can lead to differentiation of populations and could cause subsequent speciation. The two subspecies of Poecilimon veluchianus, a bushcricket endemic to central Greece, show a parapatric distribution and are partially reproductively isolated. Therefore, P. veluchianus is suitable to investigate an ongoing speciation process. We based our analysis on sequences of the internal transcribed spacer (ITS) and the mitochondrial control region (CR). The population genetic analysis based on the nuclear marker ITS revealed a barrier to gene flow within the range of Poecilimon veluchianus, which corresponds well to the described subspecies. In contrast to the results based on the nuclear ITS marker, the mitochondrial CR marker does not clearly support the separation into two subspecies with restricted gene flow and a clear contact zone. Furthermore, we could identify isolation by distance (IBD) as one important mechanism responsible for the observed genetic structure (based on the ITS marker). The population genetic analysis based on the nuclear marker ITS also suggests the existence of hybrids in the wild. Furthermore, the simultaneous lack of strong prezygotic barriers and the presence of postzygotic mating barriers, observed in previous laboratory experiments, suggest that a secondary contact after an allopatric phase is more likely than parapatric speciation.

Genetic responses to seasonal variation in altitudinal stress: whole-genome resequencing of great tit in eastern Himalayas

Species that undertake altitudinal migrations are exposed to a considerable seasonal variation in oxygen levels and temperature. How they cope with this was studied in a population of great tit (Parus major) that breeds at high elevations and winters at lower elevations in the eastern Himalayas. Comparison of population genomics of high altitudinal great tits and those living in lowlands revealed an accelerated genetic selection for carbohydrate energy metabolism (amino sugar, nucleotide sugar metabolism and insulin signaling pathways) and hypoxia response (PI3K-akt, mTOR and MAPK signaling pathways) in the high altitudinal population. The PI3K-akt, mTOR and MAPK pathways modulate the hypoxia-inducible factors, HIF-1α and VEGF protein expression thus indirectly regulate hypoxia induced angiogenesis, erythropoiesis and vasodilatation. The strategies observed in high altitudinal great tits differ from those described in a closely related species on the Tibetan Plateau, the sedentary ground tit (Parus humilis). This species has enhanced selection in lipid-specific metabolic pathways and hypoxia-inducible factor pathway (HIF-1). Comparative population genomics also revealed selection for larger body size in high altitudinal great tits.

Protein expression parallels thermal tolerance and ecologic changes in the diversification of a diving beetle species complex

Physiological changes associated with evolutionary and ecological processes such as diversification, range expansion or speciation are still incompletely understood, especially for non-model species. Here we study differences in protein expression in response to temperature in a western Mediterranean diving beetle species complex, using two-dimensional differential gel electrophoresis with one Moroccan and one Iberian population each of Agabus ramblae and Agabus brunneus. We identified proteins with significant expression differences after thermal treatments comparing them with a reference EST library generated from one of the species of the complex (A. ramblae). The colonisation during the Middle Pleistocene of the Iberian peninsula by A. ramblae, where maximum temperatures and seasonality are lower than in the ancestral north African range, was associated with changes in the response to 27 °C in proteins related to energy metabolism. The subsequent speciation of A. brunneus from within populations of Iberian A. ramblae was associated with changes in the expression of several stress-related proteins (mostly chaperons) when exposed to 4 °C. These changes are in agreement with the known tolerance to lower temperatures of A. brunneus, which occupies a larger geographical area with a wider range of climatic conditions. In both cases, protein expression changes paralleled the evolution of thermal tolerance and the climatic conditions experienced by the species. However, although the colonisation of the Iberian peninsula did not result in morphological change, the speciation process of A. brunneus within Iberia involved genetic isolation and substantial differences in male genitalia and body size and shape.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

Complex histories of repeated gene flow in Cameroon crater lake cichlids cast doubt on one of the clearest examples of sympatric speciation

One of the most celebrated examples of sympatric speciation in nature are monophyletic radiations of cichlid fishes endemic to Cameroon crater lakes. However, phylogenetic inference of monophyly may not detect complex colonization histories involving some allopatric isolation, such as double invasions obscured by genome-wide gene flow. Population genomic approaches are better suited to test hypotheses of sympatric speciation in these cases. Here, we use comprehensive sampling from all four sympatric crater lake cichlid radiations in Cameroon and outgroups across Africa combined with next-generation sequencing to genotype tens of thousands of SNPs. We find considerable evidence of gene flow between all four radiations and neighboring riverine populations after initial colonization. In a few cases, some sympatric species are more closely related to outgroups than others, consistent with secondary gene flow facilitating their speciation. Our results do not rule out sympatric speciation in Cameroon cichlids, but rather reveal a complex history of speciation with gene flow, including allopatric and sympatric phases, resulting in both reproductively isolated species and incipient species complexes. The best remaining non-cichlid examples of sympatric speciation all involve assortative mating within microhabitats. We speculate that this feature may be necessary to complete the process of sympatric speciation in nature.

Phylogenetic analysis of beak and feather disease virus across a host ring-species complex

Pathogens have been hypothesized to play a major role in host diversity and speciation. Susceptibility of hybrid hosts to pathogens is thought to be a common phenomenon that could promote host population divergence and subsequently speciation. However, few studies have tested for pathogen infection across animal hybrid zones while testing for codivergence of the pathogens in the hybridizing host complex. Over 8 y, we studied natural infection by a rapidly evolving single-strand DNA virus, beak and feather diseases virus (BFDV), which infects parrots, exploiting a host-ring species complex (Platycercus elegans) in Australia. We found that host subspecies and their hybrids varied strikingly in both BFDV prevalence and load: both hybrid and phenotypically intermediate subspecies had lower prevalence and load compared with parental subspecies, while controlling for host age, sex, longitude and latitude, as well as temporal effects. We sequenced viral isolates throughout the range, which revealed patterns of genomic variation analogous to Mayr's ring-species hypothesis, to our knowledge for the first time in any host-pathogen system. Viral phylogeny, geographic location, intraspecific host density, and parrot community diversity and composition did not explain the differences in BFDV prevalence or load between subpopulations. Overall, our analyses suggest that functional host responses to infection, or force of infection, differ between subspecies and hybrids. Our findings highlight the role of host hybridization and clines in altering host-pathogen interactions, dynamics that can have important implications for models of speciation with gene flow, and offer insights into how pathogens may adapt to diverging host populations.