Evolutionary Conservation Genomics Reveals Recent Speciation and Local Adaptation in Threatened Takins.

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan-Hengduan-Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger's rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya-Hengduan biodiversity hotspot.

Genomic Signatures of Coevolution between Nonmodel Mammals and Parasitic Roundworms.

Antagonistic coevolution between host and parasite drives species evolution. However, most of the studies only focus on parasitism adaptation and do not explore the coevolution mechanisms from the perspective of both host and parasite. Here, through the de novo sequencing and assembly of the genomes of giant panda roundworm, red panda roundworm, and lion roundworm parasitic on tiger, we investigated the genomic mechanisms of coevolution between nonmodel mammals and their parasitic roundworms and those of roundworm parasitism in general. The genome-wide phylogeny revealed that these parasitic roundworms have not phylogenetically coevolved with their hosts. The CTSZ and prolyl 4-hydroxylase subunit beta (P4HB) immunoregulatory proteins played a central role in protein interaction between mammals and parasitic roundworms. The gene tree comparison identified that seven pairs of interactive proteins had consistent phylogenetic topology, suggesting their coevolution during host-parasite interaction. These coevolutionary proteins were particularly relevant to immune response. In addition, we found that the roundworms of both pandas exhibited higher proportions of metallopeptidase genes, and some positively selected genes were highly related to their larvae's fast development. Our findings provide novel insights into the genetic mechanisms of coevolution between nonmodel mammals and parasites and offer the valuable genomic resources for scientific ascariasis prevention in both pandas.

A genome-wide scan for copy number variations using high-density single nucleotide polymorphism array in Simmental cattle.

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome-wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high-density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.

Geographical patterns and determinants of insect biodiversity in China.

Insects play important roles in the maintenance of ecosystem functioning and the provision of livelihoods for millions of people. However, compared with terrestrial vertebrates and angiosperms, such as the giant panda, crested ibis, and the metasequoia, insect conservation has not attracted enough attention, and a basic understanding of the geographical biodiversity patterns for major components of insects in China is lacking. Herein, we investigated the geographical distribution of insect biodiversity across multiple dimensions (taxonomic, genetic, and phylogenetic diversity) based on the spatial distribution and molecular DNA sequencing data of insects. Our analysis included 18 orders, 360 families, 5,275 genera, and 14,115 species of insects. The results revealed that Southwestern and Southeastern China harbored higher insect biodiversity and numerous older lineages, representing a museum, whereas regions located in Northwestern China harbored lower insect biodiversity and younger lineages, serving as an evolutionary cradle. We also observed that mean annual temperature and precipitation had significantly positive effects, whereas altitude had significantly negative effects on insect biodiversity in most cases. Moreover, cultivated vegetation harbored the highest insect taxonomic and phylogenetic diversity, and needleleaf and broadleaf mixed forests harbored the highest insect genetic diversity. These results indicated that human activities may positively contribute to insect spatial diversity on a regional scale. Our study fills a knowledge gap in insect spatial diversity in China. These findings could help guide national-level conservation plans and the post-2020 biodiversity conservation framework.

The genome of African manatee Trichechus senegalensis reveals secondary adaptation to the aquatic environment.

Sirenians exhibit unique aquatic adaptations, showcasing both convergent adaptive features shared with cetaceans and unique characteristics such as cold sensitivity and dense bones. Here, we report a chromosome-level genome of the African manatee (Trichechus senegalensis) with high continuity, completeness, and accuracy. We found that genes associated with osteopetrosis have undergone positive selection (CSF1R and LRRK1) or pseudogenized (FAM111A and IGSF23) in the African manatee, potentially contributing to the dense bone formation. The loss of KCNK18 may have increased their sensitivity to cold water temperatures. Moreover, we identified convergent evolutionary signatures in 392 genes among fully aquatic mammals, primarily enriched in skin or skeletal system development and circadian rhythm, which contributed to the transition from terrestrial to fully aquatic lifestyles. The African manatee currently possesses a small effective population size and low genome-wide heterozygosity. Overall, our study provides genetic resources for understanding the evolutionary characteristics and conservation efforts of this species.

Conservation genomics of the critically endangered Chinese pangolin.

The Chinese pangolin (Manis pentadactyla, MP) has been extensively exploited and is now on the brink of extinction, but its population structure, evolutionary history, and adaptive potential are unclear. Here, we analyzed 94 genomes from three subspecies of the Chinese pangolin and identified three distinct genetic clusters (MPA, MPB, and MPC), with MPB further divided into MPB1 and MPB2 subpopulations. The divergence of these populations was driven by past climate change. For MPB2 and MPC, recent human activities have caused dramatic population decline and small population size as well as increased inbreeding, but not decrease in genomic variation and increase in genetic load probably due to strong gene flow; therefore, it is crucial to strengthen in situ habitat management for these two populations. By contrast, although human activities have a milder impact on MPA, it is at high risk of extinction due to long-term contraction and isolation, and genetic rescue is urgently needed. MPB1 exhibited a relatively healthy population status and can potentially serve as a source population. Overall, our findings provide novel insights into the conservation of the Chinese pangolin and biogeography of the mammals of eastern Asia.

Influence of Last Glacial Maximum legacies on functional diversity and community assembly of extant Chinese terrestrial vertebrates.

Contemporary biodiversity patterns are shaped by not only modern climate but also factors such as past climate fluctuations. Investigating the relative degree of paleoclimate legacy could help us understand the formation of current biodiversity patterns. However, an assessment of this issue in China is lacking. Here, we investigated the phylogenetic structure and functional diversity patterns of Chinese terrestrial vertebrates. We found that Southern China harbored higher functional richness, while Northern and Western China were more phylogenetically clustered with higher functional divergence and evenness, indicating environmental filtering effects. Moreover, we found that drastic Last Glacial Maximum climate changes were positively related to phylogenetic clustering, lower functional richness, and higher functional divergence and evenness, although this effect varied among different taxonomic groups. We further found that mammal communities experiencing more drastic Last Glacial Maximum temperature changes were characterized by "faster" life-history trait values. Our findings provide new evidence of the paleoclimate change legacies influencing contemporary biodiversity patterns that will help guide national-level conservation plans.

Conservation priorities for global marine biodiversity across multiple dimensions.

Marine biodiversity plays important roles in ocean ecosystem services and has substantial economic value. Species diversity, genetic diversity and phylogenetic diversity, which reflect the number, evolutionary potential and evolutionary history of species in ecosystem functioning, are three important dimensions of biodiversity. Marine-protected areas have been demonstrated as an effective area-based tool for protecting marine biodiversity, but only 2.8% of the ocean has been fully protected. It is urgent to identify global conservation priority areas and percentage of the ocean across multiple dimensions of biodiversity based on Post-2020 Global Biodiversity Framework. Here, we investigate the spatial distribution of marine genetic and phylogenetic diversity using 80 075 mitochondrial DNA barcode sequences from 4316 species and a newly constructed phylogenetic tree of 8166 species. We identify that the Central Indo-Pacific Ocean, Central Pacific Ocean and Western Indian Ocean harbor high levels of biodiversity across three dimensions of biodiversity, which could be designated as conservation priority areas. We also find that strategically protecting ∼22% of the ocean would allow us to reach the target of conserving ∼95% of currently known taxonomic, genetic and phylogenetic diversity. Our study provides insights into the spatial distribution pattern of multiple marine diversities and the findings would help to design comprehensive conservation schemes for global marine biodiversity.

The rough-toothed dolphin genome provides new insights into the genetic mechanism of its rough teeth.

The rough-toothed dolphin (Steno bredanensis) is characterized by having teeth covered in finely wrinkled vertical ridges, which is a general manifestation of amelogenesis imperfecta. The rough surfaces are hypothesized to be an evolutionary morphological trait of feeding adaptation to increase the dolphin's grip on prey. Here, we assembled a rough-toothed dolphin genome and performed the comparative genomic analysis to reveal the genetic basis of the special enamel. Results showed that genes related to enamel development or dental diseases have undergone diversified adaptive changes that may shape the special enamel morphology of this dolphin species, including positive selection (CLDN19, PRKCE, SSUH2, and WDR72), rapid evolution (LAMB3), or unique amino acid substitutions (AMTN, ENAM, MMP20, and KLK4). Meanwhile, the historical demography of rough-toothed dolphin indicated several distinct population fluctuations associated with climate change. The genome-wide heterozygosity of this dolphin is in the middle of all published data for cetaceans. Although the population is considerable, there may be population or subspecies differentiation, and with the global warming and the increasing disturbance of human activities, we should pay more attention to protection in the future. Together, our study brings new insights into the genetic mechanisms that may have driven the evolution of the special enamel morphology in rough-toothed dolphins and provides the first results of genetic heterozygosity and population historical dynamics of this species, which have important guiding implications for the conservation of this dolphin species.