Search

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.

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.

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.

Live-cell imaging and genetic tools reveal a new way in which pigment cells communicate in zebrafish

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.

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.

During hibernation, animals cycle between torpor and arousal. These cycles involve dramatic but poorly understood mechanisms of dynamic physiological regulation at the level of gene expression. Each cycle, Brown Adipose Tissue (BAT) drives periodic arousal from torpor by generating essential heat. We applied digital transcriptome analysis to precisely timed samples to identify molecular pathways that underlie the intense activity cycles of hibernator BAT. A cohort of transcripts increased during torpor, paradoxical because transcription effectively ceases at these low temperatures. We show that this increase occurs not by elevated transcription but rather by enhanced stabilization associated with maintenance and/or extension of long poly(A) tails. Mathematical modeling further supports a temperature-sensitive mechanism to protect a subset of transcripts from ongoing bulk degradation instead of increased transcription. This subset was enriched in a C-rich motif and genes required for BAT activation, suggesting a model and mechanism to prioritize translation of key proteins for thermogenesis.

Human facial diversity is substantial, complex, and largely scientifically unexplained. We used spatially dense quasi-landmarks to measure face shape in population samples with mixed West African and European ancestry from three locations (United States, Brazil, and Cape Verde). Using bootstrapped response-based imputation modeling (BRIM), we uncover the relationships between facial variation and the effects of sex, genomic ancestry, and a subset of craniofacial candidate genes. The facial effects of these variables are summarized as response-based imputed predictor (RIP) variables, which are validated using self-reported sex, genomic ancestry, and observer-based facial ratings (femininity and proportional ancestry) and judgments (sex and population group). By jointly modeling sex, genomic ancestry, and genotype, the independent effects of particular alleles on facial features can be uncovered. Results on a set of 20 genes showing significant effects on facial features provide support for this approach as a novel means to identify genes affecting normal-range facial features and for approximating the appearance of a face from genetic markers.

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.

Prion diseases are rare but invariably fatal neurodegenerative disorders. They are associated with spongiform encephalopathy, a histopathology characterized by the presence of large, membrane-bound vacuolar structures in the neuropil of the brain. While the primary cause is recognized as conversion of the normal form of prion protein (PrP(C)) to a conformationally distinct, pathogenic form (PrP(Sc)), the cellular pathways and mechanisms that lead to spongiform change, neuronal dysfunction and death are not known. Mice lacking the Mahogunin Ring Finger 1 (MGRN1) E3 ubiquitin ligase develop spongiform encephalopathy by 9 months of age but do not become ill. In cell culture, PrP aberrantly present in the cytosol was reported to interact with and sequester MGRN1. This caused endo-lysosomal trafficking defects similar to those observed when Mgrn1 expression is knocked down, implicating disrupted MGRN1-dependent trafficking in the pathogenesis of prion disease. As these defects were rescued by over-expression of MGRN1, we investigated whether reduced or elevated Mgrn1 expression influences the onset, progression or pathology of disease in mice inoculated with PrP(Sc). No differences were observed, indicating that disruption of MGRN1-dependent pathways does not play a significant role in the pathogenesis of transmissible spongiform encephalopathy.

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.

The occurrence of melanism (darkening of the background coloration) is documented in 13 felid species, in some cases reaching high frequencies at the population level. Recent analyses have indicated that it arose multiple times in the Felidae, with three different species exhibiting unique mutations associated with this trait. The causative mutations in the remaining species have so far not been identified, precluding a broader assessment of the evolutionary dynamics of melanism in the Felidae. Among these, the leopard (Panthera pardus) is a particularly important target for research, given the iconic status of the 'black panther' and the extremely high frequency of melanism observed in some Asian populations. Another felid species from the same region, the Asian golden cat (Pardofelis temminckii), also exhibits frequent records of melanism in some areas. We have sequenced the coding region of the Agouti Signaling Protein (ASIP) gene in multiple leopard and Asian golden cat individuals, and identified distinct mutations strongly associated with melanism in each of them. The single nucleotide polymorphism (SNP) detected among the P. pardus individuals was caused by a nonsense mutation predicted to completely ablate ASIP function. A different SNP was identified in P. temminckii, causing a predicted amino acid change that should also induce loss of function. Our results reveal two additional cases of species-specific mutations implicated in melanism in the Felidae, and indicate that ASIP mutations may play an important role in naturally-occurring coloration polymorphism.

For most of the world, human genome structure at a population level is shaped by interplay between ancient geographic isolation and more recent demographic shifts, factors that are captured by the concepts of biogeographic ancestry and admixture, respectively. The ancestry of non-admixed individuals can often be traced to a specific population in a precise region, but current approaches for studying admixed individuals generally yield coarse information in which genome ancestry proportions are identified according to continent of origin. Here we introduce a new analytic strategy for this problem that allows fine-grained characterization of admixed individuals with respect to both geographic and genomic coordinates. Ancestry segments from different continents, identified with a probabilistic model, are used to construct and study "virtual genomes" of admixed individuals. We apply this approach to a cohort of 492 parent-offspring trios from Mexico City. The relative contributions from the three continental-level ancestral populations-Africa, Europe, and America-vary substantially between individuals, and the distribution of haplotype block length suggests an admixing time of 10-15 generations. The European and Indigenous American virtual genomes of each Mexican individual can be traced to precise regions within each continent, and they reveal a gradient of Amerindian ancestry between indigenous people of southwestern Mexico and Mayans of the Yucatan Peninsula. This contrasts sharply with the African roots of African Americans, which have been characterized by a uniform mixing of multiple West African populations. We also use the virtual European and Indigenous American genomes to search for the signatures of selection in the ancestral populations, and we identify previously known targets of selection in other populations, as well as new candidate loci. The ability to infer precise ancestral components of admixed genomes will facilitate studies of disease-related phenotypes and will allow new insight into the adaptive and demographic history of indigenous people.

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.

Agouti-related protein encodes a neuropeptide that stimulates food intake. Agrp expression in the brain is restricted to neurons in the arcuate nucleus of the hypothalamus and is elevated by states of negative energy balance. The molecular mechanisms underlying Agrp regulation, however, remain poorly defined. Using a combination of transgenic and comparative sequence analysis, we have previously identified a 760 bp conserved region upstream of Agrp which contains STAT binding elements that participate in Agrp transcriptional regulation. In this study, we attempt to improve the specificity for detecting conserved elements in this region by comparing genomic sequences from 10 mammalian species. Our analysis reveals a symmetrical organization of conserved sequences upstream of Agrp, which cluster into two inverted repeat elements. Conserved sequences within these elements suggest a role for homeodomain proteins in the regulation of Agrp and provide additional targets for functional evaluation.

Genetic variation at the melanocortin 1 receptor (MC1R) is an important risk factor for developing ultraviolet (UV) radiation-induced skin cancer, the most common form of cancer in humans. The underlying mechanisms by which the MC1R defends against UV-induced skin cancer are not known. We used neonatal mouse skin (which, like human skin, contains a mixture of melanocytes and keratinocytes) to study how pigment cells and Mc1r genotype affect the genome-level response to UV radiation. Animals without viable melanocytes (Kit(W-v)/Kit(W-v)) or animals lacking a functional Mc1r (Mc1r(e)/Mc1r(e)) were exposed to sunburn-level doses of UVB radiation, and the patterns of large-scale gene expression in the basal epidermis were compared to each other and to nonmutant animals. Our analysis revealed discrete Kit- and Mc1r-dependent UVB transcriptional responses in the basal epidermis. The Kit-dependent UVB response was characterized largely by an enrichment of oxidative and endoplasmic reticulum stress genes, highlighting a distinctive role for pigmented melanocytes in mediating antioxidant defenses against genotoxic stresses within the basal epidermal environment. By contrast, the Mc1r-dependent UVB response contained an abundance of genes associated with regulating the cell cycle and oncogenesis. To test the clinical relevance of these observations, we analyzed publicly available data sets for primary melanoma and melanoma metastases and found that the set of genes specific for the Mc1r-dependent UVB response was able to differentiate between different clinical subtypes. Our analysis also revealed that the classes of genes induced by UVB differ from those repressed by UVB with regard to their biological functions, their overall number, and their size. The findings described here offer new insights into the transcriptional nature of the UV response in the skin and provide a molecular framework for the underlying mechanisms by which melanocytes and the Mc1r independently mediate and afford protection against UV radiation.

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.

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.