High genetic load without purging in caribou, a diverse species at risk

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

Paranannizziopsis spp. infections in wild snakes and a qPCR assay for detection of the fungus

The emergence of ophidiomycosis (or snake fungal disease) in snakes has prompted increased awareness of the potential effects of fungal infections on wild reptile populations. Yet, aside from Ophidiomyces ophidiicola, little is known about other mycoses affecting wild reptiles. The closely related genus Paranannizziopsis has been associated with dermatomycosis in snakes and tuataras in captive collections, and P. australasiensis was recently identified as the cause of skin infections in non-native wild panther chameleons (Furcifer pardalis) in Florida, USA. Here we describe five cases of Paranannizziopsis spp. associated with skin lesions in wild snakes in North America and one additional case from a captive snake from Connecticut, USA. In addition to demonstrating that wild Nearctic snakes can serve as a host for these fungi, we also provide evidence that the genus Paranannizziopsis is widespread in wild snakes, with cases being identified in Louisiana (USA), Minnesota (USA), Virginia (USA), and British Columbia (Canada). Phylogenetic analyses conducted on multiple loci of the fungal strains we isolated identified P. australasiensis in Louisiana and Virginia; the remaining strains from Minnesota and British Columbia did not cluster with any of the described species of Paranannizziopsis, although the strains from British Columbia appear to represent a single lineage. Finally, we designed a pan-Paranannizziopsis real-time PCR assay targeting the internal transcribed spacer region 2. This assay successfully detected DNA of all described species of Paranannizziopsis and the two potentially novel taxa isolated in this study and did not cross-react with closely related fungi or other fungi commonly found on the skin of snakes. The assay was 100% sensitive and specific when screening clinical (skin tissue or skin swab) samples, although full determination of the assay's performance will require additional follow up due to the small number of clinical samples (n = 14 from 11 snakes) available for testing in our study. Nonetheless, the PCR assay can provide an important tool in further investigating the prevalence, distribution, and host range of Paranannizziopsis spp. and facilitate more rapid diagnosis of Paranannizziopsis spp. infections that are otherwise difficult to differentiate from other dermatomycoses.

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.

Widespread Exposure to Mosquitoborne California Serogroup Viruses in Caribou, Arctic Fox, Red Fox, and Polar Bears, Canada

Northern Canada is warming at 3 times the global rate. Thus, changing diversity and distribution of vectors and pathogens is an increasing health concern. California serogroup (CSG) viruses are mosquitoborne arboviruses; wildlife reservoirs in northern ecosystems have not been identified. We detected CSG virus antibodies in 63% (95% CI 58%-67%) of caribou (n = 517), 4% (95% CI 2%-7%) of Arctic foxes (n = 297), 12% (95% CI 6%-21%) of red foxes (n = 77), and 28% (95% CI 24%-33%) of polar bears (n = 377). Sex, age, and summer temperatures were positively associated with polar bear exposure; location, year, and ecotype were associated with caribou exposure. Exposure was highest in boreal caribou and increased from baseline in polar bears after warmer summers. CSG virus exposure of wildlife is linked to climate change in northern Canada and sustained surveillance could be used to measure human health risks.

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 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.

Incursions of rabbit haemorrhagic disease virus 2 in Canada-Clinical, molecular and epidemiological investigation

Rabbit haemorrhagic disease virus 2 (RHDV2) is a newly emerging Lagovirus belonging to the family Caliciviridae. After its first discovery in 2010 in France, this highly pathogenic virus rapidly spread to neighbouring countries and has become the dominant strain, replacing the classical RHDV strains. RHDV2 was first reported in North America in 2016 in Mont-Joli, Quebec, Canada, and it was reported again in 2018 and 2019 on Vancouver Island and the southwest mainland of British Columbia (BC). The whole genome sequence of the RHDV2 Quebec isolate resembled the 2011 RHDV2-N11 isolate from Navarra, Spain with 97% identity at the nucleotide level. The epidemiological investigation related to this outbreak involved three hobby farms and one personal residence in Quebec. In February 2018, high mortality was reported in a large colony of feral rabbits on the Vancouver Island University Campus, Nanaimo, BC. The virus identified showed only 93% identity to the Quebec RHDV2 isolate at the nucleotide level. Additional cases of RHDV2 on Vancouver Island and on the BC mainland affecting feral and captive domestic, and commercial rabbits were reported subsequently. Vaccination was recommended to control the outbreak and an inactivated bivalent vaccine was made available to the private veterinary practices. In June 2019, an isolated RHDV2 outbreak was reported in pet rabbits in an apartment building in Vancouver, BC. This virus showed only 97% identity to the RHDV2 isolates responsible for the BC outbreak in 2018 at the nucleotide level, suggesting that it was an independent incursion. The outbreak in BC killed a large number of feral European rabbits; however, there were no confirmed cases of RHD in native rabbit species in BC.

Ecology and Infection Dynamics of Multi-Host Amdoparvoviral and Protoparvoviral Carnivore Pathogens

Amdoparvovirus and Protoparvovirus are monophyletic viral genera that infect carnivores. We performed surveillance for and sequence analyses of parvoviruses in mustelids in insular British Columbia to investigate parvoviral maintenance and cross-species transmission among wildlife. Overall, 19.1% (49/256) of the tested animals were parvovirus-positive. Aleutian mink disease virus (AMDV) was more prevalent in mink (41.6%, 32/77) than martens (3.1%, 4/130), feline panleukopenia virus (FPV) was more prevalent in otters (27.3%, 6/22) than mink (5.2%, 4/77) or martens (2.3%, 3/130), and canine parvovirus 2 (CPV-2) was found in one mink, one otter, and zero ermines (N = 27). Viruses were endemic and bottleneck events, founder effects, and genetic drift generated regional lineages. We identified two local closely related AMDV lineages, one CPV-2 lineage, and five FPV lineages. Highly similar viruses were identified in different hosts, demonstrating cross-species transmission. The likelihood for cross-species transmission differed among viruses and some species likely represented dead-end spillover hosts. We suggest that there are principal maintenance hosts (otters for FPV, raccoons for CPV-2/FPV, mink for AMDV) that enable viral persistence and serve as sources for other susceptible species. In this multi-host system, viral and host factors affect viral persistence and distribution, shaping parvoviral ecology and evolution, with implications for insular carnivore conservation.

Fatal hepatic sarcocystosis in a free-ranging grizzly bear cub associated with Sarcocystis canis-like infection

We describe herein fatal hepatic sarcocystosis in a free-ranging grizzly bear ( Ursus arctos horribilis) cub with apicomplexan infection of the liver and brain, both demonstrating 100% homology for Sarcocystis canis and S. arctosi. Fatal hepatic sarcocystosis in dogs has been etiologically associated with intrahepatic schizonts of S. canis. In black and polar bears, a S. canis-like organism produces schizonts in the liver and massive hepatic necrosis. Although intramuscular sarcocysts, taxa S. arctosi and S. ursusi, have been described in healthy brown and black bears, respectively, they have not been detected in bears with hepatic sarcocystosis, to our knowledge, and it is currently unknown whether bears represent an aberrant or intermediate host.

Helminths of grizzly bears (Ursus arctos) and American black bears (Ursus americanus) in Alberta and British Columbia, Canada

Between May 2011 and June 2013, we collected the carcasses and gastrointestinal tracts of 40 American black bears (Ursus americanus Pallas, 1780) and 13 grizzly bears (Ursus arctos L., 1758) from populations of Alberta and British Columbia, Canada. Specimens were examined for helminths, which were identified to the species level by applying an integrated morphological and molecular approach. Our goal was to investigate parasite biodiversity and infection parameters in the sampled grizzly and black bears. We found seven parasite taxa: Dirofilaria ursi Yamaguti, 1941, Baylisascaris transfuga (Rudolphi, 1819), Uncinaria rauschi Olsen, 1968, Uncinaria yukonensis (Wolfgang, 1956), Taenia arctos Haukisalmi, Lavikainen, Laaksonen and Meri, 2011, Diphyllobothrium dendriticum (Nitzsch, 1824), and Diphyllobothrium nihonkaiense Yamane, Kamo, Bylund and Wikgren, 1986. The statistical significance of infection prevalence, intensity, and abundance for each helminth species was assessed relative to host species, gender, age class, sampling season, and location. This is the first unequivocal report of the potentially zoonotic tapeworms D. dendriticum and D. nihonkaiense in North American bears. Furthermore, we provide insight into the biology and ecology of the nematodes B. transfuga, D. ursi, and species of Uncinaria Frölich, 1789, and enrich the information available on the recently described tapeworm T. arctos.

Cryptosporidium and Giardia in marine-foraging river otters (Lontra canadensis) from the Puget Sound Georgia Basin ecosystem

Species of Cryptosporidium and Giardia can infect humans and wildlife and have the potential to be transmitted between these 2 groups; yet, very little is known about these protozoans in marine wildlife. Feces of river otters (Lontra canadensis), a common marine wildlife species in the Puget Sound Georgia Basin, were examined for species of Cryptosporidium and Giardia to determine their role in the epidemiology of these pathogens. Using ZnSO4 flotation and immunomagnetic separation, followed by direct immunofluorescent antibody detection (IMS/DFA), we identified Cryptosporidium sp. oocysts in 9 fecal samples from 6 locations and Giardia sp. cysts in 11 fecal samples from 7 locations. The putative risk factors of proximate human population and degree of anthropogenic shoreline modification were not associated with the detection of Cryptosporidium or Giardia spp. in river otter feces. Amplification of DNA from the IMS/DFA slide scrapings was successful for 1 sample containing > 500 Cryptosporidium sp. oocysts. Sequences from the Cryptosporidium 18S rRNA and the COWP loci were most similar to the ferret Cryptosporidium sp. genotype. River otters could serve as reservoirs for Cryptosporidium and Giardia species in marine ecosystems. More work is needed to better understand the zoonotic potential of the genotypes they carry as well as their implications for river otter health.

Cryptococcus gattii in wildlife of Vancouver Island, British Columbia, Canada

Although Cryptococcus gattii has emerged as an important pathogen of humans and domestic animals on Vancouver Island, Canada since 1999; its distribution in regional wildlife species is largely unknown. Opportunistic sampling methods were employed to obtain nasal swabs for fungal culture from wild mammal species residing within the coastal Douglas fir biogeoclimatic zone on the southeast coast of the island. Samples were collected from 91 animals representing 14 species. Cryptococcus gattii was isolated from the nasal swabs of two eastern gray squirrels (Sciurus carolinensis) trapped in Duncan, British Columbia. The relative proportion of nasal colonization in wild mammal species is consistent with findings in domestic animals, suggesting that animals may be good indicators of environmental organisms.

Population genetic structure of North American thinhorn sheep (Ovis dalli)

The thinhorn sheep (Ovis dalli ssp.) provides a rare example of a North American large mammal that occupies most of its native range and maintains close to ancestral population size. There are currently two recognized subspecies, Dall's sheep (O. d. dalli) and Stone's sheep (O. d. stonei), the validity of which remains uncertain. We investigated the spatial genetic structure of thinhorn sheep populations representing both subspecies by genotyping individuals (n = 919) from across the species range at 12 variable microsatellite loci. We found high levels of genetic diversity within (HE = 0.722) and significant genetic structure among the 24 sampled areas (FST = 0.160). Genetic distance measures and Bayesian clustering analyses revealed the presence of at least eight subpopulations that are delineated by mountain range topology. A strong overall pattern of isolation-by-distance is evident across the sampling range (r = 0.75, P < 0.001) suggesting limited dispersal and extensive philopatry. Partial Mantel tests of this relationship showed mountain range distinctions represent significant barriers to gene flow (P = 0.0001), supporting the Bayesian analyses. Genetic structure was more strongly pronounced in southern Yukon and Alaska than elsewhere. We also show evidence for genetic differences between the two currently recognized thinhorn subspecies.

Antibodies to selected pathogens in free-ranging terrestrial carnivores and marine mammals in Canada

Antibody titres to selected pathogens (canine adenovirus [CAV-2], feline herpesvirus [FHV], phocine herpesvirus [PHV-1], canine distemper virus, dolphin morbillivirus [DMV], phocine distemper virus [PDV], parainfluenza virus type 3 [PI3], rabies virus, dolphin rhabdovirus [DRV], canine coronavirus, feline coronavirus, feline leukaemia virus, Borrelia burgdorferi and Toxoplasma gondii) were determined in whole blood or serum samples from selected free-ranging terrestrial carnivores and marine mammals, including cougars (Fellis concolor), lynxes (Fellis lynx), American badgers (Taxidea taxus), fishers (Martes pennanti), wolverines (Gulo gulo), wolves (Canis lupus), black bears (Ursus americanus), grizzly bears (Ursus arctos), polar bears (Ursus maritimus), walruses (Odobenus rosmarus) and belugas (Delphinapterus leucas), which had been collected at several locations in Canada between 1984 and 2001. Antibodies to a number of viruses were detected in species in which these infections have not been reported before, for example, antibodies to CAV-2 in walruses, to PDV in black bears, grizzly bears, polar bears, lynxes and wolves, to DMV in grizzly bears, polar bears, walruses and wolves, to PI3 in black bears and fishers, and to DRV in belugas and walruses.

Potential contamination of drinking water with Toxoplasma gondii oocysts

The world's first documented toxoplasmosis outbreak associated with a municipal water supply was recognized in 1995 in Victoria, British Columbia, Canada. It was hypothesized that domestic cat (Felis catus) or cougar (Felis concolor) faeces contaminated a surface water reservoir with Toxoplasma gondii oocysts. An extensive investigation of the Victoria watershed 1 year following the outbreak documented the presence of an endemic T. gondii cycle involving the animals inhabiting the area. Cats and cougars were observed throughout the watershed. Serological evidence of T. gondii infection was demonstrated among domestic cats living in the Victoria area. Cougars were found to shed T. gondii oocysts. Serological evidence of T. gondii infection in deer mice living in the riparian environments of the watershed suggested that T. gondii oocysts were being shed near the water edge. Contamination of Victoria's water supply with T. gondii oocysts potentially occurred during the study period and future waterborne toxoplasmosis outbreaks in this and other communities are possible.