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.

Medical and psychosocial outcomes of state-funded population genomic screening

As the uptake of population screening expands, assessment of medical and psychosocial outcomes is needed. Through the Alabama Genomic Health Initiative (AGHI), a state-funded genomic research program, individuals received screening for pathogenic or likely pathogenic variants in 59 actionable genes via genotyping. Of the 3874 eligible participants that received screening results, 858 (22%) responded to an outcomes survey. The most commonly reported motivation for seeking testing through AGHI was contribution to genetic research (64%). Participants with positive results reported a higher median number of planned actions (median = 5) due to AGHI results as compared to negative results (median = 3). Interviews were conducted with survey participants with positive screening results. As determined by certified genetic counselors, 50% of interviewees took appropriate medical action based on their result. There were no negative or harmful actions taken. These findings indicate population genomic screening of an unselected adult population is feasible, is not harmful, and may have positive outcomes on participants now and in the future; however, further research is needed in order to assess clinical utility.

Parents' Perspectives on the Utility of Genomic Sequencing in the Neonatal Intensive Care Unit

BACKGROUND

It is critical to understand the wide-ranging clinical and non-clinical effects of genome sequencing (GS) for parents in the NICU context. We assessed parents' experiences with GS as a first-line diagnostic tool for infants with suspected genetic conditions in the NICU.

METHODS

Parents of newborns (N = 62) suspected of having a genetic condition were recruited across five hospitals in the southeast United States as part of the SouthSeq study. Semi-structured interviews (N = 78) were conducted after parents received their child's sequencing result (positive, negative, or variants of unknown significance). Thematic analysis was performed on all interviews.

RESULTS

Key themes included that (1) GS in infancy is important for reproductive decision making, preparing for the child's future care, ending the diagnostic odyssey, and sharing results with care providers; (2) the timing of disclosure was acceptable for most parents, although many reported the NICU environment was overwhelming; and (3) parents deny that receiving GS results during infancy exacerbated parent-infant bonding, and reported variable impact on their feelings of guilt.

CONCLUSION

Parents reported that GS during the neonatal period was useful because it provided a "backbone" for their child's care. Parents did not consistently endorse negative impacts like interference with parent-infant bonding.

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.

A gene-diet interaction controlling relative intake of dietary carbohydrates and fats

OBJECTIVE

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.

METHODS

We investigated how the choice for food rich in carbohydrate vs. fat is influenced by dietary cholesterol availability and agouti-related protein (AGRP), the orexigenic component of the central melanocortin system. We assessed how macronutrient intake and different metabolic parameters correlate with plasma AGRP in a cohort of obese humans. We also examined how both dietary cholesterol levels and inhibiting de novo cholesterol synthesis affect carbohydrate and fat intake in mice, and how dietary cholesterol deficiency during the postnatal period impacts macronutrient intake patterns in adulthood.

RESULTS

In obese human subjects, plasma levels of AGRP correlated inversely with consumption of carbohydrates over fats. Moreover, AgRP-deficient mice preferred to consume more calories from carbohydrates than fats, more so when each diet lacked cholesterol. Intriguingly, inhibiting cholesterol biosynthesis (simvastatin) promoted carbohydrate intake at the expense of fat without altering total caloric consumption, an effect that was remarkably absent in AgRP-deficient mice. Finally, feeding lactating C57BL/6 dams and pups a cholesterol-free diet prior to weaning led the offspring to prefer fats over carbohydrates as adults, indicating that altered cholesterol metabolism early in life programs adaptive changes to macronutrient intake.

CONCLUSIONS

Together, our study illustrates a specific gene-diet interaction in modulating food choice.

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.

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.

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.

A state-based approach to genomics for rare disease and population screening

PURPOSE

The Alabama Genomic Health Initiative (AGHI) is a state-funded effort to provide genomic testing. AGHI engages two distinct cohorts across the state of Alabama. One cohort includes children and adults with undiagnosed rare disease; a second includes an unselected adult population. Here we describe findings from the first 176 rare disease and 5369 population cohort AGHI participants.

METHODS

AGHI participants enroll in one of two arms of a research protocol that provides access to genomic testing results and biobank participation. Rare disease cohort participants receive genome sequencing to identify primary and secondary findings. Population cohort participants receive genotyping to identify pathogenic and likely pathogenic variants for actionable conditions.

RESULTS

Within the rare disease cohort, genome sequencing identified likely pathogenic or pathogenic variation in 20% of affected individuals. Within the population cohort, 1.5% of individuals received a positive genotyping result. The rate of genotyping results corroborated by reported personal or family history varied by gene.

CONCLUSIONS

AGHI demonstrates the ability to provide useful health information in two contexts: rare undiagnosed disease and population screening. This utility should motivate continued exploration of ways in which emerging genomic technologies might benefit broad populations.

Aberrant regulation of a poison exon caused by a non-coding variant in a mouse model of Scn1a-associated epileptic encephalopathy

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 Therapeutic Odyssey: Positioning Genomic Sequencing in the Search for a Child's Best Possible Life

Background: The desire of parents to obtain a genetic diagnosis for their child with intellectual disability and associated symptoms has long been framed as a diagnostic odyssey, an arduous and sometimes perilous journey focused on the goal of identifying a cause for the child's condition.Methods: Semi-structured interviews (N = 60) were conducted with parents of children (N = 59, aged 2-24 years) with intellectual disability and/or developmental delay (IDD) who underwent genome sequencing at a single pediatric multispecialty clinic. Interviews were conducted after parents received their child's sequencing result (positive findings, negative findings, or variants of unknown significance). Thematic analysis was performed on all interviews.Results: Parents reported that obtaining a genetic diagnosis was one important step in their overall goal of helping their child live their best life possible life. They intended to use the result as a tool to help their child by seeking the correct school placement and obtaining benefits and therapeutic services.Conclusions: For the parents of children with IDD, the search for a genetic diagnosis is best conceptualized as a part of parents' ongoing efforts to leverage various diagnoses to obtain educational and therapeutic services for their children. Cleaving parents' search for a genetic diagnosis from these broader efforts obscures the value that some parents place on a sequencing result in finding and tailoring therapies and services beyond the clinic. Interviews with parents reveal, therefore, that genomic sequencing is best understood as one important stage of an ongoing therapeutic odyssey that largely takes place outside the clinic. Findings suggest the need to expand translational research efforts to contextualize a genetic diagnosis within parents' broader efforts to obtain educational and therapeutic services outside clinical contexts.

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.

Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy

Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.

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.

Periodic patterns in Rodentia: Development and evolution

Mammalian periodic pigment patterns, such as spots and stripes, have long interested mathematicians and biologists because they arise from non-random developmental processes that are programmed to be spatially constrained, and can therefore be used as a model to understand how organized morphological structures develop. Despite such interest, the developmental and molecular processes underlying their formation remain poorly understood. Here, we argue that Arvicanthines, a clade of African rodents that naturally evolved a remarkable array of coat patterns, represent a tractable model system in which to dissect the mechanistic basis of pigment pattern formation. Indeed, we review recent insights into the process of stripe formation that were obtained using an Arvicanthine species, the African striped mouse (Rhabdomys pumilio), and discuss how these rodents can be used to probe deeply into our understanding of the factors that specify and implement positional information in the skin. By combining naturally evolved pigment pattern variation in rodents with classic and novel experimental approaches, we can substantially advance our understanding of the processes by which spatial patterns of cell differentiation are established during embryogenesis, a fundamental question in developmental biology.

Aberrant Inclusion of a Poison Exon Causes Dravet Syndrome and Related SCN1A-Associated Genetic Epilepsies

Developmental and epileptic encephalopathies (DEEs) are a group of severe epilepsies characterized by refractory seizures and developmental impairment. Sequencing approaches have identified causal genetic variants in only about 50% of individuals with DEEs.^1-3 This suggests that unknown genetic etiologies exist, potentially in the ∼98% of human genomes not covered by exome sequencing (ES). Here we describe seven likely pathogenic variants in regions outside of the annotated coding exons of the most frequently implicated epilepsy gene, SCN1A, encoding the alpha-1 sodium channel subunit. We provide evidence that five of these variants promote inclusion of a "poison" exon that leads to reduced amounts of full-length SCN1A protein. This mechanism is likely to be broadly relevant to human disease; transcriptome studies have revealed hundreds of poison exons,^4,5 including some present within genes encoding other sodium channels and in genes involved in neurodevelopment more broadly.⁶ Future research on the mechanisms that govern neuronal-specific splicing behavior might allow researchers to co-opt this system for RNA therapeutics.

De novo mutations in the GTP/GDP-binding region of RALA, a RAS-like small GTPase, cause intellectual disability and developmental delay

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⁻¹¹) 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.

While many causes of developmental disabilities have been identified, a large number of affected children cannot be diagnosed despite extensive medical testing. Previously unknown genetic factors are likely to be the culprits in many of these cases. Using DNA sequencing, and by sharing information among many doctors and researchers, we have identified a set of individuals with developmental problems who all have changes to the same gene, RALA. The affected individuals all have similar symptoms, including intellectual disability, speech delay (or no speech), and problems with motor skills like walking. In nearly all of these cases (10 of 11), the genetic change found in the child was not inherited from either parent. The locations and biological properties of these changes suggest that they are likely to disrupt the normal functions of RALA. Functional experiments also show that the genetic changes found in these individuals alter two key functions of RALA. Together, we have provided evidence that genetic changes in RALA can cause developmental disabilities. These results will allow doctors and researchers to identify additional children with the same condition, providing a clinical diagnosis to these families and leading to new research opportunities.

Systematic reanalysis of genomic data improves quality of variant interpretation

As genomic sequencing expands, so does our knowledge of the link between genetic variation and disease. Deeper catalogs of variant frequencies improve identification of benign variants, while sequencing affected individuals reveals disease-associated variation. Accumulation of human genetic data thus makes reanalysis a means to maximize the benefits of clinical sequencing. We implemented pipelines to systematically reassess sequencing data from 494 individuals with developmental disability. Reanalysis yielded pathogenic or likely pathogenic (P/LP) variants that were not initially reported in 23 individuals, 6 described here, comprising a 16% increase in P/LP yield. We also downgraded 3 LP and 6 variants of uncertain significance (VUS) due to updated population frequency data. The likelihood of identifying a new P/LP variant increased over time, as ~22% of individuals who did not receive a P/LP variant at their original analysis subsequently did after 3 years. We show here that reanalysis and data sharing increase the diagnostic yield and accuracy of clinical sequencing.

Genomic sequencing identifies secondary findings in a cohort of parent study participants

PURPOSE

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

METHODS

Exome/genome sequencing and analysis of 789 "unaffected" parents was performed.

RESULTS

Pathogenic/likely pathogenic variants were identified in 21 genes within 25 individuals (3.2%), with 11 (1.4%) participants harboring variation in a gene defined as clinically actionable by the American College of Medical Genetics and Genomics. These 25 individuals self-reported either relevant clinical diagnoses (5); relevant family history or symptoms (13); or no relevant family history, symptoms, or clinical diagnoses (7). A limited carrier screen was performed yielding 15 variants in 48 (6.1%) parents. Parents were also analyzed as mate pairs (n = 365) to identify cases in which both parents were carriers for the same recessive disease, yielding three such cases (0.8%), two of which had children with the relevant recessive disease. Four participants had two findings (one carrier and one noncarrier variant). In total, 71 of the 789 enrolled parents (9.0%) received secondary findings.

CONCLUSION

We provide an overview of the rates and types of clinically relevant secondary findings, which may be useful in the design and implementation of research and clinical sequencing efforts to identify such findings.

Inference on the Genetic Basis of Eye and Skin Color in an Admixed Population via Bayesian Linear Mixed Models

Genetic association studies in admixed populations are underrepresented in the genomics literature, with a key concern for researchers being the adequate control of spurious associations due to population structure. Linear mixed models (LMMs) are well suited for genome-wide association studies (GWAS) because they account for both population stratification and cryptic relatedness and achieve increased statistical power by jointly modeling all genotyped markers. Additionally, Bayesian LMMs allow for more flexible assumptions about the underlying distribution of genetic effects, and can concurrently estimate the proportion of phenotypic variance explained by genetic markers. Using three recently published Bayesian LMMs, Bayes R, BSLMM, and BOLT-LMM, we investigate an existing data set on eye (n = 625) and skin (n = 684) color from Cape Verde, an island nation off West Africa that is home to individuals with a broad range of phenotypic values for eye and skin color due to the mix of West African and European ancestry. We use simulations to demonstrate the utility of Bayesian LMMs for mapping loci and studying the genetic architecture of quantitative traits in admixed populations. The Bayesian LMMs provide evidence for two new pigmentation loci: one for eye color (AHRR) and one for skin color (DDB1).

Genomic diagnosis for children with intellectual disability and/or developmental delay

Background

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.

Methods

Whole-exome sequences (WES) were generated for 365 individuals (127 affected) and whole-genome sequences (WGS) were generated for 612 individuals (244 affected).

Results

Pathogenic or likely pathogenic variants were found in 100 individuals (27%), with variants of uncertain significance in an additional 42 (11.3%). We found that a family history of neurological disease, especially the presence of an affected first-degree relative, reduces the pathogenic/likely pathogenic variant identification rate, reflecting both the disease relevance and ease of interpretation of de novo variants. We also found that improvements to genetic knowledge facilitated interpretation changes in many cases. Through systematic reanalyses, we have thus far reclassified 15 variants, with 11.3% of families who initially were found to harbor a VUS and 4.7% of families with a negative result eventually found to harbor a pathogenic or likely pathogenic variant. To further such progress, the data described here are being shared through ClinVar, GeneMatcher, and dbGaP.

Conclusions

Our data strongly support the value of large-scale sequencing, especially WGS within proband-parent trios, as both an effective first-choice diagnostic tool and means to advance clinical and research progress related to pediatric neurological disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0433-1) contains supplementary material, which is available to authorized users.

Comparative tissue transcriptomics highlights dynamic differences among tissues but conserved metabolic transcript prioritization in preparation for arousal from torpor

During the hibernation season, 13-lined ground squirrels spend days to weeks in torpor with body temperatures near freezing then spontaneously rewarm. The molecular drivers of the drastic physiological changes that orchestrate and permit torpor are not well understood. Although transcription effectively ceases at the low body temperatures of torpor, previous work has demonstrated that some transcripts are protected from bulk degradation in brown adipose tissue (BAT), consistent with the importance of their protein products for metabolic heat generation during arousal from torpor. We examined the transcriptome of skeletal muscle, heart, and liver to determine the patterns of differentially expressed genes in these tissues, and whether, like BAT, a subset of these were relatively increased during torpor. EDGE-tags were quantified from five distinct physiological states representing the seasonal and torpor-arousal cycles of 13-lined ground squirrels. Supervised clustering on relative transcript abundances with Random Forest separated the two states bracketing prolonged torpor, entrance into and aroused from torpor, in all three tissues. Independent analyses identified 3347, 6784, and 2433 differentially expressed transcripts among all sampling points in heart, skeletal muscle, and liver, respectively. There were few differentially expressed genes in common across all three tissues; these were enriched in mitochondrial and apoptotic pathway components. Divisive clustering of these data revealed unique cohorts of transcripts that increased across the torpor bout in each tissue with patterns reflecting various combinations of cycling within and between seasons as well as between torpor and arousal. Transcripts that increased across the torpor bout were likewise tissue specific. These data shed new light on the biochemical pathways that alter in concert with hibernation phenotype and provide a rich resource for further hypothesis-based studies.

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.

Evolution: How the zebrafish got its stripes

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

A Decad(e) of Reasons to Contribute to a PLOS Community-Run Journal

Will be part of a Tenth Anniversary Collection

Dominant Red Coat Color in Holstein Cattle Is Associated with a Missense Mutation in the Coatomer Protein Complex, Subunit Alpha (COPA) Gene

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 (MC1R^D) and Recessive Red (MC1R^e). 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 MC1R^D. 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.

Recurrent evolution of melanism in South American felids

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.

Color polymorphism in closely related animal species provides an opportunity to study how the balance between natural selection and genetic drift shapes the evolution of appearance and form. The cat family, Felidae, is especially interesting; 13 of 37 extant species exhibit polymorphism for melanism, but evidence for any adaptive role is lacking, in part because the potential benefits of melanism to felid predators are not clear, and in part because the tools for genomic analysis of natural populations are limited. We identify the mutations responsible for melanism in three closely related South American wild felids, the pampas cat, the kodkod, and Geoffroy’s cat, then adapt a new approach for targeted genome sequencing to characterize molecular variation in the region surrounding each melanism mutation. We find that each mutation has developed independently, with strong evidence for natural selection in the black pampas cat, and reduced genetic variation in the entire population of kodkods. Our results demonstrate that some “black cats” are black not by chance, but by selection for a mutation that provides increased fitness.

Enhanced stability and polyadenylation of select mRNAs support rapid thermogenesis in the brown fat of a hibernator

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.

DOI:http://dx.doi.org/10.7554/eLife.04517.001

Many mammals hibernate to avoid food scarcity and harsh conditions during winter. Hibernation involves entering a state called torpor, which drastically reduces the amount of energy used by the body. During torpor, body temperature also decreases. This is particularly exemplified in ground squirrels, whose body temperature can hover at just above or even below the point of freezing. However, hibernating mammals cannot remain in this state continuously over the months of hibernation but instead cycle between bouts of torpor lasting for 1–3 weeks and brief periods of ‘arousal’ lasting between 12–24 hr, during which their body rapidly warms up.

The heat required to start warming up the hibernator is generated from a specialized form of fat called brown adipose tissue. Normally, the bursts of metabolic activity that are required to create this heat depend on certain proteins being produced. Making a protein involves ‘translating’ its sequence from template molecules called messenger RNA (mRNA), which are ‘transcribed’ from the gene that encodes the protein. During the low body temperatures experienced during torpor, both of these processes stop. So how is the hibernator able to quickly and efficiently heat itself up during the arousal periods of hibernation?

Grabek et al. investigated this by analyzing the relative levels of mRNA in the brown adipose tissue of hibernating 13-lined ground squirrels. Using a special technique to sample and sequence small fragments of mRNA taken from brown adipose tissue, Grabek et al. compiled a profile of the mRNA molecules present at different points in the torpor–arousal cycle and compared this with a similar profile taken from squirrels that were not hibernating.

From this analysis, Grabek et al. detected that a particular group of mRNA molecules that are required for producing heat increase in abundance during torpor, even though body temperature is low enough to stop gene transcription. This increased abundance does not occur because more of the mRNA molecules are made; instead, the mRNA molecules are modified to become more stable and long lasting. Once the animal warms up during arousal, gene transcription is reactivated and more new mRNA molecules are made.

Grabek et al. suggest that the key mRNAs required for brown adipose tissue function are selectively stabilized during torpor through a temperature-dependent protective mechanism. These mRNAs are then preferentially translated into proteins during arousal to rapidly and efficiently heat the hibernator. Most other mRNA molecules degrade throughout torpor, and so their numbers decline as replacements are not transcribed until body temperature briefly recovers during arousal. Whether this protective mechanism is also used in other tissues during torpor remains a question for future work.

DOI:http://dx.doi.org/10.7554/eLife.04517.002

Genetic architecture of skin and eye color in an African-European admixed population

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⁻⁶², SLC24A5 P = 9.6×10⁻⁹) 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⁻²⁷, TYR P = 1.1×10⁻⁹, APBA2[OCA2] P = 1.5×10⁻⁸, SLC45A2 P = 6×10⁻⁹) 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.

Differences in skin and eye color are some of the most obvious traits that underlie human diversity, yet most of our knowledge regarding the genetic basis for these traits is based on the limited range of variation represented by individuals of European ancestry. We have studied a unique population in Cape Verde, an archipelago located off the West African coast, in which extensive mixing between individuals of Portuguese and West African ancestry has given rise to a broad range of phenotypes and ancestral genome proportions. Our results help to explain how genes work together to control the full range of pigmentary phenotypic diversity, provide new insight into the evolution of these traits, and provide a model for understanding other types of quantitative variation in admixed populations.

Levels of the Mahogunin Ring Finger 1 E3 ubiquitin ligase do not influence prion disease

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.

Genetics of Pigmentation in Dogs and Cats

Color variation in companion animals has long been of interest to the breeding and scientific communities. Simple traits, like black versus brown or yellow versus black, have helped to explain principles of transmission genetics and continue to serve as models for studying gene action and interaction. We present a molecular genetic review of pigmentary variation in dogs and cats using a nomenclature and logical framework established by early leaders in the field. For most loci in which molecular variants have been identified (nine in dogs and seven in cats), homologous mutations exist in laboratory mice and/or humans. Exceptions include the K locus in dogs and the Tabby locus in cats, which give rise to alternating stripes or marks of different color, and which illustrate the continued potential of coat color genetics to provide insight into areas that transcend pigment cell biology.

Unraveling the thread of nature's tapestry: the genetics of diversity and convergence in animal pigmentation

Animals display incredibly diverse color patterns yet little is known about the underlying genetic basis of these phenotypes. However, emerging results are reshaping our view of how the process of phenotypic evolution occurs. Here, we outline recent research from three particularly active areas of investigation: melanin pigmentation in Drosophila, wing patterning in butterflies, and pigment variation in lizards. For each system, we highlight (i) the function and evolution of color variation, (ii) various approaches that have been used to explore the genetic basis of pigment variation, and (iii) conclusions regarding the genetic basis of convergent evolution which have emerged from comparative analyses. Results from these studies indicate that natural variation in pigmentation is a particularly powerful tool to examine the molecular basis of evolution, especially with regard to convergent or parallel evolution. Comparison of these systems also reveals that the molecular basis of convergent evolution is heterogeneous, sometimes involving conserved mechanisms and sometimes not. In the near future, additional work in other emerging systems will substantially expand the scope of available comparisons.

Specifying and sustaining pigmentation patterns in domestic and wild cats

Color markings among felid species display both a remarkable diversity and a common underlying periodicity. A similar range of patterns in domestic cats suggests a conserved mechanism whose appearance can be altered by selection. We identified the gene responsible for tabby pattern variation in domestic cats as Transmembrane aminopeptidase Q (Taqpep), which encodes a membrane-bound metalloprotease. Analyzing 31 other felid species, we identified Taqpep as the cause of the rare king cheetah phenotype, in which spots coalesce into blotches and stripes. Histologic, genomic expression, and transgenic mouse studies indicate that paracrine expression of Endothelin3 (Edn3) coordinates localized color differences. We propose a two-stage model in which Taqpep helps to establish a periodic pre-pattern during skin development that is later implemented by differential expression of Edn3.