Morphological variation in natural populations is a genomic test bed for studying the interface between molecular evolution and population genetics, but some of the most interesting questions involve non-model organisms that lack well annotated reference genomes. Many felid species exhibit polymorphism for melanism but the relative roles played by genetic drift, natural selection, and interspecies hybridization remain uncertain. We identify mutations of Agouti signaling protein (ASIP) or the Melanocortin 1 receptor (MC1R) as independent causes of melanism in three closely related South American species: the pampas cat (Leopardus colocolo), the kodkod (Leopardus guigna), and Geoffroy's cat (Leopardus geoffroyi). To assess population level variation in the regions surrounding the causative mutations we apply genomic resources from the domestic cat to carry out clone-based capture and targeted resequencing of 299 kb and 251 kb segments that contain ASIP and MC1R, respectively, from 54 individuals (13-21 per species), achieving enrichment of ~500-2500-fold and ~150x coverage. Our analysis points to unique evolutionary histories for each of the three species, with a strong selective sweep in the pampas cat, a distinctive but short melanism-specific haplotype in the Geoffroy's cat, and reduced nucleotide diversity for both ancestral and melanism-bearing chromosomes in the kodkod. These results reveal an important role for natural selection in a trait of longstanding interest to ecologists, geneticists, and the lay community, and provide a platform for comparative studies of morphological variation in other natural populations.

Pigmentation of the skin, hair, and eyes varies both within and between human populations. Identifying the genes and alleles underlying this variation has been the goal of many candidate gene and several genome-wide association studies (GWAS). Most GWAS for pigmentary traits to date have been based on subjective phenotypes using categorical scales. But skin, hair, and eye pigmentation vary continuously. Here, we seek to characterize quantitative variation in these traits objectively and accurately and to determine their genetic basis. Objective and quantitative measures of skin, hair, and eye color were made using reflectance or digital spectroscopy in Europeans from Ireland, Poland, Italy, and Portugal. A GWAS was conducted for the three quantitative pigmentation phenotypes in 176 women across 313,763 SNP loci, and replication of the most significant associations was attempted in a sample of 294 European men and women from the same countries. We find that the pigmentation phenotypes are highly stratified along axes of European genetic differentiation. The country of sampling explains approximately 35% of the variation in skin pigmentation, 31% of the variation in hair pigmentation, and 40% of the variation in eye pigmentation. All three quantitative phenotypes are correlated with each other. In our two-stage association study, we reproduce the association of rs1667394 at the OCA2/HERC2 locus with eye color but we do not identify new genetic determinants of skin and hair pigmentation supporting the lack of major genes affecting skin and hair color variation within Europe and suggesting that not only careful phenotyping but also larger cohorts are required to understand the genetic architecture of these complex quantitative traits. Interestingly, we also see that in each of these four populations, men are more lightly pigmented in the unexposed skin of the inner arm than women, a fact that is underappreciated and may vary across the world.

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3 × 10(-62), SLC24A5 P = 9.6 × 10(-9)) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4 × 10(-27), TYR P = 1.1 × 10(-9), APBA2[OCA2] P = 1.5 × 10(-8), SLC45A2 P = 6 × 10(-9)) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (~44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ~70% of the estimated heritability.

Insulin and leptin intracellular signaling pathways converge and act synergistically on the hypothalamic phosphatidylinositol-3-OH kinase/3-phosphoinositide-dependent protein kinase 1 (PDK1). However, little is known about whether PDK1 in agouti-related peptide (AGRP) neurons contributes to energy homeostasis. We generated AGRP neuron-specific PDK1 knockout (AGRPPdk1(-/-)) mice and mice with selective expression of transactivation-defective Foxo1 (Δ256Foxo1(AGRP)Pdk1(-/-)). The AGRPPdk1(-/-) mice showed reductions in food intake, body length, and body weight. The Δ256Foxo1(AGRP)Pdk1(-/-) mice showed increased body weight, food intake, and reduced locomotor activity. After four weeks of calorie-restricted feeding, oxygen consumption and locomotor activity were elevated in AGRPPdk1(-/-) mice and reduced in Δ256Foxo1(AGRP)Pdk1(-/-) mice. In vitro, ghrelin-induced changes in [Ca(2+)](i) and inhibition of ghrelin by leptin were significantly attenuated in AGRPPdk1(-/-) neurons compared to control neurons. However, ghrelin-induced [Ca(2+)](i) changes and leptin inhibition were restored in Δ256Foxo1(AGRP)Pdk1(-/-) mice. These results suggested that PDK1 and Foxo1 signaling pathways play important roles in the control of energy homeostasis through AGRP-independent mechanisms.

Dravet syndrome (DS) is a developmental and epileptic encephalopathy that results from mutations in the Nav1.1 sodium channel encoded by SCN1A. Most known DS-causing mutations are in coding regions of SCN1A, but we recently identified several disease-associated SCN1A mutations in intron 20 that are within or near to a cryptic and evolutionarily conserved "poison" exon, 20N, whose inclusion is predicted to lead to transcript degradation. However, it is not clear how these intron 20 variants alter SCN1A expression or DS pathophysiology in an organismal context, nor is it clear how exon 20N is regulated in a tissue-specific and developmental context. We address those questions here by generating an animal model of our index case, NM_006920.4(SCN1A):c.3969+2451G>C, using gene editing to create the orthologous mutation in laboratory mice. Scn1a heterozygous knock-in (+/KI) mice exhibited an ~50% reduction in brain Scn1a mRNA and Nav1.1 protein levels, together with characteristics observed in other DS mouse models, including premature mortality, seizures, and hyperactivity. In brain tissue from adult Scn1a +/+ animals, quantitative RT-PCR assays indicated that ~1% of Scn1a mRNA included exon 20N, while brain tissue from Scn1a +/KI mice exhibited an ~5-fold increase in the extent of exon 20N inclusion. We investigated the extent of exon 20N inclusion in brain during normal fetal development in RNA-seq data and discovered that levels of inclusion were ~70% at E14.5, declining progressively to ~10% postnatally. A similar pattern exists for the homologous sodium channel Nav1.6, encoded by Scn8a. For both genes, there is an inverse relationship between the level of functional transcript and the extent of poison exon inclusion. Taken together, our findings suggest that poison exon usage by Scn1a and Scn8a is a strategy to regulate channel expression during normal brain development, and that mutations recapitulating a fetal-like pattern of splicing cause reduced channel expression and epileptic encephalopathy.

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.

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.

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.

Melanocortin-1 receptor (MC1R) and its ligands, alpha-melanocyte stimulating hormone (alphaMSH) and agouti signaling protein (ASIP), regulate switching between eumelanin and pheomelanin synthesis in melanocytes. Here we investigated biological effects and signaling pathways of ASIP. Melan-a non agouti (a/a) mouse melanocytes produce mainly eumelanin, but ASIP combined with phenylthiourea and extra cysteine could induce over 200-fold increases in the pheomelanin to eumelanin ratio, and a tan-yellow color in pelletted cells. Moreover, ASIP-treated cells showed reduced proliferation and a melanoblast-like appearance, seen also in melanocyte lines from yellow (A(y)/a and Mc1r(e)/ Mc1r(e)) mice. However ASIP-YY, a C-terminal fragment of ASIP, induced neither biological nor pigmentary changes. As, like ASIP, ASIP-YY inhibited the cAMP rise induced by alphaMSH analog NDP-MSH, and reduced cAMP level without added MSH, the morphological changes and depigmentation seemed independent of cAMP signaling. Melanocytes genetically null for ASIP mediators attractin or mahogunin (Atrn(mg-3J/mg-3J) or Mgrn1(md-nc/md-nc)) also responded to both ASIP and ASIP-YY in cAMP level, while only ASIP altered their proliferation and (in part) shape. Thus, ASIP-MC1R signaling includes a cAMP-independent pathway through attractin and mahogunin, while the known cAMP-dependent component requires neither attractin nor mahogunin.

Nearly all neurodegenerative diseases are associated with abnormal accumulation of ubiquitin (Ub) conjugates within neuronal inclusion bodies. To directly test the hypothesis that depletion of cellular Ub is sufficient to cause neurodegeneration, we have disrupted Ubb, one of four genes that supply Ub in the mouse. Here, we report that loss of Ubb led to a progressive degenerative disorder affecting neurons within the arcuate nucleus of the hypothalamus. This neurodegenerative cytopathology was accompanied by impaired hypothalamic control of energy balance and adult-onset obesity. Ubb was highly expressed in vulnerable hypothalamic neurons and total Ub levels were selectively reduced in the hypothalamus of Ubb-null mice. These findings demonstrate that maintenance of adequate supplies of cellular Ub is essential for neuronal survival and establish that decreased Ub availability is sufficient to cause neuronal dysfunction and death.

Mutations that alter signaling of RAS/MAPK-family proteins give rise to a group of Mendelian diseases known as RASopathies. However, among RASopathies, the matrix of genotype-phenotype relationships is still incomplete, in part because there are many RAS-related proteins and in part because the phenotypic consequences may be variable and/or pleiotropic. Here, we describe a cohort of ten cases, drawn from six clinical sites and over 16,000 sequenced probands, with de novo protein-altering variation in RALA, a RAS-like small GTPase. All probands present with speech and motor delays, and most have intellectual disability, low weight, short stature, and facial dysmorphism. The observed rate of de novo RALA variants in affected probands is significantly higher (p = 4.93 x 10(-11)) than expected from the estimated random mutation rate. Further, all de novo variants described here affect residues within the GTP/GDP-binding region of RALA; in fact, six alleles arose at only two codons, Val25 and Lys128. The affected residues are highly conserved across both RAL- and RAS-family genes, are devoid of variation in large human population datasets, and several are homologous to positions at which disease-associated variants have been observed in other GTPase genes. We directly assayed GTP hydrolysis and RALA effector-protein binding of the observed variants, and found that all but one tested variant significantly reduced both activities compared to wild-type. The one exception, S157A, reduced GTP hydrolysis but significantly increased RALA-effector binding, an observation similar to that seen for oncogenic RAS variants. These results show the power of data sharing for the interpretation and analysis of rare variation, expand the spectrum of molecular causes of developmental disability to include RALA, and provide additional insight into the pathogenesis of human disease caused by mutations in small GTPases.

Coat color in Holstein dairy cattle is primarily controlled by the melanocortin 1 receptor (MC1R) gene, a central determinant of black (eumelanin) vs. red/brown pheomelanin synthesis across animal species. The major MC1R alleles in Holsteins are Dominant Black (MC1RD) and Recessive Red (MC1Re). A novel form of dominant red coat color was first observed in an animal born in 1980. The mutation underlying this phenotype was named Dominant Red and is epistatic to the constitutively activated MC1RD. Here we show that a missense mutation in the coatomer protein complex, subunit alpha (COPA), a gene with previously no known role in pigmentation synthesis, is completely associated with Dominant Red in Holstein dairy cattle. The mutation results in an arginine to cysteine substitution at an amino acid residue completely conserved across eukaryotes. Despite this high level of conservation we show that both heterozygotes and homozygotes are healthy and viable. Analysis of hair pigment composition shows that the Dominant Red phenotype is similar to the MC1R Recessive Red phenotype, although less effective at reducing eumelanin synthesis. RNA-seq data similarly show that Dominant Red animals achieve predominantly pheomelanin synthesis by downregulating genes normally required for eumelanin synthesis. COPA is a component of the coat protein I seven subunit complex that is involved with retrograde and cis-Golgi intracellular coated vesicle transport of both protein and RNA cargo. This suggests that Dominant Red may be caused by aberrant MC1R protein or mRNA trafficking within the highly compartmentalized melanocyte, mimicking the effect of the Recessive Red loss of function MC1R allele.

No abstract available

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.

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.

Clinically relevant secondary variants were identified in parents enrolled with a child with developmental delay and intellectual disability.

Many members of the animal kingdom display coat or skin color differences along their dorsoventral axis. To determine the mechanisms that control regional differences in pigmentation, we have studied how a classical mouse mutation, droopy ear (de(H)), affects dorsoventral skin characteristics, especially those under control of the Agouti gene. Mice carrying the Agouti allele black-and-tan (a(t)) normally have a sharp boundary between dorsal black hair and yellow ventral hair; the de(H) mutation raises the pigmentation boundary, producing an apparent dorsal-to-ventral transformation. We identify a 216 kb deletion in de(H) that removes all but the first exon of the Tbx15 gene, whose embryonic expression in developing mesenchyme correlates with pigmentary and skeletal malformations observed in de(H)/de(H) animals. Construction of a targeted allele of Tbx15 confirmed that the de(H) phenotype was caused by Tbx15 loss of function. Early embryonic expression of Tbx15 in dorsal mesenchyme is complementary to Agouti expression in ventral mesenchyme; in the absence of Tbx15, expression of Agouti in both embryos and postnatal animals is displaced dorsally. Transplantation experiments demonstrate that positional identity of the skin with regard to dorsoventral pigmentation differences is acquired by E12.5, which is shortly after early embryonic expression of Tbx15. Fate-mapping studies show that the dorsoventral pigmentation boundary is not in register with a previously identified dermal cell lineage boundary, but rather with the limb dorsoventral boundary. Embryonic expression of Tbx15 in dorsolateral mesenchyme provides an instructional cue required to establish the future positional identity of dorsal dermis. These findings represent a novel role for T-box gene action in embryonic development, identify a previously unappreciated aspect of dorsoventral patterning that is widely represented in furred mammals, and provide insight into the mechanisms that underlie region-specific differences in body morphology.

Developmental disabilities have diverse genetic causes that must be identified to facilitate precise diagnoses. We describe genomic data from 371 affected individuals, 309 of which were sequenced as proband-parent trios.

Preference for dietary fat vs. carbohydrate varies markedly across free-living individuals. It is recognized that food choice is under genetic and physiological regulation, and that the central melanocortin system is involved. However, how genetic and dietary factors interact to regulate relative macronutrient intake is not well understood.

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.