Low levels of genetic differentiation and structure in red fox populations in Eastern Canada and implications for Arctic fox rabies propagation potential.

Rabies is a lethal zoonosis present in most parts of the world which can be transmitted to humans through the bite from an infected mammalian reservoir host. The Arctic rabies virus variant (ARVV) persists mainly in populations of Arctic foxes (Vulpes lagopus), and to a lesser extent in red fox populations (Vulpes vulpes). Red foxes are thought to be responsible for sporadic southward movement waves of the ARVV outside the enzootic area of northern Canada. In this study, we wanted to investigate whether red foxes displayed notable levels of genetic structure across the Quebec-Labrador Peninsula, which includes portions of the provinces of Quebec and Newfoundland-Labrador in Canada, and is a region with a history of southward ARVV movement waves. We combined two datasets that were collected and genotyped using different protocols, totalling 675 red fox individuals across the whole region and genotyped across 13 microsatellite markers. We found two genetic clusters across the region, reflecting a latitudinal gradient, and characterized by low genetic differentiation. We also observed weak but significant isolation by distance, which seems to be marginally more important for females than for males. These findings suggest a general lack of resistance to movement in red fox populations across the Quebec-Labrador Peninsula, regardless of sex. Implications of these findings include additional support for the hypothesis of long-distance southward ARVV propagation through its red fox reservoir host.

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.

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.

CNVs with adaptive potential in Rangifer tarandus: genome architecture and new annotated assembly.

Rangifer tarandus has experienced recent drastic population size reductions throughout its circumpolar distribution and preserving the species implies genetic diversity conservation. To facilitate genomic studies of the species populations, we improved the genome assembly by combining long read and linked read and obtained a new highly accurate and contiguous genome assembly made of 13,994 scaffolds (L90 = 131 scaffolds). Using de novo transcriptome assembly of RNA-sequencing reads and similarity with annotated human gene sequences, 17,394 robust gene models were identified. As copy number variations (CNVs) likely play a role in adaptation, we additionally investigated these variations among 20 genomes representing three caribou ecotypes (migratory, boreal and mountain). A total of 1,698 large CNVs (length > 1 kb) showing a genome distribution including hotspots were identified. 43 large CNVs were particularly distinctive of the migratory and sedentary ecotypes and included genes annotated for functions likely related to the expected adaptations. This work includes the first publicly available annotation of the caribou genome and the first assembly allowing genome architecture analyses, including the likely adaptive CNVs reported here.

Past, present and future contributions of evolutionary biology to wildlife forensics, management and conservation.

Successfully implementing fundamental concepts into concrete applications is challenging in any given field. It requires communication, collaboration and shared will between researchers and practitioners. We argue that evolutionary biology, through research work linked to conservation, management and forensics, had a significant impact on wildlife agencies and department practices, where new frameworks and applications have been implemented over the last decades. The Quebec government's Wildlife Department (MFFP: Ministère des Forêts, de la Faune et des Parcs) has been proactive in reducing the "research-implementation" gap, thanks to prolific collaborations with many academic researchers. Among these associations, our department's outstanding partnership with Dr. Louis Bernatchez yielded significant contributions to harvest management, stocking programmes, definition of conservation units, recovery of threatened species, management of invasive species and forensic applications. We discuss key evolutionary biology concepts and resulting concrete examples of their successful implementation that derives directly or indirectly from this successful partnership. While old and new threats to wildlife are bringing new challenges, we expect recent developments in eDNA and genomics to provide innovative solutions as long as the research-implementation bridge remains open.

Natural hybrids in Atlantic eels (Anguilla anguilla, A. rostrata): evidence for successful reproduction and fluctuating abundance in space and time.

The outcome of natural hybridization is highly variable and depends on the nonexclusive effects of both pre- and post-mating reproductive barriers. The objective of this study was to address three specific questions regarding the dynamics of hybridization between the American and European eels (Anguilla rostrata and Anguilla anguilla). Using 373 AFLP loci, 1127 eels were genotyped, representing different life stages from both continents, as well as multiple Icelandic locations. We first evaluated the extent of hybridization and tested for the occurrence of hybrids beyond the first generation. Second, we tested whether hybrids were randomly distributed across continents and among Icelandic sampling sites. Third, we tested for a difference in the proportion of hybrids between glass eel and yellow eel stages in Iceland. Our results provided evidence for (i) an overall hybrid proportion of 15.5% in Iceland, with values ranging from 6.7% to 100% depending on life stages and locations; (ii) the existence of hybrids beyond the first generation; (iii) a nonrandom geographic distribution of hybrids in the North Atlantic; and (iv) a higher proportion of first and later generation hybrids in yellow eels compared to glass eels, as well as a significant latitudinal gradient in the proportion of hybrids in Icelandic freshwater. We propose that the combined effect of both differential survival of hybrids and variation in hybridization rate through time best explain these patterns. We discuss the possibility that climate change, which is impacting many environmental features in the North Atlantic, may have a determinant effect on the outcome of natural hybridization in Atlantic eels.