Conservation genetics and genomics of threatened vertebrates in China

Simplified Genomic Data Revealing the Decline of Aleuritopteris grevilleoides Population Accompanied by the Uplift of Dry-Hot Valley in Yunnan, China

Understanding the evolutionary history of endangered species is crucial for identifying the main reasons for species endangerment in the past and predicting the changing trends and evolutionary directions of their future distribution. In order to study the impact of environmental changes caused by deep valley incision after the uplift of the Qinghai-Tibet Plateau on endangered species, we collected 23 samples belonging to four populations of Aleuritopteris grevilleoides, an endangered fern endemic to the dry-hot valleys (DHV) of Yunnan. Single-nucleotide variation sites (SNPs) were obtained by the genotyping-by-sequencing (GBS) method, and approximately 8085 SNP loci were identified. Through the reconstruction and analysis of genetic diversity, population structure, population dynamics, evolution time, and ancestral geographical distribution, combined with geological historical events such as the formation of dry-hot valleys, this study explores the formation history, current situation, reasons for endangerment and scientifically sound measures for the protection of A. grevilleoides. In our study, A. grevilleoides had low genetic diversity (Obs_Het = 0.16, Exp_Het = 0.32, Pi = 0.33) and a high inbreeding coefficient (Fis = 0.45). The differentiation events were 0.18 Mya, 0.16 Mya, and 0.11 Mya in the A. grevilleoides and may have been related to the formation of terraces within the dry-hot valleys. The history of population dynamics results shows that the diversion of the river resulted in a small amount of gene flow between the two clades, accompanied by a rapid increase in the population at 0.8 Mya. After that, the effective population sizes of A. grevilleoides began to contract continuously due to topographic changes resulting from the continuous expansion of dry-hot valleys. In conclusion, we found that the environmental changes caused by geological events might be the main reason for the changing population size of A. grevilleoides.

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 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.

Chromosome-level genome assembly for giant panda provides novel insights into Carnivora chromosome evolution

BACKGROUND

Chromosome evolution is an important driver of speciation and species evolution. Previous studies have detected chromosome rearrangement events among different Carnivora species using chromosome painting strategies. However, few of these studies have focused on chromosome evolution at a nucleotide resolution due to the limited availability of chromosome-level Carnivora genomes. Although the de novo genome assembly of the giant panda is available, current short read-based assemblies are limited to moderately sized scaffolds, making the study of chromosome evolution difficult.

RESULTS

Here, we present a chromosome-level giant panda draft genome with a total size of 2.29 Gb. Based on the giant panda genome and published chromosome-level dog and cat genomes, we conduct six large-scale pairwise synteny alignments and identify evolutionary breakpoint regions. Interestingly, gene functional enrichment analysis shows that for all of the three Carnivora genomes, some genes located in evolutionary breakpoint regions are significantly enriched in pathways or terms related to sensory perception of smell. In addition, we find that the sweet receptor gene TAS1R2, which has been proven to be a pseudogene in the cat genome, is located in an evolutionary breakpoint region of the giant panda, suggesting that interchromosomal rearrangement may play a role in the cat TAS1R2 pseudogenization.

CONCLUSIONS

We show that the combined strategies employed in this study can be used to generate efficient chromosome-level genome assemblies. Moreover, our comparative genomics analyses provide novel insights into Carnivora chromosome evolution, linking chromosome evolution to functional gene evolution.