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.

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.

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.

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.

Pangolin genomes and the evolution of mammalian scales and immunity.

Pangolins, unique mammals with scales over most of their body, no teeth, poor vision, and an acute olfactory system, comprise the only placental order (Pholidota) without a whole-genome map. To investigate pangolin biology and evolution, we developed genome assemblies of the Malayan (Manis javanica) and Chinese (M. pentadactyla) pangolins. Strikingly, we found that interferon epsilon (IFNE), exclusively expressed in epithelial cells and important in skin and mucosal immunity, is pseudogenized in all African and Asian pangolin species that we examined, perhaps impacting resistance to infection. We propose that scale development was an innovation that provided protection against injuries or stress and reduced pangolin vulnerability to infection. Further evidence of specialized adaptations was evident from positively selected genes involving immunity-related pathways, inflammation, energy storage and metabolism, muscular and nervous systems, and scale/hair development. Olfactory receptor gene families are significantly expanded in pangolins, reflecting their well-developed olfaction system. This study provides insights into mammalian adaptation and functional diversification, new research tools and questions, and perhaps a new natural IFNE-deficient animal model for studying mammalian immunity.

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.

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.

A magnetic protein biocompass.

The notion that animals can detect the Earth's magnetic field was once ridiculed, but is now well established. Yet the biological nature of such magnetosensing phenomenon remains unknown. Here, we report a putative magnetic receptor (Drosophila CG8198, here named MagR) and a multimeric magnetosensing rod-like protein complex, identified by theoretical postulation and genome-wide screening, and validated with cellular, biochemical, structural and biophysical methods. The magnetosensing complex consists of the identified putative magnetoreceptor and known magnetoreception-related photoreceptor cryptochromes (Cry), has the attributes of both Cry- and iron-based systems, and exhibits spontaneous alignment in magnetic fields, including that of the Earth. Such a protein complex may form the basis of magnetoreception in animals, and may lead to applications across multiple fields.

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.

Genetic ancestry of the extinct Javan and Bali tigers.

The Bali (Panthera tigris balica) and Javan (P. t. sondaica) tigers are recognized as distinct tiger subspecies that went extinct in the 1940s and 1980s, respectively. Yet their genetic ancestry and taxonomic status remain controversial. Following ancient DNA procedures, we generated concatenated 1750bp mtDNA sequences from 23 museum samples including 11 voucher specimens from Java and Bali and compared these to diagnostic mtDNA sequences from 122 specimens of living tiger subspecies and the extinct Caspian tiger. The results revealed a close genetic affinity of the 3 groups from the Sunda Islands (Bali, Javan, and Sumatran tigers P. t. sumatrae). Bali and Javan mtDNA haplotypes differ from Sumatran haplotypes by 1-2 nucleotides, and the 3 island populations define a monophyletic assemblage distinctive and equidistant from other mainland subspecies. Despite this close phylogenetic relationship, no mtDNA haplotype was shared between Sumatran and Javan/Bali tigers, indicating little or no matrilineal gene flow among the islands after they were colonized. The close phylogenetic relationship among Sunda tiger subspecies suggests either recent colonization across the islands, or else a once continuous tiger population that had subsequently isolated into different island subspecies. This supports the hypothesis that the Sumatran tiger is the closest living relative to the extinct Javan and Bali tigers.

Sympatric Asian felid phylogeography reveals a major Indochinese-Sundaic divergence.

The dynamic geological and climatological history of Southeast Asia has spawned a complex array of ecosystems and 12 of the 37 known cat species, making it the most felid-rich region in the world. To examine the evolutionary histories of these poorly studied fauna, we compared phylogeography of six species (leopard cat Prionailurus bengalensis, fishing cat P. viverrinus, Asiatic golden cat Pardofelis temminckii, marbled cat P. marmorata, tiger Panthera tigris and leopard P. pardus) by sequencing over 5 kb of DNA each from 445 specimens at multiple loci of mtDNA, Y and X chromosomes. All species except the leopard displayed significant phylogenetic partitions between Indochina and Sundaland, with the central Thai-Malay Peninsula serving as the biogeographic boundary. Concordant mtDNA and nuclear DNA genealogies revealed deep Indochinese-Sundaic divergences around 2 MYA in both P. bengalensis and P. marmorata comparable to previously described interspecific distances within Felidae. The divergence coincided with serial sea level rises during the late Pliocene and early Pleistocene, and was probably reinforced by repeated isolation events associated with environmental changes throughout the Pleistocene. Indochinese-Sundaic differentiations within P. tigris and P. temminckii were more recent at 72-108 and 250-1570 kya, respectively. Overall, these results illuminate unexpected, deep vicariance events in Southeast Asian felids and provide compelling evidence of species-level distinction between the Indochinese and Sundaic populations in the leopard cat and marbled cat. Broader sampling and further molecular and morphometric analyses of these species will be instrumental in defining conservation units and effectively preserving Southeast Asian biodiversity.

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.

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.

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.

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.

Subspecies genetic assignments of worldwide captive tigers increase conservation value of captive populations.

Tigers (Panthera tigris) are disappearing rapidly from the wild, from over 100,000 in the 1900s to as few as 3000. Javan (P.t. sondaica), Bali (P.t. balica), and Caspian (P.t. virgata) subspecies are extinct, whereas the South China tiger (P.t. amoyensis) persists only in zoos. By contrast, captive tigers are flourishing, with 15,000-20,000 individuals worldwide, outnumbering their wild relatives five to seven times. We assessed subspecies genetic ancestry of 105 captive tigers from 14 countries and regions by using Bayesian analysis and diagnostic genetic markers defined by a prior analysis of 134 voucher tigers of significant genetic distinctiveness. We assigned 49 tigers to one of five subspecies (Bengal P.t. tigris, Sumatran P.t. sumatrae, Indochinese P.t. corbetti, Amur P.t. altaica, and Malayan P.t. jacksoni tigers) and determined 52 had admixed subspecies origins. The tested captive tigers retain appreciable genomic diversity unobserved in their wild counterparts, perhaps a consequence of large population size, century-long introduction of new founders, and managed-breeding strategies to retain genetic variability. Assessment of verified subspecies ancestry offers a powerful tool that, if applied to tigers of uncertain background, may considerably increase the number of purebred tigers suitable for conservation management.

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.

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.