Use of RNA-seq to identify cardiac genes and gene pathways differentially expressed between dogs with and without dilated cardiomyopathy

Steven G. Friedenberg Department of Clinical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607.
Center for Comparative Medicine and Translational Research, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Steven G. Friedenberg in
Current site
Google Scholar
PubMed
Close
 MBA, DVM, MS
,
Lhoucine Chdid Department of Clinical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Lhoucine Chdid in
Current site
Google Scholar
PubMed
Close
,
Bruce Keene Department of Clinical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607.
Center for Comparative Medicine and Translational Research, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Bruce Keene in
Current site
Google Scholar
PubMed
Close
 DVM
,
Barbara Sherry Department of Molecular Biomedical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607.
Center for Comparative Medicine and Translational Research, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Barbara Sherry in
Current site
Google Scholar
PubMed
Close
 PhD
,
Alison Motsinger-Reif Department of Statistics, North Carolina State University, Raleigh, NC 27607.
Center for Comparative Medicine and Translational Research, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Alison Motsinger-Reif in
Current site
Google Scholar
PubMed
Close
 PhD
, and
Kathryn M. Meurs Department of Clinical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607.
Center for Comparative Medicine and Translational Research, North Carolina State University, Raleigh, NC 27607.

Search for other papers by Kathryn M. Meurs in
Current site
Google Scholar
PubMed
Close
 PhD

Abstract

OBJECTIVE To identify cardiac tissue genes and gene pathways differentially expressed between dogs with and without dilated cardiomyopathy (DCM).

ANIMALS 8 dogs with and 5 dogs without DCM.

PROCEDURES Following euthanasia, samples of left ventricular myocardium were collected from each dog. Total RNA was extracted from tissue samples, and RNA sequencing was performed on each sample. Samples from dogs with and without DCM were grouped to identify genes that were differentially regulated between the 2 populations. Overrepresentation analysis was performed on upregulated and downregulated gene sets to identify altered molecular pathways in dogs with DCM.

RESULTS Genes involved in cellular energy metabolism, especially metabolism of carbohydrates and fats, were significantly downregulated in dogs with DCM. Expression of cardiac structural proteins was also altered in affected dogs.

CONCLUSIONS AND CLINICAL RELEVANCE Results suggested that RNA sequencing may provide important insights into the pathogenesis of DCM in dogs and highlight pathways that should be explored to identify causative mutations and develop novel therapeutic interventions.

Abstract

OBJECTIVE To identify cardiac tissue genes and gene pathways differentially expressed between dogs with and without dilated cardiomyopathy (DCM).

ANIMALS 8 dogs with and 5 dogs without DCM.

PROCEDURES Following euthanasia, samples of left ventricular myocardium were collected from each dog. Total RNA was extracted from tissue samples, and RNA sequencing was performed on each sample. Samples from dogs with and without DCM were grouped to identify genes that were differentially regulated between the 2 populations. Overrepresentation analysis was performed on upregulated and downregulated gene sets to identify altered molecular pathways in dogs with DCM.

RESULTS Genes involved in cellular energy metabolism, especially metabolism of carbohydrates and fats, were significantly downregulated in dogs with DCM. Expression of cardiac structural proteins was also altered in affected dogs.

CONCLUSIONS AND CLINICAL RELEVANCE Results suggested that RNA sequencing may provide important insights into the pathogenesis of DCM in dogs and highlight pathways that should be explored to identify causative mutations and develop novel therapeutic interventions.

Supplementary Materials

    • Supplementary Tables (PDF 108 kb)

Contributor Notes

Address correspondence to Dr. Friedenberg (steven_friedenberg@ncsu.edu).
  • 1. Meurs KM, Mauceli E, Lahmers S, et al. Genome-wide association identifies a mutation in the 3′ untranslated region of striatin, a desmosomal gene, in a canine model of arrhythmogenic right ventricular cardiomyopathy. Hum Genet 2010; 128: 315–324.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 2. Meurs KM, Lahmers S, Keene BW, et al. A splice site mutation in a gene encoding for PDK4, a mitochondrial protein, is associated with the development of dilated cardiomyopathy in the Doberman Pinscher. Hum Genet 2012; 131: 1319–1325.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 3. Meurs KM, Miller MW, Wright NA. Clinical features of dilated cardiomyopathy in Great Danes and results of a pedigree analysis: 17 cases (1990–2000). J Am Vet Med Assoc 2001; 218: 729–732.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 4. Distl O, Vollmar AC, Broschk C, et al. Complex segregation analysis of dilated cardiomyopathy (DCM) in Irish Wolfhounds. Hereditary (Edinb) 2007; 99: 460–465.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 5. Meurs KM. Genetics of cardiac disease in the small animal patient. Vet Clin North Am Small Anim Pract 2010; 40: 701–715.

  • 6. Summerfield NJ, Boswood A, O'Grady MR, et al. Efficacy of pimobendan in the prevention of congestive heart failure or sudden death in Doberman Pinschers with preclinical dilated cardiomyopathy (the PROTECT Study). J Vet Intern Med 2012; 26: 1337–1349.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 7. Khatoon Z, Figler B, Zhang H, et al. Introduction to RNA-seq and its applications to drug discovery and development. Drug Dev Res 2014; 75: 324–330.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 8. Mimura I, Kanki Y, Kodama T, et al. Revolution of nephrology research by deep sequencing: ChIP-seq and RNA-seq. Kidney Int 2014; 85: 31–38.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 9. Eizirik DL, Sammeth M, Bouckenooghe T, et al. The human pancreatic islet transcriptome: expression of candidate genes for type 1 diabetes and the impact of pro-inflammatory cytokines. PLoS Genet 2012; 8: e1002552.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 10. Costa V, Aprile M, Esposito R, et al. RNA-seq and human complex diseases: recent accomplishments and future perspectives. Eur J Hum Genet 2013; 21: 134–142.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 11. Parsch J, Ellegren H. The evolutionary causes and consequences of sex-biased gene expression. Nat Rev Genet 2013; 14: 83–87.

  • 12. Oyama MA, Chittur S. Genomic expression patterns of cardiac tissues from dogs with dilated cardiomyopathy. Am J Vet Res 2005; 66: 1140–1155.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 13. Zhao S, Fung-Leung WP, Bittner A, et al. Comparison of RNA-seq and microarray in transcriptome profiling of activated T cells. PLoS One 2014; 9: e78644.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 14. Wang C, Gong B, Bushel PR, et al. The concordance between RNA-seq and micro array data depends on chemical treatment and transcript abundance. Nat Biotechnol 2014; 32: 926–932.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 15. Hart SN, Therneau TM, Zhang Y, et al. Calculating sample size estimates for RNA sequencing data. J Comput Biol 2013; 20: 970–978.

  • 16. Dukes-McEwan J, Borgarelli M, Tidholm A, et al. Proposed guidelines for the diagnosis of canine idiopathic dilated cardiomyopathy. J Vet Cardiol 2003; 5: 7–19.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 17. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 2011; 17: 10–12.

  • 18. Lindblad-Toh K, Wade CM, Mikkelsen TS, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 2005; 438: 803–819.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 19. Kim D, Pertea G, Trapnell C, et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 2013;14:R36.

    • Search Google Scholar
    • Export Citation
  • 20. Trapnell C, Roberts A, Goff L, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc 2012; 7: 562–578.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 21. Cunningham F, Amode MR, Barrell D, et al. Ensembl 2015. Nucleic Acids Res 2015;43:D662–D669.

  • 22. Li CI, Su PF, Shyr Y. Sample size calculation based on exact test for assessing differential expression analysis in RNA-seq data. BMC Bioinformatics 2013; 14: 357.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 23. Trapnell C, Hendrickson DG, Sauvageau M, et al. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol 2013; 31: 46–53.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 24. Untergasser A, Cutcutache I, Koressaar T, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res 2012;40:e115.

  • 25. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001; 25: 402–408.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 26. Zhang B, Kirov S, Snoddy J. WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 2005;33:W741–W748.

    • Search Google Scholar
    • Export Citation
  • 27. Wang J, Duncan D, Shi Z, et al. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res 2013;41:W77–W83.

    • Search Google Scholar
    • Export Citation
  • 28. Ashburner M, Ball CA, Blake JA, et al. Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000; 25: 25–29.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 29. Kelder T, Pico AR, Hanspers K, et al. Mining biological pathways using WikiPathways web services. PLoS One 2009;4:e6447.

  • 30. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000; 28: 27–30.

  • 31. Kanehisa M, Goto S, Sato Y, et al. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res 2014;42:D199–D205.

    • Search Google Scholar
    • Export Citation
  • 32. Taniguchi CM, Emanuelli B, Kahn CR. Critical nodes in signalling pathways: insights into insulin action. Nat Rev Mol Cell Biol 2006; 7: 85–96.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 33. Thauvin-Robinet C, Auclair M, Duplomb L, et al. PIK3R1 mutations cause syndromic insulin resistance with lipoatrophy. Am J Hum Genet 2013; 93: 141–149.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 34. Puigserver P, Spiegelman BM. Peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1 alpha): transcriptional coactivator and metabolic regulator. Endocr Rev 2003; 24: 78–90.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 35. Phillips CM, Goumidi L, Bertrais S, et al. Gene-nutrient interactions with dietary fat modulate the association between genetic variation of the ACSL1 gene and metabolic syndrome. J Lipid Res 2010; 51: 1793–1800.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 36. Yu W, Chu X, Chen G, et al. Studies of human mitochondrial 2,4-dienoyl-CoA reductase. Arch Biochem Biophys 2005; 434: 195–200.

  • 37. Heinzen EL, Swoboda KJ, Hitomi Y, et al. De novo mutations in ATP1A3 cause alternating hemiplegia of childhood. Nat Genet 2012; 44: 1030–1034.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 38. Belliard A, Sottejeau Y, Duan Q, et al. Modulation of cardiac Na+,K+-ATPase cell surface abundance by simulated ischemia-reperfusion and ouabain preconditioning. Am J Physiol Heart Circ Physiol 2013;304:H94–H103.

    • Search Google Scholar
    • Export Citation
  • 39. Dewey S, Lai X, Witzmann FA, et al. Proteomic analysis of hearts from Akita mice suggests that increases in soluble epoxide hydrolase and antioxidative programming are key changes in early stages of diabetic cardiomyopathy. J Proteome Res 2013; 12: 3920–3933.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 40. Motamed M, Rajapakshe KI, Hartig SM, et al. Steroid receptor coactivator 1 is an integrator of glucose and NAD+/NADH homeostasis. Mol Endocrinol 2014; 28: 395–405.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 41. Ashrafian H, Docherty L, Leo V, et al. A mutation in the mitochondrial fission gene Dnm1l leads to cardiomyopathy. PLoS Genet 2010;6:e1001000.

    • Search Google Scholar
    • Export Citation
  • 42. Zarrouk Mahjoub S, Mehri S, Ourda F, et al. Novel m.15434C>A (p.230L>I) mitochondrial Cytb gene missense mutation associated with dilated cardiomyopathy. ISRN Cardiol 2012; 2012: 1–6.

    • Search Google Scholar
    • Export Citation
  • 43. Sibbing D, Pfeufer A, Perisic T, et al. Mutations in the mitochondrial thioredoxin reductase gene TXNRD2 cause dilated cardiomyopathy. Eur Heart J 2011; 32: 1121–1133.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 44. Tidholm A, Häggström J, Borgarelli M, et al. Canine idiopathic dilated cardiomyopathy. Part I: aetiology, clinical characteristics, epidemiology and pathology. Vet J 2001; 162: 92–107.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 45. Tidholm A, Jonsson L. Histologic characterization of canine dilated cardiomyopathy. Vet Pathol 2005; 42: 1–8.

  • 46. Xu X, Zhang Y, Williams J, et al. Parallel comparison of Illumina RNA-seq and Affymetrix microarray platforms on transcriptomic profiles generated from 5-aza-deoxy-cytidine treated HT-29 colon cancer cells and simulated datasets. BMC Bioinformatics 2013; 14 (suppl 9): S1.

    • Crossref
    • Search Google Scholar
    • Export Citation

Advertisement