Genetic diversity of captive forest musk deer (Moschus berezovskii) inferred from the mitochondrial DNA control region

Insight into mtDNA diversity of wild forest musk deer (Moschus berezovskii) in Shanxi Province mountains

Shanxi Province is the northernmost distribution range of Chinese forest musk deer (Moschus berezovskii), with wild populations scattered across Zhongtiao (ZT), Lüliang (LL), and Taiyue (TY) mountain ranges. Using line transect surveys and local guide assistance, coupled with infrared camera data, this study systematically collected fecal samples from wild forest musk deer from five nature reserves across the three mountain ranges. Genetic diversity was analyzed using the mtDNA control region, which yielded 129 effective sequences of 656 bp and 20 distinct haplotypes. Haplotype and nucleotide diversities were 0.916 and 0.01505, respectively, indicating a relatively high overall genetic diversity in the studied populations. Significant genetic differentiation was observed between the Lüliang-Taiyue (LLTY) and ZT mountain range populations, with most of the genetic variation existing within the populations. No significant correlation was detected between the geographical and genetic distances of the samples, which may have resulted from historical agricultural activities in Shanxi impeding gene flow among populations. Therefore, genetic exchange programs are recommended to improve the genetic diversity and population fitness of forest musk deer in this region to facilitate their recovery and growth.

Maslinic acid supplementation prevents di(2-ethylhexyl) phthalate-induced apoptosis via PRDX6 in peritubular myoid cells of Chinese forest musk deer

Chinese forest musk deer (FMD), an endangered species, have exhibited low reproductive rates even in captivity due to stress conditions. Investigation revealed the presence of di(2-ethylhexyl) phthalate (DEHP), an environmental endocrine disruptor, in the serum and skin of captive FMDs. Feeding FMDs with maslinic acid (MA) has been observed to alleviate the stress response and improve reproductive rates, although the precise molecular mechanisms remain unclear. Therefore, this study aims to investigate the molecular mechanisms underlying the alleviation of DEHP-induced oxidative stress and cell apoptosis in primary peritubular myoid cells (PMCs) through MA intake. Primary PMCs were isolated and exposed to DEHP in vitro. The results demonstrated that DEHP significantly suppressed antioxidant levels and promoted cell apoptosis in primary PMCs. Moreover, interfering with the expression of PRDX6 was found to induce excessive reactive oxygen species (ROS) production and cell apoptosis in primary PMCs. Supplementation with MA significantly upregulated the expression of PRDX6, thereby attenuating DEHP-induced oxidative stress and cell apoptosis in primary PMCs. These findings provide a theoretical foundation for mitigating stress levels and enhancing reproductive capacity of in captive FMDs.

mtDNA CR Evidence Indicates High Genetic Diversity of Captive Forest Musk Deer in Shaanxi Province, China

Forest musk deer (Moschus berezovskii) are endangered ruminants whose adult males secrete musk. China has been breeding forest musk deer artificially since the 1950s in an effort to restore wild populations, with Shaanxi and Sichuan provinces as the two main sites for captive breeding. Genetic diversity is a significant indicator that determines the long-term viability and status of a population, particularly for species at risk of extinction. In this study, we analyzed the current genetic makeup of seven captive forest musk deer populations in the Shaanxi province, using the mitochondrial DNA (mtDNA) control region (CR) as the molecular marker. We sequenced 604 bp of mtDNA CR, with an average content of A+T higher than G+C. We observed 111 variable sites and 39 different haplotypes from 338 sequences. The nucleotide diversity (Pi) and haplotype diversity (Hd) were 0.02887 and 0.908, respectively. Genetic differentiation between these populations was not significant, and the populations might not have experienced rapid growth. By combining our sequences with previous ones, we identified 65 unique haplotypes with 26 rare haplotypes and estimated a total of 90 haplotypes in Shaanxi province captive populations. The Shaanxi province and Sichuan province obtained 88 haplotypes, the haplotypes from the two populations were mixed together, and the two populations showed moderate genetic differentiation. Our findings suggested that captive forest musk deer populations in the Shaanxi province had high genetic diversity, with a rich founder population of about 90 maternal lines. Additionally, managers could develop genetic management plans for forest musk deer based on the haplotype database. Overall, our study will provide insights and guidelines for the conservation of genetic diversity in captive forest musk deer populations in the Shaanxi province.

Population genomics reveals moderate genetic differentiation between populations of endangered Forest Musk Deer located in Shaanxi and Sichuan

BACKGROUND

Many endangered species exist in small, genetically depauperate, or inbred populations, hence promoting genetic differentiation and reducing long-term population viability. Forest Musk Deer (Moschus berezovskii) has been subject to illegal hunting for hundreds of years due to the medical and commercial values of musk, resulting in a significant decline in population size. However, it is still unclear to what extent the genetic exchange and inbreeding levels are between geographically isolated populations. By using whole-genome data, we reconstructed the demographic history, evaluated genetic diversity, and characterized the population genetic structure of Forest Musk Deer from one wild population in Sichuan Province and two captive populations from two ex-situ centers in Shaanxi Province.

RESULTS

SNP calling by GATK resulted in a total of 44,008,662 SNPs. Principal component analysis (PCA), phylogenetic tree (NJ tree), ancestral component analysis (ADMIXTURE) and the ABBA-BABA test separated Sichuan and Shaanxi Forest Musk Deer as two genetic clusters, but no obvious genetic differentiation was observed between the two captive populations. The average pairwise FST value between the populations in Sichuan and Shaanxi ranged from 0.05-0.07, suggesting a low to moderate genetic differentiation. The mean heterozygous SNPs rate was 0.14% (0.11%-0.15%) for Forest Musk Deer at the genomic scale, and varied significantly among three populations (Chi-square = 1.22, p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest (0.11%). The nucleotide diversity of three populations varied significantly (p < 0.05, Kruskal-Wallis Test), with the Sichuan population having the lowest genetic θπ (1.69 × 10-3).

CONCLUSIONS

Genetic diversity of Forest Musk Deer was moderate at the genomic scale compared with other endangered species. Genetic differentiation between populations in Sichuan and Shaanxi may not only result from historical biogeographical factors but also be associated with contemporary human disturbances. Our findings provide scientific aid for the conservation and management of Forest Musk Deer. They can extend the proposed measures at the genomic level to apply to other musk deer species worldwide.

Genetic evidence indicates the occurrence of the Endangered Kashmir musk deer Moschus cupreus in Uttarakhand, India

The Endangered Kashmir musk deer Moschus cupreus occurs in the western Himalayan region from Nepal to Afghanistan, but there is a lack of comprehensive and reliable information on its range. The region also harbours the Endangered Himalayan musk deer Moschus leucogaster, and this range overlap may have led to misidentification of the two musk deer species and errors in the delimitation of their ranges. Here, using genetic analysis of the mitochondrial DNA control region, we examined the phylogenetic relationship among musk deer samples from three regions in India: Ganderbal District in Jammu and Kashmir, and Kedarnath Wildlife Sanctuary and Nanda Devi Biosphere Reserve, both in Uttarakhand. The Bayesian phylogenetic analysis indicated a close genetic relationship between samples from Jammu and Kashmir, Kedarnath Wildlife Sanctuary and Nanda Devi Biosphere Reserve, validated by previously published sequences of Kashmir musk deer from Nepal. Our analyses confirmed the samples from Uttarakhand to be from the Kashmir musk deer, which was not previously known from this region. Therefore, we recommend further research in this area, to validate species identification and confirm the geographical distribution of the various species of musk deer. In addition, we recommend revision of the range of M. cupreus in the IUCN Red List assessment, to facilitate effective conservation and management of this Endangered species.

Effect of probiotics on diversity and function of gut microbiota in Moschus berezovskii

The forest musk deer, Moschus berezovskii, is a nationally protected species of economic importance in China. However, in captive breeding programmes, they usually die as a result of diarrhoea. In this study, six M. berezovskii were randomly selected and divided into two groups: probiotics group (n = 3) and placebo (control) group (n = 3). The two groups were fed a basal diet that included 2 g probiotics (probiotic group) or 2 g whey powder (placebo group) for 30 days. Faecal samples were collected at day 0, 15 and 30 and evaluated for microbial diversity, species richness and metabolic function. Probiotic intervention significantly improved gut health in M. berezovskii by changing the overall community structure of the gut microbiota. Intake of probiotics reduced the relative abundance of pathogenic bacteria such as Escherichia coli and Citrobacter freundii in the intestinal flora and increased the relative abundance of beneficial Bifidobacterium species and other lactic acid bacteria. At the same time, gut microbiota in the probiotics group were involved in regulating degradation of phenylacetic acid and in dTDP-L-rhamnose synthesis; these processes have the potential to enhance immunity in M. berezovskii. This preliminary study revealed the beneficial effects of probiotics on the gut microbiota of M. berezovskii, which the potential to significantly improve the health, wellbeing and economic value of M. berezovskii.

A Genetic Evaluation System for New Zealand White Rabbit Germplasm Resources Based on SSR Markers

Simple Summary

The New Zealand white rabbit (Oryctolagus cuniculus) is one of the most important breeds of commercial and experimental rabbits in the world, and also one of the most raised rabbit breeds in China. Our goal was to develop a suite of microsatellite markers to aid future conservation genetics research for the Oryctolagus cuniculus breeds. Based on the genetic diversity of 130 New Zealand white rabbits, we obtained a set combination of 22 markers. Then, we performed a genetic analysis of 200 New Zealand white rabbits corresponding to two generations with this combination. It can be used to evaluate the breed conservation of rabbit germplasm resources.

Abstract

At present, there is an abundance of quality domestic rabbit breeds in China. However, due to the lack of technical standards for the genetic evaluation of rabbit germplasm resources, there have been a number of problems, such as poor breed conservation. By studying the genetic diversity of 130 New Zealand white rabbits (regardless of generation), we obtained the best simple sequence repeat (SSR) marker combination. We found that, when using microsatellite markers for the effective genetic evaluation of domestic rabbits, the number of records should be greater than 60 and the marker number more than 22. Through the comparative analysis of 30 combinations of 22 markers, the optimal combination of 22 markers was determined, and the 22 SSR polymorphic loci were distributed on different chromosomes. We performed a genetic analysis of 200 New Zealand white rabbits corresponding to two generations, using the best SSR polymorphic loci combination. There were no significant differences in the genetic diversity parameters between the two generations of rabbits (p > 0.05), indicating that the characteristics of this excellent rabbit germplasm have been effectively preserved. At the same time, we verified that the established method can be used to evaluate the breed conservation of rabbit germplasm resources.

Changes in the Population Genetic Structure of Captive Forest Musk Deer (Moschus Berezovskii) with the Increasing Number of Generation under Closed Breeding Conditions

We investigated the genetic diversity of the population of captive forest musk deer (Moschus berezovskii) in Barkam Musk Deer Breeding Centre using twelve microsatellite markers, and then analyzed the change in genetic structure of successive generation groups from the population. The data provide a new understanding for the evaluation and usage of the breeding management system. Microsatellite marker analysis detected 141 alleles with an average of 11.75 alleles for each marker. The average expected heterozygosity (HE) was 0.731. Performing an F-statistical analysis on the data showed that the genetic diversity of population decreased, and the inbreeding coefficient significant increased with the increase of generation, and FIS of the 1st generation is significantly lower than that of the second to fifth generation (p < 0.01). The result suggested that the captive population was facing the pressure of inbreeding (FIS = 0.115) and the subsequent loss of genetic diversity. Therefore, it is necessary to improve the breeding management system of the captive population by preventing close relatives from mating or inducing new individuals from the exotic population.

Whole-genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features

Background

Siberian musk deer, one of the seven species, is distributed in coniferous forests of Asia. Worldwide, the population size of Siberian musk deer is threatened by severe illegal poaching for commercially valuable musk and meat, habitat losses, and forest fire. At present, this species is categorized as Vulnerable on the IUCN Red List. However, the genetic information of Siberian musk deer is largely unexplored.

Results

Here, we produced 3.10 Gb draft assembly of wild Siberian musk deer with a contig N50 of 29,145 bp and a scaffold N50 of 7,955,248 bp. We annotated 19,363 protein-coding genes and estimated 44.44% of the genome to be repetitive. Our phylogenetic analysis reveals that wild Siberian musk deer is closer to Bovidae than to Cervidae. Comparative analyses showed that the genetic features of Siberian musk deer adapted in cold and high-altitude environments. We sequenced two additional genomes of Siberian musk deer constructed demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Finally, we identified several candidate genes that may play a role in the musk secretion based on transcriptome analysis.

Conclusions

Here, we present a high-quality draft genome of wild Siberian musk deer, which will provide a valuable genetic resource for further investigations of this economically important musk deer.

Analysis of Genetic Diversity and Population Structure in Three Forest Musk Deer Captive Populations with Different Origins

Musk deer (Moschidae), whose secretion is an expensive and irreplaceable component of traditional medicine, have become endangered in the wild due to habitat fragmentation and over-exploitation. In recent years, China has had success in the artificial breeding of forest musk deer, thus relieving the pressure on wild populations. However, many farmed populations are experiencing degradation, and little genetic information is available for conservation management. In this study, we selected 274 individuals from three typical captive populations (originated from the Ta-pa Mountains (Tp), the midrange of the Qinling Mountains (Ql) and the Western Sichuan Plateau (WS), respectively) to evaluate the genetic variations. A total of more than 3.15 billion high-quality clean reads and 4.37 million high-quality SNPs were generated by RAD sequencing. Based on the analysis, we found that captive forest musk deer populations exhibit a relatively low level of genetic diversity. Ql displayed a higher level of genetic diversity than the Tp and WS populations. Tp and WS had experienced population bottlenecks in the past as inferred from the values of Tajima's D. There were high levels of heterozygote deficiency caused by inbreeding within the three populations. Population structure analysis suggested that the three populations have evolved independently, and a moderate amount of genetic differentiation has developed, although there was a low level of gene flow between the Ql and Tp populations. Furthermore, the average quantities of musk secreted by musk deer in the Tp and WS populations were significantly higher than that in the Ql population. The present genetic information should be considered in management plans for the conservation and utilization of musk deer from captive breeding.

The draft genome sequence of forest musk deer (Moschus berezovskii)

Background

The forest musk deer, Moschus berezovskii, is one of seven musk deer (Moschus spp.) and is distributed in Southwest China. Akin to other musk deer, the forest musk deer has been traditionally and is currently hunted for its musk (i.e., global perfume industry). Considerable hunting pressure and habitat loss have caused significant population declines. Consequently, the Chinese government commenced captive breeding programs for musk harvesting in the 1950s. However, the prevalence of fatal diseases is considerably restricting population increases. Disease severity and extent are exacerbated by inbreeding and genetic diversity declines in captive musk deer populations. It is essential that knowledge of captive and wild forest musk deer populations' immune system and genome be gained in order to improve their physical and genetic health. We have thus sequenced the whole genome of the forest musk deer, completed the genomic assembly and annotation, and performed preliminary bioinformatic analyses.

Findings

A total of 407 Gb raw reads from whole-genome sequencing were generated using the Illumina HiSeq 4000 platform. The final genome assembly is around 2.72 Gb, with a contig N50 length of 22.6 kb and a scaffold N50 length of 2.85 Mb. We identified 24,352 genes and found that 42.05% of the genome is composed of repetitive elements. We also detected 1,236 olfactory receptor genes. The genome-wide phylogenetic tree indicated that the forest musk deer was within the order Artiodactyla, and it appeared as the sister clade of four members of Bovidae. In total, 576 genes were under positive selection in the forest musk deer lineage.

Conclusions

We provide the first genome sequence and gene annotation for the forest musk deer. The availability of these resources will be very useful for the conservation and captive breeding of this endangered and economically important species and for reconstructing the evolutionary history of the order Artiodactyla.

Next generation sequencing yields the complete mitogenome of captive forest musk deer, Moschus berezovskii (Ruminantia: Moschidae)

Moschus berezovskii is an endangered species, but its captive populations are valuable on musk secretions in traditional Chinese medicine and perfume manufacture. The mitogenome of M. berezovskii was 16,353 bp in size. Stop codons in 13 PCGs were all typical types except incomplete stop codon T for COX3, ND2 and ND4, and TA for ND3. No tandem repeat was found in control region. Phylogenetic analysis indicated that Moschidae has the closest relationship with Bovidae. We supported that M. berezovskii should be categorized into two subspecies, and suggested that the status of M. chrysogaster JQ608470 should be further investigated.

Blood transcriptomics of captive forest musk deer (Moschus berezovskii) and possible associations with the immune response to abscesses

Forest musk deer (Moschus berezovskii; FMD) are both economically valuable and highly endangered. A problem for FMD captive breeding programs has been the susceptibility of FMD to abscesses. To investigate the mechanisms of abscess development in FMD, the blood transcriptomes of three purulent and three healthy individuals were generated. A total of ~39.68 Gb bases were generated using Illumina HiSeq 4000 sequencing technology and 77,752 unigenes were identified after assembling. All the unigenes were annotated, with 63,531 (81.71%) mapping to at least one database. Based on these functional annotations, 45,798 coding sequences (CDS) were detected, along with 12,697 simple sequence repeats (SSRs) and 65,536 single nucleotide polymorphisms (SNPs). A total of 113 unigenes were found to be differentially expressed between healthy and purulent individuals. Functional annotation indicated that most of these differentially expressed genes were involved in the regulation of immune system processes, particularly those associated with parasitic and bacterial infection pathways.

Yangtze River, an insignificant genetic boundary in tufted deer (Elaphodus cephalophus): the evidence from a first population genetics study

Great rivers were generally looked at as the geographical barrier to gene flow for many taxonomic groups. The Yangtze River is the third largest river in the world, and flows across South China and into the East China Sea. Up until now, few studies have been carried out to evaluate its effect as a geographical barrier. In this study, we attempted to determine the barrier effect of the Yangtze River on the tufted deer (Elaphodus cephalophus) using the molecular ecology approach. Using mitochondrial DNA control region (CR) sequences and 13 nuclear microsatellite loci, we explored the genetic structure and gene flow in two adjacent tufted deer populations (Dabashan and Wulingshan populations), which are separated by the Yangtze River. Results indicated that there are high genetic diversity levels in the two populations, but no distinguishable haplotype group or potential genetic cluster was detected which corresponded to specific geographical population. At the same time, high gene flow was observed between Wulingshan and Dabashan populations. The tufted deer populations experienced population decrease from 0.3 to 0.09 Ma BP, then followed by a distinct population increase. A strong signal of recent population decline (T = 4,396 years) was detected in the Wulingshan population by a Markov-Switching Vector Autoregressions(MSVAR) process population demography analysis. The results indicated that the Yangtze River may not act as an effective barrier to gene flow in the tufted deer. Finally, we surmised that the population demography of the tufted deer was likely affected by Pleistocene climate fluctuations and ancient human activities.

Major histocompatibility complex class II genetic variation in forest musk deer (Moschus berezovskii) in China

The major histocompatibility complex (MHC) plays an important role in the immune system of vertebrates. We used the second exon of four MHC class II genes (DRA, DQA1, DQA2 and DRB3) to assess the overall MHC variation in forest musk deer (Moschus berezovskii). We also compared the MHC variation in captive and wild populations. We observed 22 alleles at four loci (four at DRA, four at DQA1, four at DQA2 and 10 at DRB3), 15 of which were newly identified alleles. Results suggest that forest musk deer maintain relatively high MHC variation, which may result from balancing selection. Moreover, considerable diversity was observed at the DRA locus. We found a high frequency of Mobe-DRA*02, Mobe-DQA1*01 and Mobe-DQA2*05 alleles, which may be important for pathogen resistance. A Ewens-Watterson test showed that the DRB3 locus in the wild population had experienced recent balancing selection. We detected a small divergence at the DRA locus, suggesting the effect of weak positive selection on the DRA gene. Alternatively, this locus may be young and not yet adapted a wide spectrum of alleles for pathogen resistance. The significant heterozygosity deficit observed at the DQA1 and DRB3 loci in the captive population and at all four loci in the wild population may be the result of a population bottleneck. Additionally, MHC genetic diversity was higher in the wild population than in the captive, suggesting that the wild population may have the ability to respond to a wider range of pathogens.

Major histocompatibility complex class II polymorphisms in forest musk deer (Moschus berezovskii) and their probable association with purulent disease

Genes of the major histocompatibility complex (MHC) family are crucial in immune responses because they present pathogenic peptides to T cells. In this study, we analysed the genetic variation in forest musk deer (Moschus berezovskii) MHC II genes and its potential association with musk deer purulent disease. In total, 53 purulent disease-susceptible and 46 purulent disease-resistant individuals were selected for MHC II exon 2 fragment analysis. Among them, 16 DQ alleles and four additional DR alleles were identified, with DQ exon 2 fragments displaying a low level of polymorphism. The nonsynonymous substitutions exceeded the synonymous substitutions in the peptide-binding sites of DQA2, DQB1 and DQB2. Then, 28 MHC II alleles were used to analyse the distribution patterns of purulent disease between the susceptible and resistant groups. Among them, three alleles (DQA1*01, DQA1*02 and DQA2*04) were found to be resistant, and five alleles (DRB3*07, DQA1*03, DQA1*04, DQA2*05 and DQA2*06) were found to increase susceptibility. Additionally, three haplotypes were found to be putatively associated with musk deer purulent disease. However, these three haplotypes were only found in the resistant or susceptible group, and their frequencies were low. The results from our study support a contributory role of MHC II polymorphisms in the development of purulent disease in forest musk deer.