Rapid niche expansion by selection on functional genomic variation after ecosystem recovery

Sexual selection and speciation in the Anthropocene

Anthropogenic change threatens global biodiversity by causing severe ecological disturbance and extinction. Here, we consider the effects of anthropogenic change on one process that generates biodiversity. Sexual selection (a potent evolutionary force and driver of speciation) is highly sensitive to the environment and, thus, vulnerable to anthropogenic ecological change. Anthropogenic alterations to sexual display and mate preference can make it harder to distinguish between conspecific and heterospecific mates or can weaken divergence via sexual selection, leading to higher rates of hybridization and biodiversity loss. Occasionally, anthropogenically altered sexual selection can abet diversification, but this appears less likely than biodiversity loss. In our rapidly changing world, a full understanding of sexual selection and speciation requires a global change perspective.

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Anthropogenic Change and the Process of Speciation

Anthropogenic impacts on the environment alter speciation processes by affecting both geographical contexts and selection patterns on a worldwide scale. Here we review evidence of these effects. We find that human activities often generate spatial isolation between populations and thereby promote genetic divergence but also frequently cause sudden secondary contact and hybridization between diverging lineages. Human-caused environmental changes produce new ecological niches, altering selection in diverse ways that can drive diversification; but changes also often remove niches and cause extirpations. Human impacts that alter selection regimes are widespread and strong in magnitude, ranging from local changes in biotic and abiotic conditions to direct harvesting to global climate change. Altered selection, and evolutionary responses to it, impacts early-stage divergence of lineages, but does not necessarily lead toward speciation and persistence of separate species. Altogether, humans both promote and hinder speciation, although new species would form very slowly relative to anthropogenic hybridization, which can be nearly instantaneous. Speculating about the future of speciation, we highlight two key conclusions: (1) Humans will have a large influence on extinction and "despeciation" dynamics in the short term and on early-stage lineage divergence, and thus potentially speciation in the longer term, and (2) long-term monitoring combined with easily dated anthropogenic changes will improve our understanding of the processes of speciation. We can use this knowledge to preserve and restore ecosystems in ways that promote (re-)diversification, increasing future opportunities of speciation and enhancing biodiversity.

Cryptic persistence and loss of local endemism in Lake Constance charr subject to anthropogenic disturbance

In the welcome circumstance that species believed extinct are rediscovered, it is often the case that biological knowledge acquired before the presumed extinction is limited. Efforts to address these knowledge gaps, in particular to assess the taxonomic integrity and conservation status of such species, can be hampered by a lack of genetic data and scarcity of samples in museum collections. Here, we present a proof-of-concept case study based on a multidisciplinary data evaluation approach to tackle such problems. The approach was developed after the rediscovery, 40 years after its presumed extinction, of the enigmatic Lake Constance deep-water charr Salvelinus profundus. Targeted surveys led to the capture of further species and additional sympatric normal charr, Salvelinus cf. umbla. Since the lake had been subject to massive stocking in the past, an evaluation of the genetic integrity of both extant forms was called for in order to assess possible introgression. A two-step genomic approach was developed based on restriction site associated DNA (RAD). Diagnostic population genomic (single nucleotide polymorphism [SNP]) data were harvested from contemporary samples and used for RNA bait design to perform target capture in DNA libraries of archival scale material, enabling a comparison between extant and historic samples. Furthermore, life history traits and morphological data for both extant forms were gathered and compared with historical data from the past 60-120 years. While extant deep-water charr matched historical deep-water specimens in body shape, gill raker count, and growth rates, significant differences were discovered between historical and extant normal charr. These resulted were supported by genomic analyses of contemporary samples, revealing the two extant forms to be highly divergent. The results of population assignment tests suggest that the endemic deep-water charr persisted in Lake Constance during the eutrophic phase, but not one of the historical genomic samples could be assigned to the extant normal charr taxon. Stocking with non-endemic charr seems to be the most likely reason for these changes. This proof-of-concept study presents a multidisciplinary data evaluation approach that simultaneously tests population genomic integrity and addresses some of the conservation issues arising from rediscovery of a species characterized by limited data availability.

A taxonomic revision of ten whitefish species from the lakes Lucerne, Sarnen, Sempach and Zug, Switzerland, with descriptions of seven new species (Teleostei, Coregonidae)

The taxonomy of the endemic whitefish of the lakes of the Reuss River system (Lucerne, Sarnen, Zug) and Lake Sempach, Switzerland, is reviewed and revised. Lake Lucerne harbours five species. Coregonus intermundiasp. nov. and C. suspensussp. nov., are described. Coregonus nobilis Haack, 1882, C. suidteri Fatio, 1885, and C. zugensis Nüsslin, 1882, are redescribed. Genetic studies have shown that C. suidteri and C. zugensis are composed of several distinct species endemic to different lakes. The names C. suidteri and C. zugensis are restricted to the species of lakes Sempach and Zug, respectively. The whitefish populations previously referred to as C. suidteri and C. zugensis from Lake Lucerne are described as C. litoralissp. nov. and C. muellerisp. nov., respectively. Furthermore, the whitefish from Lake Zug that were previously referred to as C. suidteri are described as C. supersumsp. nov. A holotype is designated for C. supersum that was previously one of two syntypes of C. zugensis. The other syntype is retained for C. zugensis. Coregonus obliterussp. nov. is described from Lake Zug, and C. obliterus and C. zugensis from Lake Zug are extinct. Finally, we describe C. sarnensissp. nov. from lakes Sarnen and Alpnach. Coregonus suidteri from Lake Sempach shows strong signals of introgression from deliberately translocated non-native whitefish species, which questions if the extant population still carries a genetic legacy from the original species and thus may need to be considered extinct. Coregonus suspensus is genetically partially of allochthonous origin, closely related to the radiation of Lake Constance. It is therefore compared to all known and described species of Lake Constance: C. wartmanni Bloch, 1784, C. macrophthalmus Nüsslin, 1882, C. arenicolus Kottelat,1997, and C. gutturosus Gmelin, 1818.

The role of vital dietary biomolecules in eco-evo-devo dynamics

The physiological dependence of animals on dietary intake of vitamins, amino acids, and fatty acids is ubiquitous. Sharp differences in the availability of these vital dietary biomolecules among different resources mean that consumers must adopt a range of strategies to meet their physiological needs. We review the emerging work on omega-3 long-chain polyunsaturated fatty acids, focusing predominantly on predator-prey interactions, to illustrate that trade-off between capacities to consume resources rich in vital biomolecules and internal synthesis capacity drives differences in phenotype and fitness of consumers. This can then feedback to impact ecosystem functioning. We outline how focus on vital dietary biomolecules in eco-eco-devo dynamics can improve our understanding of anthropogenic changes across multiple levels of biological organization.

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.

Genomic variation from an extinct species is retained in the extant radiation following speciation reversal

Ecosystem degradation and biodiversity loss are major global challenges. When reproductive isolation between species is contingent on the interaction of intrinsic lineage traits with features of the environment, environmental change can weaken reproductive isolation and result in extinction through hybridization. By this process called speciation reversal, extinct species can leave traces in genomes of extant species through introgressive hybridization. Using historical and contemporary samples, we sequenced all four species of an Alpine whitefish radiation before and after anthropogenic lake eutrophication and the associated loss of one species through speciation reversal. Despite the extinction of this taxon, substantial fractions of its genome, including regions shaped by positive selection before eutrophication, persist within surviving species as a consequence of introgressive hybridization during eutrophication. Given the prevalence of environmental change, studying speciation reversal and its genomic consequences provides fundamental insights into evolutionary processes and informs biodiversity conservation.

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.

Alternative splicing and gene expression play contrasting roles in the parallel phenotypic evolution of a salmonid fish

Understanding the contribution of different molecular processes to evolution and development is crucial for identifying the mechanisms of adaptation. Here, we used RNA-sequencing data to test the importance of alternative splicing and differential gene expression in a case of parallel adaptive evolution, the replicated postglacial divergence of the salmonid fish Arctic charr (Salvelinus alpinus) into sympatric benthic and pelagic ecotypes across multiple independent lakes. We found that genes differentially spliced between ecotypes were mostly not differentially expressed (<6% overlap) and were involved in different biological processes. Differentially spliced genes were primarily enriched for muscle development and functioning, while differentially expressed genes were involved in metabolism, immunity and growth. Furthermore, alternative splicing and gene expression were mostly controlled by independent cis-regulatory quantitative trait loci (<3.4% overlap). Cis-regulatory regions were associated with the parallel divergence in splicing (16.5% of intron clusters) and expression (6.7%-10.1% of differentially expressed genes), indicating shared regulatory variation across ecotype pairs. Contrary to theoretical expectation, we found that differentially spliced genes tended to be highly central in regulatory networks ("hub genes") and were annotated to significantly more gene ontology terms compared to nondifferentially spliced genes, consistent with a higher level of pleiotropy. Together, our results suggest that the concerted regulation of alternative splicing and differential gene expression through different regulatory regions leads to the divergence of complementary processes important for local adaptation. This provides novel insights into the importance of contrasting but putatively complementary molecular processes in rapid parallel adaptive evolution.

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.

Continent-wide genomic signatures of adaptation to urbanisation in a songbird across Europe

Urbanisation is increasing worldwide, and there is now ample evidence of phenotypic changes in wild organisms in response to this novel environment. Yet, the genetic changes and genomic architecture underlying these adaptations are poorly understood. Here, we genotype 192 great tits (Parus major) from nine European cities, each paired with an adjacent rural site, to address this major knowledge gap in our understanding of wildlife urban adaptation. We find that a combination of polygenic allele frequency shifts and recurrent selective sweeps are associated with the adaptation of great tits to urban environments. While haplotypes under selection are rarely shared across urban populations, selective sweeps occur within the same genes, mostly linked to neural function and development. Collectively, we show that urban adaptation in a widespread songbird occurs through unique and shared selective sweeps in a core-set of behaviour-linked genes.

A taxonomic revision of the whitefish of lakes Brienz and Thun, Switzerland, with descriptions of four new species (Teleostei, Coregonidae)

The alpha taxonomy of the endemic whitefish of lakes Brienz and Thun, Switzerland, is revised. We evaluate the status of seven known species: Coregonus steinmanni sp. nov., Coregonus profundus sp. nov. and Coregonus acrinasus sp. nov. are endemic to Lake Thun; Coregonus Brienzii sp. nov. is endemic to Lake Brienz; and C. alpinus, C. albellus, and C. fatioi from lakes Brienz and Thun are redescribed. One of these species, C. alpinus, is revised, since the lectotype for this species is incongruent with the species description given by Kottelat (1997) and Kottelat and Freyhof (2007). The name C. alpinus is thus retained for the lectotype designated by Kottelat (1997) and a new description of this taxon provided. For the species otherwise described by Kottelat (1997) and Kottelat and Freyhof (2007) as C. alpinus the new name C. profundus is designated. Coregonus acrinasus is genetically partially of allochthonous origin, closely related to the radiation of Lake Constance, and we therefore compare it to all recognized species of Lake Constance, C. wartmanni, C. macrophthalmus, C. arenicolus, and C. gutturosus.

Intraspecific variation and structuring of phenotype in a lake-dwelling species are driven by lake size and elevation

The fragmented, heterogeneous and relatively depauperate ecosystems of recently glaciated lakes present contrasting ecological opportunities for resident fish. Across a species, local adaptation may induce diverse and distinct phenotypic responses to various selection pressures. We tested for intraspecific phenotypic structuring by population in a common native lake-dwelling fish species across a medium-scale geographic region with considerable variation in lake types. We investigated potential lake-characteristic drivers of trophic morphology. Using geometric morphometric techniques, we quantified the head shapes of 759 adult brown trout (Salmo trutta L.) from 28 lakes and reservoirs across Scotland. Multivariate statistical analyses showed that almost all populations differed from one another. Trout from larger and deeper lakes had deeper, but shorter heads, and smaller eyes. Higher elevation lakes were associated with fish with shorter heads and jaws. Within-population shape variation also differed by population, and was positively correlated with lake surface area and depth. Trout within reservoirs differed subtly from those in natural lakes, having larger eyes, shorter jaws and greater variability. This study documents an extraordinary morphological variation between and within populations of brown trout, and demonstrates the role of the extrinsic environment in driving phenotypic structuring over a medium-scale and varied geographic area.

Isotopic niche size of Coregonus artedi (sensu lato) Increases in the presence of Mysis diluviana, expanded habitat use and phenotypic diversity

Post-glacial colonization of lakes in Algonquin Park, Ontario, Canada resulted in food webs with cisco (Coregonus artedi sensu lato) and either Mysis diluviana or Chaoborus spp. as the dominant diel migrator. Mysis as prey, its diel movements and benthic occupancy, are hypothesized to be key elements of ecological opportunity for cisco diversity in the Laurentian Great Lakes. If correct, the hypothesis strongly implies that lakes with Mysis would have greater trophic niche size and drive greater adaptive radiation of cisco forms relative to lakes without Mysis. The dichotomy in diel migrator in Algonquin Park lakes was an opportunity to assess the isotopic niche size of cisco (δ15N and δ13C) and determine if niche size expands with Mysis presence. We found the presence of Mysis is necessary to expand isotopic niche size in our study lakes. The use of habitats not typically associated with the ancestral form of cisco (e.g., benthic habitats) and phenotypic diversity (blackfin and cisco) also continue to expand niche size in Mysis-based food webs. Partial ecological speciation based on a large niche space appears to be present in one lake (Cauchon Lake) where use of alternative habitats is the only real difference in cisco. The presence of blackfin expands niche space in Cedar and Radiant Lakes. This was not matched in Hogan Lake where niche space was relatively smaller with similar forms. Possible reasons for this discrepancy may be related to the asymmetric basin of Hogan Lake and whether the two forms overlap during cool and cold-water periods of the annual temperature cycle. By comparing trophic niche size among lakes with and without Mysis, we conclude that Mysis provides a key ecological opportunity for cisco diversity in our study lakes and likely more widely.

Geographic hierarchical population genetic structuring in British European whitefish (Coregonus lavaretus) and its implications for conservation

The European whitefish Coregonus lavaretus complex represents one of the most diverse radiations within salmonids, with extreme morphological and genetic differentiation across its range. Such variation has led to the assignment of many populations to separate species. In Great Britain, the seven native populations of C. lavaretus (two in Scotland, four in England, one in Wales) were previously classified into three species, and recent taxonomic revision resurrected the previous nomenclature. Here we used a dataset of 15 microsatellites to: (1) investigate the genetic diversity of British populations, (2) assess the level of population structure and the relationships between British populations. Genetic diversity was highest in Welsh (HO = 0.50, AR = 5.29), intermediate in English (HO = 0.41–0.50, AR = 2.83–3.88), and lowest in Scottish populations (HO = 0.28–0.35, AR = 2.56–3.04). Population structure analyses indicated high genetic differentiation (global FST = 0.388) between all populations but for the two Scottish populations (FST = 0.063) and two English populations (FST = 0.038). Principal component analysis and molecular ANOVA revealed separation between Scottish, English, and Welsh populations, with the Scottish populations being the most diverged. We argue that the data presented here are not sufficient to support a separation of the British European whitefish populations into three separate species, but support the delineation of different ESUs for these populations.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Hybridization drives genetic erosion in sympatric desert fishes of western North America

Many species have evolved or currently coexist in sympatry due to differential adaptation in a heterogeneous environment. However, anthropogenic habitat modifications can either disrupt reproductive barriers or obscure environmental conditions which underlie fitness gradients. In this study, we evaluated the potential for an anthropogenically-mediated shift in reproductive boundaries that separate two historically sympatric fish species (Gila cypha and G. robusta) endemic to the Colorado River Basin using ddRAD sequencing of 368 individuals. We first examined the integrity of reproductive isolation while in sympatry and allopatry, then characterized hybrid ancestries using genealogical assignment tests. We tested for localized erosion of reproductive isolation by comparing site-wise genomic clines against global patterns and identified a breakdown in the drainage-wide pattern of selection against interspecific heterozygotes. This, in turn, allowed for the formation of a hybrid swarm in one tributary, and asymmetric introgression where species co-occur. We also detected a weak but significant relationship between genetic purity and degree of consumptive water removal, suggesting a role for anthropogenic habitat modifications in undermining species boundaries or expanding historically limited introgression. In addition, results from basin-wide genomic clines suggested that hybrids and parental forms are adaptively nonequivalent. If so, then a failure to manage for hybridization will exacerbate the long-term extinction risk in parental populations. These results reinforce the role of anthropogenic habitat modification in promoting interspecific introgression in sympatric species by relaxing divergent selection. This, in turn, underscores a broader role for hybridization in decreasing global biodiversity within rapidly deteriorating environments.