Nuclear and mitochondrial genetic structure in the Eurasian beaver (Castor fiber) - implications for future reintroductions

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Genetic structure of Eurasian beaver in Romania: insights after two decades from the reintroduction

Once exploited for fur, meat, and extracting the yellowish exudate called castoreum, the Eurasian beaver disappeared from Romania during the eighteenth century. After, the reintroductions carried out two decades ago are currently thriving in the Danube River basin. Using nine nSSR markers, we analysed samples from 98 individuals, and we found no genetic substructure, suggesting high dispersal and gene flow capabilities. The stepwise mutation model (SMM) indicated the existence of a recent genetic bottleneck, though the Eurasian beaver retains high levels of genetic diversity and population growth facilitated variation in nSSR loci. A fine-scale spatial correlation in females was detected, contrasting with males’ dispersal on longer distances. While the movement and establishment of individuals’ new territories were made under natural predation pressure, the mix following natural expansion improved the fitness and could contribute to a higher genetic diversity than the source population. With any reintroduction, a focus on capturing individuals from various geographic origins, as well as securing many and suitable founding individuals (adults, subadults, and juveniles) with mixed origins, could secure the post-genetic bottleneck recovery and higher genetic diversity. Beyond this conservation success, future management strategies should consider building a National Action Plan (NAP) for the species, including a permanent genetic monitoring programme for Eurasian beaver.

Eurasian beaver (Castor fiber) health surveillance in Britain: Assessing a disjunctive reintroduced population

BACKGROUND

Numerous translocations of Eurasian beavers have occurred with little implementation of standardised health screening. Pre-release health screening enables the selection of individuals with the best survival prospects and reduces potential health risks, but this is by-passed during unofficial releases. Beaver reintroduction to Britain has been haphazard and currently disjunctive populations of varying status exist.

METHODS

This observational cross section study investigated the health status of three beaver populations, with 90 live beavers tested for a range of pathogens comprising 56 from Tayside (unofficially released Scotland), nine from Knapdale (officially released Scotland) and 25 from Devon (unofficially released England). In addition, a further 32 cadavers were screened (25 from Tayside and seven from Knapdale).

RESULTS

All beavers were in good physical condition, did not harbour any non-native disease or parasites of concern and demonstrated remarkably low levels of any disease or parasite exposure.

CONCLUSION

Beavers are establishing and adapting well to British landscapes and are not acting as reservoirs of significant zoonotic diseases. Official, licensed reintroduction programmes may appear overly convoluted; however, reputational damage of unofficial releases should be considered, along with the health and welfare of the animals involved and collateral damage to other wildlife, domestic animals and humans.

First evaluation of the population structure, genetic diversity and landscape connectivity of the Endangered Arabian tahr

The Arabian tahr (Arabitragus jayakari) occurs only in the mountains of northern Oman and the United Arab Emirates. The species is classified as Endangered due to its small declining population. In this study, we combined genetic and landscape ecology techniques in order to inform landscape scale conservation and genetic management of Arabian tahr. Using 540 base pairs of mitochondrial control region in a dataset of 53 samples, we found eight haplotypes, which fell into two haplogroups. Population genetic analysis using a panel of 14 microsatellite loci also showed a weak, but significant division. Analyses of landscape connectivity supported the genetic results showing poor connectivity between populations in the far south of the study area and those in the north. The most likely location of corridors connecting Arabian tahr populations were identified. Many corridors between tahr populations are impeded by multi-lane highways and restoration of these connections is required to maintain population viability of Arabian tahr. Owing to limited genetic samples outside of Wadi Sareen, further sampling is needed to elucidate both mtDNA and the nuclear structure of Arabian tahr more fully. Our study provides a toolkit that may be used for future genetic and connectivity monitoring of the Arabian tahr population.

De novo characterization of placental transcriptome in the Eurasian beaver (Castor fiber L.)

Our pioneering data provide the first comprehensive view of placental transcriptome of the beaver during single and multiple gestation. RNA-Seq and a de novo approach allowed global pattern identification of C. fiber placental transcriptome. Non-redundant beaver transcriptome comprised 211,802,336 nt of placental transcripts, grouped into 128,459 contigs and clustered into 83,951 unigenes. An Ensembl database search revealed 14,487, 14,994, 15,004, 15,267 and 15,892 non-redundant homologs for Ictidomys tridecemlineatus, Rattus norvegicus, Mus musculus, Homo sapiens and Castor canadensis, respectively. Due to expression levels, the identified transcripts were divided into two sets: non-redundant and highly expressed (FPKM > 2 in at least three examined samples), analysed simultaneously. Among 17,009 highly expressed transcripts, 12,147 had BLASTx hits. GO annotations (175,882) were found for 4301 transcripts that were assigned to biological process (16,386), cellular component (9149) and molecular function (8338) categories; 666 unigenes were also classified into 122 KEGG pathways. Comprehensive analyses were performed for 411 and 3078 highly expressed transcripts annotated with a list of processes linked to ‘placenta’ (31 GO terms) or ‘embryo’ (324 GO terms), respectively. Among transcripts with entire CDS annotation, 281 (placenta) and 34 (embryo) alternative splicing events were identified. A total of 8499 putative SNVs (~ 6.2 SNV/transcript and 1.7 SNV/1 kb) were predicted with 0.1 minimum frequency and maximum variant quality (p value 10e⁻⁹). Our results provide a broad-based characterization of the global expression pattern of the beaver placental transcriptome. Enhancement of transcriptomic resources for C. fiber should improve understanding of crucial pathways relevant to proper placenta development and successful reproduction.

Identification of differentially expressed placental transcripts during multiple gestations in the Eurasian beaver (Castor fiber L.)

The Eurasian beaver is one of the largest rodents that, despite its high impact on the environment, is a non-model species that lacks a reference genome. Characterising genes critical for pregnancy outcome can serve as a basis for identifying mechanisms underlying effective reproduction, which is required for the success of endangered species conservation programs. In the present study, high-throughput RNA sequencing (RNA-seq) was used to analyse global changes in the Castor fiber subplacenta transcriptome during multiple pregnancy. De novo reconstruction of the C. fiber subplacenta transcriptome was used to identify genes that were differentially expressed in placentas (n=5) from two females (in advanced twin and triple pregnancy). Analyses of the expression values revealed 124 contigs with significantly different expression; of these, 55 genes were identified using MegaBLAST. Within this group of differentially expressed genes (DEGs), 18 were upregulated and 37 were downregulated in twins. Most DEGs were associated with the following gene ontology terms: cellular process, single organism process, response to stimulus, metabolic process and biological regulation. Some genes were also assigned to the developmental process, the reproductive process or reproduction. Among this group, four genes (namely keratin 19 (Krt19) and wingless-type MMTV integration site family - member 2 (Wnt2), which were downregulated in twins, and Nik-related kinase (Nrk) and gap junction protein β2 (Gjb2), which were upregulated in twins) were assigned to placental development and nine (Krt19, Wnt2 and integrin α₇ (Itga7), downregulated in twins, and Nrk, gap junction protein β6 (Gjb6), GATA binding protein 6 (Gata6), apolipoprotein A-I (ApoA1), apolipoprotein B (ApoB) and haemoglobin subunit α₁ (HbA1), upregulated in twins) were assigned to embryo development. The results of the present study indicate that the number of fetuses affects the expression profile in the C. fiber subplacental transcriptome. Enhancement of transcriptomic resources for C. fiber will improve understanding of the pathways relevant to proper placental development and successful reproduction.

Identifying source populations for the reintroduction of the Eurasian beaver, Castor fiber L. 1758, into Britain: evidence from ancient DNA

Establishing true phylogenetic relationships between populations is a critical consideration when sourcing individuals for translocation. This presents huge difficulties with threatened and endangered species that have become extirpated from large areas of their former range. We utilise ancient DNA (aDNA) to reconstruct the phylogenetic relationships of a keystone species which has become extinct in Britain, the Eurasian beaver Castor fiber. We sequenced seventeen 492 bp partial tRNAPro and control region sequences from Late Pleistocene and Holocene age beavers and included these in network, demographic and genealogy analyses. The mode of postglacial population expansion from refugia was investigated by employing tests of neutrality and a pairwise mismatch distribution analysis. We found evidence of a pre-Late Glacial Maximum ancestor for the Western C. fiber clade which experienced a rapid demographic expansion during the terminal Pleistocene to early Holocene period. Ancient British beavers were found to originate from the Western phylogroup but showed no phylogenetic affinity to any one modern relict population over another. Instead, we find that they formed part of a large, continuous, pan-Western European clade that harbored little internal substructure. Our study highlights the utility of aDNA in reconstructing population histories of extirpated species which has real-world implications for conservation planning.

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation‐by‐distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation‐by‐distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.

Haematology and Serum Biochemistry Parameters and Variations in the Eurasian Beaver (Castor fiber)

Haematology parameters (N = 24) and serum biochemistry parameters (N = 35) were determined for wild Eurasian beavers (Castor fiber), between 6 months - 12 years old. Of the population tested in this study, N = 18 Eurasian beavers were from Norway and N = 17 originating from Bavaria but now living extensively in a reserve in England. All blood samples were collected from beavers via the ventral tail vein. All beavers were chemically restrained using inhalant isoflurane in 100% oxygen prior to blood sampling. Results were determined for haematological and serum biochemical parameters for the species and were compared between the two different populations with differences in means estimated and significant differences being noted. Standard blood parameters for the Eurasian beaver were determined and their ranges characterised using percentiles. Whilst the majority of blood parameters between the two populations showed no significant variation, haemoglobin, packed cell volume, mean cell haemoglobin and white blood cell counts showed significantly greater values (p<0.01) in the Bavarian origin population than the Norwegian; neutrophil counts, alpha 2 globulins, cholesterol, sodium: potassium ratios and phosphorus levels showed significantly (p<0.05) greater values in Bavarian versus Norwegian; and potassium, bile acids, gamma globulins, urea, creatinine and total calcium values levels showed significantly (p<0.05) greater values in Norwegian versus Bavarian relict populations. No significant differences were noted between male and female beavers or between sexually immature (<3 years old) and sexually mature (≥3 years old) beavers in the animals sampled. With Eurasian beaver reintroduction encouraged by legislation throughout Europe, knowledge of baseline blood values for the species and any variations therein is essential when assessing their health and welfare and the success or failure of any reintroduction program. This is the first study to produce base-line blood values and their variations for the Eurasian beaver.

Holocene re-colonisation, central-marginal distribution and habitat specialisation shape population genetic patterns within an Atlantic European grass species

Corynephorus canescens (L.) P.Beauv. is an outbreeding, short-lived and wind-dispersed grass species, highly specialised on scattered and disturbance-dependent habitats of open sandy sites. Its distribution ranges from the Iberian Peninsula over Atlantic regions of Western and Central Europe, but excludes the two other classical European glacial refuge regions on the Apennine and Balkan Peninsulas. To investigate genetic patterns of this uncommon combination of ecological and biogeographic species characteristics, we analysed AFLP variation among 49 populations throughout the European distribution range, expecting (i) patterns of SW European glacial refugia and post-glacial expansion to the NE; (ii) decreasing genetic diversity from central to marginal populations; and (iii) interacting effects of high gene flow and disturbance-driven genetic drift. Decreasing genetic diversity from SW to NE and distinct gene pool clustering imply refugia on the Iberian Peninsula and in western France, from where range expansion originated towards the NE. High genetic diversity within and moderate genetic differentiation among populations, and a significant pattern of isolation-by-distance indicate a gene flow drift equilibrium within C. canescens, probably due to its restriction to scattered and dynamic habitats and limited dispersal distances. These features, as well as the re-colonisation history, were found to affect genetic diversity gradients from central to marginal populations. Our study emphasises the need for including the specific ecology into analyses of species (re-)colonisation histories and range centre-margin analyses. To account for discontinuous distributions, new indices of marginality were tested for their suitability in studies of centre-periphery gradients.