A high resolution RH map of the bovine major histocompatibility complex

<p>Abstract</p> <p>Background</p> <p>The cattle MHC is termed the bovine leukocyte antigen (BoLA) and, along with the MHCs of other ruminants, is unique in its genomic organization. Consequently, correct and reliable gene maps and sequence information are critical to the study of the BoLA region. The bovine genome sequencing project has produced two assemblies (Btau_3.1 and 4.0) that differ substantially from each other and from conventional gene maps in the BoLA region. To independently compare the accuracies of the different sequence assemblies, we have generated a high resolution map of BoLA using a 12,000_radradiation hybrid panel. Seventy-seven unique sequence tagged site (STS) markers chosen at approximately 50 kb intervals from the Btau 2.0 assembly and spanning the IIa-III-I and IIb regions of the bovine MHC were mapped on a 12,000_radbovine radiation hybrid (RH) panel to evaluate the different assemblies of the bovine genome sequence.</p> <p>Results</p> <p>Analysis of the data generated a high resolution RH map of BoLA that was significantly different from the Btau_3.1 assembly of the bovine genome but in good agreement with the Btau_4.0 assembly. Of the few discordancies between the RH map and Btau_4.0, most could be attributed to closely spaced markers that could not be precisely ordered in the RH panel. One probable incorrectly-assembled sequence and three missing sequences were noted in the Btau_4.0 assembly. The RH map of BoLA is also highly concordant with the sequence-based map of HLA (NCBI build 36) when reordered to account for the ancestral inversion in the ruminant MHC.</p> <p>Conclusion</p> <p>These results strongly suggest that studies using Btau_3.1 for analyses of the BoLA region should be reevaluated in light of the Btau_4.0 assembly and indicate that additional research is needed to produce a complete assembly of the BoLA genomic sequences.</p

Whole genome sequencing reveals a 7 base-pair deletion in DMD exon 42 in a dog with muscular dystrophy

Dystrophin is a key cytoskeletal protein coded by the Duchenne muscular dystrophy (DMD) gene located on the X-chromosome. Truncating mutations in the DMD gene cause loss of dystrophin and the classical DMD clinical syndrome. Spontaneous DMD gene mutations and associated phenotypes occur in several other species. The mdx mouse model and the golden retriever muscular dystrophy (GRMD) canine model have been used extensively to study DMD disease pathogenesis and show efficacy and side effects of putative treatments. Certain DMD gene mutations in high-risk, the so-called hot spot areas can be particularly helpful in modeling molecular therapies. Identification of specific mutations has been greatly enhanced by new genomic methods. Whole genome, next generation sequencing (WGS) has been recently used to define DMD patient mutations, but has not been used in dystrophic dogs. A dystrophin-deficient Cavalier King Charles Spaniel (CKCS) dog was evaluated at the functional, histopathological, biochemical, and molecular level. The affected dog’s phenotype was compared to the previously reported canine dystrophinopathies. WGS was then used to detect a 7 base pair deletion in DMD exon 42 (c.6051-6057delTCTCAAT mRNA), predicting a frameshift in gene transcription and truncation of dystrophin protein translation. The deletion was confirmed with conventional PCR and Sanger sequencing. This mutation is in a secondary DMD gene hotspot area distinct from the one identified earlier at the 5′ donor splice site of intron 50 in the CKCS breed. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00335-016-9675-2) contains supplementary material, which is available to authorized users

Assignment of chromosomal locations for unassigned SNPs/scaffolds based on pair-wise linkage disequilibrium estimates

<p>Abstract</p> <p>Background</p> <p>Recent developments of high-density SNP chips across a number of species require accurate genetic maps. Despite rapid advances in genome sequence assembly and availability of a number of tools for creating genetic maps, the exact genome location for a number of SNPs from these SNP chips still remains unknown. We have developed a locus ordering procedure based on linkage disequilibrium (LODE) which provides estimation of the chromosomal positions of unaligned SNPs and scaffolds. It also provides an alternative means for verification of genetic maps. We exemplified LODE in cattle.</p> <p>Results</p> <p>The utility of the LODE procedure was demonstrated using data from 1,943 bulls genotyped for 73,569 SNPs across three different SNP chips. First, the utility of the procedure was tested by analysing the masked positions of 1,500 randomly-chosen SNPs with known locations (50 from each chromosome), representing three classes of minor allele frequencies (MAF), namely >0.05, 0.01<MAF ≤ 0.05 and 0.001<MAF ≤ 0.01. The efficiency (percentage of masked SNPs that could be assigned a location) was 96.7%, 30.6% and 2.0%; with an accuracy (the percentage of SNPs assigned correctly) of 99.9%, 98.9% and 33.3% in the three classes of MAF, respectively. The average precision for placement of the SNPs was 914, 3,137 and 6,853 kb, respectively. Secondly, 4,688 of 5,314 SNPs unpositioned in the Btau4.0 assembly were positioned using the LODE procedure. Based on these results, the positions of 485 unordered scaffolds were determined. The procedure was also used to validate the genome positions of 53,068 SNPs placed on Btau4.0 bovine assembly, resulting in identification of problem areas in the assembly. Finally, the accuracy of the LODE procedure was independently validated by comparative mapping on the hg18 human assembly.</p> <p>Conclusion</p> <p>The LODE procedure described in this study is an efficient and accurate method for positioning SNPs (MAF>0.05), for validating and checking the quality of a genome assembly, and offers a means for positioning of unordered scaffolds containing SNPs. The LODE procedure will be helpful in refining genome sequence assemblies, especially those being created from next-generation sequencing where high-throughput SNP discovery and genotyping platforms are integrated components of genome analysis.</p

Dystrophin-deficient dogs with reduced myostatin have unequal muscle growth and greater joint contractures

Abstract Background Myostatin (Mstn) is a negative regulator of muscle growth whose inhibition promotes muscle growth and regeneration. Dystrophin-deficient mdx mice in which myostatin is knocked out or inhibited postnatally have a less severe phenotype with greater total mass and strength and less fibrosis and fatty replacement of muscles than mdx mice with wild-type myostatin expression. Dogs with golden retriever muscular dystrophy (GRMD) have previously been noted to have increased muscle mass and reduced fibrosis after systemic postnatal myostatin inhibition. Based partly on these results, myostatin inhibitors are in development for use in human muscular dystrophies. However, persisting concerns regarding the effects of long-term and profound myostatin inhibition will not be easily or imminently answered in clinical trials. Methods To address these concerns, we developed a canine (GRippet) model by crossbreeding dystrophin-deficient GRMD dogs with Mstn-heterozygous (Mstn +/−) whippets. A total of four GRippets (dystrophic and Mstn +/−), three GRMD (dystrophic and Mstn wild-type) dogs, and three non-dystrophic controls from two litters were evaluated. Results Myostatin messenger ribonucleic acid (mRNA) and protein levels were downregulated in both GRMD and GRippet dogs. GRippets had more severe postural changes and larger (more restricted) maximal joint flexion angles, apparently due to further exaggeration of disproportionate effects on muscle size. Flexors such as the cranial sartorius were more hypertrophied on magnetic resonance imaging (MRI) in the GRippets, while extensors, including the quadriceps femoris, underwent greater atrophy. Myostatin protein levels negatively correlated with relative cranial sartorius muscle cross-sectional area on MRI, supporting a role in disproportionate muscle size. Activin receptor type IIB (ActRIIB) expression was higher in dystrophic versus control dogs, consistent with physiologic feedback between myostatin and ActRIIB. However, there was no differential expression between GRMD and GRippet dogs. Satellite cell exhaustion was not observed in GRippets up to 3 years of age. Conclusions Partial myostatin loss may exaggerate selective muscle hypertrophy or atrophy/hypoplasia in GRMD dogs and worsen contractures. While muscle imbalance is not a feature of myostatin inhibition in mdx mice, findings in a larger animal model could translate to human experience with myostatin inhibitors

The evolutionary dynamics of microRNAs in domestic mammals

MiRNAs are crucial regulators of gene expression found across both the plant and animal kingdoms. While the number of annotated miRNAs deposited in miRBase has greatly increased in recent years, few studies provided comparative analyses across sets of related species, or investigated the role of miRNAs in the evolution of gene regulation. We generated small RNA libraries across 5 mammalian species (cow, dog, horse, pig and rabbit) from 4 different tissues (brain, heart, kidney and testis). We identified 1676 miRBase and 413 novel miRNAs by manually curating the set of computational predictions obtained from miRCat and miRDeep2. Our dataset spanning five species has enabled us to investigate the molecular mechanisms and selective pressures driving the evolution of miRNAs in mammals. We highlight the important contributions of intronic sequences (366 orthogroups), duplication events (135 orthogroups) and repetitive elements (37 orthogroups) in the emergence of new miRNA loci. We use this framework to estimate the patterns of gains and losses across the phylogeny, and observe high levels of miRNA turnover. Additionally, the identification of lineage-specific losses enables the characterisation of the selective constraints acting on the associated target sites. Compared to the miRBase subset, novel miRNAs tend to be more tissue specific. 20 percent of novel orthogroups are restricted to the brain, and their target repertoires appear to be enriched for neuron activity and differentiation processes. These findings may reflect an important role for young miRNAs in the evolution of brain expression plasticity. Many seed sequences appear to be specific to either the cow or the dog. Analyses on the associated targets highlight the presence of several genes under artificial positive selection, suggesting an involvement of these miRNAs in the domestication process. Altogether, we provide an overview on the evolutionary mechanisms responsible for miRNA turnover in 5 domestic species, and their possible contribution to the evolution of gene regulation

Serum Cytokines Predict Neurological Damage in Genetically Diverse Mouse Models.

Viral infections contribute to neurological and immunological dysfunction driven by complex genetic networks. Theiler's murine encephalomyelitis virus (TMEV) causes neurological dysfunction in mice and can model human outcomes to viral infections. Here, we used genetically distinct mice from five Collaborative Cross mouse strains and C57BL/6J to demonstrate how TMEV-induced immune responses in serum may predict neurological outcomes in acute infection. To test the hypothesis that serum cytokine levels can provide biomarkers for phenotypic outcomes of acute disease, we compared cytokine levels at pre-injection, 4 days post-injection (d.p.i.), and 14 d.p.i. Each strain produced unique baseline cytokine levels and had distinct immune responses to the injection procedure itself. Thus, we eliminated the baseline responses to the injection procedure itself and identified cytokines and chemokines induced specifically by TMEV infection. Then, we identified strain-specific longitudinal cytokine profiles in serum during acute disease. Using stepwise regression analysis, we identified serum immune markers predictive for TMEV-induced neurological phenotypes of the acute phase, e.g., IL-9 for limb paralysis; and TNF-α, IL-1β, and MIP-1β for limb weakness. These findings indicate how temporal differences in immune responses are influenced by host genetic background and demonstrate the potential of serum biomarkers to track the neurological effects of viral infection

Genetic and immunological contributors to virus-induced paralysis.

Infection by a single virus can evoke diverse immune responses, resulting in different neurological outcomes, depending on the host's genetic background. To study heterogenous viral response, we use Theiler's Murine Encephalomyelitis Virus (TMEV) to model virally induced neurological phenotypes and immune responses in Collaborative Cross (CC) mice. The CC resource consists of genetically distinct and reproducible mouse lines, thus providing a population model with genetic heterogeneity similar to humans. We examined different CC strains for the effect of chronic stage TMEV-induced immune responses on neurological outcomes throughout 90 days post infection (dpi), with a particular focus on limb paralysis, by measuring serum levels of 23 different cytokines and chemokines. Each CC strain demonstrated a unique set of immune responses, regardless of presence or absence of TMEV RNA. Using stepwise regression, significant associations were identified between IL-1α, RANTES, and paralysis frequency scores. To better understand these interactions, we evaluated multiple aspects of the different CC genetic backgrounds, including haplotypes of genomic regions previously linked with TMEV pathogenesis and viral clearance or persistence, individual cytokine levels, and TMEV-relevant gene expression. These results demonstrate how loci previously associated with TMEV outcomes provide incomplete information regarding TMEV-induced paralysis in the CC strains. Overall, these findings provide insight into the complex roles of immune response in the pathogenesis of virus-associated neurological diseases influenced by host genetic background

Resilience in Long-Term Viral Infection: Genetic Determinants and Interactions

Virus-induced neurological sequelae resulting from infection by Theiler’s murine encephalomyelitis virus (TMEV) are used for studying human conditions ranging from epileptic seizures to demyelinating disease. Mouse strains are typically considered susceptible or resistant to TMEV infection based on viral persistence and extreme phenotypes, such as demyelination. We have identified a broader spectrum of phenotypic outcomes by infecting strains of the genetically diverse Collaborative Cross (CC) mouse resource. We evaluated the chronic-infection gene expression profiles of hippocampi and thoracic spinal cords for 19 CC strains in relation to phenotypic severity and TMEV persistence. Strains were clustered based on similar phenotypic profiles and TMEV levels at 90 days post-infection, and we categorized distinct TMEV response profiles. The three most common profiles included “resistant” and “susceptible,” as before, as well as a “resilient” TMEV response group which experienced both TMEV persistence and mild neurological phenotypes even at 90 days post-infection. Each profile had a distinct gene expression signature, allowing the identification of pathways and networks specific to each TMEV response group. CC founder haplotypes for genes involved in these pathways/networks revealed candidate response-specific alleles. These alleles demonstrated pleiotropy and epigenetic (miRNA) regulation in long-term TMEV infection, with particular relevance for resilient mouse strains