Next Article in Journal
Intrafollicular and Systemic Dopamine, Noradrenaline and Adrenaline Concentrations in Cycling Mares
Next Article in Special Issue
Chylopericardium Effusion in a Lac Alaotra Bamboo Lemur (Hapalemur alaotrensis)
Previous Article in Journal
Who Is Pulling the Leash? Effects of Human Gender and Dog Sex on Human–Dog Dyads When Walking On-Leash
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Genomic Study of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels

by
Arianna Bionda
1,2,†,
Matteo Cortellari
1,2,†,
Mara Bagardi
1,
Stefano Frattini
1,
Alessio Negro
1,2,
Chiara Locatelli
1,*,
Paola Giuseppina Brambilla
1 and
Paola Crepaldi
1,2
1
Department of Veterinary Medicine, University of Milan, Via dell’Università 6, 26900 Lodi, Italy
2
Department of Agricultural and Environmental Sciences, University of Milan, Via Celoria 2, 20133 Milan, Italy
*
Author to whom correspondence should be addressed.
Co-first authors.
Animals 2020, 10(10), 1895; https://doi.org/10.3390/ani10101895
Submission received: 2 September 2020 / Revised: 8 October 2020 / Accepted: 12 October 2020 / Published: 16 October 2020

Abstract

:

Simple Summary

Myxomatous mitral valve disease is the most common acquired cardiopathy in dogs. The earliest onset and the highest incidence of the disease is found in cavalier King Charles spaniels. Previous studies have suggested a polygenic inheritance of the disease in this breed. Here, we aim to expand the knowledge of the genetic basis of early-onset myxomatous mitral valve disease in cavalier King Charles spaniels. The selection of cases and controls is a crucial point of the study and is based on clinical, echocardiographic, and genealogical examinations. We perform three complementary genomic analyses that have never been used for investigating this pathology to identify 10 genes differentiating the genomes of the two groups. After examining these genes, we conclude that HEPACAM2, CDK6, and FAH (which are also related to transforming growth factor β (TGF-β) pathway, probably associated with the development of the disease) are the ones that are most likely involved in the pathogenesis of myxomatous mitral valve disease.

Abstract

Cavalier King Charles spaniels (CKCSs) show the earliest onset and the highest incidence of myxomatous mitral valve disease (MMVD). Previous studies have suggested a polygenic inheritance of the disease in this breed and revealed an association with regions on canine chromosomes 13 and 14. Following clinical and echocardiographic examinations, 33 not-directly-related CKCSs were selected and classified as cases (n = 16) if MMVD was present before 5 years of age or as controls (n = 17) if no or very mild MMVD was present after 5 years of age. DNA was extracted from whole blood and genotyped with a Canine 230K SNP BeadChip instrument. Cases and controls were compared with three complementary genomic analyses (Wright’s fixation index—FST, cross-population extended haplotype homozygosity—XP-EHH, and runs of homozygosity—ROH) to identify differences in terms of heterozygosity and regions of homozygosity. The top 1% single-nucleotide polymorphisms (SNPs) were selected and mapped, and the genes were thoroughly investigated. Ten consensus genes were found localized on chromosomes 3-11-14-19, partially confirming previous studies. The HEPACAM2, CDK6, and FAH genes, related to the transforming growth factor β (TGF-β) pathway and heart development, also emerged in the ROH analysis. In conclusion, this work expands the knowledge of the genetic basis of MMVD by identifying genes involved in the early onset of MMVD in CKCSs.

1. Introduction

Myxomatous mitral valve disease (MMVD) is the most common acquired heart disease in dogs, accounting for approximately 75% of all dogs with heart disease [1], particularly in older dogs and smaller dog breeds [2]. The cavalier King Charles spaniel (CKCS) shows the earliest onset and the highest incidence of MMVD when compared with other breeds [1,3,4,5]. The primary clinical finding in dogs affected by MMVD is a systolic heart murmur, but no murmur might occur in mild cases [6]; thus, an echocardiography is considered as the gold standard for the confirmation and staging of this disease [7].
There is evidence from the literature that, at least in some breeds, such as CKCSs and dachshunds, the hereditary component plays a predominant role in the pathogenesis of MMVD [8,9,10]. This disease, in particular, has been suggested to be inherited as a polygenic trait. In fact, the proportion of offspring with heart murmurs and the intensity of these murmurs are both significantly greater with increased parental severity [8,10]. Moreover, early-onset MMVD, typically found in CKCSs, also appears to be highly heritable. In particular, the heritability is 0.67 ± 0.07 for the degree of the heart murmur and 0.33 ± 0.07 for the presence or absence of the murmur, considering dogs exclusively aged between 4 and 5 years [9].
Previous studies have demonstrated that the linkage disequilibrium in dogs is 10 to 100 times more extensive than in the human genome; therefore, the number of single-nucleotide polymorphism (SNP) markers required for genomic association studies in dogs is considerably lower than in humans [11,12,13,14]. Moreover, the peculiarities of breed structure of dogs (e.g., a small breeding population and the use of popular sires) reduce the within-breed genetic variability, generally making even a small number of cases and controls (from 20 to 100, depending on the type of trait [15]) useful to effectively detect genomic regions that are associated with a particular trait or disease when genotyping just 15,000 SNPs [11,14], which is 10 times less than those contained in the SNP chip used in the present study.
In the literature, it is possible to find several genomic studies on canine MMVD, as summarized in Table 1. Madsen et al. (2011) identified two loci on canine chromosomes (CFA) 13 and 14 that are weakly associated with the development of MMVD via a genome-wide association study (GWAS) of CKCSs [3]. French at al. (2012), in contrast, did not find any evidence for loci associated with a mitral valve murmur with a GWAS in this breed, nor regions of highly discrepant homo/heterozygosity. The authors concluded that the familial occurrence of mitral valve murmurs in the CKCS breed is not due to a single major gene effect [16]. Meurs et al. (2017) performed a whole-genome sequencing of 10 CKCSs and 10 dachshunds. They filtered the variants of canine gene orthologs of the human genes known to be associated with MMVD against a database of variants derived from whole-genome sequencing of 98 medium to large dog breeds. The latter were chosen, assuming that the prevalence of MMVD in medium to large dogs is very low, but their phenotypes were not assessed. No variant was found in any of the genes evaluated that were present in at least eight of 10 affected samples, but a single coding variant, predicted to be benign, was found in the COL5A1 gene in nine of the 10 affected CKCSs and in, at most, 5% of the other breeds examined [17].
The Maltese breed was the object of two genomic studies. Lee et al. (2018) applied a candidate gene polymorphism approach and identified six polymorphisms of the SERT gene in samples with MMVD [18]. The GWAS performed by Lee et al. (2019), instead, revealed significant SNPs in several genes associated with cardiac function, including PDZD2, ARVCF, CTNNA3, and LDLRAD4 [19].
Torres-García et al. (2016) studied the polymorphisms of the COL1A2 gene, which interestingly localizes on the region of CFA 14 (identified by Madsen et al.) in poodles, and found an association between the rs22372411 variant and susceptibility to MMVD [20].
The problem of all GWAS case-control studies concerning MMVD is the late onset of the disease, making the identification of a real control which will not develop the disease at a later age more difficult. In 2015, Stern et al. tried to solve this problem by associating their results with the whippet dog breed with a continuous variable that included the age of onset of the MMVD (considered early if under 5 years of age) and a score of severity. In this study, a genome-wide significant association was identified on the region of CFA 15, containing FSTL5, and of CFA 2, containing ARHGAP26 [21].
At present, there is no genetic test for MMVD, and the breeding plans for CKCSs put in place so far have not been as effective as expected, especially if based only on auscultatory findings [22,23,24]. Nevertheless, an accurate selection of animals for breeding is essential, since the high prevalence of this pathology in this breed makes the elimination of all the dogs diagnosed with MMVD from reproduction unfeasible. Bagardi et al. [4] clinically and echocardiographically evaluated a representative sample of the Italian population of CKCSs, identifying signs of myxomatous degeneration in 90% of them, also widely present in young dogs. Discovering the genetic basis of MMVD may increase the effectiveness of breeding protocols, allowing an early identification of subjects predisposed to a severe form of this pathology. For this reason, here, we conduct a genomic study on the Italian population analyzed by Bagardi et al. [4]. The selection of cases and controls was a crucial aspect of our study and was based on a complete echocardiographic and genealogic examination, with the latter also being checked against the genomic data. The genomic analyses we performed here have never been used for the investigation of this pathology, and the examination of previously published data about the genes and pathway we identified here also allowed us to describe their possible role in the pathogenesis of the disease. Therefore, the aim of this study was to find genomic regions associated with early-onset MMVD predisposition in the CKCS dog breed.

2. Materials and Methods

This study was carried out with 33 privately owned CKCSs. The main selection criteria were the American College of Veterinary Internal Medicine (ACVIM) classification and age. Our samples were grouped as follows:
  • Cases (n = 16): Dogs with MMVD (class B1 or more severe) diagnosed before the age of 5 or with severe disease (class C or D) before the age of 8;
  • Controls (n = 17): Dogs without MMVD (class A) or with extremely mild signs of MMVD (class B1 with a trivial mitral regurgitation characterized by a maximal ratio of the regurgitant jet area signal to left atrium area ≤ 20%) [25] over 5 years of age or those suffering from a mild form of disease (class B1) over 8 years of age.
By way of example, subjects in ACVIM class B2 that were older than 5 years or in more severe classes that were older than 8 years could not be considered as either cases or controls; thus, they were excluded from this study.
Moreover, in order to limit the rate of consanguinity among the dogs as far as possible, in cases of close kinship, the sample with more “extreme” characteristics (the youngest affected or the oldest healthy dog) was chosen.
These dogs were selected among a larger group of 90 CKCSs, examined at the Cardiology Unit of the Department of Veterinary Medicine of University of Milan between December 2018 and September 2019. Informed consent was signed by the owners, in compliance with the ethical committee statement of the University of Milan number 2/2016.
Information regarding birthdates was verified by checking each animal’s microchip number in the regional registry, while the genealogical study was derived from the consultation of the online genealogy book of the Italian Kennel Club (ENCI) (http://www.enci.it/libro-genealogico/libro-genealogico-on-line#) or the Centrale Canine site (https://www.centra-canine.fr/lofselect).
The cardiovascular system was evaluated by checking the presence/absence of a murmur and, if present, its intensity (grade I–VI/VI) and point of maximum intensity.
Since auscultation in dogs with echocardiographic evidence of this disease has often proven to be normal, in addition to clinical data, all subjects underwent a complete echocardiographic examination that was performed by three well-trained investigators using a MyLab50 Gold cardiovascular ultrasound machine (Esaote, Florence, Italy). The exams were carried out according to a standard procedure with concurrent continuous electrocardiographic monitoring [26]. Dogs were staged according to the ACVIM guidelines [7].
Peripheral venous blood sampling was performed at the end of the examination. Blood was collected from the jugular vein into 2.5 mL EDTA tubes after a 12-h fasting period.

2.1. Statistical Analysis

The statistical analysis was performed using JMP® 15.0.0 (SAS Institute Inc., Cary, NC, USA, 1989–2019). Appropriate descriptive analyses were applied. Variables were reported as the mean ± standard deviation if they were normally distributed after an Anderson–Darling test; otherwise, they were reported as median and interquartile range values. Individual and clinical data of the subjects included in different groups were compared with linear regression and χ2 tests. The statistical methods applied to the genealogical and genomic analyses are reported below.

2.2. Genealogic Analysis

Pedigree information was analyzed using Optisel [27]. The inbreeding coefficient (F) and average relatedness coefficient (AR) were both calculated.
F is the probability of an individual receiving, at one locus, two identical-by-descent alleles that are copies of a single allele carried by a common ancestor of the parents. AR is the probability that an allele randomly chosen from the whole population in the pedigree belongs to a given animal.

2.3. DNA Extraction and Genomic Analysis

The whole-blood samples, collected in EDTA tubes, were stored at −20 °C. DNA was extracted using a DNeasy Blood and Tissue Kit (QIAGEN®, Hilden, Germany) according to the manufacturer’s instructions. The concentration and quality of the DNA of each sample were both assessed using a NanoDrop 1000 spectrophotometer (Thermo Scientific®, Waltham, MA, USA). The quality and quantity of the DNA extracted were suitable for the downstream analyses.
All 33 samples were genotyped via outsourcing by Agrotis S.r.l., Laboratory of Genetics and Services, using Canine 230K SNP BeadChips (containing over 230,000 SNPs) on an iScan System (Illumina®, San Diego, CA, USA).
Raw genotype data (in the .ped and .map formats) were processed for quality control using the PLINK 1.9 software package [28] as follows: SNPs were excluded if they had a call rate of <95% or if they had a minor allele frequency (MAF) of <1%; SNPs residing on sex chromosomes were excluded because the analysis of males and females alone would have reduced their power; individuals were excluded if the genotyping rate per individual was <95%. The BEAGLE 4.1 software package was used for phasing genotype data [29].
Wright’s fixation index (FST) and single SNP cross-population extended haplotype homozygosity (XP-EHH) were calculated for the SNPs left after the quality control.
FST is the proportion of genetic diversity due to allele frequency differences between cases and controls groups. Thus, loci showing unusually large amounts of differentiation (high FST values) may identify regions of the genome that have been subject to diversifying selection in the two groups, whereas loci showing unusually small amounts of differentiation (low FST values) may identify regions that have been subject to stabilizing selection [30,31]. FST was calculated using PLINK 1.9. All markers ranking in the top 1% of the empirical distribution of FST values were considered as relevant. This threshold, which corresponds to an usual statistical threshold, was used in other similar studies [32,33] to retain only the highest signals.
XP-EHH is a linkage disequilibrium-based method that compares the lengths of haplotypes (consistent with the allele under selection and the neighbor variants in linkage disequilibrium) at each marker between two different populations, allowing the detection of strong, directional selection of one allele in one of the two populations while remaining polymorphic in the other [34,35]. XP-EHH was calculated using the SELSCAN 1.1.0 software package [36]. Similarly to the FST analysis, all markers within the top 1% of the empirical distribution of normalized XP-EHH values [37,38,39] were considered as relevant. The top 1% SNPs of both the analyses were mapped to the reference genome assembly CanFam3.1.
Since combining multiple independent tests increases power and resolution [35,40], we then compared the results of the aforementioned analyses in order to identify the SNPs found by both of them, namely, the ones that most significantly differentiated between the cases and controls.
The genes containing the relevant SNPs were examined for pathways using Enrichr [41], which looks for associations in various libraries, including KEGG Human 2019, WikiPathways Human 2019 and BioPlanet 2019. Information about the identified genes was obtained via the GeneCards database (https://www.genecards.org/) [42].
Genomic regions subjected to selection may show a reduced nucleotide diversity and increased homozygosity around the selected locus if compared to the rest of genome [43]. Therefore, on a subset of genes consisting of consensus genes and genes involved in the most relevant pathways, runs of homozygosity (ROH) were investigated using a sliding window approach in PLINK 1.9. The sliding window was 50 SNPs long, and a maximum of five missing genotypes and no heterozygous SNPs were tolerated. An ROH was called if the following criteria were fulfilled: (1) 50 or more consecutive homozygous SNPs; (2) a minimum length of 1 Mb; (3) a minimum density of one SNP per 50 kb; (4) a maximum gap between two consecutive SNPs of 100 kb. For each group, the proportion of dogs showing a ROH was calculated. The ROH-based inbreeding coefficient (FROH) was calculated for each animal, dividing the total length of all ROH in its genome by the length of the autosomal genome covered by SNPs on the chip.

3. Results

Our sample consisted of 33 subjects chosen among 90 CKCSs previously investigated by Bagardi et al. [4]. Thirty-one out of 33 (94%) dogs had signs of mitral valve degeneration. Among the controls, two were classified as A and 15 as B1; fourteen cases were included in ACVIM class B1, one in class B2, and one in class D. This subset was composed of 22 females and 11 males. Age ranged from 0.8 to 11.4 years (6.2 ± 2.6 years). Distributions in the ACVIM classes, age, and gender were not significantly different to the whole sample of 90 dogs. According to the selection criteria, age was significantly higher in controls (8.1 ± 0.5 years) than in cases (4.1 ± 0.5 years) (r2 = 0.53, p < 0.0001).
It was possible to get the pedigree information of 20 out of 33 (61%) dogs (81% of cases and 41% of controls). The average relatedness (AR) and average inbreeding (F) coefficients calculated for the dogs with a pedigree were 0.06 and 0.01, respectively. Since pedigree information was not available for all the dogs included in our study, we also used a genetic marker-based method for calculating inbreeding (FROH). F and FROH are not directly comparable, because the pedigree-based coefficient measures the mean expected autozygosity of an individual (identity by descent), whereas the latter measures the realized autozygosity (identity by state). FROH was found to be 0.24 ± 0.03 for the whole sample. No significant differences were found between FROH in cases (0.24 ± 0.04, ranging from 0.18 to 0.28) and in controls (0.23 ± 0.02, ranging from 0.19 to 0.28). These results demonstrate that the inbreeding of the two groups was similar and, therefore, could not influence the results of our genetic analyses. Moreover, the discrepancy between F and FROH was consistent with that observed in other studies on dogs [44,45,46].

3.1. Genomic Analysis

After the quality control, which excluded the SNPs with a low call rate and MAF and those localized on the sex chromosomes, 103,606 SNPs remained for downstream analyses. None of the 33 animals were excluded.
Comparing the case and control groups, there were 291 SNPs characterized by the top 1% values of FST (0.19–0.43), mapping regions containing 157 different genes (Figure 1). There were 152 SNPs characterized by the top 1% values of XP-EHH (2.73–3.82), mapping regions containing 45 different genes (Figure 2). The SNPs shared between both the analyses were mapped to genomic regions containing 10 genes (hereafter called “consensus genes” for conciseness), as reported in Table 2.
ROH analysis revealed that at least 50% of cases and no more than 9% of controls showed homozygous regions around the STEAP2, HEPACAM2, and CDK6 genes. The ROH data, including ARNT2, KIAA1024, and FAH were also interesting, because they were present in the genome of almost 80% of cases but only in approximately 40% of controls. No ROH was found in any group for the genes ADCY9, AXIN1, CACNA1H, CREBBP, PDPK1, SLC8A2, and TRAP1. All details about the proportion of cases and controls in which ROH were found are shown in the Supplementary Materials (Table S1). Figure 3 represents the most relevant genes identified by our analysis.

3.2. Pathway Analysis

To better understand the functions performed by each identified gene, single genes were evaluated in relation to the pathways they are involved in. This allowed classification of 34 of the aforementioned genes in the following most relevant pathways: the Wnt signaling pathway, apelin pathway, hippo signaling pathway, ErbB and epidermal growth factor receptor (EGFR), transforming growth factor β (TGF-β) signaling pathway, endothelins, aldosterone, renin, and body mass index. All details are reported in Table 3.

4. Discussion

Among all dog breeds, CKCSs stand out for being the breed with the highest prevalence and the earliest onset of MMVD [1,3]. Therefore, the purpose of this study was to identify selection signatures that are able to distinguish between subjects with a diagnosis of MMVD at a very young age (before 5 years) and subjects in which this pathology may appear at an older age (after 5 years) or otherwise persist at a milder stage for a long time (after 8 years).
Following the calculation of FROH in cases and controls, we verified that the degree of inbreeding of the two groups was similar. Moreover, our results were superimposable to those found in other studies on various canine breeds [44,45,46]. The identification of consensus genes was reached using two independent methods that analyze different genomic characteristics: FST highlights genomic differences between groups in terms of expected heterozygosity, whereas XP-EHH is based on the comparison of regions of homozygosity that differentiate the groups. Moreover, ROH analysis on the regions surrounding these genes showed a different presence of long homozygous portions of DNA between the cases and controls. From this analysis, we observed that ROH containing the HEPACAM2 and CDK6 genes were present in 50% of cases and only 5% of controls. It could be hypothesized that these genes may be involved in predisposition to rapidly progressing MMVD rather than its early onset. A follow-up of the dogs included in the case group may clarify if the progression of the disease in subjects in which ROH were found around the aforementioned genes is quicker than in subjects in which they were not present. Another possible reason is that the cases in which ROH were absent are heterozygous carriers of the allele that contributes to MMVD predisposition. KIAA1024 and FAH genes were also relevant, because ROH flanking them were found in the genome of 80% of cases but only in approximately 40% of controls. A possible explanation of the presence of controls showing a ROH around these genes is that they might have been misclassified due to the mitral valve pathophysiology, which makes it difficult to detect a real control.
To better describe their possible role in MMVD onset and progression, an accurate investigation about consensus genes was performed. Some of them were shown to be directly or indirectly related to mechanisms already supposed to be involved in the disease’s pathogenesis, such as the TGF-β signaling pathway, or to processes related to heart development or functionality. The most relevant consensus genes are described below.
HEPACAM2 interacts with FGFR1 (fibroblast growth factor receptor 1), which is associated with abnormal heart development. Moreover, during adult life, valves maintain a pool of mesenchymal cells responsive to FGF and producing proteoglycans, which are also increased during MMVD [47]. It should be noted that HEPACAM2 may localize on the region (CFA 14q1.3) that Madsen et al. (2011) found to be associated to MMVD [3].
FAH interacts with ADAMTSL4 (ADAMTS like 4), which is supposed to facilitate FBN1 (fibrillin 1) microfibril biogenesis [48]. FBN1 is one of candidate genes for MMVD predisposition, because it regulates TGF-β signaling and is associated with Marfan syndrome, which represents one of the syndromic forms of human mitral valve prolapse [49].
CDK6 prevents cell proliferation and negatively regulates cell differentiation but is required for the proliferation of specific cell types. Moreover, it interacts with CDKN2B (cyclin-dependent kinase inhibitor 2B), whose expression was found to be induced by TGF-β and is associated with coronary heart disease [50].
EPB41L4B promotes cellular adhesion, migration, and motility in vitro and may play a role in wound healing [51]. This gene interacts with the CASQ2 (calsequestrin 2) gene, which encodes a protein localized in cardiac muscle cells that stores calcium for muscle function [52]. Mutations in this gene cause catecholaminergic polymorphic ventricular tachycardia [53].
It is interesting to note that some genes were associated with height (CDK6 and ZRANB3) or body mass (FRRS1L, EPB41L4B, and ZRANB3) in humans. It has been well established that small dog breeds are predisposed to MMVD [54,55]. Moreover, people affected by mitral valve prolapse tend to have a low body mass index and be leaner and shorter than other individuals [56,57]. A morphometric evaluation of CKCSs could allow identification if the selection for specific physical body features is related to the predisposition to the disease. In fact, the circumference of the thorax has already been negatively correlated with mitral valve prolapse in dachshunds [8].
For the other genes able to significantly distinguish between cases and controls, reported in Table 2, no evident correlation with MMVD was found.
The main pathways associated with the genes identified by this study appear to be involved in processes related to heart development and homeostasis, as reported by several studies on humans, mice, and dogs. For example, during valvulogenesis, TGF-β is fundamental for the formation of endocardial cushions and epithelial-to-mesenchymal transition [58,59]. Cell migration and proliferation in these cushions requires the ErbB and Wnt canonical pathway, the latter being also important for maintaining the aforementioned pool of undifferentiated mesenchymal cells responsive to FGF, even in adult life [47,60,61,62]. Valvular interstitial cell (VIC) activation to myofibroblasts, one the most accredited pathogenetic mechanisms of MMVD, seems to be stimulated by TGF-β and inhibited by apelin [63,64]. Canine VIC exposition to TGF-β3, in fact, has been shown to be able to regulate myofibroblast activation and proteoglycan synthesis in an in vitro system [65]. Moreover, the knockdown of ErbB or Axin2 (which represents one of the principal Wnt canonical pathway inhibitors) determined the development of hyperplastic and myxomatous valves in mice in some studies [66,67]. The TGF-β pathway is involved in many human cardiovascular diseases, such as the Marfan, Ehlers–Danlos, Loeys–Dietz, and aneurysms–osteoarthritis syndromes, whose phenotypes often feature mitral valve prolapse [68,69,70]. It is worth mentioning that a study demonstrated that the expression of TGF-β 1 and 3 in cells and the extracellular matrix (ECM) was increased in MMVD-affected canine mitral valves [71]. It has been speculated that the endothelium can be involved in MMVD pathogenesis as well. Its damage upregulates endothelin and nitric oxide, which are involved in the production and alteration of ECMs [72,73]. Moreover, an increase in endothelin receptors has been proven to be associated with MMVD in canine mitral valves [74,75], and with age and mechanical stress in porcine ones [74,75]. Regarding renin and aldosterone, it is well documented that, when the stroke volume decreases due to MMVD progression, several neuro-endocrine compensatory mechanisms are activated, including the renin–angiotensin–aldosterone system (RAAS) [76]. Although this system contributes to the maintenance of blood pressure in heart failure, it is one of the main targets of MMVD therapy due to its role in the development of congestive heart failure itself [77,78]. Finally, the role of platelet activation has not yet been clarified. Human patients with a mitral valve prolapse seem to be more at risk of thromboembolism [79,80,81] and present increased platelet activation [82,83]. In dogs, instead, no significant association has been found between MMVD and thrombus formation [84], and, in most of the studies, platelet functionality has been decreased in canines [85,86,87].

5. Conclusions

Several mechanisms have been thought to contribute to the development of MMVD, such as the TGF-β signaling pathway [65,71], increases in serotonin receptors [88,89,90,91], alterations of ECM organization [65,71,89,92], endothelial damage [72,84,89,93,94,95], and oxidative stress [89]. The present genomic study involves a relatively small number of cases and controls. The high prevalence of MMVD in CKCSs, even at a young age, as well as the variability in its progress, makes it difficult to recruit a great number of subjects that can be clearly categorizable as either cases or controls. However, the careful selection of cases and controls among a larger population of dogs made it possible to identify genes and pathways potentially involved in the pathogenesis of early-onset MMVD in the CKCS dog breed. Particularly, our findings are consistent with the hypothesis that the role of TGF-β, ECM disruption, and/or endothelin may be relevant and deserve further investigation.
Genomic studies are not only important to identify genes associated with MMVD predisposition, but also to predict the loss of genetic diversity of the breed following the exclusion of reproduction of dogs with specific phenotypical features. Indeed, a breeding program should not exclude more than 50% of the dog population [96] and, in particular, no more than 30% when screening for a single disease [97]. For example, since the prevalence of MMVD is highly dependent on age, it would be important to choose an age limit whereby dogs with early onset of MMVD can be excluded from reproduction. On the other hand, to place the limit at an advanced age would result in the exclusion of an excessive number of animals for breeding.
The investigation of the genetic basis of canine MMVD would surely benefit from an increase in the number of samples, which would also permit the analysis of sexual chromosomes (a X-linked form of mitral valve prolapse has, in fact, been found in human [98]). Moreover, a follow-up study conducted on ACVIM B1 patients could allow us to evaluate how MMVD progresses in these dogs. This will be helpful to identify genomic haplotypes associated with early-onset and rapidly progressing MMVD predispositions that, together with morphometric, clinical, and echocardiographic characterizations, could be used as part of a screening program for CKCSs, defining early selection criteria for the exclusion of subjects from breeding. In this respect, it would be important to include collaboration with the veterinary medical community, which may inform us about dogs that should particularly be included in further genomic studies as either a case or control sample.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/10/10/1895/s1: Table S1. Proportion of cases and controls showing a run of homozygosity in regions flanking selected genes.

Author Contributions

Conceptualization, P.C., P.G.B., and A.B.; methodology, P.C.; software, M.C.; validation, P.C.; formal analysis, A.B., M.C., and A.N.; investigation, A.B., P.C., and M.B.; resources, M.B., C.L., P.G.B., and S.F.; data curation, A.B., M.B., M.C., S.F., and P.C.; writing—original draft preparation, A.B., M.B., M.C., and P.C.; writing—review and editing, A.B., M.B., C.L., P.C., P.G.B., and S.F.; visualization, A.B., M.C., and M.B.; supervision, P.C. and P.G.B.; project administration, P.C. and P.G.B.; funding acquisition, C.L. and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support for this study was provided by University of Milan (LINEA2_CLOCA_2019_AA).

Acknowledgments

The authors are grateful to the many dog owners and breeders who provided samples for this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Detweiler, D.K.; Patterson, D.F. The prevalence and types of cardiovascular disease in dogs. Ann. N. Y. Acad. Sci. 1965, 127, 481–516. [Google Scholar] [CrossRef] [PubMed]
  2. Borgarelli, M.; Crosara, S.; Lamb, K.; Savarino, P.; La Rosa, G.; Tarducci, A.; Häggström, J. Survival Characteristics and Prognostic Variables of Dogs with Preclinical Chronic Degenerative Mitral Valve Disease Attributable to Myxomatous Degeneration. J. Vet. Intern. Med. 2012, 26, 69–75. [Google Scholar] [CrossRef] [PubMed]
  3. Madsen, M.B.; Olsen, L.H.; Häggström, J.; Höglund, K.; Ljungvall, I.; Falk, T.; Wess, G.; Stephenson, H.; Dukes-McEwan, J.; Chetboul, V.; et al. Identification of 2 Loci Associated with Development of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels. J. Hered. 2011, 102, S62–S67. [Google Scholar] [CrossRef]
  4. Bagardi, M.; Bionda, A.; Locatelli, C.; Cortellari, M.; Frattini, S.; Negro, A.; Crepaldi, P.; Brambilla, P.G. Echocardiographic Evaluation of the Mitral Valve in Cavalier King Charles Spaniels. Animals 2020, 10, 1454. [Google Scholar] [CrossRef] [PubMed]
  5. 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. 2006, 20, 894–903. [Google Scholar] [CrossRef]
  6. Pedersen, H.D.; Häggström, J.; Falk, T.; Mow, T.; Olsen, L.H.; Iversen, L.; Jensen, A.L. Auscultation in mild mitral regurgitation in dogs: Observer variation, effects of physical maneuvers, and agreement with color Doppler echocardiography and phonocardiography. J. Vet. Intern. Med. 1999, 13, 56–64. [Google Scholar] [CrossRef]
  7. Keene, B.W.; Atkins, C.E.; Bonagura, J.D.; Fox, P.R.; Häggström, J.; Fuentes, V.L.; Oyama, M.A.; Rush, J.E.; Stepien, R.L.; Uechi, M. ACVIM consensus guidelines for the diagnosis and treatment of myxomatous mitral valve disease in dogs. J. Vet. Intern. Med. 2019, 1–14. [Google Scholar] [CrossRef]
  8. Olsen, L.H.; Fredholm, M.; Pedersen, H.D. Epidemiology an Inheritance of Mitral Valve Prolapse in Dachshunds. J. Vet. Intern. Med. 1999, 13, 448–456. [Google Scholar] [CrossRef]
  9. Lewis, T.W.; Swift, S.; Woolliams, J.A.; Blott, S.C. Heritability of premature mitral valve disease in Cavalier King Charles spaniels. Vet. J. 2011, 188, 73–76. [Google Scholar] [CrossRef]
  10. Swenson, L.; Häggström, J.; Kvart, C.; Juneja, R.K. Relationship between parental cardiac status in Cavalier King Charles spaniels and prevalence and severity of chronic valvular disease in offspring. J. Am. Vet. Med. Assoc. 1996, 208, 2009–2012. [Google Scholar]
  11. Lindblad-Toh, K.; Wade, C.M.; Mikkelsen, T.S.; Karlsson, E.K.; Jaffe, D.B.; Kamal, M.; Clamp, M.; Chang, J.L.; Kulbokas, E.J.; Zody, M.C.; et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 2005, 438, 803–819. [Google Scholar] [CrossRef] [PubMed]
  12. Sutter, N.B.; Eberle, M.A.; Parker, H.G.; Pullar, B.J.; Kirkness, E.F.; Kruglyak, L.; Ostrander, E.A. Extensive and breed-specific linkage disequilibrium in Canis familiaris. Genome Res. 2004, 14, 2388–2396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Parker, H.G.; Meurs, K.M.; Ostrander, E.A. Finding cardiovascular disease genes in the dog. J. Vet. Cardiol. 2006, 8, 115–127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Parker, H.G. The History and Relationships of Dog Breeds. In The Genetics of the Dog; Ostrander, E.A., Ruvinsky, A., Eds.; CAB International: Croydon, UK, 2012; pp. 38–56. [Google Scholar]
  15. Karlsson, E.K.; Baranowska, I.; Wade, C.M.; Salmon Hillbertz, N.H.; Zody, M.C.; Anderson, N.; Biagi, T.M.; Patterson, N.; Pielberg, G.R.; Kulbokas, E.J.; et al. Efficient Mapping of Mendelian Traits in Dogs Through Genome-Wide Association. Nat. Genet. 2007, 39, 1321–1328. [Google Scholar] [CrossRef]
  16. French, A.T.; Ogden, R.; Eland, C.; Hemani, G.; Pong-Wong, R.; Corcoran, B.M.; Summers, K.M. Genome-wide analysis of mitral valve disease in Cavalier King Charles Spaniels. Vet. J. 2012, 193, 283–286. [Google Scholar] [CrossRef]
  17. Meurs, K.M.; Friedenberg, S.G.; Williams, B.; Keene, B.W.; Atkins, C.E.; Adin, D.; Aona, B.; DeFrancesco, T.; Tou, S.; Mackay, T. Evaluation of genes associated with human myxomatous mitral valve disease in dogs with familial myxomatous mitral valve degeneration. Vet. J. 2018, 232, 16–19. [Google Scholar] [CrossRef]
  18. Lee, C.-M.; Han, J.-I.; Kang, M.-H.; Kim, S.-G.; Park, H.-M. Polymorphism in the serotonin transporter protein gene in Maltese dogs with degenerative mitral valve disease. J. Vet. Sci. 2018, 19, 129–135. [Google Scholar] [CrossRef]
  19. Lee, C.-M.; Song, D.-W.; Ro, W.-B.; Kang, M.-H.; Park, H.-M. Genome-wide association study of degenerative mitral valve disease in Maltese dogs. J. Vet. Sci. 2019, 20, 63–71. [Google Scholar] [CrossRef]
  20. Torres-García, O.; Rey-Buitrago, M.; Acosta-Virgüez, E.; Bernal-Rosas, Y.; Infante-González, J.; Gómez-Duarte, L. Role of COL1A2 Gene Polymorphisms in Myxomatous Mitral Valve Disease in Poodle Dogs Genetic Study of Mitral Valve Disease. J. Agric. Vet. Sci. 2016, 9, 113–118. [Google Scholar] [CrossRef]
  21. Stern, J.A.; Hsue, W.; Song, K.H.; Ontiveros, E.S.; Fuentes, V.L.; Stepien, R.L. Severity of mitral valve degeneration is associated with chromosome 15 loci in whippet dogs. PLoS ONE 2015, 10, 1–11. [Google Scholar] [CrossRef]
  22. Lundin, T.; Kvart, C. Evaluation of the Swedish breeding program for cavalier King Charles spaniels. Acta Vet. Scand. 2010, 52, 2–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Swift, S.; Baldin, A.; Cripps, P. Degenerative Valvular Disease in the Cavalier King Charles Spaniel: Results of the UK Breed Scheme 1991–2010. J. Vet. Intern. Med. 2017, 31, 9–14. [Google Scholar] [CrossRef]
  24. Birkegård, A.C.; Reimann, M.J.; Martinussen, T.; Häggström, J.; Pedersen, H.D.; Olsen, L.H. Breeding Restrictions Decrease the Prevalence of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels over an 8- to 10-Year Period. J. Vet. Intern. Med. 2016, 30, 63–68. [Google Scholar] [CrossRef] [PubMed]
  25. Chetboul, V.; Tissier, R. Echocardiographic assessment of canine degenerative mitral valve disease. J. Vet. Cardiol. 2012, 14, 127–148. [Google Scholar] [CrossRef] [PubMed]
  26. Thomas, W.P.; Gaber, C.E.; Jacobs, G.J.; Kaplan, P.M.; Lombard, C.W.; Moise, N.S.; Moses, B.L. Recommendations for Standards in Transthoracic Two-Dimensional Echocardiography in the Dog and Cat. J. Vet. Intern. Med. 1993, 7, 247–252. [Google Scholar] [CrossRef]
  27. Wellmann, R. Optimum contribution selection for animal breeding and conservation: The R package optiSel. BMC Bioinformatics 2019, 20, 20–25. [Google Scholar] [CrossRef]
  28. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A toolset for whole-genome association and population-based linkage analysis. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  29. Browning, S.R.; Browning, B.L. Rapid and accurate haplotype phasing and missing data inference for whole genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 2007, 81, 1084–1097. [Google Scholar] [CrossRef] [Green Version]
  30. Beaumont, M.A.; Balding, D.J. Identifying adaptive genetic divergence among populations from genome scans. Mol. Ecol. 2004, 13, 969–980. [Google Scholar] [CrossRef]
  31. Holsinger, K.E.; Weir, B.S. Genetics in geographically structured populations: Defining, estimating and interpreting FST. Nat. Rev. Genet. 2009, 10, 639–650. [Google Scholar] [CrossRef] [Green Version]
  32. Porto-Neto, L.R.; Sonstegard, T.S.; Liu, G.E.; Bickhart, D.M.; Da Silva, M.V.; Machado, M.A.; Utsunomiya, Y.T.; Garcia, J.F.; Gondro, C.; Van Tassell, C.P. Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping. BMC Genomics 2013, 14, 876. [Google Scholar] [CrossRef] [Green Version]
  33. Ablondi, M.; Viklund, Å.; Lindgren, G.; Eriksson, S.; Mikko, S. Signatures of selection in the genome of Swedish warmblood horses selected for sport performance. BMC Genomics 2019, 20, 717. [Google Scholar] [CrossRef]
  34. Sabeti, P.C.; Reich, D.E.; Higgins, J.M.; Levine, H.Z.P.; Richter, D.J.; Schaffner, S.F.; Gabriel, S.B.; Platko, J.V.; Patterson, N.J.; McDonald, G.J.; et al. Detecting recent positive selection in the human genome from haplotype structure. Nature 2002, 419, 832–837. [Google Scholar] [CrossRef] [PubMed]
  35. Vitti, J.J.; Grossman, S.R.; Sabeti, P.C. Detecting Natural Selection in Genomic Data. Annu. Rev. Genet. 2013, 47, 97–120. [Google Scholar] [CrossRef]
  36. Szpiech, Z.A.; Hernandez, R.D. selscan: An Efficient Multithreaded Program to Perform EHH-Based Scans for Positive Selection. Mol. Biol. Evol. 2014, 31, 2824–2827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kim, J.; Williams, F.J.; Dreger, D.L.; Plassais, J.; Davis, B.W.; Parker, H.G.; Ostrander, E.A. Genetic selection of athletic success in sport-hunting dogs. Proc. Natl. Acad. Sci. USA 2018, 115, E7212–E7221. [Google Scholar] [CrossRef] [Green Version]
  38. Liu, X.; Ong, R.T.-H.; Pillai, E.N.; Elzein, A.M.; Small, K.S.; Clark, T.G.; Kwiatkowski, D.P.; Teo, Y.-Y. Detecting and Characterizing Genomic Signatures of Positive Selection in Global Populations. Am. J. Hum. Genet. 2013, 92, 866–881. [Google Scholar] [CrossRef] [Green Version]
  39. Pitt, D.; Bruford, M.W.; Barbato, M.; Orozco-terWengel, P.; Martínez, R.; Sevane, N. Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics. Evol. Appl. 2019, 12, 105–122. [Google Scholar] [CrossRef] [PubMed]
  40. Grossman, S.R.; Shlyakhter, I.; Shylakhter, I.; Karlsson, E.K.; Byrne, E.H.; Morales, S.; Frieden, G.; Hostetter, E.; Angelino, E.; Garber, M.; et al. A composite of multiple signals distinguishes causal variants in regions of positive selection. Science 2010, 327, 883–886. [Google Scholar] [CrossRef] [Green Version]
  41. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [Green Version]
  42. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinforma. 2016, 54, 1.30.1–1.30.33. [Google Scholar] [CrossRef]
  43. Pemberton, T.J.; Absher, D.; Feldman, M.W.; Myers, R.M.; Rosenberg, N.A.; Li, J.Z. Genomic Patterns of Homozygosity in Worldwide Human Populations. Am. J. Hum. Genet. 2012, 91, 275–292. [Google Scholar] [CrossRef] [Green Version]
  44. Wiener, P.; Sánchez-Molano, E.; Clements, D.N.; Woolliams, J.A.; Haskell, M.J.; Blott, S.C. Genomic data illuminates demography, genetic structure and selection of a popular dog breed. BMC Genomics 2017, 18, 609. [Google Scholar] [CrossRef] [PubMed]
  45. Mortlock, S.-A.; Khatkar, M.S.; Williamson, P. Comparative Analysis of Genome Diversity in Bullmastiff Dogs. PLoS ONE 2016, 11, e0147941. [Google Scholar] [CrossRef] [Green Version]
  46. Sams, A.J.; Boyko, A.R. Fine-Scale Resolution of Runs of Homozygosity Reveal Patterns of Inbreeding and Substantial Overlap with Recessive Disease Genotypes in Domestic Dogs. G3 (Bethesda) 2019, 9, 117–123. [Google Scholar] [CrossRef] [Green Version]
  47. Bosada, F.M.; Devasthali, V.; Jones, K.A.; Stankunas, K. Wnt/β-catenin signaling enables developmental transitions during valvulogenesis. Development 2016, 143, 1041–1054. [Google Scholar] [CrossRef] [Green Version]
  48. Gabriel, L.A.R.; Wang, L.W.; Bader, H.; Ho, J.C.; Majors, A.K.; Hollyfield, J.G.; Traboulsi, E.I.; Apte, S.S. ADAMTSL4, a Secreted Glycoprotein Widely Distributed in the Eye, Binds Fibrillin-1 Microfibrils and Accelerates Microfibril Biogenesis. Investig. Ophthalmol. Vis. Sci. 2012, 53, 461. [Google Scholar] [CrossRef] [PubMed]
  49. Dietz, H.C.; Cutting, C.R.; Pyeritz, R.E.; Maslen, C.L.; Sakai, L.Y.; Corson, G.M.; Puffenberger, E.G.; Hamosh, A.; Nanthakumar, E.J.; Curristin, S.M.; et al. Marfan syndrome caused by a recurrent de novo missense mutation in the fibrillin gene. Nature 1991, 352, 337–339. [Google Scholar] [CrossRef]
  50. Motterle, A.; Pu, X.; Wood, H.; Xiao, Q.; Gor, S.; Ng, F.L.; Chan, K.; Cross, F.; Shohreh, B.; Poston, R.N.; et al. Functional analyses of coronary artery disease associated variation on chromosome 9p21 in vascular smooth muscle cells. Hum. Mol. Genet. 2012, 21, 4021–4029. [Google Scholar] [CrossRef] [Green Version]
  51. Bosanquet, D.C.; Ye, L.; Harding, K.G.; Jiang, W.G. Expressed in high metastatic cells (Ehm2) is a positive regulator of keratinocyte adhesion and motility: The implication for wound healing. J. Dermatol. Sci. 2013, 71, 115–121. [Google Scholar] [CrossRef] [PubMed]
  52. Otsu, K.; Fujii, J.; Periasamy, M.; Difilippantonio, M.; Uppender, M.; Ward, D.C.; MacLennan, D.H. Chromosome mapping of five human cardiac and skeletal muscle sarcoplasmic reticulum protein genes. Genomics 1993, 17, 507–509. [Google Scholar] [CrossRef] [PubMed]
  53. di Barletta, M.R.; Viatchenko-Karpinski, S.; Nori, A.; Memmi, M.; Terentyev, D.; Turcato, F.; Valle, G.; Rizzi, N.; Napolitano, C.; Gyorke, S.; et al. Clinical phenotype and functional characterization of CASQ2 mutations associated with catecholaminergic polymorphic ventricular tachycardia. Circulation 2006, 114, 1012–1019. [Google Scholar] [CrossRef] [Green Version]
  54. Serfass, P.; Chetboul, V.; Carlos Sampedrano, C.; Nicolle, A.P.; Benalloul, T.; Laforge, H.; Gau, C.; Hébert, C.; Pouchelon, J.-L.; Tissier, R. Retrospective study of 942 small-sized dogs: Prevalence of left apical systolic heart murmur and left-sided heart failure, critical effects of breed and sex. J. Vet. Cardiol. 2006, 8, 11–18. [Google Scholar] [CrossRef]
  55. Thrusfield, M.V.; Aitken, C.G.G.; Darker, P.G.G. Observations on breed and sex in relation to canine heart valve incompetence. J. Small Anim. Pract. 1985, 26, 709–717. [Google Scholar] [CrossRef]
  56. Schutte, J.E.; Gaffney, F.A.; Blend, L.; Blomqvist, C.G. Distinctive anthropometric characteristics of women with mitral valve prolapse. Am. J. Med. 1981, 71, 533–538. [Google Scholar] [CrossRef]
  57. Freed, L.A.; Levy, D.; Levine, R.A.; Larson, M.G.; Evans, J.C.; Fuller, D.L.; Lehman, B.; Benjamin, E.J. Prevalence and clinical outcome of mitral-valve prolapse. N. Engl. J. Med. 1999, 341, 1–7. [Google Scholar] [CrossRef] [PubMed]
  58. Potts, J.D.; Runyan, R.B. Epithelial-mesenchymal cell transformation in the embryonic heart can be mediated, in part, by transforming growth factor beta. Dev. Biol. 1989, 134, 392–401. [Google Scholar] [CrossRef]
  59. Azhar, M.; Runyan, R.B.; Gard, C.; Sanford, L.P.; Miller, M.L.; Andringa, A.; Pawlowski, S.; Rajan, S.; Doetschman, T. Ligand-specific function of transforming growth factor beta in epithelial-mesenchymal transition in heart development. Dev. Dyn. 2009, 238, 431–442. [Google Scholar] [CrossRef] [Green Version]
  60. Camenisch, T.D.; Schroeder, J.A.; Bradley, J.; Klewer, S.E.; McDonald, J.A. Heart-valve mesenchyme formation is dependent on hyaluronan-augmented activation of ErbB2–ErbB3 receptors. Nat. Med. 2002, 8, 850–855. [Google Scholar] [CrossRef]
  61. Chen, B.; Bronson, R.T.; Klaman, L.D.; Hampton, T.G.; Wang, J.; Green, P.J.; Magnuson, T.; Douglas, P.S.; Morgan, J.P.; Neel, B.G. Mice mutant for Egfr and Shp2 have defective cardiac semilunar valvulogenesis. Nat. Genet. 2000, 24, 296–299. [Google Scholar] [CrossRef]
  62. Lee, K.-F.; Simon, H.; Chen, H.; Bates, B.; Hung, M.-C.; Hauser, C. Requirement for neuregulin receptor erbB2 in neural and cardiac development. Nature 1995, 378, 394–398. [Google Scholar] [CrossRef]
  63. Walker, G.A.; Masters, K.S.; Shah, D.N.; Anseth, K.S.; Leinwand, L.A. Valvular myofibroblast activation by transforming growth factor-beta: Implications for pathological extracellular matrix remodeling in heart valve disease. Circ. Res. 2004, 95, 253–260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Pchejetski, D.; Foussal, C.; Alfarano, C.; Lairez, O.; Calise, D.; Guilbeau-Frugier, C.; Schaak, S.; Seguelas, M.-H.; Wanecq, E.; Valet, P.; et al. Apelin prevents cardiac fibroblast activation and collagen production through inhibition of sphingosine kinase 1. Eur. Heart J. 2012, 33, 2360–2369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Obayashi, K.; Miyagawa-Tomita, S.; Matsumoto, H.; Koyama, H.; Nakanishi, T.; Hirose, H. Effects of transforming growth factor-β3 and matrix metalloproteinase-3 on the pathogenesis of chronic mitral valvular disease in dogs. Am. J. Vet. Res. 2011, 72, 194–202. [Google Scholar] [CrossRef]
  66. Barrick, C.J.; Roberts, R.B.; Rojas, M.; Rajamannan, N.M.; Suitt, C.B.; O’Brien, K.D.; Smyth, S.S.; Threadgill, D.W. Reduced EGFR causes abnormal valvular differentiation leading to calcific aortic stenosis and left ventricular hypertrophy in C57BL/6J but not 129S1/SvImJ mice. Am. J. Physiol. Circ. Physiol. 2009, 297, H65–H75. [Google Scholar] [CrossRef]
  67. Hulin, A.; Moore, V.; James, J.M.; Yutzey, K.E. Loss of Axin2 results in impaired heart valve maturation and subsequentmyxomatous valve disease. Cardiovasc. Res. 2017, 113, 40–51. [Google Scholar] [CrossRef] [PubMed]
  68. Loeys, B.L.; Schwarze, U.; Holm, T.; Callewaert, B.L.; Thomas, G.H.; Pannu, H.; De Backer, J.F.; Oswald, G.L.; Symoens, S.; Manouvrier, S.; et al. Aneurysm Syndromes Caused by Mutations in the TGF-β Receptor. N. Engl. J. Med. 2006, 355, 788–798. [Google Scholar] [CrossRef]
  69. van de Laar, I.M.B.H.; Oldenburg, R.A.; Pals, G.; Roos-Hesselink, J.W.; de Graaf, B.M.; Verhagen, J.M.A.; Hoedemaekers, Y.M.; Willemsen, R.; Severijnen, L.-A.; Venselaar, H.; et al. Mutations in SMAD3 cause a syndromic form of aortic aneurysms and dissections with early-onset osteoarthritis. Nat. Genet. 2011, 43, 121–126. [Google Scholar] [CrossRef]
  70. Matt, P.; Schoenhoff, F.; Habashi, J.; Holm, T.; Van Erp, C.; Loch, D.; Carlson, O.D.; Griswold, B.F.; Fu, Q.; De Backer, J.F.; et al. Circulating transforming growth factor-beta in Marfan syndrome. Circulation 2009, 120, 526–532. [Google Scholar] [CrossRef] [Green Version]
  71. Aupperle, H.; März, I.; Thielebein, J.; Schoon, H.A. Expression of Transforming Growth Factor-β1, -β2 and -β3 in Normal and Diseased Canine Mitral Valves. J. Comp. Pathol. 2008, 139, 97–107. [Google Scholar] [CrossRef]
  72. Guarda, E.; Katwa, L.C.; Myers, P.R.; Tyagi, S.C.; Weber, K.T. Effects of endothelins on collagen turnover in cardiac fibroblasts. Cardiovasc. Res. 2008, 27, 2130–2134. [Google Scholar] [CrossRef]
  73. Myers, P.R.; Tanner, M.A. Vascular endothelial cell regulation of extracellular matrix collagen: Role of nitric oxide. Arterioscler. Thromb. Vasc. Biol. 1998, 18, 717–722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Pedersen, L.G.; Zhao, J.; Yang, J.; Thomsen, P.D.; Gregersen, H.; Hasenkam, J.M.; Smerup, M.; Pedersen, H.D.; Olsen, L.H. Increased expression of endothelin B receptor in static stretch exposed porcine mitral valve leaflets. Res. Vet. Sci. 2007, 82, 232–238. [Google Scholar] [CrossRef] [PubMed]
  75. Pedersen, L.G.; Offenberg, H.; Moesgaard, S.G.; Thomsen, P.D.; Pedersen, H.D.; Olsen, L.H. Transcription levels of endothelin-1 and endothelin receptors are associated with age and leaflet location in porcine mitral valves. J. Vet. Med. Ser. A Physiol. Pathol. Clin. Med. 2007, 54, 113–118. [Google Scholar] [CrossRef]
  76. Borgarelli, M.; Crosara, S. Malattia cronica mitralica. In Manuale di cardiologia del cane e del gatto; Santilli, R.A., Bussadori, C., Borgarelli, M., Eds.; Elsevier Srl: Milan, Italy, 2012; pp. 153–164. [Google Scholar]
  77. Riegger, G.A.J.; Liebau, G.; Holzschuh, M.; Witkowski, D.; Steilner, H.; Kochsiek, K. Role of the renin-angiotensin system in the development of congestive heart failure in the dog as assessed by chronic converting-enzyme blockade. Am. J. Cardiol. 1984, 53, 614–618. [Google Scholar] [CrossRef]
  78. Brilla, C.G.; Rupp, H.; Funck, R.; Maisch, B. The renin-angiotensin-aldosterone system and myocardial collagen matrix remodelling in congestive heart failure. Eur. Heart J. 1995, 16 (Suppl. O), 107–109. [Google Scholar] [CrossRef]
  79. Avierinos, J.-F.; Brown, R.D.; Foley, D.A.; Nkomo, V.; Petty, G.W.; Scott, C.; Enriquez-Sarano, M. Cerebral ischemic events after diagnosis of mitral valve prolapse: A community-based study of incidence and predictive factors. Stroke 2003, 34, 1339–1344. [Google Scholar] [CrossRef] [Green Version]
  80. Caltrider, N.D.; Irvine, A.R.; Kline, H.J.; Rosenblatt, A. Retinal Emboli in Patients with Mitral Valve Prolapse. Am. J. Ophthalmol. 1980, 90, 534–539. [Google Scholar] [CrossRef]
  81. Walsh, P.N.; Kansu, T.A.; Corbett, J.J.; Savion, P.J.; Goldburgh, W.P.; Schatz, N.J. Platelets, thromboembolism and mitral valve prolapse. Circulation 1981, 63, 552–559. [Google Scholar] [CrossRef] [Green Version]
  82. Riddle, J.M.; Stein, P.D.; Magilligan, D.J.; McElroy, H.H. Evaluation of platelet reactivity in patients with valvular heart disease. J. Am. Coll. Cardiol. 1983, 1, 1381–1384. [Google Scholar] [CrossRef] [Green Version]
  83. Tse, H.F.; Lau, C.P.; Cheng, G. Relation between mitral regurgitation and platelet activation. J. Am. Coll. Cardiol. 1997, 30, 1813–1818. [Google Scholar] [CrossRef] [Green Version]
  84. Corcoran, B.M.; Black, A.; Anderson, H.; Dukes-McEwan, J.; French, A.T.; Smith, P.; Devine, C. Identification of surface morphologic changes in the mitral valve leaflets and chordae tendineae of dogs with myxomatous degeneration. Am. J. Vet. Res. 2004, 65, 198–206. [Google Scholar] [CrossRef]
  85. Tanaka, R.; Murota, A.; Nagashima, Y.; Yamane, Y. Changes in Platelet Life Span in Dogs with Mitral Valve Regurgitation. J. Vet. Intern. Med. 2002, 16, 446–451. [Google Scholar] [CrossRef]
  86. Tarnow, I.; Kristensen, A.T.; Texel, H.; Olsen, L.H.; Pedersen, H.D. Decreased Platelet Function in Cavalier King Charles Spaniels with Mitral Valve Regurgitation. J. Vet. Intern. Med. 2003, 17, 680–686. [Google Scholar] [CrossRef]
  87. Tanaka, R.; Yamane, Y. Platelet aggregation in dogs with mitral valve regurgitation. Am. J. Vet. Res. 2000, 61, 1248–1251. [Google Scholar] [CrossRef]
  88. Lu, C.-C.; Liu, M.-M.; Culshaw, G.J.; Clinton, M.; Argyle, D.J.; Corcoran, B.M. Gene network and canonical pathway analysis in canine myxomatous mitral valve disease: A microarray study. Vet. J. 2015, 204, 23–31. [Google Scholar] [CrossRef] [Green Version]
  89. Oyama, M.A.; Chittur, S.V. Genomic expression patterns of mitral valve tissues from dogs with degenerative mitral valve disease. Am. J. Vet. Res. 2006, 67, 1307–1318. [Google Scholar] [CrossRef]
  90. Cremer, S.E.; Moesgaard, S.G.; Rasmussen, C.E.; Zois, N.E.; Falk, T.; Reimann, M.J.; Cirera, S.; Aupperle, H.; Oyama, M.A.; Olsen, L.H. Alpha-smooth muscle actin and serotonin receptors 2A and 2B in dogs with myxomatous mitral valve disease. Res. Vet. Sci. 2015, 100, 197–206. [Google Scholar] [CrossRef] [PubMed]
  91. Arndt, J.W.; Reynolds, C.A.; Singletary, G.E.; Connolly, J.M.; Levy, R.J.; Oyama, M.A. Serum Serotonin Concentrations in Dogs with Degenerative Mitral Valve Disease. J. Vet. Intern. Med. 2009, 23, 1208–1213. [Google Scholar] [CrossRef]
  92. Aupperle, H.; Disatian, S. Pathology, protein expression and signaling in myxomatous mitral valve degeneration: Comparison of dogs and humans. J. Vet. Cardiol. 2012, 14, 59–71. [Google Scholar] [CrossRef] [PubMed]
  93. Lu, C.-C.; Liu, M.-M.; Culshaw, G.J.; French, A.T.; Corcoran, B.M. Comparison of cellular changes in Cavalier King Charles spaniel and mixed breed dogs with myxomatous mitral valve disease. J. Vet. Cardiol. 2016, 18, 100–109. [Google Scholar] [CrossRef] [PubMed]
  94. Takuwa, N.; Takuwa, Y.; Yanagisawa, M.; Yamashita, K.; Masaki, T. A novel vasoactive peptide endothelin stimulates mitogenesis through inositol lipid turnover in Swiss 3T3 fibroblasts. J. Biol. Chem. 1989, 264, 7856–7861. [Google Scholar] [PubMed]
  95. Mow, T.; Pedersen, H.D. Increased Endothelin-Receptor Density in Myxomatous Canine Mitral Valve Leaflets. J. Cardiovasc. Pharmacol. 1999, 34, 254–260. [Google Scholar] [CrossRef] [PubMed]
  96. Hedhammar, Å.A.; Indrebø, A. Rules, regulations, strategies and activities within the Fédération Cynologique Internationale (FCI) to promote canine genetic health. Vet. J. 2011, 189, 141–146. [Google Scholar] [CrossRef]
  97. Häggström, J. Chronic Valvular Disease in Cavalier King Charles Spaniels—Epidemiology, Inheritance and Pathophysiology. Ph.D. Thesis, Swedish University of Agricultural Sciences, Uppsala Sweden, 1996. [Google Scholar]
  98. Kyndt, F.; Gueffet, J.-P.; Probst, V.; Jaafar, P.; Legendre, A.; Le Bouffant, F.; Toquet, C.; Roy, E.; McGregor, L.; Lynch, S.A.; et al. Mutations in the gene encoding filamin A as a cause for familial cardiac valvular dystrophy. Circulation 2007, 115, 40–49. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Manhattan plots of the Wright’s fixation index (FST) analysis comparing cases vs. controls. Each single-nucleotide polymorphism (SNP) is represented by a dot. Each chromosome is represented by a different color. The dotted red line represents the cutoff value of the top 1%, equal to 0.19. The genomic regions in which SNPs were mapped over the dotted red line were considered as significantly different between the two compared groups.
Figure 1. Manhattan plots of the Wright’s fixation index (FST) analysis comparing cases vs. controls. Each single-nucleotide polymorphism (SNP) is represented by a dot. Each chromosome is represented by a different color. The dotted red line represents the cutoff value of the top 1%, equal to 0.19. The genomic regions in which SNPs were mapped over the dotted red line were considered as significantly different between the two compared groups.
Animals 10 01895 g001
Figure 2. Manhattan plots of the cross-population extended haplotype homozygosity (XP-EHH) analysis comparing cases vs. controls. Manhattan plot of log10 XP-EHH values. Each SNP is represented by a dot. Each chromosome is represented by a different color. The dotted red line represents the cutoff value of the top 1%, equal to 2.73. The genomic regions in which SNPs were mapped over the dotted red line were considered as significantly different between the two compared groups.
Figure 2. Manhattan plots of the cross-population extended haplotype homozygosity (XP-EHH) analysis comparing cases vs. controls. Manhattan plot of log10 XP-EHH values. Each SNP is represented by a dot. Each chromosome is represented by a different color. The dotted red line represents the cutoff value of the top 1%, equal to 2.73. The genomic regions in which SNPs were mapped over the dotted red line were considered as significantly different between the two compared groups.
Animals 10 01895 g002
Figure 3. Venn diagram representing the genes identified by the FST, XP-EHH, and runs of homozygosity (ROH) analyses. The SNPs with the top 1% values of FST and XP-EHH were mapped on genomic regions containing 157 and 45 genes, respectively. The intersection between the two circles represents the consensus genes. For a subset of 44 genes, highlighted by the aforementioned analyses (see Supplementary Materials, Table S1), ROH were investigated: ROH including the STEAP2, HEPACAM2, and CDK6 genes were found in at least 50% of cases and no more than 9% of controls; ROH including the ARNT2, KIAA1024, and FAH were present in almost 80% of cases but only in approximately 40% of controls.
Figure 3. Venn diagram representing the genes identified by the FST, XP-EHH, and runs of homozygosity (ROH) analyses. The SNPs with the top 1% values of FST and XP-EHH were mapped on genomic regions containing 157 and 45 genes, respectively. The intersection between the two circles represents the consensus genes. For a subset of 44 genes, highlighted by the aforementioned analyses (see Supplementary Materials, Table S1), ROH were investigated: ROH including the STEAP2, HEPACAM2, and CDK6 genes were found in at least 50% of cases and no more than 9% of controls; ROH including the ARNT2, KIAA1024, and FAH were present in almost 80% of cases but only in approximately 40% of controls.
Animals 10 01895 g003
Table 1. Summary of previous genomic studies about myxomatous mitral valve degeneration (MMVD).
Table 1. Summary of previous genomic studies about myxomatous mitral valve degeneration (MMVD).
AuthorsResultsDog BreedSample SizeCriteria for the InclusionDiagnostic TechniquesGenomic Analysis
AgeDiagnosis
Madsen et al., 2011CFA 13q2.2.3
CFA 14q1.3
CKCS139 cases<4.5 years
<8 years
Murmur ≥1/6 and ARJ/LAA ≥20%
Heart failure symptoms
Auscultation
Echocardiography
GWAS
102 controls>8 yearsMurmur ≤2/6 and ARJ/LAA ≤50%
French et al., 2012No mutations at a single genetic locus were foundCKCS18 early-onset<5 yearsDetectable murmurAuscultationHomozygosity mapping
GWAS
18 late-onset>7 yearsDetectable murmur
Stern et al., 2015FSTL5, EEF1a1a, NAF1, NPY1R, NPY5R, TMA16, March1, ARHGAP26Whippet138 dogs5 years cut-offScoring based on age, presence and degree of mitral valve prolapse, regurgitation, and left heart enlargementAuscultation
Echocardiography
Cumulative echocardiographic score system
GWAS
Torres-García et al., 2016Allele T of the rs22372411 variant of COL1A2Poodle50 casesNo restrictions
>8 years
Diagnosis of MMVD
MMVD absent or mild
Auscultation
Echocardiography
Candidate gene polymorphisms
80 controls>8 yearsMMVD absent or mild
Meurs et al., 2018A missense mutation of COL5A1, predicted to be benign, was present in CKCSCKCS and dachshund 10 CKCSs and 10 dachshunds as casesNo restrictionsDiagnosis of MMVDAuscultation
Echocardiography
Candidate gene approach, whole genome sequencing
98 medium to large dog breeds as controlsNo restrictionsPhenotype not evaluated; low prevalence of MMVD
Lee et al., 2018SERT (SLC6A4): c.1193delT (p.Val397Gly)Maltese20 casesNo restrictionsDiagnosis of MMVDEchocardiographyCandidate gene polymorphisms
10 controlsNo restrictionsEchocardiographically healthy
Lee et al., 2019PDZD2, CTNNA3, LDLRAD4, ARVCFMaltese32 casesNo restrictionsDiagnosis of MMVDEchocardiography
Radiography
GWAS
16 controls>10 yearsEchocardiographically healthy
CFA: canine chromosome; CKCS: Cavalier King Charles spaniel; ARJ/LAA: area of regurgitant jet/left atrium area ratio; GWAS: genome-wide association study.
Table 2. Consensus genes common to the FST and XP-EHH analyses: names and chromosomal coordinates.
Table 2. Consensus genes common to the FST and XP-EHH analyses: names and chromosomal coordinates.
Gene NameCFAStartEndComplete Name
KIAA102435773974057748234KIAA1024
TBC1D1435900376659097331TBC1 domain family member 14
FAH35730058357326453Fumarylacetoacetate hydrolase
FRRS1L116428189564447650Ferric chelate reductase 1 like
EPB41L4B116431286264447436Erythrocyte membrane protein band 4.1 like 4B
CDK6141818842918420100Cyclin dependent kinase 6
HEPACAM2141869574418735539HEPACAM family member 2
RAB3GAP1193786198537957007RAB3 GTPase activating protein catalytic subunit 1
ZRANB3193800219438302593Zinc finger RANBP2-type containing 3
UBXN4193851992738568223UBX domain protein 4
Table 3. Pathways and heart diseases associated by Enrichr to the genes identified by our genomic analyses.
Table 3. Pathways and heart diseases associated by Enrichr to the genes identified by our genomic analyses.
Pathway or Diseasep-ValueAdjusted p-ValueAssociated Group of GenesLibraryAssociated Genes
Wnt signaling pathway0.0020.040Top 1% FSTKEGG 2019 humanCREBBP, TCFL1, SMAD3, AXIN1, WNT2, PLCB2
0.00030.031Top 1% FST + XP-EHHWikipathways 2019, mouseCREBBP, TCFL1, PPP2R2C, AXIN1, PRKD1, WNT2
Hippo signaling pathway0.0020.040Top 1% FSTKEGG 2019 HumanLATS1, TCF7L1, SMAD3, AXIN1, CTNNA3, WNT2
0.0010.028Top 1% FST + XP-EHHKEGG 2019 HumanLATS1, TCFL1, SMAD3, PPP2R2C, AXIN1, CTNNA3, WNNT2
Apelin 0.0010.026Top 1% FSTKEGG 2019 HumanADCY9, SMAD3, ITPR2, ADCY2, PLCB2, SLC8A2
0.0020.050Top 1% FST + XP-EHHKEGG 2019 HumanADCY9, SMAD3, ITPR2, ADCY2, PLCB2, SLC8A2
ErbB and EGFR0.00020.047Top 1% FSTBioPlanet 2019ADCY9, PDPK1, PDE1A, ITPR2, NRG1, ADCY2
TGF-β 0.0091.000Top 1% XP-EHHBioPlanet 2019ARNT2, CDK6, LPAR1, CTNNAL1, STEAP2
Body Mass Index0.00020.313ConsensusGWAS catalog 2019ZRANB3, EPB41L4B, FRRS1L, UBXN4, RAB3GAP1
Endothelins0.0030.158Top 1% FST + XP-EHHBioPlanet 2019 ADCY9, ADCY2, PLCB2, BCAR1
Aldosterone0.00010.008Top 1% FSTKEGG 2019 HumanADCY9, ITPR2, ADCY2, PRKD1, PLCB2, CACNA1H
Renin0.00020.008Top 1% FSTKEGG 2019 HumanADCYAP1R1, PDE1A, PDE3A, ITPR2, PLCB2
Platelet activation0.0030.046Top 1% FSTKEGG 2019 HumanADCY9, ITPR2, ADCY2, TLN2, PLCB2
Significant adjusted p-values and consensus genes are indicated in bold. EGFR: epidermal growth factor receptor; TGF-β: transforming growth factor β. Adjusted p-values were automatically calculated by Enrichr using the Benjamini–Hochberg method for correction for multiple hypotheses testing.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bionda, A.; Cortellari, M.; Bagardi, M.; Frattini, S.; Negro, A.; Locatelli, C.; Brambilla, P.G.; Crepaldi, P. A Genomic Study of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels. Animals 2020, 10, 1895. https://doi.org/10.3390/ani10101895

AMA Style

Bionda A, Cortellari M, Bagardi M, Frattini S, Negro A, Locatelli C, Brambilla PG, Crepaldi P. A Genomic Study of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels. Animals. 2020; 10(10):1895. https://doi.org/10.3390/ani10101895

Chicago/Turabian Style

Bionda, Arianna, Matteo Cortellari, Mara Bagardi, Stefano Frattini, Alessio Negro, Chiara Locatelli, Paola Giuseppina Brambilla, and Paola Crepaldi. 2020. "A Genomic Study of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels" Animals 10, no. 10: 1895. https://doi.org/10.3390/ani10101895

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop