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.

OBJECTIVE To characterize clinical findings for polysaccharide storage myopathy (PSSM) in warmblood horses with type 1 PSSM (PSSM1; caused by mutation of the glycogen synthase 1 gene) and type 2 PSSM (PSSM2; unknown etiology). SAMPLE Database with 3,615 clinical muscle biopsy submissions. PROCEDURES Reported clinical signs and serum creatine kinase (CK) and aspartate aminotransferase (AST) activities were retrospectively analyzed for horses with PSSM1 (16 warmblood and 430 nonwarmblood), horses with PSSM2 (188 warmblood and 646 nonwarmblood), and warmblood horses without PSSM (278). Lameness examinations were reviewed for 9 warmblood horses with PSSM2. Muscle glycogen concentrations were evaluated for horses with PSSM1 (14 warmblood and 6 nonwarmblood), warmblood horses with PSSM2 (13), and horses without PSSM (10 warmblood and 6 nonwarmblood). RESULTS Rhabdomyolysis was more common for horses with PSSM1 (12/16 [75%] warmblood and 223/303 [74%] nonwarmblood) and nonwarmblood horses with PSSM2 (221/436 [51%]) than for warmblood horses with PSSM2 (39/147 [27%]). Gait abnormality was more common in warmblood horses with PSSM2 (97/147 [66%]) than in warmblood horses with PSSM1 (1/16 [7%]), nonwarmblood horses with PSSM2 (176/436 [40%]), and warmblood horses without PSSM (106/200 [53%]). Activities of CK and AST were similar in warmblood horses with and without PSSM2. Muscle glycogen concentrations in warmblood and nonwarmblood horses with PSSM1 were significantly higher than concentrations in warmblood horses with PSSM2. CONCLUSIONS AND CLINICIAL RELEVANCE Rhabdomyolysis and elevated muscle glycogen concentration were detected in horses with PSSM1 regardless of breed. Most warmblood horses with PSSM2 had stiffness and gait abnormalities with CK and AST activities and muscle glycogen concentrations within reference limits.

Background Myofibrillar myopathy (MFM) of unknown aetiology has recently been identified in Warmblood (WB) horses. In humans, 16 genes have been implicated in various MFM‐like disorders. Objectives To identify variants in 16 MFM candidate genes and compare allele frequencies of all variants between MFM WB and non‐MFM WB and coding variants with moderate or severe predicted effects in MFM WB with publicly available data of other breeds. To compare differential gene expression and muscle fibre contractile force between MFM and non‐MFM WB. Study design Case‐control. Animals 8 MFM WB, 8 non‐MFM WB, 33 other WB, 32 Thoroughbreds, 80 Quarter Horses and 77 horses of other breeds in public databases. Methods Variants were called within transcripts of 16 candidate genes using gluteal muscle mRNA sequences aligned to EquCab3.0 and allele frequencies compared by Fisher's exact test among MFM WB, non‐MFM WB and public sequences across breeds. Candidate gene differential expression was determined between MFM and non‐MFM WB by fitting a negative binomial generalised log‐linear model per gene (false discovery rate <0.05). The maximal isometric force/cross‐sectional area generated by isolated membrane‐permeabilised muscle fibres was determined. Results None of the 426 variants identified in 16 candidate genes were associated with MFM including 26 missense variants. Breed‐specific differences existed in allele frequencies. Candidate gene differential expression and muscle fibre‐specific force did not differ between MFM WB (143.1 ± 34.7 kPa) and non‐MFM WB (140.2 ± 43.7 kPa) (P = .8). Main limitations RNA‐seq–only assays transcripts expressed in skeletal muscle. Other possible candidate genes were not evaluated. Conclusions Evidence for association of variants with a disease is essential because coding sequence variants are common in the equine genome. Variants identified in MFM candidate genes, including two coding variants offered as commercial MFM equine genetic tests, did not associate with the WB MFM phenotype.

Background Myofibrillar myopathy in humans causes protein aggregation, degeneration, and weakness of skeletal muscle. In horses, myofibrillar myopathy is a late-onset disease of unknown origin characterized by poor performance, atrophy, myofibrillar disarray, and desmin aggregation in skeletal muscle. This study evaluated molecular and ultrastructural signatures of myofibrillar myopathy in Warmblood horses through gluteal muscle tandem-mass-tag quantitative proteomics (5 affected, 4 control), mRNA-sequencing (8 affected, 8 control), amalgamated gene ontology analyses, and immunofluorescent and electron microscopy. Results We identified 93/1533 proteins and 47/27,690 genes that were significantly differentially expressed. The top significantly differentially expressed protein CSRP3 and three other differentially expressed proteins, including, PDLIM3, SYNPO2, and SYNPOL2, are integrally involved in Z-disc signaling, gene transcription and subsequently sarcomere integrity. Through immunofluorescent staining, both desmin aggregates and CSRP3 were localized to type 2A fibers. The highest differentially expressed gene CHAC1, whose protein product degrades glutathione, is associated with oxidative stress and apoptosis. Amalgamated transcriptomic and proteomic gene ontology analyses identified 3 enriched cellular locations; the sarcomere (Z-disc & I-band), mitochondrial complex I and the extracellular matrix which corresponded to ultrastructural Z-disc disruption and mitochondrial cristae alterations found with electron microscopy. Conclusions A combined proteomic and transcriptomic analysis highlighted three enriched cellular locations that correspond with MFM ultrastructural pathology in Warmblood horses. Aberrant Z-disc mechano-signaling, impaired Z-disc stability, decreased mitochondrial complex I expression, and a pro-oxidative cellular environment are hypothesized to contribute to the development of myofibrillar myopathy in Warmblood horses. These molecular signatures may provide further insight into diagnostic biomarkers, treatments, and the underlying pathophysiology of MFM.