BackgroundThe cause of immune-mediated myositis (IMM), characterized by recurrent, rapid-onset muscle atrophy in Quarter Horses (QH), is unknown. The histopathologic hallmark of IMM is lymphocytic infiltration of myofibers. The purpose of this study was to identify putative functional variants associated with equine IMM.MethodsA genome-wide association (GWA) study was performed on 36 IMM QHs and 54 breed matched unaffected QHs from the same environment using the Equine SNP50 and SNP70 genotyping arrays.ResultsA mixed model analysis identified nine SNPs within a ~ 2.87 Mb region on chr11 that were significantly (Punadjusted < 1.4 × 10− 6) associated with the IMM phenotype. Associated haplotypes within this region encompassed 38 annotated genes, including four myosin genes (MYH1, MYH2, MYH3, and MYH13). Whole genome sequencing of four IMM and four unaffected QHs identified a single segregating nonsynonymous E321G mutation in MYH1 encoding myosin heavy chain 2X. Genotyping of additional 35 IMM and 22 unaffected QHs confirmed an association (P = 2.9 × 10− 5), and the putative mutation was absent in 175 horses from 21 non-QH breeds. Lymphocytic infiltrates occurred in type 2X myofibers and the proportion of 2X fibers was decreased in the presence of inflammation. Protein modeling and contact/stability analysis identified 14 residues affected by the mutation which significantly decreased stability.ConclusionsWe conclude that a mutation in MYH1 is highly associated with susceptibility to the IMM phenotype in QH-related breeds. This is the first report of a mutation in MYH1 and the first link between a skeletal muscle myosin mutation and autoimmune disease.Electronic supplementary materialThe online version of this article (10.1186/s13395-018-0155-0) contains supplementary material, which is available to authorized users.

Background Sarcolipin (SLN), myoregulin (MRLN), and dwarf open reading frame (DWORF) are transmembrane regulators of the sarcoplasmic reticulum calcium transporting ATPase (SERCA) that we hypothesized played a role in recurrent exertional rhabdomyolysis (RER). Objectives Compare coding sequences of SLN, MRLN, DWORF across species and between RER and control horses. Compare expression of muscle Ca2+ regulatory genes between RER and control horses. Animals Twenty Thoroughbreds (TB), 5 Standardbreds (STD), 6 Quarter Horses (QH) with RER and 39 breed‐matched controls. Methods Sanger sequencing of SERCA regulatory genes with comparison of amino acid (AA) sequences among control, RER horses, human, mouse, and rabbit reference genomes. In RER and control gluteal muscle, quantitative real‐time polymerase chain reaction of SERCA regulatory peptides, the calcium release channel (RYR1), and its accessory proteins calsequestrin (CASQ1), and calstabin (FKBP1A). Results The SLN gene was the highest expressed horse SERCA regulatory gene with a uniquely truncated AA sequence (29 versus 31) versus other species. Coding sequences of SLN, MRLN, and DWORF were identical in RER and control horses. A sex‐by‐phenotype effect occurred with lower CASQ1 expression in RER males versus control males (P < .001) and RER females (P = .05) and higher FKBP1A (P = .01) expression in RER males versus control males. Conclusions and Clinical Importance The SLN gene encodes a uniquely truncated peptide in the horse versus other species. Variants in the coding sequence of SLN, MLRN, or DWORF were not associated with RER. Males with RER have differential gene expression that could reflect adaptations to stabilize RYR1.

BackgroundMarek’s disease (MD) is a commercially important neoplastic disease of chickens caused by the Marek’s disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Enhancing MD genetic resistance is desirable to augment current vaccines and other MD control measures. High throughput sequencing was used to profile splenic transcriptomes from individual F1 progeny infected with MDV at 4 days of age from both outbred broilers (meat-type) and inbred layer (egg-type) chicken lines that differed in MD genetic resistance. The resulting information was used to identify SNPs, genes, and biological pathways exhibiting allele-specific expression (ASE) in response to MDV infection in each type of chicken. In addition, we compared and contrasted the results of pathway analyses (ASE and differential expression (DE)) between chicken types to help inform on the biological response to MDV infection.ResultsWith 7 individuals per line and treatment group providing high power, we identified 6,132 single nucleotide polymorphisms (SNPs) in 4,768 genes and 4,528 SNPs in 3,718 genes in broilers and layers, respectively, that exhibited ASE in response to MDV infection. Furthermore, 548 and 434 genes in broilers and layers, respectively, were found to show DE following MDV infection. Comparing the datasets, only 72 SNPs and 850 genes for ASE and 20 genes for DE were common between the two bird types. Although the chicken types used in this study were genetically different, at the pathway level, both TLR receptor and JAK/STAT signaling pathways were enriched as well as exhibiting a high proportion of ASE genes, especially at the beginning of both above mentioned regulatory pathways.ConclusionsRNA sequencing with adequate biological replicates is a powerful approach to identify high confidence SNPs, genes, and pathways that are associated with transcriptional response to MDV infection. In addition, the SNPs exhibiting ASE in response to MDV infection provide a strong foundation for determining the extent to which variation in expression influences MD incidence plus yield genetic markers for genomic selection. However, given the paucity of overlap among ASE SNP sets (broilers vs. layers), it is likely that separate screens need to be incorporated for each population. Finally, comparison of gene lists obtained between these two diverse chicken types indicate the TLR and JAK/STAT signaling are conserved when responding to MDV infection and may be altered by selection of genes exhibiting ASE found at the start of each pathway.

BackgroundAn E321G mutation in MYH1 was recently identified in Quarter Horses (QH) with immune‐mediated myositis (IMM) defined by a phenotype of gross muscle atrophy and myofiber lymphocytic infiltrates.Hypothesis/ObjectivesWe hypothesized that the MYH1 mutation also was associated with a phenotype of nonexertional rhabdomyolysis. The objective of this study was to determine the prevalence of the MYH1 mutation in QH with exertional (ER) and nonexertional (nonER) rhabdomyolysis.AnimalsQuarter Horses: 72 healthy controls, 85 ER‐no atrophy, 56 ER‐atrophy, 167 nonER horses selected regardless of muscle atrophy.MethodsClinical and histopathologic information and DNA was obtained from a database for (1) ER > 2 years of age, with or without atrophy and (2) nonER creatine kinase (CK) ≥ 5000 U/L, <5 years of age. Horses were genotyped for E321G MYH1 by pyrosequencing.ResultsThe MYH1 mutation was present in a similar proportion of ER‐no atrophy (1/56; 2%) and in a higher proportion of ER‐atrophy (25/85; 29%) versus controls (4/72; 5%). The MYH1 mutation was present in a significantly higher proportion of nonER (113/165; 68%) than controls either in the presence (39/42; 93%) or in absence (72/123; 59%) of gross atrophy. Lymphocytes were present in <18% of muscle samples with the MYH1 mutation.Conclusions and Clinical ImportanceAlthough not associated with ER, the MYH1 mutation is associated with atrophy after ER. The MYH1 mutation is highly associated with nonER regardless of whether muscle atrophy or lymphocytic infiltrates are present. Genetic testing will enhance the ability to diagnose MYH1 myopathies (MYHM) in QH.

BackgroundMarek’s disease (MD) is a lymphoproliferative disease of poultry induced by Marek’s disease virus (MDV), a highly oncogenic alphaherpesvirus. Identifying the underlying genes conferring MD genetic resistance is desired for more efficacious control measures including genomic selection, which requires accurately identified genetic markers throughout the chicken genome.MethodsHypothesizing that variants located in transcriptional regulatory regions are the main mechanism underlying this complex trait, a genome-wide association study was conducted by genotyping a ~1,000 bird MD resource population derived from experimental inbred layers with SNPs containing 1,824 previously identified allele-specific expression (ASE) SNPs in response to MDV infection as well as 3,097 random SNPs equally spaced throughout the chicken genome. Based on the calculated associations, genomic predictions were determined for 200 roosters and selected sires had their progeny tested for Marek’s disease incidence.ResultsOur analyses indicate that these ASE SNPs account for more than 83 % of the genetic variance and exhibit nearly all the highest associations. To validate these findings, 200 roosters had their genetic merit predicted from the ASE SNPs only, and the top 30 and bottom 30 ranked roosters were reciprocally mated to random hens. The resulting progeny showed that after only one generation of bidirectional selection, there was a 22 % difference in MD incidence and this approach gave a 125 % increase in accuracy compared to current pedigree-based estimates.ConclusionsWe conclude that variation in transcriptional regulation is the major driving cause for genetic resistance to MD, and ASE SNPs identify the underlying genes and are sufficiently linked to the causative polymorphisms that they can be used for accurate genomic prediction as well as help define the underlying molecular basis. Furthermore, this approach should be applicable to other complex traits.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2016-0) contains supplementary material, which is available to authorized users.

Background Immune‐mediated myositis (IMM) in American Quarter Horses (QHs) causes acute muscle atrophy and lymphocytic infiltration of myofibers. Recently, an E321G mutation in a highly conserved region of the myosin heavy chain 1 (MYH1) gene was associated with susceptibility to IMM and nonexertional rhabdomyolysis. Objectives To estimate prevalence of the E321G MYH1 variant in the QH breed and performance subgroups. Animals Three‐hundred seven elite performance QHs and 146 random registered QH controls. Methods Prospective genetic survey. Elite QHs from barrel racing, cutting, halter, racing, reining, Western Pleasure, and working cow disciplines and randomly selected registered QHs were genotyped for the E321G MYH1 variant and allele frequencies were calculated. Results The E321G MYH1 variant allele frequency was 0.034 ± 0.011 in the general QH population (6.8% of individuals in the breed) and the highest among the reining (0.135 ± 0.040; 24.3% of reiners), working cow (0.085 ± 0.031), and halter (0.080 ± 0.027) performance subgroups. The E321G MYH1 variant was present in cutting (0.044 ± 0.022) and Western Pleasure (0.021 ± 0.015) QHs at lower frequency and was not observed in barrel racing or racing QHs. Conclusions and Clinical Importance Knowing that reining and working cow QHs have the highest prevalence of the E321G MYH1 variant and that the variant is more prevalent than the alleles for hereditary equine regional dermal asthenia and hyperkalemic periodic paralysis in the general QH population will guide the use of genetic testing for diagnostic and breeding purposes.

Marek's disease (MD) is an important neoplastic disease of chickens caused by the Marek's disease virus (MDV), an oncogenic alphaherpesvirus. In this study, dysbiosis induced by MDV on the core gut flora of chicken was assessed using next generation sequence (NGS) analysis. Total fecal and cecum-derived samples from individual birds were used to estimate the influence of MDV infection on the gut microbiome of chicken. Our analysis shows that MDV infection alters the core gut flora in the total fecal samples relatively early after infection (2-7 days) and in the late phase of viral infection (28-35 days) in cecal samples, corresponding well with the life cycle of MDV. Principle component analyses of total fecal and cecal samples showed clustering at the early and late time points, respectively. The genus Lactobacillus was exclusively present in the infected samples in both total fecal and cecal bird samples. The community colonization of core gut flora was altered by viral infection, which manifested in the enrichment of several genera during the early and late phases of MDV replication. The results suggest a relationship between viral infection and microbial composition of the intestinal tract that may influence inflammation and immunosuppression of T and B cells in the host.

Equine myofibrillar myopathy (MFM) causes exertional muscle pain and is characterized by myofibrillar disarray and ectopic desmin aggregates of unknown origin. To investigate the pathophysiology of MFM, we compared resting and 3 h postexercise transcriptomes of gluteal muscle and the resting skeletal muscle proteome of MFM and control Arabian horses with RNA sequencing and isobaric tags for relative and absolute quantitation analyses. Three hours after exercise, 191 genes were identified as differentially expressed (DE) in MFM vs. control muscle with >1 log2 fold change (FC) in genes involved in sulfur compound/cysteine metabolism such as cystathionine-beta-synthase (CBS, ↓4.51), a cysteine and neutral amino acid membrane transporter (SLC7A10, ↓1.80 MFM), and a cationic transporter (SLC24A1, ↓1.11 MFM). In MFM vs. control at rest, 284 genes were DE with >1 log2 FC in pathways for structure morphogenesis, fiber organization, tissue development, and cell differentiation including > 1 log2 FC in cardiac alpha actin (ACTC1 ↑2.5 MFM), cytoskeletal desmoplakin (DSP ↑2.4 MFM), and basement membrane usherin (USH2A ↓2.9 MFM). Proteome analysis revealed significantly lower antioxidant peroxiredoxin 6 content (PRDX6, ↓4.14 log2 FC MFM), higher fatty acid transport enzyme carnitine palmitoyl transferase (CPT1B, ↑3.49 MFM), and lower sarcomere protein tropomyosin (TPM2, ↓3.24 MFM) in MFM vs. control muscle at rest. We propose that in MFM horses, altered cysteine metabolism and a deficiency of cysteine-containing antioxidants combined with a high capacity to oxidize fatty acids and generate ROS during aerobic exercise causes chronic oxidation and aggregation of key proteins such as desmin.