Whole-Genome sequencing and genetic variant analysis of a quarter Horse mare

BACKGROUND: The catalog of genetic variants in the horse genome originates from a few select animals, the majority originating from the Thoroughbred mare used for the equine genome sequencing project. The purpose of this study was to identify genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs) in the genome of an individual Quarter Horse mare sequenced by next-generation sequencing. RESULTS: Using massively parallel paired-end sequencing, we generated 59.6 Gb of DNA sequence from a Quarter Horse mare resulting in an average of 24.7X sequence coverage. Reads were mapped to approximately 97% of the reference Thoroughbred genome. Unmapped reads were de novo assembled resulting in 19.1 Mb of new genomic sequence in the horse. Using a stringent filtering method, we identified 3.1 million SNPs, 193 thousand INDELs, and 282 CNVs. Genetic variants were annotated to determine their impact on gene structure and function. Additionally, we genotyped this Quarter Horse for mutations of known diseases and for variants associated with particular traits. Functional clustering analysis of genetic variants revealed that most of the genetic variation in the horse's genome was enriched in sensory perception, signal transduction, and immunity and defense pathways. CONCLUSIONS: This is the first sequencing of a horse genome by next-generation sequencing and the first genomic sequence of an individual Quarter Horse mare. We have increased the catalog of genetic variants for use in equine genomics by the addition of novel SNPs, INDELs, and CNVs. The genetic variants described here will be a useful resource for future studies of genetic variation regulating performance traits and diseases in equids

Identification of Genomic Loci Associated with Rhodococcus equi Susceptibility in Foals

Pneumonia caused by Rhodococcus equi is a common cause of disease and death in foals. Although agent and environmental factors contribute to the incidence of this disease, the genetic factors influencing the clinical outcomes of R. equi pneumonia are ill-defined. Here, we performed independent single nucleotide polymorphism (SNP)- and copy number variant (CNV)-based genome-wide association studies to identify genomic loci associated with R. equi pneumonia in foals. Foals at a large Quarter Horse breeding farm were categorized into 3 groups: 1) foals with R. equi pneumonia (clinical group [N = 43]); 2) foals with ultrasonographic evidence of pulmonary lesions that never developed clinical signs of pneumonia (subclinical group [N = 156]); and, 3) foals without clinical signs or ultrasonographic evidence of pneumonia (unaffected group [N = 49]). From each group, 24 foals were randomly selected and used for independent SNP- and CNV-based genome-wide association studies (GWAS). The SNP-based GWAS identified a region on chromosome 26 that had moderate evidence of association with R. equi pneumonia when comparing clinical and subclinical foals. A joint analysis including all study foals revealed a 3- to 4-fold increase in odds of disease for a homozygous SNP within the associated region when comparing the clinical group with either of the other 2 groups of foals or their combination. The region contains the transient receptor potential cation channel, subfamily M, member 2 (TRPM2) gene, which is involved in neutrophil function. No associations were identified in the CNV-based GWAS. Collectively, these data identify a region on chromosome 26 associated with R. equi pneumonia in foals, providing evidence that genetic factors may indeed contribute to this important disease of foals.The open access fee for this work was funded through the Texas A&M University Open Access to Knowledge (OAK) Fund

Epigenetic profiling at mouse imprinted gene clusters reveals novel epigenetic and genetic features at differentially methylated regions

Genomic imprinting arises from allele-specific epigenetic modifications that are established during gametogenesis and that are maintained throughout somatic development. These parental-specific modifications include DNA methylation and post-translational modifications to histones, which create allele-specific active and repressive domains at imprinted regions. Through the use of a high-density genomic tiling array, we generated DNA and histone methylation profiles at 11 imprinted gene clusters in the mouse from DNA and from chromatin immunoprecipitated from sperm, heart, and cerebellum. Our analysis revealed that despite high levels of differential DNA methylation at non-CpG islands within these regions, imprinting control regions (ICRs) and secondary differentially methylated regions (DMRs) were identified by an overlapping pattern of H3K4 trimethylation (active chromatin) and H3K9 trimethylation (repressive chromatin) modifications in somatic tissue, and a sperm differentially methylated region (sDMR; sperm ≠ somatic tissue). Using these features as a common signature of DMRs, we identified 11 unique regions that mapped to known imprinted genes, to uncharacterized genes, and to intergenic regions flanking known imprinted genes. A common feature among these regions was the presence of a CpG island and an array of tandem repeats. Collectively, this study provides a comprehensive analysis of DNA methylation and histone H3K4me3 and H3K9me3 modifications at imprinted gene clusters, and identifies common epigenetic and genetic features of regions regulating genomic imprinting

Epigenetic and genomic imprinting analysis in nuclear transfer derived Bos gaurus/Bos taurus hybrid fetuses

Somatic cell nuclear transfer (NT) in cattle is an inefficient process, whereby the production of calves is hindered by low pregnancy rates as well as fetal and placental abnormalities. Interspecies models have been previously used to facilitate the identification of single nucleotide polymorphisms (SNPs) within coding regions of genes to discriminate between parental alleles in the offspring. Here we report the use of a bovine interspecies model (Bos gaurus x Bos taurus) for the assessment and characterization of epigenetic modifications and genomic imprinting in Day 40-old female NT-derived fetuses and placenta. Analysis of NT and control pregnancies indicated disruption of genomic imprinting at the X inactivation-specific transcript (XIST) locus in the chorion, but not the fetus of clones, whereas proper allelic expression of the insulin-like growth factor II (IGF2) and gene trap locus 2 (GTL2) loci was maintained in both the fetus and placenta. Analysis of the XIST differentially methylated region (DMR) in clones indicated normal patterns of methylation; however, bisulfite sequencing of the satellite I repeat element and epidermal cytokeratin promoter indicated hypermethylation in the chorion of clones when compared with controls. No differences were detected in methylation levels in the fetus proper. These results indicate that the nuclear transfer process affects gene expression patterns in the trophectoderm- and inner cell mass-derived tissues to different extents

Unique Features of the Gut Microbiome Characterized in Animal Models of Angelman Syndrome

A large subset of patients with Angelman syndrome (AS) suffer from concurrent gastrointestinal (GI) issues, including constipation, poor feeding, and reflux. AS is caused by the loss of ubiquitin ligase E3A (UBE3A) gene expression in the brain. Clinical features of AS, which include developmental delays, intellectual disability, microcephaly, and seizures, are primarily due to the deficient expression or function of the maternally inherited UBE3A allele. The association between neurodevelopmental delay and GI disorders is part of the increasing evidence suggesting a link between the brain and the gut microbiome via the microbiota-gut-brain axis. To investigate the associations between colonization of the gut microbiota in AS, we characterized the fecal microbiome in three animal models of AS involving maternal deletions of Ube3A, including mouse, rat, and pig, using 16S rRNA amplicon sequencing. Overall, we identified changes in bacterial abundance across all three animal models of AS. Specific bacterial groups were significantly increased across all animal models, including Lachnospiraceae Incertae sedis, Desulfovibrios sp., and Odoribacter, which have been correlated with neuropsychiatric disorders. Taken together, these findings suggest that specific changes to the local environment in the gut are driven by a Ube3a maternal deletion, unaffected by varying housing conditions, and are prominent and detectable across multiple small and large animal model species. These findings begin to uncover the underlying mechanistic causes of GI disorders in AS patients and provide future therapeutic options for AS patients. IMPORTANCE Angelman syndrome (AS)-associated gastrointestinal (GI) symptoms significantly impact quality of life in patients. In AS models in mouse, rat, and pig, AS animals showed impaired colonization of the gut microbiota compared to wild-type (healthy) control animals. Common changes in AS microbiomes across all three animal models may play a causal effect for GI symptoms and may help to identify ways to treat these comorbidities in patients in the future

Early Developmental EEG and Seizure Phenotypes in a Full Gene Deletion of Ubiquitin Protein Ligase E3A Rat Model of Angelman Syndrome.

Angelman syndrome (AS) is a neurodevelopmental disorder with unique behavioral phenotypes, seizures, and distinctive electroencephalographic (EEG) patterns. Recent studies identified motor, social communication, and learning and memory deficits in a CRISPR engineered rat model with a complete maternal deletion of the Ube3a gene. It is unknown whether this model recapitulates other aspects of the clinical disorder. We report here the effect of Ube3a maternal deletion in the rat on epileptiform activity, seizure threshold, and quantitative EEG. Using video-synchronized EEG (vEEG) monitoring, we assessed spectral power and epileptiform activity early postnatally through adulthood. While EEG power was similar to wild-type (WT) at 1.5 weeks postnatally, at all other ages analyzed, our findings were similar to the AS phenotype in mice and humans with significantly increased δ power. Analysis of epileptiform activity in juvenile and adult rats showed increased time spent in epileptiform activity in AS compared with WT rats. We evaluated seizure threshold using pentylenetetrazol (PTZ), audiogenic stimulus, and hyperthermia to provoke febrile seizures (FSs). Behavioral seizure scoring following PTZ induction revealed no difference in seizure threshold in AS rats, however behavioral recovery from the PTZ-induced seizure was longer in the adult group with significantly increased hippocampal epileptiform activity during this phase. When exposed to hyperthermia, AS rat pups showed a significantly lower temperature threshold to first seizure than WT. Our findings highlight an age-dependence for the EEG and epileptiform phenotypes in a preclinical model of AS, and support the use of quantitative EEG and increased δ power as a potential biomarker of AS

Manhattan plots of standard chi-squared significance values for the 3 genome-wide association studies.

<p>Manhattan plots for (A) comparison 1, (B) comparison 2, and (C) comparison 3; (D) Mixed effects-model analysis of comparison 3.</p