Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals.

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

An Enhanced Method for the Use of Reptile Skin Sheds as a High-Quality DNA Source for Genome Sequencing.

With the emergence of high-throughput sequencing technology, a number of non-avian reptile species have been sequenced at the genome scale, shedding light on various scientific inquiries related to reptile ecology and evolution. However, the routine requirement of tissue or blood samples for genome sequencing often poses challenges in many elusive reptiles, hence limiting the application of high-throughput sequencing technologies to reptile studies. An alternative reptilian DNA resource suitable for genome sequencing is in urgent need. Here, we used the corn snake (Pantherophis guttatus) as a reptile model species to demonstrate that the shed skin is a high-quality DNA source for genome sequencing. Skin sheds provide a noninvasive type of sample that can be easily collected without restraining or harming the animal. Our findings suggest that shed skin from corn snakes yields DNA of sufficient quantity and quality that are comparable to tissue DNA extracts. Genome sequencing data analysis revealed that shed skin DNA is subject to bacteria contamination at variable levels, which is a major issue related to shed skin DNA and may be addressed by a modified DNA extraction method through introduction of a 30 min pre-digestion step. This study provides an enhanced method for the use of reptile shed skins as a high-quality DNA source for whole genome sequencing. Utilizing shed skin DNA enables researchers to overcome the limitations generally associated with obtaining traditional tissue or blood samples and promises to facilitate the application of genome sequencing in reptilian research.

Genomic analysis reveals a cryptic pangolin species.

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Ancient DNA reveals genetic admixture in China during tiger evolution.

The tiger (Panthera tigris) is a charismatic megafauna species that originated and diversified in Asia and probably experienced population contraction and expansion during the Pleistocene, resulting in low genetic diversity of modern tigers. However, little is known about patterns of genomic diversity in ancient populations. Here we generated whole-genome sequences from ancient or historical (100-10,000 yr old) specimens collected across mainland Asia, including a 10,600-yr-old Russian Far East specimen (RUSA21, 8× coverage) plus six ancient mitogenomes, 14 South China tigers (0.1-12×) and three Caspian tigers (4-8×). Admixture analysis showed that RUSA21 clustered within modern Northeast Asian phylogroups and partially derived from an extinct Late Pleistocene lineage. While some of the 8,000-10,000-yr-old Russian Far East mitogenomes are basal to all tigers, one 2,000-yr-old specimen resembles present Amur tigers. Phylogenomic analyses suggested that the Caspian tiger probably dispersed from an ancestral Northeast Asian population and experienced gene flow from southern Bengal tigers. Lastly, genome-wide monophyly supported the South China tiger as a distinct subspecies, albeit with mitochondrial paraphyly, hence resolving its longstanding taxonomic controversy. The distribution of mitochondrial haplogroups corroborated by biogeographical modelling suggested that Southwest China was a Late Pleistocene refugium for a relic basal lineage. As suitable habitat returned, admixture between divergent lineages of South China tigers took place in Eastern China, promoting the evolution of other northern subspecies. Altogether, our analysis of ancient genomes sheds light on the evolutionary history of tigers and supports the existence of nine modern subspecies.

Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

Genetic mapping and molecular mechanism behind color variation in the Asian vine snake.

BACKGROUND: Reptiles exhibit a wide variety of skin colors, which serve essential roles in survival and reproduction. However, the molecular basis of these conspicuous colors remains unresolved. RESULTS: We investigate color morph-enriched Asian vine snakes (Ahaetulla prasina), to explore the mechanism underpinning color variations. Transmission electron microscopy imaging and metabolomics analysis indicates that chromatophore morphology (mainly iridophores) is the main basis for differences in skin color. Additionally, we assemble a 1.77-Gb high-quality chromosome-anchored genome of the snake. Genome-wide association study and RNA sequencing reveal a conservative amino acid substitution (p.P20S) in SMARCE1, which may be involved in the regulation of chromatophore development initiated from neural crest cells. SMARCE1 knockdown in zebrafish and immunofluorescence verify the interactions among SMARCE1, iridophores, and tfec, which may determine color variations in the Asian vine snake. CONCLUSIONS: This study reveals the genetic associations of color variation in Asian vine snakes, providing insights and important resources for a deeper understanding of the molecular and genetic mechanisms related to reptilian coloration.

Phylogeography of the endangered "eyed" turtles (genus Sacalia) and the discovery of a lineage derived from natural interspecific hybridization.

The herpetofauna of the Indomalayan bioregion of Asia suffers from severe habitat loss, unsustainable harvesting, and lack of research and conservation. Here, we investigated the range-wide phylogeography of the endangered "eyed" turtles (genus Sacalia, including the Beale's Eyed Turtle S. bealei and the Four-eyed Turtle S. quadriocellata) and discovered a natural interspecific hybrid turtle population in China. Based on phylogeny of the mitochondrial Cytochrome b gene of 101 samples in this study and public data, three major clades and six subclades were identified: S. bealei (SBE) in eastern-southern China, east S. quadriocellata in South China (northern east [SQUen] and southern east [SQUes] subclades), and west S. quadriocellata mainly in Vietnam (northern west [SQUwn], central west [SQUwc], and southern west [SQUws] subclades). We sequenced 16 nuclear DNA loci of 87 samples from SBE, SQUen, SQUes, and SQUwn subclades. Population genetic clustering analysis suggested a structure similar to the mitochondrial phylogeny, where most samples were classified into four genetic clusters corresponding to the four mtDNA subclades. However, a proportion of samples carrying SQUen mtDNA haplotypes formed an additional distinct cluster SHY. Those samples are found in the contact zone of the two species bearing mosaic and intermediate morphological characteristics. We detected an admixed ancestry in SHY from SBE and SQUen that conformed to an intrapopulation breeding scenario for at least hundreds of generations after the initial hybrid event, leading to a conclusion that SHY is a distinct and near-panmictic population derived from natural interspecific hybridization. In addition, SQUes (Hainan Island endemic) is of special concern due to significant isolation and low genetic diversity. We suggest that seven evolutionarily significant units should be recognized to facilitate appropriate conservation actions. These findings also highlight the urgent need for further herpetological research and conservation in this region.

The genetics and evolution of eye color in domestic pigeons (Columba livia).

The eye color of birds, generally referring to the color of the iris, results from both pigmentation and structural coloration. Avian iris colors exhibit striking interspecific and intraspecific variations that correspond to unique evolutionary and ecological histories. Here, we identified the genetic basis of pearl (white) iris color in domestic pigeons (Columba livia) to explore the largely unknown genetic mechanism underlying the evolution of avian iris coloration. Using a genome-wide association study (GWAS) approach in 92 pigeons, we mapped the pearl iris trait to a 9 kb region containing the facilitative glucose transporter gene SLC2A11B. A nonsense mutation (W49X) leading to a premature stop codon in SLC2A11B was identified as the causal variant. Transcriptome analysis suggested that SLC2A11B loss of function may downregulate the xanthophore-differentiation gene CSF1R and the key pteridine biosynthesis gene GCH1, thus resulting in the pearl iris phenotype. Coalescence and phylogenetic analyses indicated that the mutation originated approximately 5,400 years ago, coinciding with the onset of pigeon domestication, while positive selection was likely associated with artificial breeding. Within Aves, potentially impaired SLC2A11B was found in six species from six distinct lineages, four of which associated with their signature brown or blue eyes and lack of pteridine. Analysis of vertebrate SLC2A11B orthologs revealed relaxed selection in the avian clade, consistent with the scenario that during and after avian divergence from the reptilian ancestor, the SLC2A11B-involved development of dermal chromatophores likely degenerated in the presence of feather coverage. Our findings provide new insight into the mechanism of avian iris color variations and the evolution of pigmentation in vertebrates.

Genomic evidence for the Chinese mountain cat as a wildcat conspecific (Felis silvestris bieti) and its introgression to domestic cats.

The Qinghai-Tibet Plateau endemic Chinese mountain cat has a controversial taxonomic status, whether it is a true species or a wildcat (Felis silvestris) subspecies and whether it has contributed to cat (F. s. catus) domestication in East Asia. Here, we sampled F. silvestris lineages across China and sequenced 51 nuclear genomes, 55 mitogenomes, and multilocus regions from 270 modern or museum specimens. Genome-wide analyses classified the Chinese mountain cat as a wildcat conspecific F. s. bieti, which was not involved in cat domestication of China, thus supporting a single domestication origin arising from the African wildcat (F. s. lybica). A complex hybridization scenario including ancient introgression from the Asiatic wildcat (F. s. ornata) to F. s. bieti, and contemporary gene flow between F. s. bieti and sympatric domestic cats that are likely recent Plateau arrivals, raises the prospect of disrupted wildcat genetic integrity, an issue with profound conservation implications.

Sorting Out the Genetic Background of the Last Surviving South China Tigers.

The South China tiger (Panthera tigris amoyensis) is endemic to China and also the most critically endangered subspecies of living tigers. It is considered extinct in the wild and only about 150 individuals survive in captivity to date, whose genetic heritage, however, is ambiguous and controversial. Here, we conducted an explicit genetic assessment of 92 studbook-registered South China tigers from 14 captive facilities using a subspecies-diagnostic system in the context of comparison with other voucher specimens to evaluate the genetic ancestry and level of distinctiveness of the last surviving P. t. amoyensis. Three mtDNA haplotypes were identified from South China tigers sampled in this study, including a unique P. t. amoyensis AMO1 haplotype not found in other subspecies, a COR1 haplotype that is widespread in Indochinese tigers (P. t. corbetti), and an ALT haplotype that is characteristic of Amur tigers (P. t. altaica). Bayesian STRUCTURE analysis and parentage verification confirmed the verified subspecies ancestry (VSA) as the South China tiger in 74 individuals. Genetic introgression from other tigers was detected in 18 tigers, and subsequent exclusion of these and their offspring from the breeding program is recommended. Both STRUCTURE clustering and microsatellite-based phylogenetic analyses demonstrated a close genetic association of the VSA South China tigers to Indochinese tigers, an issue that could only be elucidated by analysis of historical South China tiger specimens with wild origin. Our results also indicated a moderate level of genetic diversity in the captive South China tiger population, suggesting a potential for genetic restoration.

Tigers of the World: Genomics and Conservation.

Of all the big cats, or perhaps of all the endangered wildlife, the tiger may be both the most charismatic and most well-recognized flagship species in the world. The rapidly changing field of molecular genetics, particularly advances in genome sequencing technologies, has provided new tools to reconstruct what characterizes a tiger. Here we review how applications of molecular genomic tools have been used to depict the tiger's ancestral roots, phylogenetic hierarchy, demographic history, morphological diversity, and genetic patterns of diversification on both temporal and geographical scales. Tiger conservation, stabilization, and management are important areas that benefit from use of these genome resources for developing survival strategies for this charismatic megafauna both in situ and ex situ.

Genome-Wide Evolutionary Analysis of Natural History and Adaptation in the World's Tigers.

No other species attracts more international resources, public attention, and protracted controversies over its intraspecific taxonomy than the tiger (Panthera tigris) [1, 2]. Today, fewer than 4,000 free-ranging tigers survive, covering only 7% of their historical range, and debates persist over whether they comprise six, five, or two subspecies [3-6]. The lack of consensus over the number of tiger subspecies has partially hindered the global effort to recover the species from the brink of extinction, as both captive breeding and landscape intervention of wild populations increasingly require an explicit delineation of the conservation management units [7]. The recent coalescence to a late Pleistocene bottleneck (circa 110 kya) [5, 8, 9] poses challenges for detecting tiger subspecific morphological traits, suggesting that elucidating intraspecific evolution in the tiger requires analyses at the genomic scale. Here, we present whole-genome sequencing analyses from 32 voucher specimens that resolve six statistically robust monophyletic clades corresponding to extant subspecies, including the recently recognized Malayan tiger (P. tigris jacksoni). The intersubspecies gene flow is very low, corroborating the recognized phylogeographic units. We identified multiple genomic regions that are candidates for identifying the adaptive divergence of subspecies. The body-size-related gene ADH7 appears to have been strongly selected in the Sumatran tiger, perhaps in association with adaptation to the tropical Sunda Islands. The identified genomic signatures provide a solid basis for recognizing appropriate conservation management units in the tiger and can benefit global conservation strategic planning for this charismatic megafauna icon.

The Complete Phylogeny of Pangolins: Scaling Up Resources for the Molecular Tracing of the Most Trafficked Mammals on Earth.

Pangolins, considered the most-trafficked mammals on Earth, are rapidly heading to extinction. Eight extant species of these African and Asian scale-bodied anteaters are commonly recognized, but their evolutionary relationships remain largely unexplored. Here, we present the most comprehensive phylogenetic assessment of pangolins, based on genetic variation of complete mitogenomes and 9 nuclear genes. We confirm deep divergence among Asian and African pangolins occurring not later than the Oligocene-Miocene boundary ca. 23 million years ago (Ma) (95% HPD = 18.7-27.2), limited fossil evidence suggesting dispersals from Europe. We recognize 3 genera including Manis (Asian pangolins), Smutsia (large African pangolins), and Phataginus (small African pangolins), which first diversified in the Middle-Upper Miocene (9.8-13.3 Ma) through a period of gradual cooling coinciding with a worldwide taxonomic diversification among mammals. Based on large mitogenomic distances among the 3 genera (18.3-22.8%) and numerous (18) morphological traits unique to Phataginus, we propose the subfamily Phatagininae subfam. nov. to designate small African pangolins. In contrast with the morphological-based literature, our results establish that the thick-tailed pangolin (Manis crassicaudata) is sister-species of the Sunda (Manis javanica) and Palawan (Manis culionensis) pangolins. Mitogenomic phylogenetic delineations supported additional pangolin species subdivisions (n = 13), including 6 African common pangolin (Phataginus tricuspis) lineages, but these patterns were not fully supported by our multi-locus approach. Finally, we identified more than 5000 informative mitogenomic sites and diagnostic variation from 5 nuclear genes among all species and lineages of pangolins, providing an important resource for further research and for effectively tracing the worldwide pangolin trade.

Effects of the Qinghai-Tibet Railway on the Landscape Genetics of the Endangered Przewalski's Gazelle (Procapra przewalskii).

The Przewalski's gazelle (Procapra przewalskii) is one of the most endangered ungulates in the world, with fewer than 2,000 individuals surviving in nine habitat fragments on the Qinghai-Tibet Plateau and isolated by human settlements and infrastructure. In particular, the Qinghai-Tibet railway, which crosses the largest part of the gazelle's distribution, remains a major concern because of its potential to intensify landscape genetic differentiation. Here, using mtDNA sequencing and microsatellite genotyping to analyze 275 Przewalski's gazelle samples collected throughout the range, we observed low level of genetic diversity (mtDNA π = 0.0033) and strong phylogeographic structure. Overall, the nine patches of gazelles can be further clustered into five populations, with a strong division between the eastern vs. western side of Qinghai Lake. Our study provides the first evidence of the genetic divergence between the Haergai North and Haergai South gazelle populations, corresponding to the recent construction of a wired enclosure along the Qinghai-Tibet railway less than ten years ago, an equivalent of five generations. Well-designed wildlife corridors across the railway along with long-term monitoring of the anthropogenic effects are therefore recommended to alleviate further habitat fragmentation and loss of genetic diversity in Przewalski's gazelle.

Origin of the Reflectin Gene and Hierarchical Assembly of Its Protein.

Cephalopods, the group of animals including octopus, squid, and cuttlefish, have remarkable ability to instantly modulate body coloration and patterns so as to blend into surrounding environments [1, 2] or send warning signals to other animals [3]. Reflectin is expressed exclusively in cephalopods, filling the lamellae of intracellular Bragg reflectors that exhibit dynamic iridescence and structural color change [4]. Here, we trace the possible origin of the reflectin gene back to a transposon from the symbiotic bioluminescent bacterium Vibrio fischeri and report the hierarchical structural architecture of reflectin protein. Intrinsic self-assembly, and higher-order assembly tightly modulated by aromatic compounds, provide insights into the formation of multilayer reflectors in iridophores and spherical microparticles in leucophores and may form the basis of structural color change in cephalopods. Self-assembly and higher-order assembly in reflectin originated from a core repeating octapeptide (here named protopeptide), which may be from the same symbiotic bacteria. The origin of the reflectin gene and assembly features of reflectin protein are of considerable biological interest. The hierarchical structural architecture of reflectin and its domain and protopeptide not only provide insights for bioinspired photonic materials but also serve as unique "assembly tags" and feasible molecular platforms in biotechnology.

Phylogeography of the heavily poached African common pangolin (Pholidota, Manis tricuspis) reveals six cryptic lineages as traceable signatures of Pleistocene diversification.

Knowledge on faunal diversification in African rainforests remains scarce. We used phylogeography to assess (i) the role of Pleistocene climatic oscillations in the diversification of the African common pangolin (Manis tricuspis) and (ii) the utility of our multilocus approach for taxonomic delineation and trade tracing of this heavily poached species. We sequenced 101 individuals for two mitochondrial DNA (mtDNA), two nuclear DNA and one Y-borne gene fragments (totalizing 2602 bp). We used a time-calibrated, Bayesian inference phylogenetic framework and conducted character-based, genetic and phylogenetic delineation of species hypotheses within African common pangolins. We identified six geographic lineages partitioned into western Africa, Ghana, the Dahomey Gap, western central Africa, Gabon and central Africa, all diverging during the Middle to Late Pleistocene. MtDNA (cytochrome b + control region) was the sole locus to provide diagnostic characters for each of the six lineages. Tree-based Bayesian delimitation methods using single- and multilocus approaches gave high support for 'species' level recognition of the six African common pangolin lineages. Although the diversification of African common pangolins occurred during Pleistocene cyclical glaciations, causative correlation with traditional rainforest refugia and riverine barriers in Africa was not straightforward. We conclude on the existence of six cryptic lineages within African common pangolins, which might be of major relevance for future conservation strategies. The high discriminative power of the mtDNA markers used in this study should allow an efficient molecular tracing of the regional origin of African common pangolin seizures.