Characterizing the demographic history and prion protein variation to infer susceptibility to chronic wasting disease in a naïve population of white-tailed deer (Odocoileus virginianus)

Large-scale assessment of genetic structure to assess risk of populations of a large herbivore to disease

Chronic wasting disease (CWD) can spread among cervids by direct and indirect transmission, the former being more likely in emerging areas. Identifying subpopulations allows the delineation of focal areas to target for intervention. We aimed to assess the population structure of white-tailed deer (Odocoileus virginianus) in the northeastern United States at a regional scale to inform managers regarding gene flow throughout the region. We genotyped 10 microsatellites in 5701 wild deer samples from Maryland, New York, Ohio, Pennsylvania, and Virginia. We evaluated the distribution of genetic variability through spatial principal component analysis and inferred genetic structure using non-spatial and spatial Bayesian clustering algorithms (BCAs). We simulated populations representing each inferred wild cluster, wild deer in each state and each physiographic province, total wild population, and a captive population. We conducted genetic assignment tests using these potential sources, calculating the probability of samples being correctly assigned to their origin. Non-spatial BCA identified two clusters across the region, while spatial BCA suggested a maximum of nine clusters. Assignment tests correctly placed deer into captive or wild origin in most cases (94%), as previously reported, but performance varied when assigning wild deer to more specific origins. Assignments to clusters inferred via non-spatial BCA performed well, but efficiency was greatly reduced when assigning samples to clusters inferred via spatial BCA. Differences between spatial BCA clusters are not strong enough to make assignment tests a reliable method for inferring the geographic origin of deer using 10 microsatellites. However, the genetic distinction between clusters may indicate natural and anthropogenic barriers of interest for management.

Estimating sequence diversity of prion protein gene (PRNP) in Portuguese populations of two cervid species: red deer and fallow deer

Among the transmissible spongiform encephalopathies (TSEs), chronic wasting disease (CWD) in cervids is now a rising concern in wildlife within Europe, after the detection of the first case in Norway in 2016, in a wild reindeer and until June 2022 a total of 34 cases were described in Norway, Sweden and Finland. The definite diagnosis is post-mortem, performed in target areas of the brain and lymph nodes. Samples are first screened using a rapid test and, if positive, confirmed by immunohistochemistry and Western immunoblotting. The study of the genetics of the prion protein gene, PRNP, has been proved to be a valuable tool for determining the relative susceptibility to TSEs. In the present study, the exon 3 of PRNP gene of 143 samples from red deer (Cervus elaphus) and fallow deer (Dama dama) of Portugal was analysed. Three single nucleotide polymorphisms (SNPs) were found in red deer - codon A136A, codon T98A, codon Q226E - and no sequence variation was detected in fallow deer. The low genetic diversity found in our samples is compatible with previous studies in Europe. The comparison with results from North America suggests that the free-ranging deer from our study may present susceptibility to CWD, although lack of experimental data and the necessity of continuous survey are necessary to evaluate these populations.

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 (F_ROH ) 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 F_ROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low N_e 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 F_ROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

Genomic architecture of phenotypic extremes in a wild cervid

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses.