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.

A palaeogenomic investigation of overharvest implications in an endemic wild reindeer subspecies.

Overharvest can severely reduce the abundance and distribution of a species and thereby impact its genetic diversity and threaten its future viability. Overharvest remains an ongoing issue for Arctic mammals, which due to climate change now also confront one of the fastest changing environments on Earth. The high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus), endemic to Svalbard, experienced a harvest-induced demographic bottleneck that occurred during the 17-20th centuries. Here, we investigate changes in genetic diversity, population structure, and gene-specific differentiation during and after this overharvesting event. Using whole-genome shotgun sequencing, we generated the first ancient and historical nuclear (n = 11) and mitochondrial (n = 18) genomes from Svalbard reindeer (up to 4000 BP) and integrated these data with a large collection of modern genome sequences (n = 90) to infer temporal changes. We show that hunting resulted in major genetic changes and restructuring in reindeer populations. Near-extirpation followed by pronounced genetic drift has altered the allele frequencies of important genes contributing to diverse biological functions. Median heterozygosity was reduced by 26%, while the mitochondrial genetic diversity was reduced only to a limited extent, likely due to already low pre-harvest diversity and a complex post-harvest recolonization process. Such genomic erosion and genetic isolation of populations due to past anthropogenic disturbance will likely play a major role in metapopulation dynamics (i.e., extirpation, recolonization) under further climate change. Our results from a high-arctic case study therefore emphasize the need to understand the long-term interplay of past, current, and future stressors in wildlife conservation.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights.

BACKGROUND: Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining. RESULTS: We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids. CONCLUSIONS: Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants.

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

Genomic health is dependent on long-term population demographic history.

Current genetic methods of population assessment in conservation biology have been challenged by genome-scale analyses due to their quantitatively novel insights. These analyses include assessments of runs-of-homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white-tailed deer, caribou, and mountain goat. The white-tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome-wide signatures of inbreeding using the inbreeding coefficient F and %ROH (FROH ) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and FROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low Ne of the mountain goat, this supports the idea that the genome-wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white-tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased FROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

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.

Telomere shortening rate predicts species life span.

Telomere shortening to a critical length can trigger aging and shorter life spans in mice and humans by a mechanism that involves induction of a persistent DNA damage response at chromosome ends and loss of cellular viability. However, whether telomere length is a universal determinant of species longevity is not known. To determine whether telomere shortening can be a single parameter to predict species longevities, here we measured in parallel the telomere length of a wide variety of species (birds and mammals) with very different life spans and body sizes, including mouse (Mus musculus), goat (Capra hircus), Audouin's gull (Larus audouinii), reindeer (Rangifer tarandus), griffon vulture (Gyps fulvus), bottlenose dolphin (Tursiops truncatus), American flamingo (Phoenicopterus ruber), and Sumatran elephant (Elephas maximus sumatranus). We found that the telomere shortening rate, but not the initial telomere length alone, is a powerful predictor of species life span. These results support the notion that critical telomere shortening and the consequent onset of telomeric DNA damage and cellular senescence are a general determinant of species life span.

Detecting genetic signals of selection in heavily bottlenecked reindeer populations by comparing parallel founder events.

Founder populations are of special interest to both evolutionary and conservation biologists, but the detection of genetic signals of selection in these populations is challenging due to their demographic history. Geographically separated founder populations likely to have been subjected to similar selection pressures provide an ideal but rare opportunity to overcome these challenges. Here we take advantage of such a situation generated when small, isolated founder populations of reindeer were established on the island of South Georgia, and using this system we look for empirical evidence of selection overcoming strong genetic drift. We generated a 70 k ddRADseq single nucleotide polymorphism database for the two parallel reindeer founder populations and screened for signatures of soft sweeps. We find evidence for a genomic region under selection shared among the two populations, and support our findings with Wright-Fisher model simulations to assess the power and specificity of interpopulation selection scans-namely Bayescan, OutFLANK, PCadapt and a newly developed scan called Genome Wide Differentiation Scan (GWDS)-in the context of pairwise source-founder comparisons. Our simulations indicate that loci under selection in small founder populations are most probably detected by GWDS, and strengthen the hypothesis that the outlier region represents a true locus under selection. We explore possible, relevant functional roles for genes in linkage with the detected outlier loci.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum.

Pleistocene glacial cycles influenced the diversification of high-latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger-scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20-30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre-LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high-latitude species.

Large-scale prion protein genotyping in Canadian caribou populations and potential impact on chronic wasting disease susceptibility.

Polymorphisms within the prion protein gene (Prnp) are an intrinsic factor that can modulate chronic wasting disease (CWD) pathogenesis in cervids. Although wild European reindeer (Rangifer tarandus tarandus) were infected with CWD, as yet there have been no reports of the disease in North American caribou (R. tarandus spp.). Previous Prnp genotyping studies on approximately 200 caribou revealed single nucleotide polymorphisms (SNPs) at codons 2 (V/M), 129 (G/S), 138 (S/N), 146 (N/n) and 169 (V/M). The impact of these polymorphisms on CWD transmission is mostly unknown, except for codon 138. Reindeer carrying at least one allele encoding for asparagine (138NN or 138SN) are less susceptible to clinical CWD upon infection by natural routes, with the majority of prions limited to extraneural tissues. We sequenced the Prnp coding region of two caribou subspecies (n = 986) from British Columbia, Saskatchewan, Yukon, Nunavut and the Northwest Territories, to identify SNPs and their frequencies. Genotype frequencies at codon 138 differed significantly between barren-ground (R. t. groenlandicus) and woodland (R. t. caribou) caribou when we excluded the Chinchaga herd (p < .05). We also found new variants at codons 153 (Y/F) and 242 (P/L). Our findings show that the 138N allele is rare among caribou in areas with higher risk of contact with CWD-infected species. As both subspecies are classified as Threatened and play significant roles in North American Indigenous culture, history, food security and the economy, determining frequencies of Prnp genotypes associated with susceptibility to CWD is important for future wildlife management measures.

Adult survival in migratory caribou is negatively associated with MHC functional diversity.

Genes of the major histocompatibility complex (MHC) are involved in acquired immunity in vertebrates. Only a few studies have investigated the fitness consequences of MHC gene diversity in wild populations. Here, we looked at the association between annual survival and body mass and MHC-DRB exon 2 (MHC-DRB) genetic diversity, obtained from high-throughput sequencing, in two declining migratory caribou (Rangifer tarandus) herds. To disentangle the potential direct and general effects of MHC-DRB genetic diversity, we compared different indices of diversity that were either based on DNA-sequence variation or on physicochemical divergence of the translated peptides, thereby covering a gradient of allelic-to-functional diversity. We found that (1) body mass was not related to MHC-DRB diversity or genotype, and (2) adult survival probability was negatively associated with point accepted mutation distance, a corrected distance that considers the likelihood of each amino acid substitution to be accepted by natural selection. In addition, we found no evidence of fluctuating selection over time on MHC-DRB diversity. We concluded that direct effects were involved in the negative relationship between MHC functional diversity and survival, although the mechanism underlying this result remains unclear. A possible explanation could be that individuals with higher MHC diversity suffer higher costs of immunity (immunopathology). Our results suggest that genetic diversity is not always beneficial even in genes that are likely to be strongly shaped by balancing selection.

Genetic Specificity of the Siberian Forest Reindeer (Rangifer tarandus valentinae Flerov, 1932) of the Kuznetsk Alatau.

This is the first study to show the genetic identity of the Altai-Sayan population of the forest reindeer of the Kuznetsk Alatau (Rangifer tarandus valentinae). The population is characterized by the existence of unique mitochondrial lines, the absence of signs of introgression of domestic rein deer mtDNA, as well as a low level of genetic diversity. In the sample studied, only two nucleotide substitutions (both of them transitions) were revealed, the nucleotide diversity (0.0015 ± 0.00136) was almost ten times lower than in most populations of wild reindeer in Russia and was comparable only with that of some wild reindeer populations of Norway and Svalbard. The haplotype diversity (h) was also relatively low (0.615 ± 0.102).

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

Selection forces that favour different phenotypes in different environments can change frequencies of genes between populations along environmental clines. Clines are also compatible with balancing forces, such as negative frequency-dependent selection (NFDS), which maintains phenotypic polymorphisms within populations. For example, NFDS is hypothesized to maintain partial migration, a dimorphic behavioural trait prominent in species where only a fraction of the population seasonally migrates. Overall, NFDS is believed to be a common phenomenon in nature, yet a scarcity of studies were published linking naturally occurring allelic variation with bimodal or multimodal phenotypes and balancing selection. We applied a Pool-seq approach and detected selection on alleles associated with environmental variables along a North-South gradient in western North American caribou, a species displaying partially migratory behaviour. On 51 loci, we found a signature of balancing selection, which could be related to NFDS and ultimately the maintenance of the phenotypic polymorphisms known within these populations. Yet, remarkably, we detected directional selection on a locus when our sample was divided into two behaviourally distinctive groups regardless of geographic provenance (a subset of GPS-collared migratory or sedentary individuals), indicating that, within populations, phenotypically homogeneous groups were genetically distinctive. Loci under selection were linked to functional genes involved in oxidative stress response, body development and taste perception. Overall, results indicated genetic differentiation along an environmental gradient of caribou populations, which we found characterized by genes potentially undergoing balancing selection. We suggest that the underlining balancing force, NFDS, plays a strong role within populations harbouring multiple haplotypes and phenotypes, as it is the norm in animals, plants and humans too.

No evidence of inbreeding depression in fast declining herds of migratory caribou.

Identifying inbreeding depression early in small and declining populations is essential for management and conservation decisions. Correlations between heterozygosity and fitness (HFCs) provide a way to identify inbreeding depression without prior knowledge of kinship among individuals. In Northern Quebec and Labrador, the size of two herds of migratory caribou (Rivière-George, RG and Rivière-aux-Feuilles, RAF) has declined by one to two orders of magnitude in the last three decades. This raises the question of a possible increase in inbreeding depression originating from, and possibly contributing to, the demographic decline in those populations. Here, we tested for the association of genomic inbreeding indices (estimated with 22,073 SNPs) with body mass and survival in 400 caribou sampled in RG and RAF herds between 1996 and 2016. We found no association of individual heterozygosity or inbreeding coefficient with body mass or annual survival. Furthermore, those genomic inbreeding indices remained stable over the period monitored. These results suggest that the rapid and intense demographic decline of the herds did not cause inbreeding depression in those populations. Although we found no evidence for HFCs, if demographic decline continues, it is possible that such inbreeding depression would be triggered.

Post-glacial recolonization of insular Newfoundland across the Strait of Belle Isle gave rise to an endemic subspecies of woodland caribou, Rangifer tarandus terranovae (Bangs, 1896): evidence from mtDNA haplotypes.

Post-glacial origins of woodland caribou (Rangifer tarandus subsp.) on the island of Newfoundland and their relationship to mainland populations have been uncertain. Sequence analysis of 2223 bp of the mitochondrial DNA control region and cytochrome b gene from 233 Newfoundland caribou identified 32 haplotypes in four major clades. Comparison with other Nearctic caribou confirms a closer affinity of the basal Clade A with animals from the mainland, and as an outgroup to Clades B, C, and D that are endemic to the island. This indicates re-entry of caribou to post-glacial Newfoundland across the Strait of Belle Isle from Labrador, rather than from southern coastal refugia. Newfoundland caribou are a distinct subspecies, Rangifer tarandus terranovae (Bangs, 1896). Hierarchical AMOVA shows significant clinal differentiation of the major clades from northwest to southeast across the island. The isolated Avalon Peninsula population in the extreme southeast is genetically depauperate. Founder effects are evident in herds introduced to previously unoccupied areas by wildlife managers over the past 40-50 years. Reindeer introduced in the early 20th century have not contributed to mtDNA diversity in Newfoundland caribou.

Molecular characterization and expression patterns of ghrelin in the reindeer (Rangifer tarandus).

Ghrelin is a novel growth hormone (GH)-releasing peptide, which has been identified as an endogenous ligand for the GH-sretagogue receptor. The sequence and expression of ghrelin has been determined in many species. In this study, to reveal the molecular characterization and expression patterns of ghrelin in the reindeer (Rangifer tarandus), the full-length DNA and cDNA encoding ghrelin were cloned from reindeer stomach using genome walking and rapid amplification of complementary deoxyribonucleic acid ends (RACE). The expression of ghrelin in almost all tissues was examined by real-time quantitative PCR (RT-qPCR). The 4076 bp amplicon of the ghrelin gene consisting of 4 exons and 3 introns was cloned from reindeer. Results of cDNA cloning and sequence analysis revealed that the full-length ghrelin cDNA was composed of 539 bp that included a 5'-untranslated region (46 bp), an open reading frame (ORF) (351 bp), and a 3'-untranslated region (142 bp). In addition, ghrelin was expressed in the all tissues examined, with the expression in the abomasum significantly higher than that in other tissues (p<0.05), followed by the pancreas, duodenum, testis and oesophagus. The results show that the expression of ghrelin in the reindeer gastrointestinal tract is extensive, suggesting its may have a role in regulating the digestive function.

Loss of connectivity among island-dwelling Peary caribou following sea ice decline.

Global warming threatens to reduce population connectivity for terrestrial wildlife through significant and rapid changes to sea ice. Using genetic fingerprinting, we contrasted extant connectivity in island-dwelling Peary caribou in northern Canada with continental-migratory caribou. We next examined if sea-ice contractions in the last decades modulated population connectivity and explored the possible impact of future climate change on long-term connectivity among island caribou. We found a strong correlation between genetic and geodesic distances for both continental and Peary caribou, even after accounting for the possible effect of sea surface. Sea ice has thus been an effective corridor for Peary caribou, promoting inter-island connectivity and population mixing. Using a time series of remote sensing sea-ice data, we show that landscape resistance in the Canadian Arctic Archipelago has increased by approximately 15% since 1979 and may further increase by 20-77% by 2086 under a high-emission scenario (RCP8.5). Under the persistent increase in greenhouse gas concentrations, reduced connectivity may isolate island-dwelling caribou with potentially significant consequences for population viability.

Cross-Species Application of SNP Chips is Not Suitable for Identifying Runs of Homozygosity.

Cross-species application of single-nucleotide polymorphism (SNP) chips is a valid, relatively cost-effective alternative to the high-throughput sequencing methods generally required to obtain a genome-wide sampling of polymorphisms. Kharzinova et al. (2015) examined the applicability of SNP chips developed in domestic bovids (cattle and sheep) to a semi-wild cervid (reindeer). The ancestors of bovids and cervids diverged between 20 and 30 million years ago (Hassanin and Douzery 2003; Bibi et al. 2013). Empirical work has shown that for a SNP chip developed in a bovid and applied to a cervid species, approximately 50% genotype success with 1% of the loci being polymorphic is expected (Miller et al. 2012). The genotyping of Kharzinova et al. (2015) follows this pattern; however, these data are not appropriate for identifying runs of homozygosity (ROH) and can be problematic for estimating linkage disequilibrium (LD) and we caution readers in this regard.