Model-based demographic inference of introgression history in European whitefish species pairs'

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Phylogenetic analysis of museum specimens of houting Coregonus oxyrinchus shows the need for a revision of its extinct status

According to the IUCN Red List the anadromous houting Coregonus oxyrinchus is categorized as 'extinct'. However, this extinct status might be incorrect because taxonomic difference between C. oxyrinchus and the closely related C. lavaretus is based on a disputable morphological comparison. Also, phylogenetic studies on mtDNA only focused on recent obtained coregonids. We are the first to perform a mtDNA analysis on both historic and recent specimens, including the syntype specimen which was used for species description by Linnaeus originally. Two primer pairs for mitochondrial CytB and ND3 were used to extract sequences for phylogenetic analysis. Sequences from 14 out of 21 C. oxyrinchus museum specimens were successfully obtained and compared with sequences from recent obtained C. lavaretus. The sequences were combined with GenBank data from a previous phylogenetic study on houting to create a phylogenetic tree and two minimum spanning haplotype networks. Results show that C. oxyrinchus and C. lavaretus form a clade with limited genetic variation. Low bootstrap values also show weak support for geographical patterns in distribution of mitochondrial haplotypes. Statistical analysis of the haplotype networks also shows that historic and recent specimens are similar species. Our results suggest that C. oxyrinchus is a junior synonym of C. lavaretus. A definitive taxonomic revision could not be made because only CytB sequencing was successful for the syntype specimen. We discuss taxonomic consequences and the species-specific focus in nature conservation. We propose a shift in nature conservation to a more functional approach based on traits rather than species.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Genetic Causes and Consequences of Sympatric Morph Divergence in Salmonidae: A Search for Mechanisms

Repeatedly and recently evolved sympatric morphs exhibiting consistent phenotypic differences provide natural experimental replicates of speciation. Because such morphs are observed frequently in Salmonidae, this clade provides a rare opportunity to uncover the genomic mechanisms underpinning speciation. Such insight is also critical for conserving salmonid diversity, the loss of which could have significant ecological and economic consequences. Our review suggests that genetic differentiation among sympatric morphs is largely nonparallel apart from a few key genes that may be critical for consistently driving morph differentiation. We discuss alternative levels of parallelism likely underlying consistent morph differentiation and identify several factors that may temper this incipient speciation between sympatric morphs, including glacial history and contemporary selective pressures. Our synthesis demonstrates that salmonids are useful for studying speciation and poses additional research questions to be answered by future study of this family. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Highly Replicated Evolution of Parapatric Ecotypes

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

Complex and divergent histories gave rise to genome-wide divergence patterns amongst European whitefish (Coregonus lavaretus)

Pleistocene glaciations dramatically affected species distribution in regions that were impacted by ice cover and subsequent postglacial range expansion impacted contemporary biodiversity in complex ways. The European whitefish, Coregonus lavaretus, is a widely distributed salmonid fish species on mainland Europe, but in Britain it has only seven native populations, all of which are found on the western extremes of the island. The origins and colonization routes of the species into Britain are unknown but likely contributed to contemporary genetic patterns and regional uniqueness. Here, we used up to 25,751 genome‐wide polymorphic loci to reconstruct the history and to discern the demographic and evolutionary forces underpinning divergence between British populations. Overall, we found lower genetic diversity in Scottish populations but high differentiation (F_ST = 0.433–0.712) from the English/Welsh and other European populations. Differentiation was elevated genome‐wide rather than in particular genomic regions. Demographic modelling supported a postglacial colonization into western Scotland from northern refugia and a separate colonization route for the English/Welsh populations from southern refugia, with these two groups having been separated for more than ca. 50 Ky. We found cyto‐nuclear discordance at a European scale, with the Scottish populations clustering closely with Baltic population in the mtDNA analysis but not in the nuclear data, and with the Norwegian and Alpine populations displaying the same mtDNA haplotype but being distantly related in the nuclear tree. These findings suggest that neutral processes, primarily drift and regionally distinct pre‐glacial evolutionary histories, are important drivers of genomic divergence in British populations of European whitefish. This sheds new light on the establishment of the native British freshwater fauna after the last glacial maximum.

Genetic relationships between sympatric and allopatric Coregonus ciscoes in North and Central Europe

BACKGROUND

Sympatric speciation along ecological gradients has been studied repeatedly, in particular in freshwater fishes. Rapid post-glacial ecological divergence has resulted in numerous endemic species or ecologically distinct populations in lakes of the temperate zones. Here, we focus on the Baltic cisco (Coregonus albula) complex, to study the genetic similarity among two pairs of sympatric autumn- and spring-spawning populations from post-glacial German Lakes Stechlin and Breiter Luzin. For comparison, we included a similar pair of sympatric populations from the Swedish Lake Fegen. We wanted to explore potential genetic similarities between the three sympatric cisco population pairs in the three lakes, to evaluate whether the pairs may have emerged independently in the three lakes, or whether two different species may have colonized all three lakes independently. Furthermore, we considered allopatric C. albula populations from three Polish, three Finnish, and four Swedish locations, added one Siberian population of the sister species C. sardinella and a Swedish C. maraena (whitefish) population. By genotyping nine microsatellite markers in 655 individuals from these 18 populations, we wanted to elucidate how strongly the cisco populations differ across a larger geographical area within Europe. Finally, we compared the genetic differences between the spring- and autumn-spawning populations of ciscoes in the two German lakes to infer the potentially deteriorating effect of strong anthropogenic pressure on the lakes.

RESULTS

Dendrogram, Principal Coordinate Analysis and admixture analysis all indicated strong correspondence between population differentiation and geographical location for most cisco populations in Europe, including the Siberian population of C. sardinella. However, populations from some Swedish lakes deviated from this general pattern, by showing a distinct genetic structure. We found evidence for independent evolution of the three sympatric population pairs, because the populations co-occurring in the same lake were always most closely related. However, genetic differentiation was weak in the two German population pairs, but strong in the Swedish Lake Fegen, indicating that the weak differentiation in the German pairs reported earlier has eroded further.

CONCLUSIONS

Our results suggest that the genetic differentiation at neutral genetic markers among populations of the Baltic cisco complex has evolved (and is maintained) by random genetic drift in isolated populations. However, earlier studies on the Swedish populations combining mitochondrial DNA and microsatellite data indicate that also post-glacial immigration from separate glacial refugia has shaped the present genetic population structure. The low neutral differentiation of the German sympatric pairs in contrast to the Swedish pair suggests that recent anthropogenic effects on the lakes in Germany may put the endemic spring-spawners at risk to extinction.

Rapid adaptation through genomic and epigenomic responses following translocations in an endangered salmonid

Identifying the molecular mechanisms facilitating adaptation to new environments is a key question in evolutionary biology, especially in the face of current rapid and human-induced changes. Translocations have become an important tool for species conservation, but the attendant small population sizes and new ecological pressures might affect phenotypic and genotypic variation and trajectories dramatically and in unknown ways. In Scotland, the European whitefish (Coregonus lavaretus) is native to only two lakes and vulnerable to extirpation. Six new refuge populations were established over the last 30 years as a conservation measure. In this study, we examined whether there is a predictable ecological and evolutionary response of these fishes to translocation. We found eco-morphological differences, as functional traits relating to body shape differed between source and refuge populations. Dual isotopic analyses suggested some ecological release, with the diets in refuge populations being more diverse than in source populations. Analyses of up to 9117 genome-mapped SNPs showed that refuge populations had reduced genetic diversity and elevated inbreeding and relatedness relative to source populations, though genomic differentiation was low (F ST = 0.002-0.030). We identified 14 genomic SNPs that showed shared signals of a selective response to translocations, including some located near or within genes involved in the immune system, nervous system and hepatic functions. Analysis of up to 120,897 epigenomic loci identified a component of consistent differential methylation between source and refuge populations. We found that epigenomic variation and genomic variation were associated with morphological variation, but we were not able to infer an effect of population age because the patterns were also linked with the methodology of the translocations. These results show that conservation-driven translocations affect evolutionary potential by impacting eco-morphological, genomic and epigenomic components of diversity, shedding light on acclimation and adaptation process in these contexts.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Biases in Demographic Modeling Affect Our Understanding of Recent Divergence

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (N_e) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in N_e as well as the effects of linked selection and barrier loci. Failure to account for changes in N_e resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for N_e changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in N_e on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling.

Speciation-by-depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second-generation sequence data. Here, we revisit the case of two depth-segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coral Eunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth-segregated populations (Agaricia fragilis and Pocillopora damicornis) suggest divergence leading to depth-segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence-with-gene flow may account for much of the differentiation that separates closely related, depth-segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.