Endangered Père David's deer genome provides insights into population recovering

Increased incidences of cervical ribs in deer indicate extinction risk

Mammals as a rule have seven cervical vertebrae, a number which remains remarkably conserved. Occasional deviations of this number are usually due to the presence of cervical ribs on the seventh vertebra, indicating a homeotic transformation from a cervical rib-less vertebra into a thoracic rib-bearing vertebra. These transformations are often associated with major congenital abnormalities or pediatric cancers (pleiotropic effects) that are, at least in humans, strongly selected against. Based on data from Late Pleistocene mammoths (Mammuthus primigenius) and woolly rhinoceroses (Coelodonta antiquitatis) from the North Sea, we hypothesized that high incidences of cervical ribs in declining populations are due to inbreeding and/or adverse conditions impacting early pregnancies. In this study, we investigated the incidence of cervical ribs in an extinct Late Pleistocene megaherbivore, giant deer (Megaloceros giganteus) from Ireland and in the extant highly inbred Père David deer (Elaphurus davidianus) and in twenty other extant species. We show that the incidence of cervical ribs is exceptionally high in both the Irish giant deer and the Père David deer and much higher than in extant outbred deer. Our data support the hypothesis that inbreeding and genetic drift increase the frequencies of maladaptive alleles in populations at risk of extinction. The high incidence of cervical ribs indicates a vulnerable condition, which may have contributed to the extinction of megaherbivore species in the Late Pleistocene. We argue that cervical rib frequency may be a good proxy for extinction risk in inbred populations.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Genetic Differentiation of Reintroduced Père David's Deer (Elaphurus davidianus) Based on Population Genomics Analysis

The reintroduction is an important conservation tool to restore a species in its historically distribution area, but the rate of reintroduction success varies across species or regions due to different reasons. Genetic evaluation is important to the conservation management of reintroduced species. Conservation concerns relate to genetic threats for species with a small population size or severely historically bottle-necked species, such as negative consequences associated with loss of genetic diversity and inbreeding. The last 40years have seen a rapid increasing of population size for Père David's deer (Elaphurus davidianus), which originated from a limited founder population. However, the genetic structure of reintroduced Père David's deer has not been investigated in terms of population genomics, and it is still not clear about the evolutionary history of Père David's deer and to what extent the inbreeding level is. Conservation genomics methods were used to reconstruct the demographic history of Père David's deer, evaluate genetic diversity, and characterize genetic structure among 18 individuals from the captive, free-ranging and wild populations. The results showed that 1,456,457 single nucleotide polymorphisms (SNPs) were obtained for Père David's deer, and low levels of genome-wide genetic diversity were observed in Père David's deer compared with Red deer (Cervus elaphus) and Sika deer (Cervus nippon). A moderate population genetic differentiation was detected among three populations of Père David's deer, especially between the captive population in Beijing Père David's deer park and the free-ranging population in Jiangsu Dafeng National Nature Reserve. The effective population size of Père David's deer started to decline ~25.8ka, and the similar levels of three populations' LD reflected the genetic impacts of long-term population bottlenecks in the Père David's deer. The findings of this study could highlight the necessity of individual exchange between different facilities, and genetic management should generally be integrated into conservation planning with other management considerations.

Host Bias in Diet-Source Microbiome Transmission in Wild Cohabitating Herbivores: New Knowledge for the Evolution of Herbivory and Plant Defense

It is commonly understood that dietary nutrition will influence the composition and function of the animal gut microbiome. However, the transmission of organisms from the diet-source microbiome to the animal gut microbiome in the natural environment remains poorly understood, and elucidating this process may help in understanding the evolution of herbivores and plant defenses. Here, we investigated diet-source microbiome transmission across a range of herbivores (insects and mammals) living in both captive and wild environments. We discovered a host bias among cohabitating herbivores (leaf-eating insects and deer), where a significant portion of the herbivorous insect gut microbiome may originate from the diet, while in deer, only a tiny fraction of the gut microbiome is of dietary origin. We speculated that the putative difference in the oxygenation level in the host digestion systems would lead to these host biases in plant-source (diet) microbiome transmission due to the oxygenation living condition of the dietary plant's symbiotic microbiome. IMPORTANCE We discovered a host bias among cohabitating herbivores (leaf-eating insects and deer), where a significant portion of the herbivorous insect gut microbiome may originate from the diet, while in deer, only a tiny fraction of the gut microbiome is of dietary origin. We speculated that the putative difference in the oxygenation level in the host digestion systems would lead to these host biases in plant-source (diet) microbiome transmission due to the oxygenation living condition of the dietary plant's symbiotic microbiome. This study shed new light on the coevolution of herbivory and plant defense.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.

Whole-genome sequencing of Tarim red deer (Cervus elaphus yarkandensis) reveals demographic history and adaptations to an arid-desert environment

Background

The initiation of desert conditions in the Tarim Basin in China since the late Miocene has led to the significant genetic structuring of local organisms. Tarim Red Deer (Cervus elaphus yarkandensis, TRD) have adapted to the harsh environmental conditions in this basin, including high solar radiation and temperature, aridity, and poor nutritional conditions. However, the underlying genetic basis of this adaptation is poorly understood.

Results

We sequenced the whole genomes of 13 TRD individuals, conducted comparative genomic analyses, and estimated demographic fluctuation. The ∂a∂i model estimated that the TRD and Tule elk (Cervus canadensis nannodes) populations diverged approximately 0.98 Mya. Analyses revealed a substantial influence of the Earth’s climate on the effective population size of TRD, associated with glacial advances and retreat, and human activities likely underlie a recent serious decline in population. A marked bottleneck may have profoundly affected the genetic diversity of TRD populations. We detected a set of candidate genes, pathways, and GO categories related to oxidative stress, water reabsorption, immune regulation, energy metabolism, eye protection, heat stress, respiratory system adaptation, prevention of high blood pressure, and DNA damage and repair that may directly or indirectly be involved in the adaptation of TRD to an arid-desert environment.

Conclusions

Our analyses highlight the role of historical global climates in the population dynamics of TRD. In light of ongoing global warming and the increasing incidence of droughts, our study offers insights into the genomic adaptations of animals, especially TRD, to extreme arid-desert environments and provides a valuable resource for future research on conservation design and biological adaptations to environmental change.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

Chromosome-level genome assembly of Tarim red deer, Cervus elaphus yarkandensis

Tarim red deer (Cervus elaphus yarkandensis) is the only subspecies of red deer (of 22 subspecies) from Central Asia. This species is a desert dweller of the Tarim Basin of southern Xinjiang, China, and exhibits some unique adaptations to the dry and extreme hot climate. We report here the assembly of a Tarim red deer genome employing a 10X Genomics library, termed CEY_v1. Our genome consisted of 2.6 Gb with contig N50 and scaffold N50 of 275.5 Kb and 31.7 Mb, respectively. Around 96% of the assembled sequences were anchored onto 34 chromosomes based on the published high-quality red deer genetic linkage map. More than 94% BUSCOs complete genes (including 90.5% single and 3.6% duplicated ones) were detected in the CEY_v1 and 20,653 genes were annotated. The CEY_v1 is expected to contribute to comparative analysis of genome biology, to evolutionary studies within Cervidae, and to facilitating investigation of mechanisms underlying adaptation of this species to the extreme dry and hot climate.

Measurement(s)	DNA • genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Cervus elaphus

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12410564

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (N_e ), despite very large census sizes in modern C. capreolus populations and show that low N_e has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.

Lessons from bamboo-eating pandas and their gut microbiome: Gut microbiome flow and applications

The giant panda is one of the most endangered mammals in the world, and many studies have revealed their evolutionary adaptation to the local environment (e.g., dietary cellulose and cyanide) on the evidences from population genetics and their gut microbiome. Here, based on the results of our analysis of the giant panda gut microbiome, we concluded that instability and resilience are the two primary characteristics of the giant panda gut microbiome. This basic information may have an impact on giant panda conservation, as well the management of other animal species.

Evaluation of the Function of Wild Animal Gut Microbiomes Using Next-Generation Sequencing and Bioinformatics and its Relevance to Animal Conservation

The relationship between animal conservation and the animal gut microbiome is a hot topic in current microbial ecology research. Our group has recently revealed that the occurrence of diverse combinations of gut microbial compositions and functions (metagenomics) in Père David's deer (Elaphurus davidianus) populations is likely to lead to increased evolutionary potential and resilience in response to environmental changes. Thus, considering the effects of diet on the gut microbiome and the importance of a stable gut microbial community to host health, we suggest that a transitional buffer period (with feeding on a regular diet and a diet from the translocation habitat) is needed before animal translocation. When the gut microbiome enters into relatively stable stages and adapts to the new diet from the translocation site, the time is suitable for translocation. Long-term monitoring of the gut microbiomes of translocated animals (by collecting fresh feces and carrying out next-generation sequencing) is still necessary after their translocation.

Père David's deer gut microbiome changes across captive and translocated populations: Implications for conservation

The gut microbial composition and function are shaped by different factors (e.g., host diet and phylogeny). Gut microbes play an important role in host nutrition and development. The gut microbiome may be used to evaluate the host potential environmental adaptation. In this study, we focused on the coevolution of the gut microbiome of captive and translocated Père David's deer populations (Elaphurus davidianus; Chinese: Père David's deer). To address this, we used several different macro- and micro-ecological approaches (landscape ecology, nutritional methods, microscopy, isotopic analysis, and metagenomics). In this long-term study (2011-2014), we observed some dissimilarities in gut microbiome community and function between the captive and wild/translocated Dafeng Père David's deer populations. These differences might link microbiome composition with deer diet within a given season. The proportion of genes coding for putative enzymes (endoglucanase, beta-glucosidase, and cellulose 1,4-beta-cellobiosidase) involved in cellulose digestion in the gut microbiome of the captive populations was higher than that of the translocated population, possibly because of the high proportion of cellulose, hemicellulose, and lignin in the plants most consumed by the captive populations. However, the two enzymes (natA and natB) involved in sodium transport system were enriched in the gut microbiome in translocated population, possibly because of their high salt diet (e.g., Spartina alterniflora). Thus, our results suggested that Père David's deer gut microorganisms potentially coevolved with host diet, and reflected the local adaptation of translocated population in the new environment (e.g., new dietary plants: Spartina alterniflora). A current problem for Père David's deer conservation is the saturation of captive populations. Given that the putative evolutionary adaptation of Père David's deer gut microbiome and its possible applications in conservation, the large area of wetlands along the Yellow Sea dominated by S. alterniflora might be the major translocation region in the future.