A defect in dystrophin causes a novel porcine stress syndrome
© Nonneman et al.; licensee BioMed Central Ltd. 2012
Received: 22 September 2011
Accepted: 25 May 2012
Published: 12 June 2012
Skip to main content
© Nonneman et al.; licensee BioMed Central Ltd. 2012
Received: 22 September 2011
Accepted: 25 May 2012
Published: 12 June 2012
Losses of slaughter-weight pigs due to transport stress are both welfare and economic concerns to pork producers. Historically, the HAL-1843 mutation in ryanodine receptor 1 was considered responsible for most of the losses; however, DNA testing has effectively eliminated this mutation from commercial herds. We identified two sibling barrows in the USMARC swine herd that died from apparent symptoms of a stress syndrome after transport at 12 weeks of age. The symptoms included open-mouth breathing, skin discoloration, vocalization and loss of mobility.
We repeated the original mating along with sire-daughter matings to produce additional offspring. At 8 weeks of age, heart rate and electrocardiographs (ECG) were monitored during isoflurane anesthesia challenge (3% for 3 min). Four males from the original sire-dam mating and two males from a sire-daughter mating died after one minute of anesthesia. Animals from additional litters were identified as having a stress response, sometimes resulting in death, during regular processing and weighing. Affected animals had elevated plasma creatine phosphokinase (CPK) levels before and immediately after isoflurane challenge and cardiac arrhythmias. A pedigree containing 250 pigs, including 49 affected animals, was genotyped with the Illumina PorcineSNP60 Beadchip and only one chromosomal region, SSCX at 25.1-27.7 Mb over the dystrophin gene (DMD), was significantly associated with the syndrome. An arginine to tryptophan (R1958W) polymorphism in exon 41 of DMD was the most significant marker associated with stress susceptibility. Immunoblots of affected heart and skeletal muscle showed a dramatic reduction of dystrophin protein and histopathology of affected hearts indicated muscle fiber degeneration.
A novel stress syndrome was characterized in pigs and the causative genetic factor most likely resides within DMD that results in less dystrophin protein and cardiac abnormalities that can lead to death under stressful conditions. The identification of predictive markers will allow us to determine the prevalence of this disease in commercial swine populations. This defect also provides a unique biomedical model for human cardiomyopathy associated with muscular dystrophy that may be superior to those available because of the similarities in anatomy and physiology and allow advances in gene therapies for human disease.
Losses of market-weight pigs present economic, legal and animal welfare issues to U.S. swine producers and include dead, non-ambulatory, fatigued and injured pigs . Although the incidence occurs at low frequency (less than 1%) these losses are substantial due to added labor and disposal costs, along with loss of full-value product and have been estimated to be over $50,000,000 per year in the U.S. Fatigued, non-ambulatory pigs exhibit acute stress symptoms including open-mouth breathing, discoloration (red to purple) and blotching of skin, muscle tremors, abnormal vocalization and refusal to move, similar to what is observed in pigs with malignant hyperthermia due to a mutation in the ryanodine receptor 1 gene (RYR1) . A DNA-based test has been available since the discovery of the mutation and producers have been able eliminate the unfavorable RYR1 allele from their herds . Over a decade ago it was found that half of the pigs arriving dead or dying at the packing plants carried at least one copy of the RYR1 mutation . By 2006 the number of dead and non-ambulatory pigs that carried the RYR1 mutation was about 5% ; thus, the mutation was still present at a low frequency at that time . However, a high proportion of RYR1 normal pigs show a sensitivity to halothane anesthesia  and are more prone to becoming non-ambulatory after handling . Pigs that are more sensitive to halothane exposure may also have inferior pork quality . These relationships are believed to have a genetic basis; however, the specific cause has yet to be identified.
In the U.S. Meat Animal Research Center’s swine research population, a novel stress syndrome was detected. Test matings were made to characterize the syndrome’s physiological effects and its genetic basis. This study reports the identification of a novel porcine stress syndrome and mapping of the defect to dystrophin.
Plasma creatine phosphokinase activity in normal and stress-responsive pigs
Normal (n = 103)b
Carrier (n = 37)
Affected (n = 35)c
482.42 ± 41.30
573.93 ± 89.14
1635.27 ± 148.22
644.37 ± 68.80
469.26 ± 51.08
1928.46 ± 136.37
The exons and flanking intron sequences of DMD were obtained from normal and affected pigs by PCR amplification and direct sequencing. No polymorphisms were identified near splice sites. The cDNA from normal and affected animals was also amplified and sequenced and there was no evidence of alternative splicing or deletion of exons in the affected animals. The promoter regions P1, P2 and the muscle-specific promoter were sequenced and no obvious mutations were found. Forty-seven SNPs were found in DMD by sequencing 12 members of the affected family; of these, five non-synonymous polymorphisms were found in the coding region (Additional file 2: Table S2). One SNP in particular, 85890_783, that causes the amino acid change arginine to tryptophan at amino acid 1958 (R1958W) in exon 41 was predicted by PolyPhen-2  (http://genetics.bwh.harvard.edu/pph2/) to be damaging with a probability score of 0.983.
Associations of SNP in dystrophin with the stress response
We have identified a porcine stress syndrome that can be induced by transport, handling or isoflurane anesthesia. Affected animals quickly become non-ambulatory, have difficulty breathing and rarely recover. Isoflurane anesthesia was tested to elicit a stress response because we thought the stress syndrome was a malignant hyperthermia, similar to the classic porcine stress syndrome. Isoflurane is more readily available, less irritating to the respiratory system and has replaced halothane as a common anesthesia for veterinary applications. Although the response is delayed, compared to halothane, isoflurane can also induce malignant hyperthermia in swine . Using the Illumina PorcineSNP60 Beadchip we were able to map the defect to dystrophin on the X chromosome and show that affected animals have reduced dystrophin protein in heart and skeletal muscle and elevated CPK levels in blood. These animals, however, have no apparent muscle impairment later in life. Dystrophin is the largest gene in the mammalian genome covering 2.4 megabases of DNA, and contains 79 exons and multiple tissue-specific promoters and transcripts . Mutations in dystrophin cause Duchennes Muscular Dystrophy (DMD), Becker Muscular Dystrophy (BMD) and X-linked dilated cardiomyopathy (XLCM) in humans . Most cases of DMD or BMD involve deletions/duplications, followed by nonsense mutations and/or microdeletions/insertions, most of which cause a frameshift or exon skipping; missense mutations are extremely rare [13–16]. DMD is usually caused by deletions, duplications or nonsense mutations that disrupt or truncate the reading frame, while BMD mutations usually preserve the reading frame [13, 15]. These allelic forms are generally characterized by the severity and age of onset of the disease with DMD manifesting within the first few years of life as a progressive weakness and wasting of muscle, while BMD occurs in later decades and is more variable in severity. Patients with either DMD or BMD usually develop dilated cardiomyopathy independent of the degree of skeletal muscle involvement  and usually have higher than normal circulating levels of CPK [15, 18]. Likewise, patients with XLCM [19, 20], also have elevated CPK, probably due to increased membrane permeability , because the absence of dystrophin makes the membrane more susceptible to tears under the force of contraction. Cardiomyopathy that occurs without skeletal muscle involvement is frequently caused by missense mutations in the dystrophin gene . Dystrophinopathies can also present with arrhythmia  or cardiorespiratory arrest after isoflurane anesthesia [23, 24]. The pathophysiological reasons for anesthesia-induced cardiac arrest in patients with dystrophinopathies are unknown, however it is postulated that damaged, leaky muscle fibers (rhabdomyolysis) can lead to increased intracellular calcium and potassium released into the blood (hyperkalemia) and cardiac arrest [21, 24]. It has been reported that isoflurane anesthesia or stress induced rhabdomyolosis, hyperkalemia and death in dystrophin-deficient cats . This is the only report of stress-induced death in a dystrophic animal.
Dystrophin protein is usually absent by Western blot in DMD patients, while BMD patients usually have a decreased amount of normal or abnormal size protein , which is similar to the phenotype seen in the affected pigs. While it is not clear why missense mutations cause a reduction in dystrophin protein [26, 27] it may be due to damage of the protein by contractile stress.
One of the nonsynonymous SNPs, (85890_783; ss410758971), which causes an amino acid change of arginine to tryptophan at amino acid 1958, was the most significant marker for the stress response and could possibly be the causative mutation. The SNP is in exon 41 of dystrophin, which codes for spectrin repeat 15 of the central coiled-coil rod domain . This basic repeat region is in the second actin-binding domain of dystrophin (ABD2)  and has also been shown to bind lipid . Although in vitro mutation studies have not been reported for ABD2, missense mutations in actin-binding domain 1 (ABD1) can cause protein instability, misfolding and aggregation [31, 32], leading to a reduction of expressed protein.
The variant T allele of 85890_783 (Trp) could be tracked to some of the founding sires and dams in the affected families. The frequency of this allele was higher than expected, based on the rate of spontaneous deaths of feeder-aged pigs in the population. Thus, if this is the causative mutation, the penetrance of this defect is not extremely high. Further studies are required to determine the causal genetic variation, the penetrance of this phenotype and if this locus accounts for losses during transportation of market weight animals.
Genetic mapping indicated the causative variation for a novel porcine stress syndrome is likely in the dystrophin gene. In addition to muscular dystrophies, mutations in human dystrophin can cause dilated cardiomyopathy, rhabdomyolosis and a malignant hyperthermia-like reaction in response to inhaled anesthesia, supporting the conclusion that this locus causes the observed phenotypes in pigs. The identification of the causative mutation in these families will allow investigation of the prevalence of this disease in commercial populations and its pleiotropic effects.
The animals were from a multigenerational Landrace-Duroc-Yorkshire composite population developed at the U.S. Meat Animal Research Center (USMARC) by mating Yorkshire-Landrace composite females (n = 220) to either Duroc or Landrace boars selected from the industry (n = 12 boars of each breed). This population was further developed by mating female descendants (n = 10/boar) to either a Duroc or Landrace-sired boar; Duroc-sired animals were mated to Landrace-sired animals and vice versa. Subsequent matings were random except that half-sib matings were avoided. The first affected animals identified were from the 7th generation and the sire and dam were from the 6th generation of this population. Litters were evaluated through generation 11. All animal procedures were reviewed and approved by the U.S. Meat Animal Research Center Animal Care and Use Committee and procedures for handling pigs complied with those specified in the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching .
We identified 2 barrows in the USMARC swine herd that died after a 1 km transport to a research location at 12 weeks of age. We reproduced the original mating along with sire-daughter matings to produce additional offspring. To induce a reproducible stress response, pigs were challenged with 3% isoflurane anesthesia for 3 min at 8 weeks of age. During the course of this study, pigs from other litters were identified from their reactions to regular processing, weighing or transport to the sale barn at 5–8 weeks of age. Their littermates were also challenged with anesthesia and additional affected animals were identified. A stress response was usually indicated when the pigs stopped breathing after 1–2 minutes of anesthetic, when they were slow to awaken after anesthesia, or displayed rigidity and blotching during handling. Matings were continued with these related offspring to produce additional animals. This produced a final population of 250 animals containing 49 affected pigs for genetic analysis.
Pigs were challenged at 8 weeks of age with isoflurane anesthesia (3% for 3 min) and body temperature, heart rate and ECG were monitored during anesthesia. During the 3 minute isoflurane challenge, ECG wave forms were obtained using a 3-wiring system recorded at 100 Hz using a 12-bit A/D. The signal was conditioned prior to recording using both a differential amplifier and a low pass filter to ensure proper signal. Blood was collected one week before challenge and immediately after isoflurane administration on most animals (n = 192). Blood was drawn into Li-heparin tubes for the measurement of creatine phosphokinase. Plasma creatine phosphokinase was measured using a 2-part reagent system (Pointe Scientific Inc., Canton, MI) in a SpectraMax M5 microplate plate reader (Molecular Devices, Sunnyvale, CA). Twenty-five μl samples were measured in duplicate and the rate of NADH formation was monitored at 340 nm and 37°C. Samples that were greater than 2500 U/L were diluted in PBS and reassayed.
Initially, a pedigree of 58 pigs, including 15 affected animals, was genotyped with the Illumina PorcineSNP60 Beadchip containing 64,232 SNPs and for candidate gene SNP using the Sequenom MASSARRAY® system (Sequenom, San Diego, CA). Additional pigs were phenotyped as they were generated and subsequently genotyped. In total, 250 pigs, including 49 affected animals, their parents and unaffected siblings, were genotyped with the Illumina PorcineSNP60 Beadchip (Illumina, San Diego, CA) . Genotyping of additional SNPs identified in the dystrophin gene and for candidate gene SNP was done using the Sequenom MASSARRAY® system (Sequenom, San Diego, CA). Multiplex assays with about 30 SNP/group were designed using massarray® Assay Design software. Amplicon lengths were approximately 100 bp. Reaction conditions were performed as suggested by Sequenom iPLEX chemistry. Genome-wide associations were performed using the case–control option in PLINK .
Genomic DNA was extracted from blood buffy coats or from tails collected shortly after birth. Exons, flanking intronic regions and promoter regions of DMD were amplified and sequenced from twelve animals with primers (Integrated DNA Technologies, Coralville, IA) designed using Primer 3  in a Dyad DNA engine (Bio-Rad, Hercules, CA) using 0.5 U Hot Star Taq polymerase (Qiagen, Valencia, CA); 1x supplied buffer; 1.5 mM MgCl2; 200 μM dNTPs; 0.8 μM each primer; and 100 ng genomic DNA in 25 μl reactions. The complete mRNA of the dystrophin muscle isoform was amplified and sequenced from cDNA using overlapping fragments and exon primers. Promoter regions were identified using Cister  and the cis-element matrices of TRANSFAC . PCR reactions were treated with 0.1 U exonuclease (USB, Cleveland, OH), and prepared for sequencing on an ABI3730 capillary sequencer (Applied Biosystems, Foster City, CA).
Left ventricular and longissimus lumborum muscle samples from at least twelve affected and twelve age-matched unaffected pigs were frozen in liquid nitrogen within 10 minutes of the time of death. Protein for Western blotting was extracted as reported before  with the addition of protease inhibitor cocktail (Sigma Chemical Co., St. Louis, MO). Dystrophin protein was detected using 15 μg protein from cardiac and skeletal muscle run in 5% polyacrylamide gels using MANDYS8 and MANDRA1 anti-dystrophin monoclonal antibodies (Sigma Chemical Co., St. Louis, MO) diluted 1:400, goat anti-mouse IgG-HRP and SuperSignal® Western Blotting kit (Pierce, Rockford, IL). MANDYS8 recognizes residues 1431–1505  and MANDRA1 recognizes residues 3667–3671 . Equal loading of protein in the wells was confirmed by Coomassie Blue staining.
Total RNA was extracted from 50 mg of porcine heart and longissimus lumborum muscle from 8-week old or 6-month old pigs by homogenization with TriZol reagent (Invitrogen, Carlsbad, CA) following the manufacturer’s protocol. The RNA pellets were resuspended in 100 μL DEPC treated water and absorbance was measured by spectrophotometry at 260 nm. Total RNA (2 μg) was reverse transcribed using oligo dT18 primer and M-MLV reverse transcriptase (Promega, Madison, WI). Primers for DMD (GenBank NM_001012408.1) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH, GenBank AF017079) were designed for real-time RT-PCR using Primer 3. Primers used to amplify DMD were located in exons 4 and 6. Dystrophin forward primer sequence (nucleotides 423 to 444) was 5′- tgcaggtcttgcagaaaaataa -3′ and the reverse primer sequence was 5′- gcaatccagccatgatattctt -3′ (nucleotides 554 to 575). GAPDH forward primer sequence (nucleotides 739 to 757) was 5′-gcgtgaaccatgagaagtatga-3′ and the reverse primer sequence (nucleotides 947 to 967) was 5′-ggtagaagcagggatgatgttc-3′. Real-time PCR was performed using 1X iQ™ SYBR® Green Supermix (Bio-Rad, Hercules, CA), 0.4 ng of cDNA and 0.3 μM of each primer. The PCR reaction was performed at 95°C for 3 min followed by 40 cycles at 95°C for 20s, 58°C for 30s, and 70°C for 1 s on a MJ PTC-200 with a Chromo-4 detector (MJ Research, Watertown, MA). A quantity of 0.4 ng of cDNA was in the linear range of amplification and used for the real-time assays. The threshold cycle (Ct) for DMD and GAPDH of each sample was determined and used to calculate the fold difference between samples (n = 12 affected and 12 normal pigs) using the 2-ΔΔCt method .
The authors thank Sue Hauver, Dale Janssen, Karen Shuck and Kris Simmerman for technical assistance, Sherry Kluver for manuscript preparation and the National Swine Registry for providing samples for genotyping. USDA is an equal opportunity provider and employer. Mention of trade names or commercial products is solely for the purpose of providing information and does not imply recommendation, endorsement or exclusion of other suitable products by the U.S. Department of Agriculture.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.