Phenotypic Characterization of Recombinant Marek's Disease Virus in Live Birds Validates Polymorphisms Associated with Virulence.

Marek's disease (MD) is a highly infectious lymphoproliferative disease in chickens with a significant economic impact. Mardivirus gallidalpha 2, also known as Marek's disease virus (MDV), is the causative pathogen and has been categorized based on its virulence rank into four pathotypes: mild (m), virulent (v), very virulent (vv), and very virulent plus (vv+). A prior comparative genomics study suggested that several single-nucleotide polymorphisms (SNPs) and genes in the MDV genome are associated with virulence, including nonsynonymous (ns) SNPs in eight open reading frames (ORF): UL22, UL36, UL37, UL41, UL43, R-LORF8, R-LORF7, and ICP4. To validate the contribution of these nsSNPs to virulence, the vv+MDV strain 686 genome was modified by replacing nucleotides with those observed in the vMDV strains. Pathogenicity studies indicated that these substitutions reduced the MD incidence and increased the survival of challenged birds. Furthermore, using the best-fit pathotyping method to rank the virulence, the modified vv+MDV 686 viruses resulted in a pathotype similar to the vvMDV Md5 strain. Thus, these results support our hypothesis that SNPs in one or more of these ORFs are associated with virulence but, as a group, are not sufficient to result in a vMDV pathotype, suggesting that there are additional variants in the MDV genome associated with virulence, which is not surprising given this complex phenotype and our previous finding of additional variants and SNPs associated with virulence.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants.

BACKGROUND: The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS: We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS: We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS: We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Novel Expression of GLUT3, GLUT6 and GLUT10 in Equine Gluteal Muscle Following Glycogen-Depleting Exercise: Impact of Dietary Starch and Fat.

Horses have a slow rate of muscle glycogen repletion relative to other species for unknown reasons. Our aim was to determine the expression of glucose transporters (GLUT) and genes impacting GLUT4 expression and translocation in the gluteal muscle. Five fit Thoroughbred horses performed glycogen-depleting exercises on high-starch (HS, 2869 g starch/day) and low-starch, high-fat diets (LS-HF, 358 g starch/d) with gluteal muscle biopsies obtained before and after depletion and during repletion. Muscle glycogen declined by ≈30% on both diets with little increase during repletion on LS-HF. Transcriptomic analysis identified differential expression (DE) of only 2/12 genes impacting GLUT4 translocation (two subunits of AMP protein kinase) and only at depletion on LS-HF. Only 1/13 genes encoding proteins that promote GLUT4 transcription had increased DE (PPARGC1A at depletion LS-HF). GLUT4 comprised ≈30% of total GLUT mRNA expression at rest. Remarkably, by 72 h of repletion expression of GLUT3, GLUT6 and GLUT10 increased to ≈25% of total GLUT mRNA. Expression of GLUT6 and GLUT10 lagged from 24 h of repletion on HS to 72 h on LS-HF. Lacking an increase in GLUT4 gene expression in response to glycogen-depleting exercise, equine muscle increases GLUT3, GLUT6 and GLUT10 expression potentially to enhance glucose transport, resembling responses observed in resistance trained GLUT4-null mice.

Mapping splice QTLs reveals distinct transcriptional and post-transcriptional regulatory variation of gene expression and identifies putative alternative splicing variation mediating complex trait variation in pigs.

BACKGROUND: Alternative splicing is an important step in gene expression, generating multiple isoforms for the same genes and greatly expanding the diversity of proteomes. Genetic variation in alternative splicing contributes to phenotypic diversity in natural populations. However, the genetic basis of variation in alternative splicing in livestock including pigs remains poorly understood. RESULTS: In this study, using a Duroc x Pietrain F2 pig population, we performed genome-wide analysis of alternative splicing estimated from stranded RNA-Seq data in skeletal muscle. We characterized the genetic architecture of alternative splicing and compared its basic features with those of overall gene expression. We detected a large number of novel alternative splicing events that were not previously annotated. We found heritability of quantitative alternative splicing scores (percent spliced in or PSI) to be lower than that of overall gene expression. In addition, heritabilities showed little correlation between alternative splicing and overall gene expression. We mapped expression QTLs (eQTLs) and splice QTLs (sQTLs) and found them to be largely non-overlapping. Finally, we integrated sQTL mapping with phenotype QTL (pQTL mapping to identify potential mediator of pQTL effect by alternative splicing. CONCLUSIONS: Our results suggest that regulatory variation exists at multiple levels and that their genetic controls are distinct, offering opportunities for genetic improvement.

Absence of myofibrillar myopathy in Quarter Horses with a histopathological diagnosis of type 2 polysaccharide storage myopathy and lack of association with commercial genetic tests.

BACKGROUND: Genetic tests for variants in MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323 or P3b; rs1142918816) and MYOZ3 (P4; rs1142544043) genes are offered commercially to diagnose myofibrillar myopathy (MFM) and type 2 polysaccharide storage myopathy (PSSM2) in Quarter Horses (QH). OBJECTIVES: To determine if PSSM2-QH has histopathological features of MFM. To compare genotype and allele frequencies of variants P2, P3, P4 between control-QH and PSSM2-QH diagnosed by histopathology. STUDY DESIGN: Retrospective cross-sectional. METHODS: The study includes a total of 229 healthy control-QH, 163 PSSM2-QH GYS1 mutation negative. Desmin stains of gluteal/semimembranosus muscle were evaluated. Purported disease alleles P2, P3a, P3b, P4 were genotyped by pyrosequencing. Genotype, allele frequency and total number of variant alleles or loci were compared between phenotypes using additive/genotypic and dominant models and quantitative effects evaluated by multivariable logistic regression. RESULTS: Histopathological features of MFM were absent in all QH. A P variant allele at any locus was not associated (P > .05) with a histopathological diagnosis of PSSM2 and one or more P variants were common in control-QH (57%) and PSSM2-QH (61%). Allele frequencies (control/PSSM2) were: 0.24/0.21 (P2), 0.07/0.12 (P3a), 0.07/0.11 (P3b) and 0.06/0.08 (P4). P3a and P3b loci were not independent (r2  = 0.894); and not associated with PSSM2 histopathology comparing the haplotype of both P3a and P3b variants to other haplotypes. A receiver operator curve did not accurately predict the PSSM2 phenotype (AUC = 0.67, 95% CI 0.62-0.72), and there was no difference in the total number of variant loci or total variant allele count between control-QH and PSSM2-QH. MAIN LIMITATIONS: P3a and P3b were not in complete linkage disequilibrium. CONCLUSIONS: The P2, P3 and P4 variants in genes associated with human MFM were not associated with PSSM2 in 392 QH. Their use would improperly diagnose PSSM2/MFM in 57% of healthy QH and fail to diagnose PSSM2 in 40% of QH with histopathological evidence of PSSM2.

CONTEXTO: Testes genéticos para detecção das mutações MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323 ou P3b; rs1142918816) e MYOZ3 (P4; rs1142544043) são oferecidos comercialmente para diagnosticar miopatia miofibrilar (MMF) e miopatia por acúmulo de polissacarídeo tipo 2 (PSSM2) em cavalos Quarto de Milha (QM). HIPÓTESES/OBJETIVOS: Determinar se PSSM2-QM tem características similares à MMF. Comparar o genótipo e a frequência dos alelos variantes P2, P3, e P4 entre cavalos QM controle, e PSSM2-QM diagnosticados por histologia. MÉTODOS: 229 cavalos QH saudáveis como controle, e 163 PSSM2-QM positivos na histologia e negativos para a mutação GYS1. METODOLOGIA: Amostras dos músculos glúteo/semimembranoso foram avaliadas após coloração com desmina. Os pretensos genes alelos P2, P3a, P3b e P4 foram genotipados por pirosequenciamento. Genótipo, frequência alélica, e número total de variância alélica ou loci foram comparados entre os fenótipos usado aditivo/genotípico e modelos dominantes e efeitos quantitativos através de regressão logística multivariável. RESULTADOS: Características histopatológicas de MMF não foram encontradas em nenhum QM. Uma variante alélica P em qualquer uma dos loci não foi associada (P > .05) com o diagnóstico histopatológicos de PSSM2 e uma ou mais variante P foram comuns em QM controles (57%) e PSSM2-QM (61%). Frequência alélica (controle/PSSM2) foram: 0.24/0.21 (P2), 0.07/0.12 (P3a), 0.07/0.11 (P3b), e 0.06/0.08 (P4). P3a e P3b loci não foram independentes (r2  = 0.894); e não foram associados com achados histopatológicos de PSSM quando comparando o haplótipo de ambas as variantes P3a e P3b com os outros haplótipos. A curva característica de operação do receptor não previu acuradamente o fenótipo PSSM2 (AUC = 0.67, 95% IC 0.62-0.72), e não houve diferença no número dotal de variantes no loci ou na contagem de variantes alélicas total entre QM controles e PSSM2-QM. PRINCIPAIS LIMITAÇÕES: P3a e P3b não estavam em desequilíbrio de ligação. CONCLUSÕES: As variantes P2, P3 e P4 em genes associados com MMF em humanos não foram associadas com PSSM em 392 QM. O seu uso diagnosticaria impropriamente PSSM2 e MMF em 57% dos cavalos saudáveis utilizados como controle e não diagnosticaria PSSM2 em 40% dos QM com evidência histológica de PSSM2.

Type 2 polysaccharide storage myopathy in Quarter Horses is a novel glycogen storage disease causing exertional rhabdomyolysis.

BACKGROUND: Both type 1 (PSSM1) and type 2 polysaccharide storage myopathy (PSSM2) are characterised by aggregates of abnormal polysaccharide in skeletal muscle. Whereas the genetic basis for PSSM1 is known (R309H GYS1), the cause of PSSM2 in Quarter Horses (PSSM2-QH) is unknown and glycogen concentrations not defined. OBJECTIVES: To characterise the histopathological and biochemical features of PSSM2-QH and determine if an associated monogenic variant exists in genes known to cause glycogenosis. STUDY DESIGN: Retrospective case control. METHODS: Sixty-four PSSM2-QH, 30 PSSM1-QH and 185 control-QH were identified from a biopsy repository and clinical data, histopathology scores (0-3), glycogen concentrations and selected glycolytic enzyme activities compared. Coding sequences of 12 genes associated with muscle glycogenoses were identified from whole genome sequences and compared between seven PSSM2-QH and five control-QH. RESULTS: Exertional rhabdomyolysis in PSSM2-QH occurred predominantly in barrel racing and working cow/roping performance types and improved with regular exercise and a low starch/fat-supplemented diet. Histopathological scores, including the amount of amylase-resistant polysaccharide (PSSM2-QH 1.4 ± 0.6, PSSM1-QH 2.1 ± 0.3, control-QH 0 ± 0, p < 0.001), and glycogen concentrations (PSSM2-QH 129 ± 62, PSSM1-QH 175 ± 9, control-QH 80 ± 27 mmol/kg, p < 0.0001) were intermediate in PSSM2-QH with significant differences among groups. In PSSM2-QH, abnormal polysaccharide had a less filamentous ultrastructure than PSSM1-QH and phosphorylase and phosphofructokinase activities were normal. Seventeen of 30 PSSM2-QH with available pedigrees descended from one of three stallions within four generations. Of the 29 predicted high or moderate impact genetic variants identified in candidate genes, none were present in only PSSM2-QH and absent in control-QH. MAIN LIMITATIONS: Analyses of PSSM2-QH and PSSM1-QH were performed on shipped samples, controls on frozen samples. CONCLUSIONS: PSSM2-QH is a novel glycogen storage disorder that is not the result of a mutation in genes currently known to cause muscle glycogenoses in other species.

CONTEXTO: Ambos os tipos 1 e 2 de miopatia por acúmulo de polissacarídeo (PSSM) são caracterizados por agregados de polissacarídeos anormais no músculo esquelético. Enquanto a base genética do PSSM 1 é conhecida (R309H GYS1), a causa do PSSM2 em cavalos Quarto de Milha (PSSM2-QH) é desconhecida, e a concentração de glicogênio não é definida. OBJETIVOS: Identificar as características histopatológicas e bioquímicas do PSSM-QH e determinar se há uma variante monogênica em genes conhecidos por causar glicogenose. DELINEAMENTO DO ESTUDO: Caso controlado retrospectivo. METODOLOGIA: 64 PSSM2-QH, 30 PSSM1-QH e 185 QH controles foram identificados em um arquivo de dados. Informação clínica, achados histológicos (escala 0-3), concentração de glicogênio e atividade enzimática de algumas enzimas glicolíticas foram comparadas. Sequências codificadas de 12 genes associados com glicogenose muscular foram identificados nas sequências genômicas completas, e comparadas entre 7 PSSM2-QH e 5 QH controles. RESULTADOS: Rabdomiólise por exercício em PSSM2-QH ocorreu predominantemente em cavalos de corrida de tambor e cavalos de team roping/trabalho com gado, e melhorou com exercício regular e uma dieta com baixo amido e alta gordura. A escala histopatológica, incluindo a quantidade de polissacarídeos resistentes à amilase (PSSM2-QH 1.4 ± 0.6, PSSM1-QH 2.1 ± 0.3, controle-QH 0 ± 0, P < 0.001), e concentrações de glicogênio (PSSM2-QH 129 ± 62, PSSM1-QH 175 ± 9, controle-QH 80 ± 27 mmol/kg, P < 0.0001) foram intermediárias em PSSM2-QH com diferença significante entre grupos. Em PSSM2-QH, polissacarídeo anormal teve uma ultraestrutura menos filamentosa do que PSSM1-QH e as atividades de fosforilase e fosfofrutoquinase foram normais. Dezessete dos 30 PSSM2-QH com pedigree disponível descendiam de 1 de 3 garanhões dentro de 4 gerações. Das 29 variações genéticas preditas a terem impacto moderado ou alto como genes candidatos, nenhuma estava presente apenas em PSSM2-QH e ausente no grupo controle-QH. PRINCIPAIS LIMITAÇÕES: As análises feitas nas amostras de PSSM2-QH e PSSM1-QH foram realizadas em amostras enviadas por correio, e as amostras dos animais controles eram amostras congeladas. CONCLUSÕES: PSSM2-QH é uma nova doença por acúmulo de glicogênio que não é o resultado de uma mutação nos genes conhecidos por causarem glicogenose muscular em outras espécies.

Enriched Pathways of Calcium Regulation, Cellular/Oxidative Stress, Inflammation, and Cell Proliferation Characterize Gluteal Muscle of Standardbred Horses between Episodes of Recurrent Exertional Rhabdomyolysis.

Certain Standardbred racehorses develop recurrent exertional rhabdomyolysis (RER-STD) for unknown reasons. We compared gluteal muscle histopathology and gene/protein expression between Standardbreds with a history of, but not currently experiencing rhabdomyolysis (N = 9), and race-trained controls (N = 7). Eight RER-STD had a few mature fibers with small internalized myonuclei, one out of nine had histologic evidence of regeneration and zero out of nine degeneration. However, RER-STD versus controls had 791/13,531 differentially expressed genes (DEG). The top three gene ontology (GO) enriched pathways for upregulated DEG (N = 433) were inflammation/immune response (62 GO terms), cell proliferation (31 GO terms), and hypoxia/oxidative stress (31 GO terms). Calcium ion regulation (39 GO terms), purine nucleotide metabolism (32 GO terms), and electron transport (29 GO terms) were the top three enriched GO pathways for down-regulated DEG (N = 305). DEG regulated RYR1 and sarcoplasmic reticulum calcium stores. Differentially expressed proteins (DEP ↑N = 50, ↓N = 12) involved the sarcomere (24% of DEP), electron transport (23%), metabolism (20%), inflammation (6%), cell/oxidative stress (7%), and other (17%). DEP included ↑superoxide dismutase, ↑catalase, and DEP/DEG included several cysteine-based antioxidants. In conclusion, gluteal muscle of RER-susceptible Standardbreds is characterized by perturbation of pathways for calcium regulation, cellular/oxidative stress, inflammation, and cellular regeneration weeks after an episode of rhabdomyolysis that could represent therapeutic targets.

The Impact of N-Acetyl Cysteine and Coenzyme Q10 Supplementation on Skeletal Muscle Antioxidants and Proteome in Fit Thoroughbred Horses.

Horses have one of the highest skeletal muscle oxidative capacities amongst mammals, which, combined with a high glycolytic capacity, could perturb redox status during maximal exercise. We determined the effect of 30 d of oral coenzyme Q10 and N-acetyl-cysteine supplementation (NACQ) on muscle glutathione (GSH), cysteine, ROS, and coenzyme Q10 concentrations, and the muscle proteome, in seven maximally exercising Thoroughbred horses using a placebo and randomized cross-over design. Gluteal muscle biopsies were obtained the day before and 1 h after maximal exercise. Concentrations of GSH, cysteine, coenzyme Q10, and ROS were measured, and citrate synthase, glutathione peroxidase, and superoxide dismutase activities analyzed. GSH increased significantly 1 h post-exercise in the NACQ group (p = 0.022), whereas other antioxidant concentrations/activities were unchanged. TMT proteomic analysis revealed 40 differentially expressed proteins with NACQ out of 387 identified, including upregulation of 13 mitochondrial proteins (TCA cycle and NADPH production), 4 Z-disc proteins, and down regulation of 9 glycolytic proteins. NACQ supplementation significantly impacted muscle redox capacity after intense exercise by enhancing muscle glutathione concentrations and increasing expression of proteins involved in the uptake of glutathione into mitochondria and the NAPDH-associated reduction of oxidized glutathione, without any evident detrimental effects on performance.

Integrated proteomic and transcriptomic profiling identifies aberrant gene and protein expression in the sarcomere, mitochondrial complex I, and the extracellular matrix in Warmblood horses with myofibrillar myopathy.

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.

Integrated Genome-Wide Analysis of MicroRNA Expression Quantitative Trait Loci in Pig Longissimus Dorsi Muscle.

Determining mechanisms regulating complex traits in pigs is essential to improve the production efficiency of this globally important protein source. MicroRNAs (miRNAs) are a class of non-coding RNAs known to post-transcriptionally regulate gene expression affecting numerous phenotypes, including those important to the pig industry. To facilitate a more comprehensive understanding of the regulatory mechanisms controlling growth, carcass composition, and meat quality phenotypes in pigs, we integrated miRNA and gene expression data from longissimus dorsi muscle samples with genotypic and phenotypic data from the same animals. We identified 23 miRNA expression Quantitative Trait Loci (miR-eQTL) at the genome-wide level and examined their potential effects on these important production phenotypes through miRNA target prediction, correlation, and colocalization analyses. One miR-eQTL miRNA, miR-874, has target genes that colocalize with phenotypic QTL for 12 production traits across the genome including backfat thickness, dressing percentage, muscle pH at 24 h post-mortem, and cook yield. The results of our study reveal genomic regions underlying variation in miRNA expression and identify miRNAs and genes for future validation of their regulatory effects on traits of economic importance to the global pig industry.

Pathways of calcium regulation, electron transport, and mitochondrial protein translation are molecular signatures of susceptibility to recurrent exertional rhabdomyolysis in Thoroughbred racehorses.

Recurrent exertional rhabdomyolysis (RER) is a chronic muscle disorder of unknown etiology in racehorses. A potential role of intramuscular calcium (Ca2+) dysregulation in RER has led to the use of dantrolene to prevent episodes of rhabdomyolysis. We examined differentially expressed proteins (DEP) and gene transcripts (DEG) in gluteal muscle of Thoroughbred race-trained mares after exercise among three groups of 5 horses each; 1) horses susceptible to, but not currently experiencing rhabdomyolysis, 2) healthy horses with no history of RER (control), 3) RER-susceptible horses treated with dantrolene pre-exercise (RER-D). Tandem mass tag LC/MS/MS quantitative proteomics and RNA-seq analysis (FDR <0.05) was followed by gene ontology (GO) and semantic similarity of enrichment terms. Of the 375 proteins expressed, 125 were DEP in RER-susceptible versus control, with 52 ↑DEP mainly involving Ca2+ regulation (N = 11) (e.g. RYR1, calmodulin, calsequestrin, calpain), protein degradation (N = 6), antioxidants (N = 4), plasma membranes (N = 3), glyco(geno)lysis (N = 3) and 21 DEP being blood-borne. ↓DEP (N = 73) were largely mitochondrial (N = 45) impacting the electron transport system (28), enzymes (6), heat shock proteins (4), and contractile proteins (12) including Ca2+ binding proteins. There were 812 DEG in RER-susceptible versus control involving the electron transfer system, the mitochondrial transcription/translational response and notably the pro-apoptotic Ca2+-activated mitochondrial membrane transition pore (SLC25A27, BAX, ATP5 subunits). Upregulated mitochondrial DEG frequently had downregulation of their encoded DEP with semantic similarities highlighting signaling mechanisms regulating mitochondrial protein translation. RER-susceptible horses treated with dantrolene, which slows sarcoplasmic reticulum Ca2+ release, showed no DEG compared to control horses. We conclude that RER-susceptibility is associated with alterations in proteins, genes and pathways impacting myoplasmic Ca2+ regulation, the mitochondrion and protein degradation with opposing effects on mitochondrial transcriptional/translational responses and mitochondrial protein content. RER could potentially arise from excessive sarcoplasmic reticulum Ca2+ release and subsequent mitochondrial buffering of excessive myoplasmic Ca2+.

Commercial genetic testing for type 2 polysaccharide storage myopathy and myofibrillar myopathy does not correspond to a histopathological diagnosis.

BACKGROUND: Commercial genetic tests for type 2 polysaccharide storage myopathy (PSSM2) and myofibrillar myopathy (MFM) have not been validated by peer-review, and formal regulation of veterinary genetic testing is lacking. OBJECTIVES: To compare genotype and allele frequencies of commercial test variants (P variants) in MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323), FLNC (P3b; rs1142918816) and MYOZ3 (P4; rs1142544043) between Warmblood (WB) and Arabian (AR) horses diagnosed with PSSM2/MFM by muscle histopathology, and phenotyped breed-matched controls. To quantify variant frequency in public repositories of ancient and modern horse breeds. STUDY DESIGN: Cross sectional using archived clinical material and publicly available data. METHODS: We studied 54 control-WB, 68 PSSM2/MFM-WB, 30 control-AR, 30 PSSM2/MFM-AR and 205 public genotypes. Variants were genotyped by pyrosequencing archived DNA. Genotype and allele frequency, and number of variant alleles or loci were compared within breed between controls, PSSM2/MFM combined and MFM or PSSM2 horses considered separately using additive/genotypic and dominant models (Fisher's exact tests). Variant frequencies in modern, early domestic and Przewalski horses were determined from a public data repository. RESULTS: There was no significant association between any P locus and a histopathological diagnosis of PSSM2/MFM, and no difference between control and myopathic horses in total loci with alternative alleles, or total alternate alleles when PSSM2/MFM was considered combined or separately as PSSM2 or MFM. For all tests, sensitivity was <0.33. Allele frequencies in WB (controls/cases) were: 8%/15% (P2), 5%/6% (P3a/b) and 9%/13% (P4); in AR, frequencies were: 12%/17% (P2), 2%/2% (P3a/b) and 7%/12% (P4). All P variants were present in early domestic (400- to 5500-year-old) horses and P2 present in the Przewalski. CONCLUSIONS: Because of the lack of significant association between a histopathological diagnosis of PSSM2 or MFM and the commercial genetic test variants P2, P3 and P4 in WB and AR, we cannot recommend the use of these variant genotypes for selection and breeding, prepurchase examination or diagnosis of a myopathy.

Candidate gene expression and coding sequence variants in Warmblood horses with myofibrillar myopathy.

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.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs.

BACKGROUND: Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs. RESULTS: Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat. CONCLUSION: Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Profiling and characterization of a longissimus dorsi muscle microRNA dataset from an F2 Duroc × Pietrain pig resource population.

To elucidate the effects of microRNA (miRNA) regulation in skeletal muscle of adult pigs, miRNA expression profiling was performed with RNA extracted from longissimus dorsi (LD) muscle samples from 174 F2 pigs (~ 5.5 months of age) from a Duroc × Pietrain resource population. Total RNA was extracted from LD samples, and libraries were sequenced on an Illumina HiSeq 2500 platform in 1 × 50 bp format. After processing, 232,826,977 total reads were aligned to the Sus scrofa reference genome (v10.2.79), with 74.8% of total reads mapping successfully. The miRDeep2 software package was utilized to quantify annotated Sus scrofa mature miRNAs from miRBase (Release 21) and to predict candidate novel miRNA precursors. Among the retained 295 normalized mature miRNA expression profiles ssc-miR-1, ssc-miR-133a-3p, ssc-miR-378, ssc-miR-206, and ssc-miR-10b were the most abundant, all of which have previously been shown to be expressed in pig skeletal muscle. Additionally, 27 unique candidate novel miRNA precursors were identified exhibiting homologous sequence to annotated human miRNAs. The composition of classes of small RNA present in this dataset was also characterized; while the majority of unique expressed sequence tags were not annotated in any of the queried databases, the most abundantly expressed class of small RNA in this dataset was miRNAs. This data provides a resource to evaluate miRNA regulation of gene expression and effects on complex trait phenotypes in adult pig skeletal muscle. The raw sequencing data were deposited in the Sequence Read Archive, BioProject PRJNA363073.