Abstract
Dilated cardiomyopathy (DCM) is characterized by decreased systolic function and dilation of one or both ventricles, often leading to heart failure or sudden death. Two 10-month-old sibling Nova Scotia Duck Tolling Retrievers (NSDTR) died acutely with evidence of dilated cardiomyopathy with myocardial fibrosis. Association analysis using two cases and 35 controls identified three candidate regions homozygous in the two cases. Whole genome sequencing identified a frameshift deletion in the LMNA gene (NC_049228.1:g.41688530del, NP_001274080:p.(Asp576ThrfsTer124)). Three retrospectively identified NSDTRs with sudden death before 2 years of age and severe myocardial fibrosis were also homozygous for the deletion. One 5 year old with sudden death and myocardial fibrosis was heterozygous for the deletion. This variant was not identified in 722 dogs of other breeds, nor was it identified to be homozygous in 784 NSDTR. LMNA codes for lamin A/C proteins, which are type V intermediate filaments that provide structural support to the nuclear membrane. In humans, LMNA variants can cause DCM with sudden death as well as diseases of striated muscles, lipodystrophy, neuropathies, and accelerated aging disorders. This frameshift deletion is predicted to affect processing of prelamin A into lamin A. Pedigree analysis in the NSDTR and functional evaluation of heterozygotes is consistent with a predominantly recessive mode of inheritance and possibly low penetrance in heterozygotes in contrast to people, where most pathogenic LMNA variants are dominantly inherited.
Similar content being viewed by others
Introduction
Cardiovascular disease is the fourth most common cause of death in dogs1 and one of the most common causes of death in humans2. Dilated cardiomyopathy (DCM) is the second most prevalent form of heart disease in dogs, accounting for 10% of cardiac diagnoses3. Characteristic findings of adult onset DCM in dogs include progressive dilation of the left or both ventricles with concurrent systolic dysfunction and eventual development of congestive heart failure or sudden death4.
While a number of phenocopies for DCM in dogs are reported, such as nutritionally mediated5,6, infectious disease-associated7,8, or metabolically mediated9, idiopathic canine DCM is considered a heritable disease10. Large breed dogs are more susceptible to adult onset DCM and several breeds, such as the Doberman Pinscher, Great Dane, Newfoundland, Boxer, Scottish Deerhound, and Irish Wolfhound, are over-represented. Mendelian early onset DCM has been identified in a number of breeds. A highly penetrant dominant missense variant in phospholamban was identified in a family of Welsh Springer Spaniels with young adult onset of DCM and sudden death11. Another highly penetrant recessive variant in the RBM20 gene was identified in the Standard and Giant Schnauzer breeds, leading to significantly reduced lifespan12. A juvenile form of DCM has also been documented in the Portuguese water dog and Toy Manchester Terrier, although a causative variant has not been identified13,14. One challenge in identifying inherited forms of DCM leading to sudden death is that necropsies are not routinely performed.
Many different forms of DCM have been identified in humans. Nineteen genes with high evidence for causation come from different ontologies, underscoring the genetic complexity of the disease process15. One of the genes with strong evidence for causation is the LMNA gene which codes for lamin A and lamin C, two alternatively spliced isoforms that are type V intermediate filament proteins localized to the nuclear lamina16. Autosomal dominant dilated cardiomyopathy and conduction system disease may result from variants in LMNA17. Patients with LMNA variants are at risk for sudden death and end-stage heart failure18.
Over 450 disease causing variants in the LMNA gene lead to metabolic, nervous system, and skeletal and cardiac muscle diseases (UMD-LMNA at http://www.umd.be)19,20. The genotype–phenotype correlation is unclear and the same variants may be found in patients with distinct disease phenotypes19,20. The nuclear lamina has important roles in nuclear shape and structure, as well as functional roles in transcriptional regulation and heterochromatin organization. Disruption of these processes contributes to the diverse pathogenesis of laminopathies21.
A 10-month-old Nova Scotia Duck Tolling Retriever (NSDTR) was identified with a ventricular tachyarrhythmia and DCM after presentation for evaluation of gastrointestinal disease. Within 2 weeks of presentation the dog died suddenly. Phenotypic evaluation of family members, GWAS and whole genome sequencing identified a frame shift variant in the LMNA gene segregating with DCM and/or sudden death and myocardial fibrosis.
Methods
Phenotype
The proband dog had an echocardiogram and thoracic radiographs. The dog died suddenly and had a complete necropsy. An affected littermate had an echocardiogram and thoracic radiographs. The sire, dam and grand dam of these dogs as well as two siblings each had an echocardiogram. Normal echocardiographic values were determined using previously published reference ranges22,23. An echocardiographic diagnosis of DCM was made based on the presence of both reduced left ventricle (LV) systolic function (%FS < 25) and systolic LV enlargement. Pedigree analysis was performed by using an online database for the NSDTR (https://www.k9data.com/). Life span for NSDTR was also determined when available from this online database. A second related litter with three siblings with sudden death at 15 months of age was born 11 years previously. Necropsy findings on one of these dogs were available for review. An additional affected dog (not closely related to the first two litters) was identified by screening the Bannasch Laboratory DNA Database for any cases of sudden death within this breed. A full necropsy report was available for this dog. Upon release of a genetic test for dog breeders, two additional affected dogs that had died 2 years prior were identified and had undergone necropsy. DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissues to determine the genotypes of these two dogs (Supplementary Table S1). A full sibling to the sire of the proband had sudden death at 5 years of age and a complete necropsy was performed. Comparison of lifespan between dogs heterozygous for the LMNA variant and wildtype dogs was performed using an unpaired t test as implemented in Graphpad Prism v0.97 (Dotmatics, Boston, MA, USA).
Samples used for genetic analyses
Three hundred and sixty-seven North American NSDTR DNA samples were used for this study. Thirty-five related NSDTR were used as controls for GWAS with two full sibling cases. One case, one obligate carrier and grandam of the proband were used for whole genome sequencing and variant calling. 300 NSDTR unrelated within 2 generations were genotyped for the LMNA variant identified. Additional related NSDTR were also genotyped to confirm the mode of inheritance (N = 32). An online database (https://www.k9data.com/) was used to identify related individuals. DNA samples were collected under Institutional Animal Care and Use Committees (IACUC) #s15356, 16892, 18561 and 22035. Four hundred and twenty-two DNA samples from NSDTR of mostly Swiss origin (351 Swiss, 57 other European countries, 14 unknown) from the Vetsuisse Biobank were used to estimate the allele frequency in Europe. None of those dogs had known reports of cardiac phenotypes. Blood samples were collected with the approval of the Cantonal Committee for Animal Experiments (Canton of Bern; permit BE 71/19). All sample collections were performed in accordance with the guidelines of the University of California Davis and the Canton of Bern respectively. The authors have complied with the ARRIVE guidelines in the design and execution of this project.
DNA extraction
Genomic DNA was extracted from buccal swabs and EDTA blood samples using the Gentra Puregene Kit (Qiagen, Valencia, CA, USA). FFPE kidney or liver scrolls were extracted using Zymo Research’s (Orange, CA, USA) Quick-DNA FFPE Miniprep kit. For European samples DNA was extracted from EDTA blood samples according to standard methods using the Maxwell RSX Whole Blood DNA kit in combination with the Maxwell RSC instrument (Promega, Dübendorf, Switzerland).
Genome-wide association study
The Illumina Canine HD BeadChip was used to generate genotype calls for 37 NSDTR, including two cases and 35 related NSDTR controls. PLINK v1.9 software was used to prune the data for minor allele frequency below 5% and individuals and single nucleotide variants (SNVs) with more than 10% missing genotypes24. A kinship matrix was implemented using GEMMA to correct for population stratification and a modified form of homozygosity mapping was performed. Only variants homozygous in the two cases were evaluated further 3636 variants were homozygous in the two cases. Bonferroni threshold was calculated for the 3636 variants at P = 0.05 and only SNV with significant associations were considered. A Manhattan plot was generated for SNV that were homozygous in the cases using the qqman package in RStudio25. UU_Cfam_GSD_1.0 dog reference genome assembly coordinates are used.
Whole genome sequencing and variant calling
One case, one obligate carrier and the grandam of the case were whole genome sequenced. PCR free Library preparation and 2 × 150 bp Illumina paired-end sequencing was performed at the UC Davis DNA Technologies Core. The raw read data was filtered using HTStream (version 1.3.3)26 which included screening for contaminants (such as PhiX), overlapping reads, quality-based trimming, and adapter trimming. BWA MEM (version 0.7.17)27 and Samtools (version 1.14)28 were used to align the processed data to the canine genome ((UU_Cfam_GSD_1.0). Sequencing duplicates were marked using Picard tools (version 2.26.11) (broadinstitute.github.io/picard/). Variant calling was accomplished using the GATK (version 4.2.5.0)29 pipeline which included the following steps: 1. HaplotypeCaller to call variants in the gVCF format, 2. GenomicsDB to create a datastore to store the variant call data, 3. GenotypeGVCFs to perform joint genotyping on all samples. Finally, SnpEff (version 5.1)30 was used to add effect prediction to the final VCFs using GSD_1.0 annotation31. Whole genome sequencing data from an affected NSDTR, one obligate carrier and one likely carrier, and 128 unaffected control dogs (Bioprojects PRJNA776905, PRJNA377155 and PRJNA961733) were analyzed for protein-coding variants unique to the affected and carrier dogs. WebGQT32 was utilized to filter based on a recessive pedigree model for variants.
Gene analysis
UU_Cfam_GSD_1.0 dog reference genome assembly was used to investigate the LMNA gene. IGV was used to visually confirm the variant33. NCBI RefSeq accession numbers NM_001287151.1 (mRNA) and NP_001274080.1 (protein) corresponding to the canine LMNA gene were used.
Database searches and functional predictions
The Ensembl Variant Effect Predictor release 107 was used to predict the biological outcome of the identified candidate variant. Plink v1.9 software was used to evaluate the LMNA genotype from a VCF file comprising 722 canine genomess34.
PCR and Sanger sequencing
PCR primers 5′-GAAGAGCCAGAGGAGATGGA and 5′-GTGCGAGCAGGAGTACAGAG were used to generate an amplicon containing the LMNA variant NM_001287151.1:c.1726del. PCR products from North American samples were directly Sanger sequenced on an Applied Biosystems 3500 Genetic Analyzer after enzymatic cleanup. Sanger sequences were analyzed using SnapGene® software (Dotmatics, Boston, MA, USA). European samples were sequenced on an ABI 3730 DNA Analyzer (Thermo Fisher Scientific, Reinach, Switzerland) and viewed in Sequencher 5.1 software (Gene Codes, Ann Arbor, MI, USA).
Quantification of collagen in cardiac muscle
FFPE blocks from dogs with no known cardiac conditions euthanized between 17 months and 3.5 years of age were used as references to examine fibrosis levels in diseased heart tissue. Masson’s trichrome stained slides of left ventricle from three control (one German shepherd, one Labrador retriever and one Rottweiler mix) and one homozygous case were prepared by the Anatomic Pathology Service at the UC Davis Veterinary Medical Teaching Hospital. Slides from a second and third homozygous case were prepared by Pathobiology at University of Guelph’s Ontario Veterinary College and the Veterinary Diagnostic Laboratory at the University of Arkansas respectively. Ten randomly selected fields of cross-sectional myocardium were imaged from each slide at 20 × magnification on a Revolve hybrid digital microscope. Fibrosis, quantified by collagen tissue proportion, was measured using Masson Trichrome vectors in Fiji’s Colour Deconvolution 2 plug-in35,36. Blue-stained collagen pixels were quantified in deconvolution Colour_1 according to Yen’s thresholding method and divided by total tissue pixels counted according to Intermodes thresholding37,38. Normality was tested using the Anderson–Darling test and an unpaired two-tailed t test was used to compare the collagen-to-total tissue proportion in LMNA mutant dogs to wildtype controls. Statistical tests and graphs were prepared using Prism Version 9.2.0
Verification of the mutant RNA
RNA was extracted from the white blood cells of a carrier and reverse transcribed to cDNA to verify the sequence of the mutant transcript. The C-terminal domain of LMNA was amplified with primers 5′-GAGGACGATGAGGATGAGGA (Exon 10) and 5′-AGCAGGGGAGATTGACATAGA (3′UTR). RT-PCR products were TOPO TA cloned and plasmids purified by column miniprep. Inserts were sequenced by the UCDNA Sequencing Facility at UC Davis.
Results
A 10-month old male NSDTR presented for diarrhea and vomiting. On physical exam an abnormal rhythm was auscultated and the dog was referred to a cardiologist. Paroxysmal ventricular tachycardia, severe DCM and mild suspected mitral valve dysplasia were diagnosed. Two weeks later the dog died suddenly. Concomitantly, a female littermate was reported to be experiencing inappetence and lethargy. The littermate was also diagnosed with DCM and sinus tachycardia with a left bundle branch block and died 2 weeks later. The parents, paternal granddam and two apparently unaffected siblings had echocardiograms and showed no evidence of DCM or heart disease at the time of screening (Table 1). Both affected puppies showed significant reduction in fractional shortening (< 25%) along with systolic and diastolic LV and LA dilation consistent with a diagnosis of advanced DCM22,23.
Pedigree analysis of the litter revealed that there was a recent inbreeding loop (Litter A, Fig. 1). A litter related to the proband had three puppies that died suddenly at ~ 15 months of age within a few weeks of one another (Litter B, Fig. 1). The proband and one puppy from the second litter had necropsies performed and both showed evidence of significant myocardial fibrosis. DNA samples were not available from the affected dogs from the second litter (Litter B, (Supplementary Table S1 and Fig. 1).
Pedigree of affected dogs. Affected dogs are outlined in red. Squares represent male dogs and circles represent female dogs. The asterisk indicates the proband case. Numbers below the symbols indicate lifespan of dogs when available. Genotype data for the LMNA variant is shown inside each symbol for the dogs with DNA available. A M indicates the mutant allele (NC_049228.1:g.41688530del) and a + indicates the normal allele. The three litters that produced affected puppies are indicated with letters A, B and C. The three litters have inbreeding loops. Age at death is noted under the pedigree symbols in months (m) or years (y).
Although DNA was initially only available from two cases, a genome-wide association study was performed using the two cases and 35 related, clinically unaffected NSDTR. Four candidate chromosomal regions homozygous in the cases were identified: chr7 (18693546–74786398), chr19 (33841298–34756758), chr24 (40956242–41731762), and a single SNP on chr3 (54049478)(Fig. 2).
Identification of candidate causative variants
Whole genome sequence was obtained from the affected sibling of the proband, the sire and paternal granddam. Candidate causal variants for DCM were identified by applying a hard filtering strategy on protein coding variants. The affected dog was required to be homozygous for the non-reference allele, the two presumed carriers were required to be heterozygous, and 128 previously sequenced dogs from other breeds were required to be homozygous for the reference allele.
Using this strategy two variants were identified. One was a 9 base pair deletion in the 5′-UTR of POLI (NC_049222.1:g.21302953-21302961del) on CFA1 (UU_Cfam_GSD_1.0). This variant did not fall within our candidate regions. The second variant was a single base deletion in the LMNA gene (NC_049228.1:g.41688530del) (UU_Cfam_GSD_1.0). LMNA is located on CFA 7 within the largest region of homozygosity identified in the association analysis. Based on its function and location, the LMNA variant was a very likely candidate causal variant. The LMNA variant causes a frameshift mutation, NM_001287151.1:c.1726del, NP_001274080:p.(Asp576ThrfsTer124).
In order to verify splicing of exons flanking this insertion, Sanger sequencing of RT-PCR products generated from RNA isolated from white blood cells from a heterozygous dog confirmed the intron–exon boundaries in the mutant and wildtype alleles for exons 10, 11 (which contains the deletion) and 12. The frameshift does not introduce a premature stop codon and is predicted to result in a mutant protein with 124 altered amino acids at its C-terminus (34 amino acids longer than the wildtype protein). Cleavage of lamin C is not predicted to be affected since residues 567 and 568 are unaltered. Due to the altered C-terminus, prelamin A processing is predicted to be disrupted. The first step of processing is farnesylation of the carboxy-terminal CaaX motif which is no longer present in the predicted mutant protein. This step is usually followed by cleavage at AA 646 which is also altered (Fig. 3). Expression from the mutant allele is predicted to lead to normal lamin C, aberrant C-terminal non-farnesylated prelamin A, and absent mature lamin A.
LMNA protein alignment of mutant dog, wildtype dog and human. The deletion results in a frameshift starting at codon 576 (red box). The predicted mutant protein contains 34 amino acids beyond the normal carboxy-terminus (red underline). The CaaX motif, CSIM, (blue box) is required for initial processing of prelamin to mature lamin A. The last 13 amino acids of prelamin are cleaved between a tyrosine and leucine residue (blue dotted box) to complete the processing of mature lamin A. Lamin C is cleaved between amino acids 566 and 567 (red dashed line) for processing and is predicted to occur normally.
Genotyping for the LMNA variant in the pedigree of the proband demonstrated that the variant segregated with disease and was consistent with a recessive mode of inheritance (Fig. 1). An additional dog with sudden death at 12 months of age and necropsy findings of severe myocardial fibrosis was identified in our collection of DNA samples from this breed. It was confirmed to be homozygous for the LMNA variant. Genotypes from relatives of this dog were consistent with a recessive mode of inheritance (Litter C, Fig. 1). Four dogs heterozygous for the LMNA variant were evaluated by echocardiography and did not show evidence of cardiac disease (Table 1 and Supplementary Table S1). One 8-year-old heterozygote also had a 24-h portable electrocardiogram monitor that showed no evidence of arrhythmias (Supplementary Table S1). Five dogs in the pedigree that were heterozygous or obligate carriers lived from 11 to 17 years. A 5-year old heterozygote sibling to the sire of Litter A had sudden death and a necropsy which showed myocardial fibrosis.
We genotyped 300 additional NSDTR from North America unrelated to 2 generations to the proband; none of the dogs were homozygous for the LMNA variant and the carrier frequency was 8.7%. The carrier frequency in an additional 422 NSDTRs from a European cohort was 0.2%. None of the dogs were homozygous for the variant. The LMNA variant also was not present in 722 whole genome sequenced dogs from various sources in the public domain34. Evaluation of the age of death of heterozygotes (N = 20) compared to wildtype (N = 44) dogs did not reveal a statistically significant difference in age (Het—mean age 12.12 years, WT—mean age 13.12 years).
After the NSDTR breeder community was informed about the presence of a genetic variant that could cause sudden death, two additional dogs were identified that died suddenly at a young age and had necropsies performed. Genotyping these dogs from FFPE samples demonstrated that both were homozygous for the LMNA variant (Supplementary Table S1). In order to quantify myocardial fibrosis, Masson’s trichrome stained left ventricle for these two cases and the proband were compared to three controls (Fig. 4a–d). Significant differences in blue staining consistent with cardiomyocyte fibrosis was observed between cases and controls (Fig. 4e). Cases were also observed to have a variable but subjective increase of adipocytes infiltrating the myocardial sections.
Severe myocardial Fibrosis. Sections are stained with Masson’s trichrome. Panel (a) is from a normal dog for comparison. Panels (b–d) are from three cases homozygous for the LMNA variant who died suddenly. Panel (e) showed the quantitative difference of the ratio of fibrosis/tissue between cases (M/M) and controls (+/+). The black calibration bar in panel (d) is 90 µm in length and the P value was < 0.0001.
Discussion
Young NSDTRs were identified that died suddenly and were diagnosed with DCM prior to death or severe myocardial fibrosis upon necropsy after sudden cardiac death. Genome-wide association analyses and whole genome sequencing identified a deletion resulting in a frameshift in the LMNA gene. Five young affected individuals were homozygous and one older affected dog was heterozygous for the deletion. Four heterozygotes were phenotypically normal and did not show abnormalities by echocardiography or ambulatory electrocardiogram when available. The frameshift variant identified in LMNA is predicted to result in normal production of lamin C. However, any expressed prelamin A would have a modified C-terminal tail, most likely interfering with the processing into mature lamin A due to disruption of key residues required for post-translational processing.
Nuclear lamins are an integral part of the inner nuclear membrane where they form an interlinked protein lattice structure that provides nuclear structural integrity and acts as a scaffold for chromatin organization39,40. Variants in LMNA result in more than 15 distinct tissue-specific diseases and over 450 disease causing variants have been identified in this gene in humans19,20. Dilated cardiomyopathy and conduction system disturbances are a common sequela of LMNA variants and can occur in combination with other phenotypes or in isolation17,41,42,43,44,45,46. In humans, 6–33% of cases presenting with both DCM and cardiac conduction defects are due to variants in LMNA41,47. LMNA-related DCM are inherited in an autosomal dominant manner either due to missense variants or truncating frameshift variants which are predicted to alter both lamin A and lamin C proteins41,47.
The LMNA variant identified in NSDTR is predicted to produce exclusively non-farnesylated prelamin A and normal lamin C. Mice that produce only non-farnesylated prelamin A develop DCM with mild to moderate fibrosis and die prematurely; however, in this mouse model no lamin C is produced48. Mice without lamin A or prelamin A have also been shown to be phenotypically normal49. The accumulation of non-farnesylated prelamin A is a gain of function and would be predicted to result in a dominant phenotype, however in the NSDTR, there is evidence for incomplete dominance with low penetrance in heterozygotes. The severe phenotype in homozygous mutant NDSTRs suggests loss of function of mature lamin A as a major contributor to the disease phenotype.
Similar to the mouse mutants which accumulate non-farnesylated prelamin A, affected NSDTR cardiac tissue showed significant fibrosis. Pathological myocardial fibrosis occurs at the first signs of conduction disturbances and prior to echocardiographic changes in humans with dominant variants in LMNA50. Cardiovascular magnetic resonance (CMR), a technique not used in veterinary medicine, has demonstrated 88% of asymptomatic or mildly symptomatic carriers of LMNA variants have myocardial fibrosis51,52. Recently, activation of the DNA damage response (DDR) pathway has been implicated in the induction of inflammation in laminopathies53,54. In fact, blocking the DDR pathway by disrupting the Mb21d1 gene lead to increased survival and decreased apoptosis and fibrosis in Lmna-deficient mice55.
A human variant very similar to the one identified in the NSDTR is LMNA p.Thr655AsnTer49, which leads to an altered and extended C-terminus, although the frameshift occurs 80 codons later than in the NSDTR variant. This variant causes accumulation of non-farnesylated prelamin A56. People heterozygous and homozygous for this variant have lipodystrophy. Homozygotes can also have a cardiac phenotype characterized by DCM and/or rhythm/conduction disturbances and sometimes sudden death57,58. Interestingly there is no evidence for lipodystrophy in the NSDTR. One significant disease manifestation of lipodystrophy in humans is an increased risk of type 2 diabetes. Unlike humans, insulin resistant diabetes is not recognized clinically in dogs since their beta cells do not appear to be sensitive to hyperglycemia59. Even when obesity is induced in dogs, type 2 diabetes only develops with the addition of streptozotocin60.
One distinct difference between human cardiac laminopathies and the NSDTR and mouse mutants is the mode of inheritance. Virtually all human variants that lead to DCM and cardiac fibrosis are dominantly inherited with the exception of LMNA:p.Thr655AsnTer49 where some homozygotes showed cardiac defects61. Mouse Lmna mutants (both null and missense) also show a recessive mode of inheritance48,62,63. Four dogs that were carriers of the LMNA p.Asp576ThrTer124 variant received echocardiograms which showed no evidence of DCM (Table 1, Supplementary Table S1). In addition, other heterozygotes in this pedigree lived normal lifespans (N = 4; 11, 12, 15 and 17 years). However, the identification of a heterozygote with cardiac fibrosis and sudden death at 5 years of age suggests that the mode of inheritance might not be a simple recessive. On the other hand, It is possible that there are other risk factors that contributed to this dog’s death.
The clinical course of the NSDTR homozygous for the LMNA variant was similar in all six affected dogs. They were all young (10–15 months) and died suddenly. The proband was diagnosed with heart disease prior to sudden death but his history and clinical signs were gastrointestinal in nature. His littermate was subsequently diagnosed with DCM. The other affected homozygotes and the single heterozygote had no history of illness prior to sudden death. There is no evidence of incomplete penetrance when having two copies of the LMNA variant, however only a small number of homozygous dogs have been identified. In the two litters with multiple affected dogs, the age of onset was slightly different (Litter A 10 months and Litter B 15 months) opening the possibility of environmental or genetic modification of the phenotype. Importantly five cases were identified because their owners opted to have a complete necropsy performed. The importance of a necropsy for any young animal that dies suddenly cannot be overemphasized.
Familial DCM has not been described in the NSDTR dog breed prior to this study. Genotyping of NSDTR dogs from North America revealed a carrier frequency of 8.7% whereas European origin NSDTR dogs had a carrier frequency of 0.2%. By increasing awareness of this genetic defect within the NSDTR community we hope breeders are able to identify carriers and avoid pairings that would produce affected puppies.
In purebred dogs there are long regions of linkage disequilibrium which can challenge the identification of the causative variant underlying a phenotype64. The evidence that NM_001287151.1:c.1726del is a pathogenic variant is very strong65. The variant segregates with disease in distantly related families of NSDTR. The variant is not found in other breeds and was not found to be homozygous in 784 unaffected NSDTR. There is also a large body of evidence from human and mouse studies that LMNA variants can cause myocardial fibrosis and conduction system disturbances similar to what was seen in the NSDTR. Lastly the nature of the change to the LMNA isoforms is very similar but more severe than one seen in humans which can cause myocardial disease. The identification of LMNA p.(Asp576ThrfsTer124) as a pathogenic variant associated with sudden death and severe myocardial fibrosis provides a large animal model of a cardiac laminopathy. The underlying severe myocardial fibrosis seen in these young dogs also underscores the importance of inflammation and fibrosis to conduction disturbances of LMNA variants.
Data availability
The whole genome sequence data generated during the current study are available in the SRA repository, https://www.ncbi.nlm.nih.gov/sra, Bioproject PRJNA961733.
References
Fleming, J. M., Creevy, K. E. & Promislow, D. E. Mortality in north American dogs from 1984 to 2004: An investigation into age-, size-, and breed-related causes of death. J. Vet. Intern. Med. 25, 187–198. https://doi.org/10.1111/j.1939-1676.2011.0695.x (2011).
Mathers, C. D. & Loncar, D. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med. 3, e442. https://doi.org/10.1371/journal.pmed.0030442 (2006).
Egenvall, A., Bonnett, B. N. & Häggström, J. Heart disease as a cause of death in insured Swedish dogs younger than 10 years of age. J. Vet. Intern. Med. 20, 894–903. https://doi.org/10.1892/0891-6640(2006)20[894:hdaaco]2.0.co;2 (2006).
Bonagura, J. D. & Visser, L. C. Echocardiographic assessment of dilated cardiomyopathy in dogs. J. Vet. Cardiol. 40, 15–50. https://doi.org/10.1016/j.jvc.2021.08.004 (2022).
Freid, K. J. et al. Retrospective study of dilated cardiomyopathy in dogs. J. Vet. Intern. Med. 35, 58–67. https://doi.org/10.1111/jvim.15972 (2021).
Ontiveros, E. S. et al. Development of plasma and whole blood taurine reference ranges and identification of dietary features associated with taurine deficiency and dilated cardiomyopathy in golden retrievers: A prospective, observational study. PLoS One 15, e0233206. https://doi.org/10.1371/journal.pone.0233206 (2020).
Ford, J., McEndaffer, L., Renshaw, R., Molesan, A. & Kelly, K. Parvovirus infection is associated with myocarditis and myocardial brosis in young dogs. Vet. Pathol. 54, 964–971 (2017).
Santilli, R. A. et al. Bartonella-associated inflammatory cardiomyopathy in a dog. J. Vet. Cardiol. 19, 74–81 (2017).
Phillips, D. E. & Harkin, K. R. Hypothyroidism and myocardial failure in two Great Danes. J. Am. Anim. Hosp. Assoc. 39, 133–137. https://doi.org/10.5326/0390133 (2003).
Gaar-Humphreys, K. R. et al. Genetic basis of dilated cardiomyopathy in dogs and its potential as a bidirectional model. Animals (Basel) https://doi.org/10.3390/ani12131679 (2022).
Yost, O., Friedenberg, S. G., Jesty, S. A., Olby, N. J. & Meurs, K. M. The R9H phospholamban mutation is associated with highly penetrant dilated cardiomyopathy and sudden death in a spontaneous canine model. Gene 697, 118–122. https://doi.org/10.1016/j.gene.2019.02.022 (2019).
Leach, S. B., Briggs, M., Hansen, L. & Johnson, G. S. Prevalence, geographic distribution, and impact on lifespan of a dilated cardiomyopathy-associated RNA-binding motif protein 20 variant in genotyped dogs. J. Vet. Cardiol. 40, 119–125. https://doi.org/10.1016/j.jvc.2021.05.002 (2022).
Werner, P. et al. A novel locus for dilated cardiomyopathy maps to canine chromosome 8. Genomics 91, 517–521. https://doi.org/10.1016/j.ygeno.2008.03.007 (2008).
Legge, C. H. et al. Histological characterization of dilated cardiomyopathy in the juvenile toy Manchester terrier. Vet. Pathol. 50, 1043–1052. https://doi.org/10.1177/0300985813480509 (2013).
Jordan, E. et al. Evidence-based assessment of genes in dilated cardiomyopathy. Circulation 144, 7–19. https://doi.org/10.1161/circulationaha.120.053033 (2021).
Fisher, D. Z., Chaudhary, N. & Blobel, G. cDNA sequencing of nuclear lamins A and C reveals primary and secondary structural homology to intermediate filament proteins. Proc. Natl. Acad. Sci. USA 83, 6450–6454. https://doi.org/10.1073/pnas.83.17.6450 (1986).
Fatkin, D. et al. Missense mutations in the rod domain of the lamin A/C gene as causes of dilated cardiomyopathy and conduction-system disease. N. Engl. J. Med. 341, 1715–1724. https://doi.org/10.1056/nejm199912023412302 (1999).
Taylor, M. R. et al. Natural history of dilated cardiomyopathy due to lamin A/C gene mutations. J. Am. Coll. Cardiol. 41, 771–780. https://doi.org/10.1016/s0735-1097(02)02954-6 (2003).
Crasto, S., My, I. & Di Pasquale, E. The broad spectrum of LMNA cardiac diseases: From molecular mechanisms to clinical phenotype. Front. Physiol. 11, 761. https://doi.org/10.3389/fphys.2020.00761 (2020).
Storey, E. C. & Fuller, H. R. Genotype–phenotype correlations in human diseases caused by mutations of LINC complex-associated genes: A systematic review and meta-summary. Cells https://doi.org/10.3390/cells11244065 (2022).
Capell, B. C. & Collins, F. S. Human laminopathies: Nuclei gone genetically awry. Nat. Rev. Genet. 7, 940–952. https://doi.org/10.1038/nrg1906 (2006).
Cornell, C. C. et al. Allometric scaling of M-mode cardiac measurements in normal adult dogs. J. Vet. Intern. Med. 18, 311–321. https://doi.org/10.1892/0891-6640(2004)18%3c311:asomcm%3e2.0.co;2 (2004).
Rishniw, M. & Erb, H. N. Evaluation of four 2-dimensional echocardiographic methods of assessing left atrial size in dogs. J. Vet. Intern. Med. 14, 429–435. https://doi.org/10.1892/0891-6640(2000)014%3c0429:eofemo%3e2.3.co;2 (2000).
Purcell, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. https://doi.org/10.1086/519795 (2007).
qqman: An R package for visualizing GWAS results using Q-Q and Manhattan plots. J. Open Source Soft. (2018).
Petersen, K. R., Streett, D. A., Gerritsen, A. T., Hunter, S. S. & Settles, M. L. Super deduper, fast PCR duplicate detection in fastq files. In Proceedings of the 6th ACM Conference on Bioinformatics, Computational Biology and Health Informatics (2015).
Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 26, 589–595. https://doi.org/10.1093/bioinformatics/btp698 (2010).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079. https://doi.org/10.1093/bioinformatics/btp352 (2009).
McKenna, A. et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. https://doi.org/10.1101/gr.107524.110 (2010).
Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6, 80–92. https://doi.org/10.4161/fly.19695 (2012).
Wang, C. et al. A novel canine reference genome resolves genomic architecture and uncovers transcript complexity. Commun. Biol. 4, 185. https://doi.org/10.1038/s42003-021-01698-x (2021).
Arumilli, M., Layer, R. M., Hytönen, M. K. & Lohi, H. webGQT: A shiny server for genotype query tools for model-based variant filtering. Front. Genet. 11, 152. https://doi.org/10.3389/fgene.2020.00152 (2020).
Robinson, J. T. et al. Integrative genomics viewer. Nat. Biotechnol. 29, 24–26. https://doi.org/10.1038/nbt.1754 (2011).
Plassais, J. et al. Whole genome sequencing of canids reveals genomic regions under selection and variants influencing morphology. Nat. Commun. 10, 1489. https://doi.org/10.1038/s41467-019-09373-w (2019).
Landini, G., Martinelli, G. & Piccinini, F. Colour deconvolution—Stain unmixing in histological imaging. Bioinformatics 37, 1485–1487. https://doi.org/10.1093/bioinformatics/btaa847 (2020).
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682. https://doi.org/10.1038/nmeth.2019 (2012).
Prewitt, J. M. S. & Mendelsohn, M. L. The analysis of cell images. Ann. N. Y. Acad. Sci. 128, 1035–1053. https://doi.org/10.1111/j.1749-6632.1965.tb11715.x (2006).
Yen, J.-C., Chang, F.-J. & Chang, S. A new criterion for automatic multilevel thresholding. IEEE Trans. Image Process. 4, 370–378. https://doi.org/10.1109/83.366472 (1995).
Patil, S. & Sengupta, K. Role of A- and B-type lamins in nuclear structure-function relationships. Biol. Cell 113, 295–310. https://doi.org/10.1111/boc.202000160 (2021).
Shevelyov, Y. Y. & Ulianov, S. V. The nuclear lamina as an organizer of chromosome architecture. Cells https://doi.org/10.3390/cells8020136 (2019).
van Tintelen, J. P. et al. High yield of LMNA mutations in patients with dilated cardiomyopathy and/or conduction disease referred to cardiogenetics outpatient clinics. Am. Heart J. 154, 1130–1139. https://doi.org/10.1016/j.ahj.2007.07.038 (2007).
Kane, M. S. et al. LMNA-associated cardiocutaneous progeria: An inherited autosomal dominant premature aging syndrome with late onset. Am. J. Med. Genet. A 161a, 1599–1611. https://doi.org/10.1002/ajmg.a.35971 (2013).
Renou, L. et al. Heart-hand syndrome of Slovenian type: A new kind of laminopathy. J. Med. Genet. 45, 666–671. https://doi.org/10.1136/jmg.2008.060020 (2008).
Chen, L. et al. LMNA mutations in atypical Werner’s syndrome. Lancet 362, 440–445. https://doi.org/10.1016/s0140-6736(03)14069-x (2003).
Nguyen, D. et al. Collagen expression in fibroblasts with a novel LMNA mutation. Biochem. Biophys. Res. Commun. 352, 603–608. https://doi.org/10.1016/j.bbrc.2006.11.070 (2007).
van der Kooi, A. J. et al. Lamin A/C mutations with lipodystrophy, cardiac abnormalities, and muscular dystrophy. Neurology 59, 620–623. https://doi.org/10.1212/wnl.59.4.620 (2002).
Parks, S. B. et al. Lamin A/C mutation analysis in a cohort of 324 unrelated patients with idiopathic or familial dilated cardiomyopathy. Am. Heart J. 156, 161–169. https://doi.org/10.1016/j.ahj.2008.01.026 (2008).
Davies, B. S. et al. An accumulation of non-farnesylated prelamin A causes cardiomyopathy but not progeria. Hum. Mol. Genet. 19, 2682–2694. https://doi.org/10.1093/hmg/ddq158 (2010).
Fong, L. G. et al. Prelamin A and lamin A appear to be dispensable in the nuclear lamina. J. Clin. Investig. 116, 743–752. https://doi.org/10.1172/jci27125 (2006).
van Tintelen, J. P. et al. Severe myocardial fibrosis caused by a deletion of the 5′ end of the lamin A/C gene. J. Am. Coll. Cardiol. 49, 2430–2439. https://doi.org/10.1016/j.jacc.2007.02.063 (2007).
Holmström, M. et al. Late gadolinium enhanced cardiovascular magnetic resonance of lamin A/C gene mutation related dilated cardiomyopathy. J. Cardiovasc. Magn. Reson. 13, 30. https://doi.org/10.1186/1532-429x-13-30 (2011).
Fontana, M. et al. CMR-verified interstitial myocardial fibrosis as a marker of subclinical cardiac involvement in LMNA mutation carriers. JACC Cardiovasc. Imaging 6, 124–126. https://doi.org/10.1016/j.jcmg.2012.06.013 (2013).
Cho, S. et al. Mechanosensing by the lamina protects against nuclear rupture, DNA damage, and cell-cycle arrest. Dev. Cell 49, 920-935.e925. https://doi.org/10.1016/j.devcel.2019.04.020 (2019).
Chen, S. N. et al. DNA damage response/TP53 pathway is activated and contributes to the pathogenesis of dilated cardiomyopathy associated with LMNA (lamin A/C) mutations. Circ. Res. 124, 856–873. https://doi.org/10.1161/CIRCRESAHA.118.314238 (2019).
Cheedipudi, S. M., Asghar, S. & Marian, A. J. Genetic ablation of the DNA damage response pathway attenuates lamin-associated dilated cardiomyopathy in mice. JACC Basic Transl. Sci. 7, 1232–1245. https://doi.org/10.1016/j.jacbts.2022.06.015 (2022).
Le Dour, C. et al. A homozygous mutation of prelamin-A preventing its farnesylation and maturation leads to a severe lipodystrophic phenotype: New insights into the pathogenicity of nonfarnesylated prelamin-A. J. Clin. Endocrinol. Metab. 96, E856–E862. https://doi.org/10.1210/jc.2010-2234 (2011).
Andre, P. et al. Metabolic and cardiac phenotype characterization in 37 atypical Dunnigan patients with nonfarnesylated mutated prelamin A. Am. Heart J. 169, 587–593. https://doi.org/10.1016/j.ahj.2014.12.021 (2015).
Treiber, G. et al. A recurrent familial partial lipodystrophy due to a monoallelic or biallelic LMNA founder variant highlights the multifaceted cardiac manifestations of metabolic laminopathies. Eur. J. Endocrinol. 185, 453–462. https://doi.org/10.1530/eje-21-0282 (2021).
Verkest, K. R., Rand, J. S., Fleeman, L. M. & Morton, J. M. Spontaneously obese dogs exhibit greater postprandial glucose, triglyceride, and insulin concentrations than lean dogs. Domest. Anim. Endocrinol. 42, 103–112. https://doi.org/10.1016/j.domaniend.2011.10.002 (2012).
Ionut, V. et al. Novel canine models of obese prediabetes and mild type 2 diabetes. Am. J. Physiol. Endocrinol. Metab. 298, E38–E48. https://doi.org/10.1152/ajpendo.00466.2009 (2010).
Amberger, J. S., Bocchini, C. A., Schiettecatte, F., Scott, A. F. & Hamosh, A. OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 43, D789–D798. https://doi.org/10.1093/nar/gku1205 (2015).
Mounkes, L. C., Kozlov, S. V., Rottman, J. N. & Stewart, C. L. Expression of an LMNA-N195K variant of A-type lamins results in cardiac conduction defects and death in mice. Hum. Mol. Genet. 14, 2167–2180. https://doi.org/10.1093/hmg/ddi221 (2005).
Nikolova, V. et al. Defects in nuclear structure and function promote dilated cardiomyopathy in lamin A/C-deficient mice. J. Clin. Investig. 113, 357–369. https://doi.org/10.1172/jci19448 (2004).
Lindblad-Toh, K. et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 438, 803–819. https://doi.org/10.1038/nature04338 (2005).
Richards, S. et al. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 17, 405–424. https://doi.org/10.1038/gim.2015.30 (2015).
Acknowledgements
We thank the Drs. Madi Garcia, Rebecca Egan, Don Stanton and Van Wettere for their assistance and the owners of the dogs whose samples were used for this study. This work was supported by the Center for Companion Animal Health CCAH 2021-88-F and the Maxine Adler Endowed Chair funds. D.T.O. was supported by the Students Training in Advanced Research (STAR) Program NIH 5T35OD010956-23. Support for J.L.K. was provided in part by NIH T32 OD011147 and Morris Animal Foundation D22FE-027.
Author information
Authors and Affiliations
Contributions
D.B. conceived the study. D.B., D.O., K.B., J.V., J.S. and R.L. contributed to the study design. Data collection, analysis and figures were performed by D.B., D.O. and J.V. Data collection was performed by J.K., J.S., I.M., M.C., T.L., and N.J. The first draft of the manuscript was written by D.O. and D.B. and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bannasch, D.L., Oertle, D.T., Vo, J. et al. Naturally occurring canine laminopathy leading to a dilated and fibrosing cardiomyopathy in the Nova Scotia Duck Tolling Retriever. Sci Rep 13, 19077 (2023). https://doi.org/10.1038/s41598-023-46601-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-46601-2