Identifying designatable units for intraspecific conservation prioritization: a hierarchical approach applied to the lake whitefish species complex (Coregonus spp.)

Pelvic spine reduction affects diet but not gill raker morphology in two polymorphic brook stickleback (Culaea inconstans) populations

Pelvic spine polymorphism occurs in several species in the stickleback family (Gasterosteidae). Given the similar phenotypic polymorphisms in multiple stickleback species, we sought to determine the extent of parallelism in the ecological correlates of pelvic spine reduction. Based on a metabarcoding analysis of brook stickleback gut contents in two polymorphic populations, we found that significant diet differences were associated with pelvic spine reduction, but we found no clear or consistent trend supporting a tendency for benthic feeding in pelvic-reduced brook sticklebacks. These results contrast with those found in threespine sticklebacks where pelvic spine reduction is often associated with a benthic diet. Hence, we found non-parallel consequences of spine polymorphism across species. Furthermore, a difference in gill raker morphology has been frequently observed between ecomorphs with different diets in many fish species. However, we found no evidence of any difference in gill raker morphology associated with pelvic spine polymorphism in brook sticklebacks.

Multifaceted framework for defining conservation units: An example from Atlantic salmon (Salmo salar) in Canada

Conservation units represent important components of intraspecific diversity that can aid in prioritizing and protecting at-risk populations, while also safeguarding unique diversity that can contribute to species resilience. In Canada, identification and assessments of conservation units is done by the Committee on the Status of Endangered Wildlife in Canada (COSEWIC). COSEWIC can recognize conservation units below the species level (termed "designatable units"; DUs) if the unit has attributes that make it both discrete and evolutionarily significant. There are various ways in which a DU can meet criteria of discreteness and significance, and increasing access to "big data" is providing unprecedented information that can directly inform both criteria. Specifically, the incorporation of genomic data for an increasing number of non-model species is informing more COSEWIC assessments; thus, a repeatable, robust framework is needed for integrating these data into DU characterization. Here, we develop a framework that uses a multifaceted, weight of evidence approach to incorporate multiple data types, including genetic and genomic data, to inform COSEWIC DUs. We apply this framework to delineate DUs of Atlantic salmon (Salmo salar, L.), an economically, culturally, and ecologically significant species, that is also characterized by complex hierarchical population structure. Specifically, we focus on an in-depth example of how our approach was applied to a previously data limited region of northern Canada that was defined by a single large DU. Application of our framework with newly available genetic and genomic data led to subdividing this DU into three new DUs. Although our approach was developed to meet criteria of COSEWIC, it is widely applicable given similarities in the definitions of a conservation unit.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Population structure of lake whitefish (Coregonus clupeaformis) from the Mississippian lineage in North America

The lake whitefish ( Coregonus clupeaformis) is a commercially valuable freshwater species with a broad distribution in North America. Some phylogeographic work has been done on this species, but little is known about genetic population subdivision among populations of the widely dispersed Mississippian lineage. We used 3,173 single nucleotide polymorphisms in 508 lake whitefish from 22 different lakes to examine population structure across central Canada and the United States. Bayesian clustering, ordination, and fixation indices identified population subdivision that largely reflected geographic distance and hydrological connectivity, with greater differentiation between lakes that are farther apart. Population subdivision was hierarchical, with greater differentiation between Canadian provinces and less differentiation based on river basins within provincial boundaries. Interestingly, isolation by distance alone was not sufficient to account for all of the observed genetic differentiation among populations. We conclude that important components of lake whitefish genetic diversity are present at different spatial scales, and that populations within the Mississippian lineage have differentiated widely across their range.

Split it up and see: using proxies to highlight divergent inter-populational performances in aquaculture standardised conditions

Background

Considering wild inter-populational phenotypic differentiation can facilitate domestication and subsequent production of new species. However, comparing all populations across a species range to identify those exhibiting suitable key traits for aquaculture (KTA; i.e. important for domestication and subsequent production) expressions is not feasible. Therefore, proxies highlighting inter-populational divergences in KTA are needed. The use of such proxies would allow to identify, prior to bioassays, the wild population pairs which are likely to present differentiations in KTA expressions in aquaculture conditions. Here, we assessed the relevance of three alternative proxies: (i) genetic distance, (ii) habitat divergence, and (iii) geographic/hydrologic distances. We performed this evaluation on seven allopatric populations of Perca fluviatilis for which divergences in KTA had already been shown.

Results

We showed differences in the correlation degree between the alternative proxy-based and KTA-based distance matrices, with the genetic proxy being correlated to the highest number of KTA. However, no proxy was correlated to all inter-populational divergences in KTA.

Conclusion

For future domestication trials, we suggest using a multi-proxy assessment along with a prioritisation strategy to identify population pairs which are of interest for further evaluation in bioassays.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01937-z.

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.

Reversing extinction trends: new uses of (old) herbarium specimens to accelerate conservation action on threatened species

Although often not collected specifically for the purposes of conservation, herbarium specimens offer sufficient information to reconstruct parameters that are needed to designate a species as 'at-risk' of extinction. While such designations should prompt quick and efficient legal action towards species recovery, such action often lags far behind and is mired in bureaucratic procedure. The increase in online digitization of natural history collections has now led to a surge in the number new studies on the uses of machine learning. These repositories of species occurrences are now equipped with advances that allow for the identification of rare species. The increase in attention devoted to estimating the scope and severity of the threats that lead to the decline of such species will increase our ability to mitigate these threats and reverse the declines, overcoming a current barrier to the recovery of many threatened plant species. Thus far, collected specimens have been used to fill gaps in systematics, range extent, and past genetic diversity. We find that they also offer material with which it is possible to foster species recovery, ecosystem restoration, and de-extinction, and these elements should be used in conjunction with machine learning and citizen science initiatives to mobilize as large a force as possible to counter current extinction trends.

Limited long-distance dispersal success in a Western European fairy shrimp evidenced by nuclear and mitochondrial lineage structuring

Anostraca are known by their ability for long-distance dispersal, but the existence in several species of deep, geographically structured mtDNA lineages suggests their populations are subjected to allopatric differentiation, isolation, and prevalence of local scale dispersion. Tanymastix stagnalis is one of the most widespread species of Anostraca and previous studies revealed an unclear geographical pattern of mtDNA genetic diversity. Here, we analyze populations from the Iberian and Italian Peninsulas, Central Europe, and Scandinavia, with the aim to characterize the patterns of genetic diversity in a spatio-temporal framework using mtDNA and nuclear markers to test gene flow among close populations. For these aims we built a time-calibrated phylogeny and carried out Bayesian phylogeographic analyses using a continuous diffusion model. Our results indicated that T. stagnalis presents a deeply structured genetic diversity, including 7 ancient lineages, some of them even predating the Pleistocene. The Iberian Peninsula harbors high diversity of lineages, with strong isolation and recent absence of gene flow between populations. Dispersal at local scale seems to be the prevailing dispersal mode of T. stagnalis, which exhibits a pattern of isolation-by-distance in the Iberian Peninsula. We remark the vulnerability of most of these lineages, given the limited known geographic distribution of some of them, and the high risk of losing important evolutionary potential for the species.

Genotyping-by-sequencing illuminates high levels of divergence among sympatric forms of coregonines in the Laurentian Great Lakes

Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the "coregonine problem," where recent postglacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site-associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially misidentified specimens. Population assignment rates for these forms using our RAD data were 93%-100% with the only mis-assignments arising from putative hybrids, an improvement from 62% to 77% using microsatellites. Estimates of pairwise differentiation (F ST: 0.045-0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long-term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well-informed management strategies and stress the importance of re-evaluating conservation units with genomic tools to ensure they accurately reflect species diversity.

Northern European Salmo trutta (L.) populations are genetically divergent across geographical regions and environmental gradients

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large-scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene-environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub-set of regions, suggesting that locus-specific variation was associated with local processes. A paired-population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection (pcadapt and bayescan) showed low overlap in loci identified as statistical outliers across sub-sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units.

How "simple" methodological decisions affect interpretation of population structure based on reduced representation library DNA sequencing: A case study using the lake whitefish

Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these “simple” methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (F_ST, DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high -m values for F_ST and ADMIXTURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the north-west, Search Bay, North Point and Hammond Bay, showing slight differentiation. Overall, we show that apparently simple decisions about library construction and bioinformatics parameters can have important impacts on the interpretation of population subdivision. Although potentially more costly on a per-locus basis, early investment in striking a balance between the number of loci and sequencing effort is well worth the reduced genomic coverage for population genetics studies. More conservative stringency settings on STACKS parameters lead to a final dataset that was more consistent and robust when examining both weak and strong population differentiation. Overall, we recommend that researchers approach “simple” methodological decisions with caution, especially when working on non-model species for the first time.

Evolution and diversity of two cisco forms in an outlet of glacial Lake Algonquin

The diversity of Laurentian Great Lakes ciscoes (Coregonus artedi, sensu lato) arose via repeated local adaptive divergence including deepwater ciscoes that are now extirpated or threatened. The nigripinnis form, or Blackfin Cisco, is extirpated from the Great Lakes and remains only in Lake Nipigon. Putative nigripinnis populations were recently discovered in sympatry with artedi in a historical drainage system of glacial Lake Algonquin, the precursor of lakes Michigan and Huron. Given the apparent convergence on Great Lakes form, we labeled this form blackfin. Here, we test the hypothesis that nigripinnis may have colonized this area from the Great Lakes as a distinct lineage. It would then represent a relict occurrence of the historical diversity of Great Lakes ciscoes. Alternatively, blackfin could have evolved in situ in several lakes. We captured more than 600 individuals in the benthic or pelagic habitat in 14 lakes in or near Algonquin Provincial Park (Ontario, Canada). Fish were compared based on habitat, morphology, and genetic variation at 6,676 SNPs. Contrary to our expectations, both cisco and blackfin belonged to an Atlantic lineage that colonized the area from the east, not from the Great Lakes. Sympatric cisco and blackfin were closely related while fish from different lakes were genetically differentiated, strongly suggesting the repeated in situ origin of each form. Across lakes, there was a continuum of ecological, morphological, and genetic differentiation that could be associated with alternative resources and lake characteristics. This study uncovers a new component of cisco diversity in inland lakes of Canada that evolved independently from ciscoes of the Laurentian Great lakes. The diversity of cisco revealed in this study and across their Canadian range presents a challenge for designating conservation units at the intraspecific level within the framework of the Committee on the Status of Endangered Wildlife in Canada (COSEWIC).

Cisco diversity in a historical drainage of glacial Lake Algonquin

Cisco (Coregonus artedi (sensu lato) Lesueur, 1818) forms matching in appearance to Blackfin Cisco from the Laurentian Great Lakes occur in four lakes in Algonquin Park, Ontario, Canada, a historical drainage of glacial Lake Algonquin (precursor of lakes Michigan and Huron). Their occurrence may represent colonization from glacial Lake Algonquin drainage patterns 13 000 calibrated years BP or independent evolution within each lake. Gill-raker numbers, temperature at capture depth during lake stratification, and hurdle models of habitat distribution are summarized. Blackfin (nigripinnis-like) in the four lakes had higher gill-raker numbers than artedi-like cisco captured in nearby lakes or within the same lake. Two lakes have a bimodal gill-raker distribution that indicate co-occurrence of two forms. Blackfin occupied the hypolimnion with a peak depth distribution at 20–25 m. Maximum depth for blackfin was 35–40 m. The presence of the opossum shrimp (Mysis diluviana Audzijonyte and Väinölä, 2005) appears necessary for the occurrence of cisco diversity in lakes but not sufficient in all cases. The presence of two forms of cisco in at least two lakes points to the possibility of the colonization hypothesis or the ecological speciation hypothesis as accounting for this phenomenon. Genetic analysis is needed to determine which of these hypotheses best accounts for the occurrence of blackfin in Algonquin Park.

Are we accurately estimating the potential role of pollution in the decline of species at risk in Canada?

Pollution is a pervasive, albeit often invisible, threat to biodiversity in Canada. Currently, the Committee on the Status of Endangered Wildlife in Canada (COSEWIC) relies on expert opinion to assess the scope (i.e., the proportion of a species’ population that may be affected) of pollution to species at risk. Here, we describe a spatially explicit, quantitative method for assessing the scope of pollution as a threat to species at risk in Canada. Using this method, we quantified the geographic co-occurrence of 488 terrestrial and freshwater species and pollution sources and determined that, on average, 57% of the mapped occurrences of each species at risk co-occurred with at least one pollution source. Furthermore, we found a weak correlation between the scope of the threat of pollution as assessed by COSEWIC expert panels and the geographic overlap of species occurrences and pollution sources that we determined with our quantitative method. Experts frequently identified scope of pollution as absent or negligible even for species with extensive co-occurrence with pollution sources, especially vascular plants. Clearly, a quantitative approach is needed to make accurate estimates of the scope of pollution as a threat to species at risk in Canada.

Fine-Scale Ecological and Genetic Population Structure of Two Whitefish (Coregoninae) Species in the Vicinity of Industrial Thermal Emissions

Thermal pollution from industrial processes can have negative impacts on the spawning and development of cold-water fish. Point sources of thermal effluent may need to be managed to avoid affecting discrete populations. Correspondingly, we examined fine-scale ecological and genetic population structure of two whitefish species (Coregonus clupeaformis and Prosopium cylindraceum) on Lake Huron, Canada, in the immediate vicinity of thermal effluent from nuclear power generation. Niche metrics using δ¹³C and δ¹⁵N stable isotopes showed high levels of overlap (48.6 to 94.5%) in resource use by adult fish captured in areas affected by thermal effluent compared to nearby reference locations. Isotopic niche size, a metric of resource use diversity, was 1.3- to 2.8-fold higher than reference values in some thermally affected areas, indicative of fish mixing. Microsatellite analyses of genetic population structure (F_st, STRUCTURE and DAPC) indicated that fish captured at all locations in the vicinity of the power plant were part of a larger population extending beyond the study area. In concert, ecological and genetic markers do not support the presence of an evolutionarily significant unit in the vicinity of the power plant. Thus, future research should focus on the potential impacts of thermal emissions on development and recruitment.