Comparative transcriptomic profiling of myxomatous mitral valve disease in the Cavalier King Charles Spaniel

doi:10.21203/rs.3.rs-23508/v1

Download PDF

Research article

Comparative transcriptomic profiling of myxomatous mitral valve disease in the Cavalier King Charles Spaniel

https://doi.org/10.21203/rs.3.rs-23508/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 23 Sep, 2020

Read the published version in BMC Veterinary Research →

You are reading this older preprint version

Read the latest preprint version →

Background.

All dogs develop myxomatous mitral valve disease by the end of their life, but the Cavalier King Charles Spaniel (CKCS) has a heightened susceptibility, frequently resulting in death at a young age and suggesting that there is a genetic component to the condition in this breed. Transcriptional profiling can reveal the impact of genetic variation through differences in gene expression levels. The aim of this study was to determine whether expression patterns were different in CKCS diseased mitral valves compared to valves from non-CKCS dogs.

Results.

Gene expression patterns in three groups of canine valves resulted in distinct separation of normal valves, CKCS diseased valves and diseased valves from other breeds; the latter were more similar to the normal valves than were the CKCS valves. CKCS valve gene expression patterns were quite different from those in the other dogs, both affected and normal. Patterns in all diseased valves (CKCS and other breeds) were also somewhat different from normal non-diseased samples. Analysis of differentially expressed genes showed enrichment in GO terms relating to cardiac development and function and to calcium signalling canonical pathway in the genes down-regulated in CKCS, compared to normal valves and to diseased valves from other breeds. F2 (prothrombin) (CKCS valves compared to normal) and MEF2C pathway activation (CKCS valves compared to non-CKCS) had the strongest association with the gene changes. A large number of DEGs in the CKCSs were associated with cardiomyocytes including CASQ2, TNNI3 and RYR2.

Conclusion.

Transcriptomic profiling identified gene expression changes in CKCS diseased valves that were not present in age and disease severity-matched non-CKCS valves. These genes are associated with cardiomyocytes, coagulation and extra-cellular matrix remodelling. Identification of genes that vary in the CKCS will allow exploration of genetic variation to understand the aetiology of the disease in this breed, and ultimately development of breeding strategies to eliminate this disease from the breed.

Small Animal Medicine

Genes expression

gene clustering

gene networks

myxomatous mitral valve disease

Cavalier King Charles Spaniel

Myxomatous mitral valve disease (MMVD) is very common in elderly dogs of all genetic backgrounds [1, 2]. The Cavalier King Charles Spaniel (CKCS) breed has a particularly high prevalence of myxomatous mitral valve disease (MMVD), with earlier development and a more rapid progression to severe end-stage disease compared to most other dog breeds [3-6]. This breed association suggests that there is a genetic basis for MMVD in the CKCS. Interestingly, the CKCS breed exhibits unique physiological differences from other breeds such as defects in platelet function, macro-thrombocytopenia and elevated circulating 5-hydroxytryptamine (5HT, also known as serotonin) levels [8-12]. A rare familial macro-thrombocytopenia has been found to be associated with cardiac valvulopathies in children, and there is evidence to suggest dysregulation of serotonergic signalling might contribute to the aetiopathogenesis of MMVD in both the dog and human [13]. For example, rare instances of acquired valvulopathies have been reported in human patients with carcinoid syndrome (excessive circulating 5HT) and chronic usage of appetite suppressants and anti-Parkinsonian drugs that target 5HT receptors [14-16]. Furthermore, there is evidence of enhanced 5HT signalling in both human and canine MMVD valve tissue and valve interstitial cells (VIC) exposed to tensile strain, and genetic variation in the exons of the 5HT transporter gene (SLC6A4, also known as SERT) may be associated with MMVD development in the Maltese terrier [18, 19]. Changes in expression of the 5HT receptor gene are reported for the canine mitral valve transcriptome, and in the dog HTR2B expression is associated with disease progression (Markby et al, in press; doi: 10.3389/fgene.2020.00372) [20-22]. However, 5HT does not induce disease phenotype in cultured VICs and expression of 5HT receptor genes is controlled by TGFβ1 [23]. Nevertheless, taken together, the evidence would suggest that the CKCS has a variant of MMVD that differs somewhat from other dogs, despite the end-stage valve pathology and outcome being the same [24].

For some time it has also been recognised that MMVD has a high degree of heritability in the CKCS and the success of breeding programmes increasing the age of onset support this [5]. One of the only two reported genome-wide association studies (GWAS) identified two loci with a total of thirty one protein coding genes within these regions that could be affected, but follow up on these proteins has not been published [25, 26]. Attempts to find similar genetic mutations in CKCSs and Dachshunds to those reported for human familial MMVD have been unsuccessful [7]. While features of the condition show high heritability in certain breeds, it is unlikely that a single or small number of genetic variants are implicated.

There are limited tissue gene expression data for canine MMVD, but one study reported the transcriptomic profile of aged CKCS with advanced disease [20]. The analysis found that the most affected biologically relevant functions were inflammatory/immune response, cellular movement, cardiovascular development, extracellular matrix organization and epithelial-to-mesenchymal transition (EMT), not dissimilar to the one other canine report not involving CKCSs [20, 21]. Of particular note were changes in expression of a range of cathepsins, matrix metalloproteinases and ADAMTS (a disintegrin and metalloproteinase with thrombospondin motif) family members, important in extracellular matrix remodelling, and a range of genes involved in endothelial-to-mesenchymal transition (EndoMT) [20, 27]. More recently we have shown the canine MMVD valve transcriptome changes progressively as the disease worsens in a range of dog breeds, but also have identified by sample-to-sample network analysis that the Cavalier King Charles Spaniel (CKCS) valve transcriptome was distinct from that seen with non-CKCSs with similar disease severity (Markby et al, in press; doi: 10.3389/fgene.2020.00372).

While global transcriptomic profiling of all dogs can give insight into the possible signalling pathways involved in MMVD pathogenesis, it does not clarify any breed-specific changes. Considering the prevalence of MMVD in the CKCS, the physiological differences that might impact on MMVD pathogenesis and the apparent differential clustering of gene changes compared to non-CKCSs, the aim of this study was to compare the mitral valve transcriptome profile of the CKCS to non-CKCS with severe to very severe (Whitney grade 3-4) MMVD, and against the profile of normal non-CKCS young adult dogs.

Characteristics of the dogs in the study

Details of the animals used in the study are presented in Table 1. There was no significant difference in age between the two diseased groups, but both diseased groups were significantly older than the normal group (P<0.001). Both sexes were represented in all groups and gender did not have an effect on gene expresison. Six Grade 3 and 4 valves were from CKCS and five were from a range of other breeds. The normal dogs were Beagles and mixed breeds.

Network analysis of samples

To compare normal and diseased valves we initially performed a sample-to-sample analysis using Graphia Pro network analysis software, with the Pearson correlation coefficient threshold set at 0.98. This showed that nodes representing grade 3 and 4 samples were separated from the normal samples, and there was also a clear separation of CKCS samples away from the non-CKCS diseased samples, with the non-CKCS diseased valve samples more similar to the normal valves than were the CKCS diseased valves (Figure 1). This suggested that the diseased valves from non-CKCS and CKCS could be separated on the basis of gene expression. We therefore constructed a gene co-expression network (GCN) [28] using a high stringency threshold correlation coefficient of 0.9. 7,247 genes making 20,612 edges were included in the analysis. As we have seen in other studies [29, 30] there was considerable individual dog-specific variation and many clusters were driven by high expression of a subset of genes in a single or small number of individuals, independent of disease status. However, there were also clusters showing overall high or low expression in the majority of diseased valves and a small number where expression was high or low only in CKCS. The largest component of the network is shown in Figure 2A and clusters that demonstrate differences between diseased and normal valves are shown in Figure 2B. Histograms indicate the average expression of genes in the clusters.

We searched the GCN for proliferation genes and found that there was no association of cell division markers [31] with disease or breed status, indicating that perturbed patterns of cell division are not a consistent feature of MMVD in either CKCS or other breeds. In addition, mitochondrial genes [32, 33] were spread amongst a number of clusters and showed no relationship to breed or disease status suggesting that mitochondrial dysfunction is also not a feature of MMVD. Since immune response had been identified as altered in previous studies [20, 21] we also looked for clusters associated with changes in gene expression of immune/inflammatory genes [34]. A small number of macrophage and T-cell genes were in cluster011 (higher in all diseased valves; Figure 2B) but most of the genes associated with immune function were found in different clusters, suggesting that there is no co-ordinated immune gene expression in the diseased valves.

The genes within the disease associated clusters shown in Figure 2B were subject to enrichment analysis using DAVID (see Methods). Cluster002, which was a large group of genes apparently down-regulated in CKCS was analysed separately, while the genes in other down-regulated clusters and the up-regulated clusters were each pooled because they were close to each other in the network. This ensured that the groups were of sufficient size for meaningful analysis. The gene lists for the three analyses are given in Additional file 1, Table S1 which also shows the histograms of average gene expression for each cluster.

For Cluster002, which contained genes that were lower in CKCS valves than those from other diseased dogs or normal dogs, there was enrichment of Biological Process (BP) terms related to muscle structure and activity, cardiac conduction and calcium ion release. For Cellular Component (CC), terms related to sarcolemma were enriched. The top ten GO terms are presented in Additional file 1, Table S2. For the group of clusters in which genes were generally down-regulated in all diseased valves, DAVID GO enrichment analysis found enrichment for extracellular matrix and cell surface terms. For the clusters which contained genes that were up-regulated in diseased valves, there was enrichment of GO terms related to inflammatory response, monocyte chemotaxis and extracellular exosome.

Comparison of CKCS diseased valve transcriptome with normal valve transcriptome

The Graphia Pro analysis identified groups of genes that appeared to be up- or down-regulated in diseased valves, and a cluster of genes that were down-regulated only in CKCS valves. The analysis revealed that there were many genes showing idiosyncratic sample-specific expression patterns, which may have concealed some meaningful differences and resulted in the paucity of significant GO term enhancement. We therefore generated lists of differentially expressed genes (DEG) with the Affymetrix Transcription Analysis Console, using a one-way between-subject ANOVA (unpaired), which takes into account variance among samples in calculating the p-value and determines a false discovery rate (FDR) to allow for multiple comparisons.

To confirm the findings of the GCN analysis, an initial comparison was made between CKCS diseased valves and normal valves. Using a fold change of at least 1.5 in either direction and an FDR q-value < 0.05, transcripts detected by 755 probesets were differentially expressed, representing 599 annotated genes and a number of unannotated probesets. 271 annotated genes were higher in CKCS (up-regulated) and 328 genes were lower in CKCS (down-regulated) (Figure 3A, full list of genes in Additional file 2, Table S3). These 599 genes were analysed using DAVID GO enrichment analysis. For the down-regulated genes, there was enrichment of similar GO terms to those found for Cluster002 of the GCN analysis, as shown for the top ten enriched GO terms in Additional file 2, Table S4. The enriched GO terms were related to cardiac muscle cell function and structure and calcium channel activity. For up-regulated genes terms related to immune response and ERK1/ERK2 activity were listed.

Comparison of CKCS diseased valve transcriptome with transcriptome of diseased valves from other breeds

The sample-to-sample analysis suggested that gene expression in CKCS diseased valves was different from that in diseased valves from other breeds. Therefore we next compared the two sets of diseased valves. 161 annotated genes were differentially expressed, 27 with higher expression in CKCS valves than other diseased valves and 134 with lower expression in CKCS valves (Figure 3B, Additional file 3, Table S5). For the genes that were lower in CKCS diseased valves than other breed diseased valves, the GO terms were similar to those distinguishing CKCS from normal valves, with an emphasis on cardiac muscle structure and function (Additional file 3, Table S6). A single GO term was enriched for the genes that were higher in CKCS valves than other breed diseased valves.

Eighty-one genes were downregulated in CKCS valves compared both to normal valves and to diseased valves from other breeds, suggesting that these represent a breed-specific effect. Similar GO terms were enriched in this set of down-regulated genes, primarily terms associated with cardiac development and function (Additional file 4, Table S7 and Table S8). This may indicate that CKCS have abnormalities of expression of genes involved in cardiac development and function which may be related to the early onset of MMVD. Nineteen genes were up-regulated in CKCS compared with both normal valves and diseased valves from other breeds. The genes are shown in Additional file 4, Table S7. A single GO term was enriched, as for the comparison between CKCS and other breed diseased valves (Additional file 4, Table S8).

Comparison of diseased valve transcriptome with normal valve transcriptome

Since the GCN analysis showed that there were a number of genes where the majority of diseased valves were different from the normal valves, we generated a list of DEGs comparing all diseased valves with normal valves (using the same stringency) (Additional file 5, Table S9). 116 genes were differentially expressed, 50 with lower expression and 66 with higher expression in diseased valves than normal valves. Enrichment of GO terms for the down- and up-regulated genes is shown in Additional file 5, Table S10. For the up-regulated genes, terms related to skeletal system and mesenchyme migration were found, consistent with the increased EndoMT in the diseased valve. For the downregulated genes, the term calcium ion binding was enriched, supporting the idea that calcium homeostasis is perturbed in MMVD.

The Graphia Pro sample-to-sample network showed that the non-CKCS diseased valves were close in gene expression pattern to the normal valves, and there were no significant DEGs when comparing these two groups at the stringency used for the other comparisons (FDR q-value < 0.05). Consistent with this, the volcano plot for this comparison (Figure 3C) showed that the fold changes were lower, the p-values were higher and the genes were less scattered than in the other plots, supporting the observation that the normal diseased and non-CKCS diseased valves were closer in gene expression pattern than either were to the CKCS valves. Relaxing the stringency of the analysis (unadjusted p < 0.05, no FDR correction) showed that there were 278 differentially expressed genes, of which 166 were higher in the diseased valves and 112 were lower. For this lower stringency set, there were several enriched GO terms that overlapped with the set produced comparing CKCS and normal valves .

Pathway analysis of differentially expressed genes

IPA identified 77 canonical pathways when comparing the transcriptome of normal valves to CKCS valves, 28 when comparing CKCS valves to non-CKCS diseased valves and 56 when comparing all diseased valves with normal valves. The top three pathways for each comparison are shown in Table 3 and the top four upstream regulators for each dataset comparison are shown in Table 4. Of particular note are changes in calcium signaling in the CKCS compared to both normal and non-CKCSs datasets (Tables 4 and 5), whereas hepatic fibrosis/hepatic stellate cell activation was changed in CKCS and all diseased valves compared to normal. This association was also found in the lower stringency analysis of normal valves with non-CKCS diseased valves. For calcium signaling there were 23 DEGs (17 lower in CKCS) comparing CKCS with normal, 10 DEGs (all lower in CKCS) comparing CKCS and non-CKCS and 5 DEGs (1 lower in diseased valves) comparing all diseased valves to normal. Shared DEGs included genes associated with calcium homeostasis, non-canonical TGFβ signalling pathways (ERK1/2, IP3, RhoGTPase), cytoskeleton, muscle contraction (including cardiomyocytes), and the calcineurin/NFAT (nuclear factor of activated T cells) pathway. In the hepatic fibrosis/hepatic stellate cell activation gene dataset DEGs included genes associated with collagen homeostasis, cytoskeleton and cell growth and differentiation.

In general, comparing CKCS to normal, canonical pathways tended to be activated, but the pathways tended to be inhibited in CKCS relative to non-CKCS. For CKCSs compared to normal there were 1831 molecules with F2 (prothrombin) signalling pathway having the strongest association, but also TNF and TGFβ1 signalling featured in the top four upstream regulators. Comparing the two disease datasets, 377 molecules were associated, with the top upstream regulator being myocyte-specific factor 2c (MEF2C). The downstream effects, as predicted by IPA, for TGFβ1, MEF2C and F2 (prothrombin) are shown in Additional file 6, Figure S1 A-C.

The top five disease and function annotations assigned to the datasets by IPA and an illustrative graphical representation are shown in the Additional file 7, Table S11 and Figure S2. Comparing the datasets, these annotations generally matched to the same general themes, but identifying the subtle differences required examination of the graphic representations, with the illustrative example for Skeletal and Muscular Disorders and Developmental Disorder and Hereditary Disorder, comparing CKCS and non-CKCS (Additional file 7, Figure S2).

We have recently shown that the transcriptome of the mitral valve changes as MMVD advances from the normal valve to the severely diseased Whitney grade 4 valve (Markby et al, in press; doi: 10.3389/fgene.2020.00372). We found that TGFβ was the dominant signalling pathway controlling pathogenesis, consistent with findings in cultured valve interstitial cells [23]. The current study has identified gene expression differences in mitral valves of CKCS with severe valve pathology compared to valves from age and disease severity-matched non-CKCS dogs. These differences may be associated with the earlier onset of severe MMVD in the CKCS or may be an unrelated feature of the breed. CKCS and non-CKCS have been shown to have closely matched valve pathology [24], and the gene expression differences that are found in all diseased valves compared to normal valves may account for this. However, there are several limitations to the study. Firstly, all diseased dogs were older than the normal dogs and differences could simply be due to an age effect unrelated to MMVD status. Secondly, the other breed diseased valves and the normal valves came from several breeds of dogs, including some mixed breed animals, and the differences from CKCS valves could reflect the greater genetic similarity of the CKCS in contrast to the greater admixture of the other dogs [35]. Examining gene expression in normal and diseased age-matched CKCS valves would be ideal, but obtaining such valves is not feasible due to the ubiquity of the disease in all elderly dogs. Gene expression could also be examined in age matched valves of another MMVD-predisposed breed with a similar level of genetic diversity. However the early onset and high prevalence of MMVD in CKCS means that it would be difficult to find a matched breed. Similarly, comparison of normal valves from young CKCS with diseased valves from older CKCS would be interesting, but young CKCS with healthy valves are rarely euthanased. Nonetheless, it is possible to draw some conclusions from our study.

Among the genes that had lower expression in CKCS valves than non-CKCS diseased valves or normal dog valves were a range of cardiomyocyte-related genes including CASQ2 (calsequestrin), TNNI3 (troponin I type3) and RYR2 (ryanodine receptor 2), and various myosin genes. A striking feature of the GO term analysis was the reduced expression of genes involved in cardiac development and function in the CKCS valves. There may be a more general effect on the cardiovascular system of the CKCS. There were clear similarities in gene enrichment analysis and GO terms comparing both diseased groups to normal, including genes recognised as hallmarks of canine MMVD, such as ACTA2 (encoding α-SMA) and HTR2B (encoding 5HT2B receptor) [20, 27]. This likely reflects the global gene signature typical of MMVD where there is aberrant extra-cellular matrix remodelling as a consequence of changes in TGFβ signalling (Markby et al, in press; doi: 10.3389/fgene.2020.00372) [20, 21, 36]. The other shared DEGs, in particular those coding for the myosin heavy chains and growth factors, suggest the involvement of TGFβ non-canonical signalling pathways affecting down-stream signals including MEK, ERK1/2, IP3 and RhoGTPase [37]. All these pathways have regulatory roles in stress fibre formation, EndoMT and apoptosis, and can be induced by the TGFβs, 5HT and Ca⁺⁺signalling. TGFβ1 was one of the top upstream regulators in the CKCS valves, which would suggest TGFβ1 is an important driver of MMVD pathogenesis in dogs.

In the previously published study comparing diseased CKCS valves with normal valves using the Affymetrix Canine Gene 1.0ST Array, similar changes in gene ontology and gene networks were found, including regulation of EndoMT and caveolar-mediated endocytosis, but calcium signalling as a dysregulated canonical pathway was not reported [20]. In that study IPA identified only 33 canonical pathways compared to 77 in the current study.

The identification of the hepatic fibrosis/hepatic stellate cell activation pathway in diseased valves fits with what is already known about MMVD and the role of activated myofibroblasts in ECM remodelling and disease pathogenesis [24, 38-40]. Of interest was the identification of GO terms and canonical pathways and upstream regulators exclusive to the CKCS group. These included the canonical pathway calcium signalling and upstream regulator F2 (prothrombin), and various cardiovascular and immune-related GO terms. The potential role of immunity in MMVD might be difficult to explain as there is currently no evidence of inflammatory cell involvement in canine MMVD [24, 38]. There is however, evidence for inflammatory changes in the human MMVD and the CKCS breed does appear to be pre-disposed to a range of inflammatory conditions [22, 24, 38, 41, 42]. Network analysis using IPA has identified genes associated with “inflammation”, but this likely refers to pathways shared with other functions such as cell signalling, cell migration and EndoMT, since many of the same genes can be associated with immune responses [20, 21, 27]. Furthermore, expression of interleukins, chemokines and growth factors can be affected by the same up-steam regulators that control TGFβ signalling pathways.

F2 (prothrombin) was identified as the top upstream regulator in diseased CKCS valves, compared to normal valves. While typically associated with blood coagulation, prothrombin/thrombin also has a pro-inflammatory role and down-stream regulatory effects on endothelial, smooth muscle and fibroblast (interstitial) cells, and can both induce and inhibit TGFβ1 signalling in a context specific manner [43]. However, the reported CKCS predisposition to platelet dysfunction and macro-thrombocytopenia is well recognised, and macro-thrombocytopenia has been associated with early MMVD development in humans [8, 11, 13]. Considering these effects, and what is known about MMVD pathology, a potential role for prothrombin/thrombin in CKCS MMVD pathogenesis can be postulated, possibly through platelet interaction and enhanced 5HT signalling [16, 44, 45].

Changes in expression of cardiomyocyte-related genes were found in the CKCS dataset, including MB (myoglobin), CASQ2 (calsequestrin2), NEBL (nebulette) and various actin and myosin genes. While the changes in calcium signalling and MEF2C expression might reflect global changes in cell signalling and ECM remodelling, they might also be due to changes in the cardiomyocyte content of the tissue samples. Care was taken to ensure the same dissection protocol, removing atrial myocardium and the annular attachment, was used for each valve. Nevertheless, in whole valve samples atrial cardiomyocytes can be found extending a variable distance into the valve, and this extent declines in diseased leaflets [46-49]. Morphological studies are needed to determine if there are differences in CKCS valve myocardium content compared to non-CKCS that may explain some of these transcriptomic differences. Indeed, there is emerging interest in valve contraction mechanics and how this is affected by or contributes to the diseased state, and in this context the relative contribution of valve cardiomyocytes needs to be considered [50].

The myocardial gene MEF2C (myocyte-specific factor 2c) was the top upstream regulator when CKCS were compared to non-CKCS diseased valves. MEF2C is important in myogenesis and cardiomyocyte development, regulating cardiac alpha-actin and alpha myosin heavy chain [51]. However, as a transcription factor, it also regulates the expression of the ECM protein cartilage link protein 1 (HAPLN1 gene) in VICs during valve development and has a wide range of other effects, including in the EGF/EGFR and the apelin signalling pathways, that contribute to non-canonical TGFβ1 signalling, proper development of megakaryocytes and platelets, and cell migration and differentiation [51].

Calcium signalling was identified as the top canonical pathway in the CKCS dataset, and transcription factors in the top upstream regulator list and GO terms related to cardiac muscle contraction, control and differentiation, were down-regulated in CKCS diseased valves. Understanding how down-regulation of cardiac muscle, smooth muscle-related and calcium binding pathways impact on CKCS MMVD may provide insight into the early-onset pathogenesis in CKCS. For example, 5HT signalling through Gq protein-coupled receptors (such as 5HTR2B) causes an influx of calcium into the cell as part of 5HT downstream signalling pathways, and there is consistent and significant increased 5HTR2B expression in the canine mitral valve transcriptome [16, 20, 21, 27]. In the light of the data presented here, the recent identification of the beneficial effect of the calcium sensitizer and positive inotrope pimobendan in protecting against progression to heart failure is an interesting finding [52].

The notable limitations of this study include the small sample size, the age disparity between normal and diseased animals and the genetic dissimilarity of the CKCS to the other dogs. Difficulties in age-matching are an accepted confounding factor for MMVD research considering the ubiquity of valve pathology changes in all elderly dogs, some of which may be associated with healthy aging. CKCS have a distinct genetic profile [35] and this may account, in part, for the disparity in DEG numbers comparing CKCS with the two other groups. Studies are needed to determine whether this differential gene expression underlies the susceptibility of the CKCS to develop the disease sooner and to progress more rapidly than other dogs, bearing in mind that they might not have life-long differential expression [25]. Lastly, the more stringent FDR analysis could not be applied to the normal and non-CKCS analysis. This likely reflects their genetic heterogeneity and difference in baseline expression of genes, and in that circumstance is an accepted limitation of genomic profiling, with similar problems found when examining human populations [22, 53]. Relaxing the strict FDR threshold revealed candidate genes and pathways that were enriched in this comparison, consistent with the observations of the higher stringency comparisons and supporting the need for further analysis of these candidates.

This study has identified interesting differences in the transcriptomic profile of CKCS mitral valves compared with non-CKCS valves, all with severe MMVD. Studies are needed to determine if these gene expression changes simply reflect differences in valve cardiomyocyte density, which might then affect valve mechanics, valve coaptation and pathology development. Alternatively, the gene expression changes may impact in a CKCS-specific manner on valve cell function and ECM synthesis and remodelling. These results provide the basis for further studies to examine the specific gene and signalling pathway changes and their contribution to disease pathogenesis, and any breed-specific susceptibility.

Aims, design and setting of the study

The aim of this study was to determine whether expression patterns were different in CKCS diseased mitral valves compared to non-CKCS dogs. Valve samples were collected from dogs were presented to the Hospital for Small Animals, Royal (Dick) School of Veterinary Studies, the University of Edinburgh. Euthanasia had been requested by the owners because of terminal conditions or intractable conditions which impacted severely on their pets welfare. No dogs were euthanased for the purpose of the study and all tissue was collected and used with full informed owner consent, and with institutional ethical approval (Veterinary Ethics in Research Committee).

Tissue samples

Details of the dogs used in this study are shown in Table 1. All valves were scored by two of the authors (GRM & BMC). Diseased valves had changes typical of Whitney grade 3 or 4 (severe to very severe disease) [2]. Normal dogs had no evidence of any disease, and mitral valves were considered to be normal on gross inspection. Valves were collected shortly after euthanasia as previously described (intra-venous pentobarbitone overdose) [20], washed gently in warm phosphate buffered saline (PBS), immediately placed in RNAlater (Invitrogen USA) and stored for future RNA extraction [20]. Any attached atrial myocardium, annulus and chordae were removed, and the remaining whole valve was used for RNA extraction. Dogs with advanced disease were on a combination of standard medication for the treatment of congestive heart failure and any effect this might have on valve gene expression is unknown.

RNA extraction, quantification, quality control and transcriptomic profiling

RNA was extracted from the whole valve (anterior and posterior leaflet) following a standard protocol as previously described [20]. RNA extraction and DNA digestion were performed using the Qiagen RNeasy mini kit (Qiagen, Germany) according to the manufacturer’s instructions. Total RNA was eluted in 30μl nuclease free water and stored at -70°C. Quantification of RNA was performed by spectrophotometry in a NanoDrop™1000 (Thermo Scientific), measuring absorbance at 260nm wavelength in 1µl of extracted RNA solution, and absorbance ratios of 260/280 and 260/230 were analysed to check for impurities. The ratio of ribosomal 28S to 18S RNA was measured to assess for degradation using the Agilent RNA Screentape system and Agilent 2200 tapestation analyser (Agilent Technologies, USA) according to the manufacturer’s instructions. RNA integrity number (RIN) was then calculated with RIN ≥7 being taken as optimal for transcriptomic analysis. The Affymetrix GeneChip™ Canine Gene 1.1 ST Array plate was used for transcriptomic profiling. Arrays were run by Edinburgh Genomics, University of Edinburgh, UK. The Affymetrix Expression Console (Build 1.4.1.46) was used to normalise the data from the generated .cel files and for quality control. An annotated file of expression results was generated and filtered for all genes with a maximum relative expression level of less than 15, consistent with the background level detected during the quality control process. This filtered set of genes was used for subsequent expression and enhancement analysis.

Network analysis and functional clustering of canine valve samples

The network analysis tool Graphia Pro 1.4 (formerly BioLayout Express^3D; https://kajeka.com) [54] was used to examine expression of genes across the grades of valve disease. Graphia Pro 1.4 clusters data based on similarity of gene expression pattern with nodes representing a data point and edges the relation between nodes. In a sample-to-sample analysis (similar to a principal components analysis) nodes represent samples and the network layout shows the similarity of samples based on the expression of all genes in the sample. For the sample-to-sample comparison of mitral valve transcriptomic data using the filtered set of genes (relative expression ≥ 15 in at least one sample), a Pearson correlation coefficient of r ≥ 0.98 was used as it was the highest that included all samples. Meta-data (grade of disease, age, sex, breed) associated with the dataset were then examined to see if any parameter accounted for the network layout.

Gene-to-gene analysis with Graphia Pro was used to generate a gene co-expression network (GCN) where nodes represent genes and edges the correlation between them at or above the chosen threshold. The network layout shows the similarity of gene expression patterns across all samples. For the network layout a high stringency threshold r value of 0.9 was used. Subsequent analysis using the Markov clustering algorithm (MCL) [55] identified groups of highly connected genes within the elements of the network. The inflation value was set at 1.7 to control granularity of the clusters.

Differentially expressed genes

The Affymetrix transcriptome analysis console (TAC, version 3.1.0.5) was used to perform unpaired one-way analysis of variance and detect differentially expressed genes (DEG). DEG lists were created for genes with a fold change of >1.5 or <-1.5. A Benjamini-Hochberg false discovery rate (FDR) correction (Q-value <0.05) was applied. Annotation information for un-annotated transcript probes was found in TAC through an interface with the Affymetrix online browser. Gene lists were then used for gene enrichment analysis.

Reverse transcriptase quantitative PCR (RT-qPCR) was undertaken using the Takyon 2X low Rox SYBR green mastermix dTTP blue (Eurogentec, Belgium) to validate the microarray data and included the following genes; ACTA2, HTR2B, TAGLN, ACTG2, SLIT3, CDKN2A, SLC10A6, CILP, MMP12, ADAMTS5 and ADAMTS19. Primer sequences are shown in Table 2. RT-qPCR showed equivalent direction and magnitude to the results with the microarray for these genes.

Gene enrichment analysis

Gene lists from the network analysis and the differentially expressed gene lists were analysed with the Database for Annotation, Visualisation and Integrated Discovery (DAVID v6.8; http://www.david.ncifcrf.gov) [56, 57]. DAVID collates the biological processes that are associated with genes in a list using the gene ontology (GO) terms. Gene lists deriving from differential expression analysis or Graphia Pro 1.4 clustering were uploaded to DAVID for analysis. GO terms were selected for biological processes (GOTERM_BP_DIRECT), cellular components (GOTERM_CC_DIRECT) and molecular function (GOTERM_MF_DIRECT). The reference list was from Canis lupus familiaris and enrichment was assessed using Fisher’s exact test.

In addition, differentially expressed gene lists were uploaded, with related fold change, to the online Ingenuity Pathway Analysis (IPA) server (Qiagen, Germany) where core analysis was performed. This compared the submitted genes to all published literature and assigned reported attributes and pathways to them. From these canonical pathways, upstream regulators, disease and biological functions, as well as other factors, were inferred by the software. These results were generated based on the number of genes in the submitted list that matched the genes reported in the literature to be involved in a certain process, with a statistical association calculated. Core analysis was performed with default settings on filtered gene sets.

FDR, false discovery rate

IPA, Ingenuity Pathway Analysis

MMVD, myxomatous mitral valve disease

Ethical approval and consent to participate

The study was performed with ethical approval of the Veterinary Ethics Research Committee, University of Edinburgh and tissues were collected with full consent of the owners.

Consent for publication

Not applicable

Availability of data and materials

The datasets generated during the current study are available in the University of Edinburgh Datashare repository, https://doi.org/10.7488/ds/2754.

Competing interests

The authors declare they have no competing interests.

Funding

GRM was funded by the Dogs Trust. VEM is supported by funding from the Biotechnology and Biological Sciences Research Council (BBSRC) in the form of an Institute Strategic Programme Grant (BB/J004316/1). KMS is supported by the Mater Foundation, Brisbane, Australia. The Translational Research Institute receives funding from the Australian Government.

Authors’ contributions

BMC and KMS conceived the project and secured the funding. GRM and BMC removed and scored the mitral valves from all dogs. GRM extracted and prepared the RNA for microarray analysis. GRM and KMS performed the bioinformatic analyses. KMS, VEM and BMC supervised the project. GRM wrote the first draft of the manuscript. KMS, BMC and VEM edited and modified the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We wish to thank the owners for allowing us to collect valve tissue from their pets at an emotionally difficult time, and without whose support this study would not have been possible.

Additional file 1: .pdf: Gene lists (Table S1) and GO term enrichment analysis (Table S2) for clusters detected by GCN analysis with Graphia Pro.

Additional file 2: .pdf: Lists of differentially expressed genes (Table S3) and GO term enrichment analysis (Table S4) for comparisons between CKCS and normal valves.

Additional file 3: .pdf: Lists of differentially expressed genes (Table S5) and GO term enrichment analysis (Table S6) for comparisons between CKCS and other breed diseased valves.

Additional file 4: .pdf: List of genes (Table S7) that were down-regulated in CKCS compared with both other breed diseased valves and normal valves and GO enrichment analysis (Table S8).

Additional file 5: .pdf: Lists of differentially expressed genes (Table S9) and GO term enrichment analysis (Table S10) for comparisons between all diseased valves and normal valves.

Additional file 6: .pdf: Graphical representations of predicted down-stream regulator effects identified by IPA (Figure S1).

Additional file 7: .pdf: Disease and function networks identified by IPA (Table S11) and illustrative graphical representation (Figure S2).

Table 1

Metadata for valve samples analysed. Whitney gross pathology grade (0 normal to 4 very severe) was assigned independently by two of the authors (GRM and BMC). There was no significant difference in age between the CKCS and non-CKCS group, but there was for both compared to normal group (P<0.001)

Breed	Gender	Age	Whitney Grade
Normal
Cross-terrier	Male	2yrs	0
Cross-terrier	Male	3yrs	0
Beagle	Male	3yrs	0
Cross-Staffordshire bull terrier	Male	3yrs	0
Beagle	Female	4yrs	0
Beagle	Male	4yrs	0
Mean Age +/-S.E.	3.1yrs +/-0.35
Non-CKCS
Cross-English bull terrier	Female	10yrs	3
West Highland white terrier	Male	10yrs	3
Jack Russel terrier	Female	11yrs	4
Border collie	Male	13yrs	4
Border collie	Male	13yrs	4
Mean Age +/-S.E.	11.4yrs +/-0.27
CKCS
	Male	12yrs	3
	Male	11yrs	3
	Male	12yrs	3
	Female	10yrs	3
	Male	16yrs	4
	Female	12yrs	4
Mean Age +/-S.E.	12.5yrs+/-0.33

Table 2

Primer sequences for selected genes used in RT-qPCR to validate the microarray data.

Gene Symbol	Forward Primer Sequence	Reverse Primer Sequence
ACTA2	5'CGGCTACTCCTTTGTGACG3'	5'CGTGGCCATCTCGTTCTC3'
HTR2B	5'CCAATCCAGGCCAATCAAAG3'	5'CAGGTGATGTTGCTTGGGTT3'
TAGLN	5'GACATGTTCCAGACCGTCGA3'	5'CAATGACGTGCTTTCCCTCC3'
ACTG2	5'TGCCAACAATGTCCTTTCCG3'	5'GCCTCCAATCCAGACTGAGT3'
SLIT3	5'CTGACAAGGACAACGGCATC3'	5'CCCATCATTCACCGTCTCCA3'
CDKN2A	5'CATGTTGGCTCAGAATCGGG3'	5'CTCACGTCCAAGGCACAAAA3'
SLC10A6	5'GCTGTTGGATGGGTTTCTCA3'	5'TCCAAGAAAGCACCAGTCTCT3'
CILP	5'TGCTCCAATTATACCGTGCG3'	5'CAGAACACTTGCTCCAGGGA3'
MMP12	5'GACACAATTCATGGACCCTGG3'	5'TCAAATACGTCAGGTCCTTGGA3'
ADAMTS5	5'GTTCCCAAATATGCAGGCGT3'	5'AGCTTCGAACCAATGATGCC3'
ADAMTS19	5' GGACGGTGAGGTGTACTAAC 3'	5'ACTGCATTCCTTTACCACAGG 3'

Table 3

Top three canonical pathways associated with each dataset. The number of genes altered in each pathway as well as the total number of genes changed in each pathway is shown. Of note is change in calcium signalling comparing CKCS to the other two data sets highlighted in bold. The P-value score shows the strength of association of the gene list to the pathway.

Analysis	Canonical Pathway	Up	Down	Gene changes in pathway	p-value
CKCS vs Normal	Calcium signaling	6	17	23/179	7.943E-11
	Hepatic fibrosis/Hepatic stellate cell activation	8	12	20/183	2.988E-08
	Axonal guidance signaling	11	16	27/452	1.995E-05
CKCS vs non-CKCS	Calcium signaling	0	10	10/179	6.456E-06
	LPS/IL-1 mediated inhibition of RXR function	2	5	7/168	0.0009
	Gluconeogenesis I	0	3	3/22	0.001
All diseased vs normal	Paxillin signaling	3	1	4/108	0.0006
	Calcium signaling	4	1	5/179	0.0009
	STAT3 pathway	4	0	4/135	0.001

Table 4

The top four upstream regulators associated with the differentially expressed genes lists for each dataset. For each upstream regulator, the molecule type, Z-score and P-value are given. The activation Z-score is used to infer likely activation states of upstream regulators based on comparison with a model that assigns random regulation directions. https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/

Analysis	Upstream regulator	Molecule type	Activation Z-score	P-value
CKCS vs Normal	F2	Peptidase	1.503	2.2E-12
	TNF	Cytokine	1.892	2.18E-11
	AGT	Growth factor	2.018	1.15E-10
	TGFB1	Growth factor	1.486	9.03E-10
CKCS vs non-CKCS	MEF2C	Transcription regulator	-3.087	1.11E-09
	MYOCD	Transcription regulator	-2.768	1.33E-07
	2,3 butanedione monoxime	Chemical drug	-1.4	3.71E-07
	DNMT3A	Enzyme	1.667	6.27E-07
All diseased vs Normal	NOTCH4	Transcription regulator	1.777	6.33E-10
	MED28	Other	-1.957	6.32E-09
	HEY1	Transcription regulator	-1.547	2.61E-08
	MYCOD	Transcription regulator	2.571	4.6E-08

Table 5

Gene expression changes associated with calcium signalling and hepatic fibrosis/hepatic stellate cell activation canonical pathways. Down-regulated genes are shown in bold.

Datasets	Gene name
Calcium Signalling
CKCS vs Normal	*ACTA1, ACTA2, AKAP5, ATP2A2, CASQ, GRIN2A, GRIA3,* HDAC9, MYH7, MYH7B, MYH1, MYL4, NFATC1, RCAN2, RYR2, SLC8B1, TNNI3, TPM2, TRDN
CKCS vs non-CKCS	*ACTA1, ATP2A2, CACNA1G, CACNA1H, CASQ2, MYH7, RCAN2, SLC8B1, TNNI3, TRDN*
All diseased vs Normal	ACTA2, HDAC9, MYH11, TPM2, TRPC5
Hepatic fibrosis/hepatic stellate cell activation
CKCS vs Normal	ACTA2, CCL5, COL6A3, EDNRA, IL1RL1, KDR, LAMA1, MYH7, MYH11, MYH7B, MYL4, NGFR, PDGFR, SERPINE1,
All diseased vs Normal	ACTA2, COL11A2,* MYH11, SERPINE1*

Whitney JC: Cardiovascular pathology. J Small Anim Pract 1967, 8(8):459-465.
Whitney JC: Observations on the effect of age on the severity of heart valve lesions in the dog. J Small Anim Pract 1974, 15(8):511-522.
Beardow AW, Buchanan JW: Chronic mitral valve disease in cavalier King Charles spaniels: 95 cases (1987-1991). J Am Vet Med Assoc 1993, 203(7):1023-1029.
Darke PG: Valvular incompetence in cavalier King Charles spaniels. Vet Rec 1987, 120(15):365-366.
Lewis T, Swift S, Woolliams JA, Blott S: Heritability of premature mitral valve disease in Cavalier King Charles spaniels. Vet J 2011, 188(1):73-76.
Swift S: The problem of inherited diseases. 5: Valvular disease in Cavalier King Charles spaniels. J Small Anim Pract 1996, 37(10):505-506.
Meurs KM, Friedenberg SG, Williams B, Keene BW, Atkins CE, 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.
Cowan SM, Bartges JW, Gompf RE, Hayes JR, Moyers TD, Snider CC, Gerard DA, Craft RM, Muenchen RA, Carroll RC: Giant platelet disorder in the Cavalier King Charles Spaniel. Exp Hematol 2004, 32(4):344-350.
Cremer SE, Singletary GE, Olsen LH, Wallace K, Haggstrom J, Ljungvall I, Hoglund K, Reynolds CA, Pizzinat N, Oyama MA: Serotonin concentrations in platelets, plasma, mitral valve leaflet, and left ventricular myocardial tissue in dogs with myxomatous mitral valve disease. J Vet Intern Med 2014, 28(5):1534-1540.
Cremer SE, Zois NE, Moesgaard SG, Ravn N, Cirera S, Honge JL, Smerup MH, Hasenkam JM, Sloth E, Leifsson PS et al: Serotonin markers show altered transcription levels in an experimental pig model of mitral regurgitation. Vet J 2015, 203(2):192-198.
Olsen LH, Kristensen AT, Haggstrom J, Jensen AL, Klitgaard B, Hansson H, Pedersen HD: Increased platelet aggregation response in Cavalier King Charles Spaniels with mitral valve prolapse. J Vet Intern Med 2001, 15(3):209-216.
Tarnow I, Kristensen AT, Texel H, Olsen LH, Pedersen HD: Decreased platelet function in Cavalier King Charles Spaniels with mitral valve regurgitation. J Vet Intern Med 2003, 17(5):680-686.
Becker PS, Clavell LA, Beardsley DS: Giant platelets with abnormal surface glycoproteins: a new familial disorder associated with mitral valve insufficiency. J Pediatr Hematol Oncol 1998, 20(1):69-73.
Soga J, Yakuwa Y, Osaka M: Carcinoid syndrome: a statistical evaluation of 748 reported cases. J Exp Clin Cancer Res 1999, 18(2):133-141.
Connolly JM, Bakay MA, Fulmer JT, Gorman RC, Gorman JH, 3rd, Oyama MA, Levy RJ: Fenfluramine disrupts the mitral valve interstitial cell response to serotonin. Am J Pathol 2009, 175(3):988-997.
Driesbaugh KH, Branchetti E, Grau JB, Keeney SJ, Glass K, Oyama MA, Rioux N, Ayoub S, Sacks MS, Quackenbush J et al: Serotonin receptor 2B signaling with interstitial cell activation and leaflet remodeling in degenerative mitral regurgitation. J Mol Cell Cardiol 2018, 115:94-103.
Fitzgerald LW, Burn TC, Brown BS, Patterson JP, Corjay MH, Valentine PA, Sun JH, Link JR, Abbaszade I, Hollis JM et al: Possible role of valvular serotonin 5-HT(2B) receptors in the cardiopathy associated with fenfluramine. Mol Pharmacol 2000, 57(1):75-81.
Lacerda CM, Maclea HB, Kisiday JD, Orton EC: Static and cyclic tensile strain induce myxomatous effector proteins and serotonin in canine mitral valves. J Vet Cardiol 2012, 14(1):223-230.
Disatian S, Orton EC: Autocrine serotonin and transforming growth factor beta 1 signaling mediates spontaneous myxomatous mitral valve disease. J Heart Valve Dis 2009, 18(1):44-51.
Lu CC, Liu MM, Culshaw G, Clinton M, Argyle DJ, Corcoran BM: Gene network and canonical pathway analysis in canine myxomatous mitral valve disease: a microarray study. Vet J 2015, 204(1):23-31.
Oyama MA, Chittur SV: Genomic expression patterns of mitral valve tissues from dogs with degenerative mitral valve disease. Am J Vet Res 2006, 67(8):1307-1318.
Thalji NM, Hagler MA, Zhang H, Casaclang-Verzosa G, Nair AA, Suri RM, Miller JD: Nonbiased Molecular Screening Identifies Novel Molecular Regulators of Fibrogenic and Proliferative Signaling in Myxomatous Mitral Valve Disease. Circulation Cardiovascular genetics 2015, 8(3):516-528.
Tan K, Markby G, Muirhead R, Blake R, Bergeron L, Fici G, Summers K, Macrae V, Corcoran B: Evaluation of canine 2D cell cultures as models of myxomatous mitral valve degeneration. PLoS One 2019, 14(8):e0221126.
Lu CC, Liu MM, Culshaw G, French A, Corcoran B: Comparison of cellular changes in Cavalier King Charles spaniel and mixed breed dogs with myxomatous mitral valve disease. J Vet Cardiol 2016, 18(2):100-109.
French AT, Ogden R, Eland C, Hemani G, Pong-Wong R, Corcoran B, Summers KM: Genome-wide analysis of mitral valve disease in Cavalier King Charles Spaniels. Vet J 2012, 193(1):283-286.
Madsen MB, Olsen LH, Haggstrom J, Hoglund 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 Suppl 1:S62-67.
Lu CC, Liu MM, Clinton M, Culshaw G, Argyle DJ, Corcoran BM: Developmental pathways and endothelial to mesenchymal transition in canine myxomatous mitral valve disease. Vet J 2015.
Wolfe CJ, Kohane IS, Butte AJ: Systematic survey reveals general applicability of "guilt-by-association" within gene coexpression networks. BMC Bioinformatics 2005, 6:227.
Young R, Lefevre L, Bush SJ, Joshi A, Singh SH, Jadhav SK, Dhanikachalam V, Lisowski ZM, Iamartino D, Summers KM et al: A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis). Front Genet 2019, 10:668.
Freem L, Summers KM, Gheyas AA, Psifidi A, Boulton K, MacCallum A, Harne R, O'Dell J, Bush SJ, Hume DA: Analysis of the Progeny of Sibling Matings Reveals Regulatory Variation Impacting the Transcriptome of Immune Cells in Commercial Chickens. Front Genet 2019, 10:1032.
Giotti B, Chen SH, Barnett MW, Regan T, Ly T, Wiemann S, Hume DA, Freeman TC: Assembly of a Parts List of the Human Mitotic Cell Cycle Machinery. J Mol Cell Biol 2018.
Hume DA, Summers KM, Raza S, Baillie JK, Freeman TC: Functional clustering and lineage markers: insights into cellular differentiation and gene function from large-scale microarray studies of purified primary cell populations. Genomics 2010, 95(6):328-338.
Freeman TC, Ivens A, Baillie JK, Beraldi D, Barnett MW, Dorward D, Downing A, Fairbairn L, Kapetanovic R, Raza S et al: A gene expression atlas of the domestic pig. BMC Biol 2012, 10:90.
Nirmal AJ, Regan T, Shih BB, Hume DA, Sims AH, Freeman TC: Immune Cell Gene Signatures for Profiling the Microenvironment of Solid Tumors. Cancer Immunol Res 2018, 6(11):1388-1400.
Mellanby RJ, Ogden R, Clements DN, French AT, Gow AG, Powell R, Corcoran B, Schoeman JP, Summers KM: Population structure and genetic heterogeneity in popular dog breeds in the UK. Vet J 2013, 196(1):92-97.
Aupperle H, Marz I, Thielebein J, Schoon HA: Expression of transforming growth factor-beta1, -beta2 and -beta3 in normal and diseased canine mitral valves. J Comp Pathol 2008, 139(2-3):97-107.
Zhang YE: Mechanistic insight into contextual TGF-beta signaling. Curr Opin Cell Biol 2018, 51:1-7.
Disatian S, Ehrhart EJ, 3rd, Zimmerman S, Orton EC: Interstitial cells from dogs with naturally occurring myxomatous mitral valve disease undergo phenotype transformation. J Heart Valve Dis 2008, 17(4):402-411; discussion 412.
Han RI, Black A, Culshaw GJ, French AT, Else RW, Corcoran BM: Distribution of myofibroblasts, smooth muscle-like cells, macrophages, and mast cells in mitral valve leaflets of dogs with myxomatous mitral valve disease. Am J Vet Res 2008, 69(6):763-769.
Liu MM, Flanagan TC, Lu CC, French AT, Argyle DJ, Corcoran BM: Culture and characterisation of canine mitral valve interstitial and endothelial cells. Vet J 2015, 204(1):32-39.
Summers JF, O'Neill DG, Church DB, Thomson PC, McGreevy PD, Brodbelt DC: Prevalence of disorders recorded in Cavalier King Charles Spaniels attending primary-care veterinary practices in England. Canine Genet Epidemiol 2015, 2:4.
Hulin A, Anstine LJ, Kim AJ, Potter SJ, DeFalco T, Lincoln J, Yutzey KE: Macrophage Transitions in Heart Valve Development and Myxomatous Valve Disease. Arteriosclerosis, thrombosis, and vascular biology 2018, 38(3):636-644.
Bae JS, Kim IS, Rezaie AR: Thrombin down-regulates the TGF-beta-mediated synthesis of collagen and fibronectin by human proximal tubule epithelial cells through the EPCR-dependent activation of PAR-1. J Cell Physiol 2010, 225(1):233-239.
Hoglund K, Haggstrom J, Hanas S, Merveille AC, Gouni V, Wiberg M, Lundgren Willesen J, Entee KM, Mejer Sorensen L, Tiret L et al: Interbreed variation in serum serotonin (5-hydroxytryptamine) concentration in healthy dogs. J Vet Cardiol 2018, 20(4):244-253.
Oyama MA, Levy RJ: Insights into serotonin signaling mechanisms associated with canine degenerative mitral valve disease. J Vet Intern Med 2010, 24(1):27-36.
Culshaw GJ, French AT, Han RI, Black A, Pearson GT, Corcoran BM: Evaluation of innervation of the mitral valves and the effects of myxomatous degeneration in dogs. Am J Vet Res 2010, 71(2):194-202.
Buchanan JW: Chronic valvular disease (endocardiosis) in dogs. Adv Vet Sci Comp Med 1977, 21:75-106.
Fenoglio JJ, Jr., Tuan Duc P, Wit AL, Bassett AL, Wagner BM: Canine mitral complex. Ultrastructure and electromechanical properties. Circ Res 1972, 31(3):417-430.
Fox PR: Pathology of myxomatous mitral valve disease in the dog. J Vet Cardiol 2012, 14(1):103-126.
Easley TF, Bloodworth CHt, Bhal V, Yoganathan AP: Effects of annular contraction on anterior leaflet strain using an in vitro simulator with a dynamically contracting mitral annulus. J Biomech 2018, 66:51-56.
Lockhart MM, Wirrig EE, Phelps AL, Ghatnekar AV, Barth JL, Norris RA, Wessels A: Mef2c regulates transcription of the extracellular matrix protein cartilage link protein 1 in the developing murine heart. PLoS One 2013, 8(2):e57073.
Boswood A, Haggstrom J, Gordon SG, Wess G, Stepien RL, Oyama MA, Keene BW, Bonagura J, MacDonald KA, Patteson M et al: Effect of Pimobendan in Dogs with Preclinical Myxomatous Mitral Valve Disease and Cardiomegaly: The EPIC Study-A Randomized Clinical Trial. J Vet Intern Med 2016, 30(6):1765-1779.
Hulin A, Deroanne CF, Lambert CA, Dumont B, Castronovo V, Defraigne JO, Nusgens BV, Radermecker MA, Colige AC: Metallothionein-dependent up-regulation of TGF-beta2 participates in the remodelling of the myxomatous mitral valve. Cardiovasc Res 2012, 93(3):480-489.
Theocharidis A, van Dongen S, Enright AJ, Freeman TC: Network visualization and analysis of gene expression data using BioLayout Express(3D). Nat Protoc 2009, 4(10):1535-1550.
van Dongen S, Abreu-Goodger C: Using MCL to extract clusters from networks. In: Bacterial molecular networks: Mtehods and protocols. edn. Edited by Van Helden J, Toussaint A, Theiffry D. New York, NY, USA: Springer; 2012.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009, 4(1):44-57.
Huang da W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 2009, 37(1):1-13.

Download PDF

Journal Publication

published 23 Sep, 2020

Read the published version in BMC Veterinary Research →

Editorial decision: Major revision
30 May, 2020
Review #2 received at journal
27 May, 2020
Review #1 received at journal
25 May, 2020
Reviewer #2 agreed at journal
06 May, 2020
Reviewer #1 agreed at journal
05 May, 2020
Editor assigned by journal
23 Apr, 2020
Reviewers invited by journal
23 Apr, 2020
Submission checks completed at journal
22 Apr, 2020
Editor invited by journal
22 Apr, 2020
First submitted to journal
17 Apr, 2020

You are reading this older preprint version

Read the latest preprint version →

Comparative transcriptomic profiling of myxomatous mitral valve disease in the Cavalier King Charles Spaniel

Status:

Journal Publication

Version 1

Abstract

Background.

Results.

Conclusion.

Figures

Introduction

Results

Discussion

Conclusions

Materials And Methods

List of Abbreviations

Declarations

Supplementary Information

Tables

References

Supplementary Files

Status:

Journal Publication

Version 1