Genomics, environment and balancing selection in behaviourally bimodal populations: The caribou case

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Reconstructing past migratory behaviour of reindeer (Rangifer tarandus): Insights from geometric morphometric analysis of proximal phalanx morphology from extant caribou populations

Reindeer mobility patterns vary widely in modern ecosystems, notably between more open or more wooded environments. This renders the reconstruction of past reindeer mobility patterns challenging, while being at the same time key if archaeologists want to better understand the role that reindeer herds played in the subsistence and territorial organisation of Prehistoric hunter-gatherer societies. Studying the morphology associated with different habitats and mobility patterns can be a useful method for understanding past reindeer behaviour. To access paleoecological information, the relationship between locomotor anatomy and substrate type can be explored in modern animals and transposed to the past, as previous studies have proven that an animal´s environment affects bone morphology. In this study, 3D Geometric Morphometrics are used to explore the impact of extant reindeer habitat type and mobility pattern on phalanx morphology. Results obtained reflects on the potential archaeological application of such an approach for paleoecological reconstructions. Size and shape of phalanx vary significantly, yet complex to interpret in light of interplaying factors such as subspecies, sexual dimorphism and the influence of migration costs, snow cover and substrate type. If direct application to the archaeological record remains preliminary, this first study permits to highlight promising avenues for future research.

Genomics of founders for conservation breeding: the Jasper caribou case

Conservation breeding programs are increasingly used as recovery actions for wild animals; bringing founders into captivity to rear captive populations for future reintroduction into the wild. The International Union for the Conservation of Nature recommends that founders should come from genetically close populations and should have sufficient genetic diversity to avoid mating among relatives. Genomic data are highly informative for evaluating founders due to their high resolution and ability to capture adaptive divergence, yet, their application in that context remains limited. Woodland caribou are federally listed as a Species at Risk in Canada, with several populations facing extirpation, such as those in the Rocky Mountains of Alberta and British Columbia (BC). To prevent local extirpation, Jasper National Park (JNP) is proposing a conservation breeding program. We examined single nucleotide polymorphisms for 144 caribou from 11 populations encompassing a 200,0002 km area surrounding JNP to provide information useful for identifying appropriate founders for this program. We found that this area likely hosts a caribou metapopulation historically characterized by high levels of gene flow, which indicates that multiple sources of founders would be appropriate for initiating a breeding program. However, population structure and adaptive divergence analyses indicate that JNP caribou are closest to populations in the BC Columbia range, which also have suitable genetic diversity for conservation breeding. We suggest that collaboration among jurisdictions would be beneficial to implement the program to promote recovery of JNP caribou and possibly other caribou populations in the surrounding area, which is strategically at the periphery of the distribution of this endangered species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01540-3.

Negative frequency dependent selection unites ecology and evolution

From genes to communities, understanding how diversity is maintained remains a fundamental question in biology. One challenging to identify, yet potentially ubiquitous, mechanism for the maintenance of diversity is negative frequency dependent selection (NFDS), which occurs when entities (e.g., genotypes, life history strategies, species) experience a per capita reduction in fitness with increases in relative abundance. Because NFDS allows rare entities to increase in frequency while preventing abundant entities from excluding others, we posit that negative frequency dependent selection plays a central role in the maintenance of diversity. In this review, we relate NFDS to coexistence, identify mechanisms of NFDS (e.g., mutualism, predation, parasitism), review strategies for identifying NFDS, and distinguish NFDS from other mechanisms of coexistence (e.g., storage effects, fluctuating selection). We also emphasize that NFDS is a key place where ecology and evolution intersect. Specifically, there are many examples of frequency dependent processes in ecology, but fewer cases that link this process to selection. Similarly, there are many examples of selection in evolution, but fewer cases that link changes in trait values to negative frequency dependence. Bridging these two well-developed fields of ecology and evolution will allow for mechanistic insights into the maintenance of diversity at multiple levels.

Population transcriptomics uncover the relative roles of positive selection and differential expression in Batrachium bungei adaptation to the Qinghai-Tibetan plateau

KEY MESSAGE

Positive selection genes are related to metabolism, while differentially expressed genes are related to photosynthesis, suggesting that genetic adaptation and expression regulation may play independent roles in different gene classes. Genome-wide investigation of the molecular mechanisms for high-altitude adaptation is an intriguing topic in evolutionary biology. The Qinghai-Tibet Plateau (QTP) with its extremely variable environments is an ideal site for studying high-altitude adaptation. Here, we used transcriptome data of 100 individuals from 20 populations collected from various altitudes on the QTP to investigate the adaptive mechanisms of the aquatic plant Batrachium bungei at both the genetic and transcriptional level. To explore genes and biological pathways that may contribute to QTP adaptation, we employed a two-step approach, in which we identified positively selected genes and differentially expressed genes using the landscape genomic and differential expression approaches. The positive selection analysis showed that genes involved in metabolic regulation played a crucial role in B. bungei adaptation to the extreme environments of the QTP, especially intense ultraviolet radiation. Altitude-based differential expression analysis suggested that B. bungei could increase the rate of energy dissipation or reduce the efficiency of light energy absorption by down regulating the expression of photosynthesis-related genes to adapt to the strong ultraviolet radiation. Weighted gene co-expression network analysis identified ribosomal genes as hubs of altitude adaptation in B. bungei. Only a small part of genes (about 10%) overlapped between positively selected genes and differentially expressed genes in B. bungei, suggesting that genetic adaptation and gene expression regulation might play relatively independent roles in different categories of functional genes. Taken together, this study enriches our understanding of the high-altitude adaptation mechanism of B. bungei on the QTP.

Genomic correlates for migratory direction in a free-ranging cervid

Animal migrations are some of the most ubiquitous and one of the most threatened ecological processes globally. A wide range of migratory behaviours occur in nature, and this behaviour is not uniform among and within species, where even individuals in the same population can exhibit differences. While the environment largely drives migratory behaviour, it is necessary to understand the genetic mechanisms influencing migration to elucidate the potential of migratory species to cope with novel conditions and adapt to environmental change. In this study, we identified genes associated with a migratory trait by undertaking pooled genome-wide scans on a natural population of migrating mule deer. We identified genomic regions associated with variation in migratory direction, including FITM1, a gene linked to the formation of lipids, and DPPA3, a gene linked to epigenetic modifications of the maternal line. Such a genetic basis for a migratory trait contributes to the adaptive potential of the species and might affect the flexibility of individuals to change their behaviour in the face of changes in their environment.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus)

BACKGROUND

Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species.

RESULTS

From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip.

CONCLUSION

This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Seasonal movements in caribou ecotypes of Western Canada

BACKGROUND

Several migratory ungulates, including caribou, are dramatically declining. Caribou of the Barren-ground ecotype, which forms its own subspecies, are known to be mainly migratory. By contrast, within the Woodland subspecies, animals of the Boreal ecotype are known to be mainly sedentary, while those within the Northern and Central Mountain ecotypes to be partially migratory, with only some individuals migrating. Promotion of conservation actions (e.g., habitat protection) that are specific to both residents and migrants, as well as to the areas they frequent seasonally (which may be separate for migrants), requires distinguishing migration from other movement behaviours, which might be a challenge.

METHODS

We aimed at assessing seasonal movement behaviours, including migratory, resident, dispersing, and nomadic, for caribou belonging to the Barren-ground and Woodland subspecies and ecotypes. We examined seasonal displacement, both planar and altitudinal, and seasonal ranges overlap for 366 individuals that were GPS-collared in Northern and Western Canada. Lastly, we assessed the ability of caribou individuals to switch between migratory and non-migratory movement behaviours between years.

RESULTS

We detected migratory behaviour within each of the studied subspecies and ecotypes. However, seasonal ranges overlap (an index of sedentary behaviour) varied, with proportions of clear migrants (0 overlap) of 40.94% for Barren-ground caribou and 23.34% for Woodland caribou, and of 32.95%, 54.87%, and 8.86% for its Northern Mountain, Central Mountain, and Boreal ecotype, respectively. Plastic switches of individuals were also detected between migratory, resident, dispersing, and nomadic seasonal movements performed across years.

CONCLUSIONS

Our unexpected findings of marked seasonal movement plasticity in caribou indicate that this phenomenon should be better studied to understand the resilience of this endangered species to habitat and climatic changes. Our results that a substantial proportion of individuals engaged in seasonal migration in all studied ecotypes indicate that caribou conservation plans should account for critical habitat in both summer and winter ranges. Accordingly, conservation strategies are being devised for the Woodland subspecies and its ecotypes, which were found to be at least partially migratory in this study. Our findings that migration is detectable with both planar and altitudinal analyses of seasonal displacement provide a tool to better define seasonal ranges, also in mountainous and hilly environments, and protect habitat there.

Selection of both habitat and genes in specialized and endangered caribou

Genetic mechanisms determining habitat selection and specialization of individuals within species have been hypothesized, but not tested at the appropriate individual level in nature. In this work, we analyzed habitat selection for 139 GPS-collared caribou belonging to three declining ecotypes sampled throughout Northwestern Canada. We used Resource Selection Functions (RSFs) comparing resources at used and available locations. We found that the three caribou ecotypes differed in their use of habitat suggesting specialization. On expected grounds, we also found differences in habitat selection between summer and winter, but also, originally, among the individuals within an ecotype. We next obtained Single Nucleotide Polymorphisms (SNPs) for the same caribou individuals, we detected those associated to habitat selection, and then identified genes linked to these SNPs. These genes had functions related in other organisms to habitat and dietary specializations, and climatic adaptations. We therefore suggest that individual variation in habitat selection was based on genotypic variation in the SNPs of individual caribou, indicating that genetic forces underlie habitat and diet selection in the species. We also suggest that the associations between habitat and genes that we detected may lead to lack of resilience in the species, thus contributing to caribou endangerment. Our work emphasizes that similar mechanisms may exist for other specialized, endangered species. This article is protected by copyright. All rights reserved.

Genomic legacy of migration in endangered caribou

Wide-ranging animals, including migratory species, are significantly threatened by the effects of habitat fragmentation and habitat loss. In the case of terrestrial mammals, this results in nearly a quarter of species being at risk of extinction. Caribou are one such example of a wide-ranging, migratory, terrestrial, and endangered mammal. In populations of caribou, the proportion of individuals considered as "migrants" can vary dramatically. There is therefore a possibility that, under the condition that migratory behavior is genetically determined, those individuals or populations that are migratory will be further impacted by humans, and this impact could result in the permanent loss of the migratory trait in some populations. However, genetic determination of migration has not previously been studied in an endangered terrestrial mammal. We examined migratory behavior of 139 GPS-collared endangered caribou in western North America and carried out genomic scans for the same individuals. Here we determine a genetic subdivision of caribou into a Northern and a Southern genetic cluster. We also detect >50 SNPs associated with migratory behavior, which are in genes with hypothesized roles in determining migration in other organisms. Furthermore, we determine that propensity to migrate depends upon the proportion of ancestry in individual caribou, and thus on the evolutionary history of its migratory and sedentary subspecies. If, as we report, migratory behavior is influenced by genes, caribou could be further impacted by the loss of the migratory trait in some isolated populations already at low numbers. Our results indicating an ancestral genetic component also suggest that the migratory trait and their associated genetic mutations could not be easily re-established when lost in a population.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.