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.

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.

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.

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.

Genome sequencing as a first-line diagnostic test for hospitalized infants.

PURPOSE: SouthSeq is a translational research study that undertook genome sequencing (GS) for infants with symptoms suggestive of a genetic disorder. Recruitment targeted racial/ethnic minorities and rural, medically underserved areas in the Southeastern United States, which are historically underrepresented in genomic medicine research. METHODS: GS and analysis were performed for 367 infants to detect disease-causal variation concurrent with standard of care evaluation and testing. RESULTS: Definitive diagnostic (DD) or likely diagnostic (LD) genetic findings were identified in 30% of infants, and 14% of infants harbored an uncertain result. Only 43% of DD/LD findings were identified via concurrent clinical genetic testing, suggesting that GS testing is better for obtaining early genetic diagnosis. We also identified phenotypes that correlate with the likelihood of receiving a DD/LD finding, such as craniofacial, ophthalmologic, auditory, skin, and hair abnormalities. We did not observe any differences in diagnostic rates between racial/ethnic groups. CONCLUSION: We describe one of the largest-to-date GS cohorts of ill infants, enriched for African American and rural patients. Our results show the utility of GS because it provides early-in-life detection of clinically relevant genetic variations not detected by current clinical genetic testing, particularly for infants exhibiting certain phenotypic features.

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.

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.

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.

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

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

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.

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.

Identifying rare, medically relevant variation via population-based genomic screening in Alabama: opportunities and pitfalls.

PURPOSE: To evaluate the effectiveness and specificity of population-based genomic screening in Alabama. METHODS: The Alabama Genomic Health Initiative (AGHI) has enrolled and evaluated 5369 participants for the presence of pathogenic/likely pathogenic (P/LP) variants using the Illumina Global Screening Array (GSA), with validation of all P/LP variants via Sanger sequencing in a CLIA-certified laboratory before return of results. RESULTS: Among 131 variants identified by the GSA that were evaluated by Sanger sequencing, 67 (51%) were false positives (FP). For 39 of the 67 FP variants, a benign/likely benign variant was present at or near the targeted P/LP variant. Variants detected within African American individuals were significantly enriched for FPs, likely due to a higher rate of nontargeted alternative alleles close to array-targeted P/LP variants. CONCLUSION: In AGHI, we have implemented an array-based process to screen for highly penetrant genetic variants in actionable disease genes. We demonstrate the need for clinical validation of array-identified variants in direct-to-consumer or population testing, especially for diverse populations.

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.

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.

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

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