Unraveling the genomic diversity and admixture history of captive tigers in the United States.

Genomic studies of endangered species have primarily focused on describing diversity patterns and resolving phylogenetic relationships, with the overarching goal of informing conservation efforts. However, few studies have investigated genomic diversity housed in captive populations. For tigers (Panthera tigris), captive individuals vastly outnumber those in the wild, but their diversity remains largely unexplored. Privately owned captive tiger populations have remained an enigma in the conservation community, with some believing that these individuals are severely inbred, while others believe they may be a source of now-extinct diversity. Here, we present a large-scale genetic study of the private (non-zoo) captive tiger population in the United States, also known as "Generic" tigers. We find that the Generic tiger population has an admixture fingerprint comprising all six extant wild tiger subspecies. Of the 138 Generic individuals sequenced for the purpose of this study, no individual had ancestry from only one subspecies. We show that the Generic tiger population has a comparable amount of genetic diversity relative to most wild subspecies, few private variants, and fewer deleterious mutations. We observe inbreeding coefficients similar to wild populations, although there are some individuals within both the Generic and wild populations that are substantially inbred. Additionally, we develop a reference panel for tigers that can be used with imputation to accurately distinguish individuals and assign ancestry with ultralow coverage (0.25×) data. By providing a cost-effective alternative to whole-genome sequencing (WGS), the reference panel provides a resource to assist in tiger conservation efforts for both ex- and in situ populations.

High frequency of an otherwise rare phenotype in a small and isolated tiger population.

Most endangered species exist today in small populations, many of which are isolated. Evolution in such populations is largely governed by genetic drift. Empirical evidence for drift affecting striking phenotypes based on substantial genetic data are rare. Approximately 37% of tigers (Panthera tigris) in the Similipal Tiger Reserve (in eastern India) are pseudomelanistic, characterized by wide, merged stripes. Camera trap data across the tiger range revealed the presence of pseudomelanistic tigers only in Similipal. We investigated the genetic basis for pseudomelanism and examined the role of drift in driving this phenotype's frequency. Whole-genome data and pedigree-based association analyses from captive tigers revealed that pseudomelanism cosegregates with a conserved and functionally important coding alteration in Transmembrane Aminopeptidase Q (Taqpep), a gene responsible for similar traits in other felid species. Noninvasive sampling of tigers revealed a high frequency of the Taqpep p.H454Y mutation in Similipal (12 individuals, allele frequency = 0.58) and absence from all other tiger populations (395 individuals). Population genetic analyses confirmed few (minimal number) tigers in Similipal, and its genetic isolation, with poor geneflow. Pairwise FST (0.33) at the mutation site was high but not an outlier. Similipal tigers had low diversity at 81 single nucleotide polymorphisms (mean heterozygosity = 0.28, SD = 0.27). Simulations were consistent with founding events and drift as possible drivers for the observed stark difference of allele frequency. Our results highlight the role of stochastic processes in the evolution of rare phenotypes. We highlight an unusual evolutionary trajectory in a small and isolated population of an endangered species.

PEA15 loss of function and defective cerebral development in the domestic cat.

Cerebral cortical size and organization are critical features of neurodevelopment and human evolution, for which genetic investigation in model organisms can provide insight into developmental mechanisms and the causes of cerebral malformations. However, some abnormalities in cerebral cortical proliferation and folding are challenging to study in laboratory mice due to the absence of gyri and sulci in rodents. We report an autosomal recessive allele in domestic cats associated with impaired cerebral cortical expansion and folding, giving rise to a smooth, lissencephalic brain, and that appears to be caused by homozygosity for a frameshift in PEA15 (phosphoprotein expressed in astrocytes-15). Notably, previous studies of a Pea15 targeted mutation in mice did not reveal structural brain abnormalities. Affected cats, however, present with a non-progressive hypermetric gait and tremors, develop dissociative behavioral defects and aggression with age, and exhibit profound malformation of the cerebrum, with a 45% average decrease in overall brain weight, and reduction or absence of the ectosylvian, sylvian and anterior cingulate gyrus. Histologically, the cerebral cortical layers are disorganized, there is substantial loss of white matter in tracts such as the corona radiata and internal capsule, but the cerebellum is relatively spared. RNA-seq and immunohistochemical analysis reveal astrocytosis. Fibroblasts cultured from affected cats exhibit increased TNFα-mediated apoptosis, and increased FGFb-induced proliferation, consistent with previous studies implicating PEA15 as an intracellular adapter protein, and suggesting an underlying pathophysiology in which increased death of neurons accompanied by increased proliferation of astrocytes gives rise to abnormal organization of neuronal layers and loss of white matter. Taken together, our work points to a new role for PEA15 in development of a complex cerebral cortex that is only apparent in gyrencephalic species.

Recent Evolutionary History of Tigers Highlights Contrasting Roles of Genetic Drift and Selection.

Species conservation can be improved by knowledge of evolutionary and genetic history. Tigers are among the most charismatic of endangered species and garner significant conservation attention. However, their evolutionary history and genomic variation remain poorly known, especially for Indian tigers. With 70% of the world's wild tigers living in India, such knowledge is critical. We re-sequenced 65 individual tiger genomes representing most extant subspecies with a specific focus on tigers from India. As suggested by earlier studies, we found strong genetic differentiation between the putative tiger subspecies. Despite high total genomic diversity in India, individual tigers host longer runs of homozygosity, potentially suggesting recent inbreeding or founding events, possibly due to small and fragmented protected areas. We suggest the impacts of ongoing connectivity loss on inbreeding and persistence of Indian tigers be closely monitored. Surprisingly, demographic models suggest recent divergence (within the last 20,000 years) between subspecies and strong population bottlenecks. Amur tiger genomes revealed the strongest signals of selection related to metabolic adaptation to cold, whereas Sumatran tigers show evidence of weak selection for genes involved in body size regulation. We recommend detailed investigation of local adaptation in Amur and Sumatran tigers prior to initiating genetic rescue.

Population structure, inbreeding and stripe pattern abnormalities in plains zebras.

One of the most iconic wild equids, the plains zebra occupies a broad region of sub-Saharan Africa and exhibits a wide range of phenotypic diversity in stripe patterns that have been used to classify multiple subspecies. After decades of relative stability, albeit with a loss of at least one recognized subspecies, the total population of plains zebras has undergone an approximate 25% decline since 2002. Individuals with abnormal stripe patterns have been recognized in recent years but the extent to which their appearance is related to demography and/or genetics is unclear. Investigating population genetic health and genetic structure are essential for developing effective strategies for plains zebra conservation. We collected DNA from 140 plains zebra, including seven with abnormal stripe patterns, from nine locations across the range of plains zebra, and analyzed data from restriction site-associated and whole genome sequencing (RAD-seq, WGS) libraries to better understand the relationships between population structure, genetic diversity, inbreeding, and abnormal phenotypes. We found that genetic structure did not coincide with described subspecific variation, but did distinguish geographic regions in which anthropogenic habitat fragmentation is associated with reduced gene flow and increased evidence of inbreeding, especially in certain parts of East Africa. Further, zebras with abnormal striping exhibited increased levels of inbreeding relative to normally striped individuals from the same populations. Our results point to a genetic cause of stripe pattern abnormalities, and dramatic evidence of the consequences of habitat fragmentation.

Long live the king: chromosome-level assembly of the lion (Panthera leo) using linked-read, Hi-C, and long-read data.

BACKGROUND: The lion (Panthera leo) is one of the most popular and iconic feline species on the planet, yet in spite of its popularity, the last century has seen massive declines for lion populations worldwide. Genomic resources for endangered species represent an important way forward for the field of conservation, enabling high-resolution studies of demography, disease, and population dynamics. Here, we present a chromosome-level assembly from a captive African lion from the Exotic Feline Rescue Center (Center Point, IN) as a resource for current and subsequent genetic work of the sole social species of the Panthera clade. RESULTS: Our assembly is composed of 10x Genomics Chromium data, Dovetail Hi-C, and Oxford Nanopore long-read data. Synteny is highly conserved between the lion, other Panthera genomes, and the domestic cat. We find variability in the length of runs of homozygosity across lion genomes, indicating contrasting histories of recent and possibly intense inbreeding and bottleneck events. Demographic analyses reveal similar ancient histories across all individuals during the Pleistocene except the Asiatic lion, which shows a more rapid decline in population size. We show a substantial influence on the reference genome choice in the inference of demographic history and heterozygosity. CONCLUSIONS: We demonstrate that the choice of reference genome is important when comparing heterozygosity estimates across species and those inferred from different references should not be compared to each other. In addition, estimates of heterozygosity or the amount or length of runs of homozygosity should not be taken as reflective of a species, as these can differ substantially among individuals. This high-quality genome will greatly aid in the continuing research and conservation efforts for the lion, which is rapidly moving towards becoming a species in danger of extinction.

Electrostatic Similarity Analysis of Human ß-Defensin Binding in the Melanocortin System.

The ß-defensins are a class of small cationic proteins that serve as components of numerous systems in vertebrate biology, including the immune and melanocortin systems. Human ß-defensin 3 (HBD3), which is produced in the skin, has been found to bind to melanocortin receptors 1 and 4 through complementary electrostatics, a unique mechanism of ligand-receptor interaction. This finding indicates that electrostatics alone, and not specific amino acid contact points, could be sufficient for function in this ligand-receptor system, and further suggests that other small peptide ligands could interact with these receptors in a similar fashion. Here, we conducted molecular-similarity analyses and functional studies of additional members of the human ß-defensin family, examining their potential as ligands of melanocortin-1 receptor, through selection based on their electrostatic similarity to HBD3. Using Poisson-Boltzmann electrostatic calculations and molecular-similarity analysis, we identified members of the human ß-defensin family that are both similar and dissimilar to HBD3 in terms of electrostatic potential. Synthesis and functional testing of a subset of these ß-defensins showed that peptides with an HBD3-like electrostatic character bound to melanocortin receptors with high affinity, whereas those that were anticorrelated to HBD3 showed no binding affinity. These findings expand on the central role of electrostatics in the control of this ligand-receptor system and further demonstrate the utility of employing molecular-similarity analysis. Additionally, we identified several new potential ligands of melanocortin-1 receptor, which may have implications for our understanding of the role defensins play in melanocortin physiology.

Coordinated regulation of hepatic energy stores by leptin and hypothalamic agouti-related protein.

Like obesity, prolonged food deprivation induces severe hepatic steatosis; however, the functional significance of this phenomenon is not well understood. In this study, we show that the fall in plasma leptin concentration during fasting is required for the development of hepatic steatosis in mice. Removal of leptin receptors from AGRP neurons diminishes fasting-induced hepatic steatosis. Furthermore, the suppressive effects of leptin on fasting-induced hepatic steatosis are absent in mice lacking the gene encoding agouti-related protein (Agrp), suggesting that this function of leptin is mediated by AGRP. Prolonged fasting leads to suppression of hepatic sympathetic activity, increased expression of acyl CoA:diacylglycerol acyltransferase-2 in the liver, and elevation of hepatic triglyceride content and all of these effects are blunted in the absence of AGRP. AGRP deficiency, despite having no effects on feeding or body adiposity in the free-fed state, impairs triglyceride and ketone body release from the liver during prolonged fasting. Furthermore, reducing CNS Agrp expression in wild-type mice by RNAi protected against the development of hepatic steatosis not only during starvation, but also in response to consumption of a high-fat diet. These findings identify the leptin-AGRP circuit as a critical modulator of hepatic triglyceride stores in starvation and suggest a vital role for this circuit in sustaining the supply of energy from the liver to extrahepatic tissues during periods of prolonged food deprivation.

Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

Genetic architecture and evolution of color variation in American black bears.

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Developmental genetics of color pattern establishment in cats.

Intricate color patterns are a defining aspect of morphological diversity in the Felidae. We applied morphological and single-cell gene expression analysis to fetal skin of domestic cats to identify when, where, and how, during fetal development, felid color patterns are established. Early in development, we identify stripe-like alterations in epidermal thickness preceded by a gene expression pre-pattern. The secreted Wnt inhibitor encoded by Dickkopf 4 plays a central role in this process, and is mutated in cats with the Ticked pattern type. Our results bring molecular understanding to how the leopard got its spots, suggest that similar mechanisms underlie periodic color pattern and periodic hair follicle spacing, and identify targets for diverse pattern variation in other mammals.

Epigenetic models developed for plains zebras predict age in domestic horses and endangered equids.

Effective conservation and management of threatened wildlife populations require an accurate assessment of age structure to estimate demographic trends and population viability. Epigenetic aging models are promising developments because they estimate individual age with high accuracy, accurately predict age in related species, and do not require invasive sampling or intensive long-term studies. Using blood and biopsy samples from known age plains zebras (Equus quagga), we model epigenetic aging using two approaches: the epigenetic clock (EC) and the epigenetic pacemaker (EPM). The plains zebra EC has the potential for broad application within the genus Equus given that five of the seven extant wild species of the genus are threatened. We test the EC's ability to predict age in sister taxa, including two endangered species and the more distantly related domestic horse, demonstrating high accuracy in all cases. By comparing chronological and estimated age in plains zebras, we investigate age acceleration as a proxy of health status. An interaction between chronological age and inbreeding is associated with age acceleration estimated by the EPM, suggesting a cumulative effect of inbreeding on biological aging throughout life.

Dog colour patterns explained by modular promoters of ancient canid origin.

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Genetics of Sex-linked yellow in the Syrian hamster.

Alternating patches of black and yellow pigment are a ubiquitous feature of mammalian color variation that contributes to camouflage, species recognition, and morphologic diversity. X-linked determinants of this pattern--recognized by variegation in females but not in males--have been described in the domestic cat as Orange, and in the Syrian hamster as Sex-linked yellow (Sly), but are curiously absent from other vertebrate species. Using a comparative genomic approach, we develop molecular markers and a linkage map for the euchromatic region of the Syrian hamster X chromosome that places Sly in a region homologous to the centromere-proximal region of human Xp. Comparison to analogous work carried out for Orange in domestic cats indicates, surprisingly, that the cat and hamster mutations lie in nonhomologous regions of the X chromosome. We also identify the molecular cause of recessively inherited black coat color in hamsters (historically referred to as nonagouti) as a Cys115Tyr mutation in the Agouti gene. Animals doubly mutant for Sly and nonagouti exhibit a Sly phenotype. Our results indicate that Sly represents a melanocortin pathway component that acts similarly to, but is genetically distinct from, Mc1r and that has implications for understanding both the evolutionary history and the mutational mechanisms of pigment-type switching.

PI3K integrates the action of insulin and leptin on hypothalamic neurons.

Central control of energy balance depends on the ability of proopiomelanocortin (POMC) or agouti-related protein (Agrp) hypothalamic neurons to sense and respond to changes in peripheral energy stores. Leptin and insulin have been implicated as circulating indicators of adiposity, but it is not clear how changes in their levels are perceived or integrated by individual neuronal subtypes. We developed mice in which a fluorescent reporter for PI3K activity is targeted to either Agrp or POMC neurons and used 2-photon microscopy to measure dynamic regulation of PI3K by insulin and leptin in brain slices. We show that leptin and insulin act in parallel to stimulate PI3K in POMC neurons but in opposite ways on Agrp neurons. These results suggest a new view of hypothalamic circuitry, in which the effects of leptin and insulin are integrated by anorexigenic but not by orexigenic neurons.

Effects of hypothalamic neurodegeneration on energy balance.

Normal aging in humans and rodents is accompanied by a progressive increase in adiposity. To investigate the role of hypothalamic neuronal circuits in this process, we used a Cre-lox strategy to create mice with specific and progressive degeneration of hypothalamic neurons that express agouti-related protein (Agrp) or proopiomelanocortin (Pomc), neuropeptides that promote positive or negative energy balance, respectively, through their opposing effects on melanocortin receptor signaling. In previous studies, Pomc mutant mice became obese, but Agrp mutant mice were surprisingly normal, suggesting potential compensation by neuronal circuits or genetic redundancy. Here we find that Pomc-ablation mice develop obesity similar to that described for Pomc knockout mice, but also exhibit defects in compensatory hyperphagia similar to what occurs during normal aging. Agrp-ablation female mice exhibit reduced adiposity with normal compensatory hyperphagia, while animals ablated for both Pomc and Agrp neurons exhibit an additive interaction phenotype. These findings provide new insight into the roles of hypothalamic neurons in energy balance regulation, and provide a model for understanding defects in human energy balance associated with neurodegeneration and aging.

Inactivation of signal transducer and activator of transcription 3 in proopiomelanocortin (Pomc) neurons causes decreased pomc expression, mild obesity, and defects in compensatory refeeding.

Leptin is an adipocyte-derived hormone that signals body energy status to the brain by acting on multiple neuronal subgroups in the hypothalamus, including those that express proopiomelanocortin (Pomc) and agouti-related protein (Agrp). Signal transducer and activator of transcription 3 (Stat3) is an important intracellular signaling molecule activated by leptin, and previous studies have shown that mice carrying a mutated leptin receptor that abolished Stat3 binding are grossly obese. To determine the extent to which Stat3 signaling in Pomc neurons was responsible for these effects, we constructed Pomc-specific Stat3 mutants using a Cre recombinase transgene driven by the Pomc promoter. We find that Pomc expression is diminished in the mutant mice, suggesting that Stat3 is required for Pomc transcription. Pomc-specific Stat3 female mutant mice exhibit a 2-fold increase in fat pad mass but only a slight increase in total body weight. Mutant mice remain responsive to leptin-induced hypophagia and are not hypersensitive to a high-fat diet; however, mutant mice fail to mount a normal compensatory refeeding response. These results demonstrate a requirement for Stat3 in transcriptional regulation of Pomc but indicate that this circuit is only one of several components that underlie the neuronal response to leptin and the role of Stat3 in that response.

Signal transducer and activator of transcription (stat) binding sites but not stat3 are required for fasting-induced transcription of agouti-related protein messenger ribonucleic acid.

Energy homeostasis depends on the regulation of hypothalamic neurons by leptin, an adipocyte hormone whose circulating levels communicate body energy stores. Leptin activates the transcription factor signal transducer and activator of transcription 3 (Stat3) in hypothalamic neurons, including neuronal subtypes producing Agouti-related protein (Agrp), a neuropeptide that stimulates feeding. Previous studies have suggested a model in which high levels of Agrp transcription during fasting represent a default state that is actively repressed by phospho-Stat3 induced by leptin signaling in the fed state. We identify putative Stat3 binding elements in the Agrp promoter that have been highly conserved during vertebrate evolution. Using a reporter assay in transgenic mice that faithfully recapitulates normal regulation of Agrp, we show that these sites are required, but in a way opposite to that predicted by the existing model: mutation of the sites leads to a default state characterized by a low level of Agrp transcription and insensitivity to fasting. We also find that removing activatable Stat3 from Agrp neurons has no detectable effect on steady-state levels of Agrp mRNA in the fed or fasted state. These results suggest a new model for transcriptional regulation of orexigenic neuropeptides in which the default level of expression is low in the fed state, and transcriptional activation in response to fasting is mediated by factors other than Stat3.