Prospective pilot study on the predictive significance of plasma miR-30b-5p through the study of echocardiographic modifications in Cavalier King Charles Spaniels affected by different stages of myxomatous mitral valve disease: The PRIME study

Abstract

Specific microRNAs expressions may accurately characterize different stages of canine myxomatous mitral valve disease. This preliminary pilot study aimed to (1) describe the clinical and echocardiographic parameters of Cavalier King Charles Spaniels affected by myxomatous mitral valve disease at different American College of Veterinary Internal Medicine (ACVIM) stages (B1, B2 and C) and healthy control group (ACVIM A), comparing the parameters collected during the first examination (T0) and the end of the follow-up (T1); (2) assess the association between the values of echocardiographic parameters at T1 and the expression profile of miR-30b-5p at T0. Thirty-five Cavalier King Charles Spaniels (median age 4.29 years and median weight 9 Kg) in different ACVIM stages were included (7 A, 19 B1, 6 B2 and 3 C). Inverse probability weighting analysis was performed to estimate the association of the exposure variable (miR-30b-5p) with the outcome variables (clinical and echocardiographic variables). Time was included as variable. The results pointed out that high levels of plasma miR-30b-5p corresponded to lower values of left ventricular end-diastolic diameter normalized for body weight, end-diastolic and end-systolic volumes indexed for body weight, and left atrium-to aortic root ratio. Hence, higher miR-30b-5p expressions were associated with milder forms of mitral valve disease in our study population. In contrast, the results obtained for the intensity of heart murmur, the mitral regurgitation severity, and the Mitral INsufficiency Echocardiographic score) were not statistically significant. A relationship between high abundance of miR-30b-5p and myxomatous mitral valve disease that appear echocardiographically more stable over time has been demonstrated. In conclusion, Cavalier King Charles Spaniels affected by myxomatous mitral valve disease that at the first cardiologic evaluation showed an upregulation of miR-30b-5p are expected to experience lesser variations on their echocardiographic examination between T0 and T1.

Introduction

Myxomatous mitral valve disease (MMVD) is the most common cardiovascular disease in dogs and a frequent cause of heart failure [1,2]. Despite its large diffusion, the underlying causes and mechanisms involved in the development of valve degeneration and disease progression remain to be conclusively elucidated [1,2]. Although a causative mutation has not been identified yet, canine MMVD is reported to be hereditable in some breeds, including the Cavalier King Charles Spaniel (CKCS). Approximately 100% of CKCSs develop MMVD at the age of 10 years, and a polygenic mode of inheritance has been suggested [3–5]. The application of advanced genomic technologies has potential to provide information that may impact treatment, prevalence, or severity of MMVD through the elucidation of pathogenic mechanisms and the detection of predisposing genetic loci of major effect [6].

Transthoracic echocardiography is currently considered the gold standard for the diagnosis of MMVD and the evaluation of the follow-up of the affected dogs [2]. However, this test needs specialized equipment and well-trained veterinarians to reduce interobserver variability, especially in the early stages of the disease as valves affected by mild changes may be misinterpreted as normal.

The identification of reliable biomarkers may enhance the robustness of clinical decision-making and may pave the way for the development of alternative/complementary diagnostic strategies [7]. In recent years, the study of microRNAs (miRNAs)as promising diagnostic and prognostic biomarkers for the detection of disease increased [8]. They are short non-coding RNA molecules (21–25 nucleotides in size) found in plants, animals and some viruses, that regulate post-transcriptional silencing of target genes involved in pathophysiologic conditions in many mammalian species [9–15], including cardiovascular diseases [16,17]. In human medicine, plasma and serum miRNAs are currently being intensively investigated and specific miRNA expression patterns have been reported for various cardiovascular diseases [17]. In contrast, in veterinary medicine, the list of studies on circulating miRNAs is currently brief and few of them concern canine MMVD [6,18–22]. Specifically, Hulanicka et al., in 2014 [18], analysed the expression of 9 specific miRNAs in the plasma of Dachshunds suffering from MMVD. Enrolled dogs were classified adapting the American College of Veterinary Internal Medicine (ACVIM) classification as stage A (i.e., unaffected control dogs; n. 8), B (asymptomatic dogs, considering together the ACVIM stages B1 [i.e., dogs without cardiac enlargement] and B2 [dogs with cardiac enlargement]; n. 8) and C (considering together the ACVIM stages C [i.e., dogs that have experienced an episode of congestive heart failure] and D [i.e., dogs affected by refractory congestive heart failure]; n. 8). Expression analysis using the Real-Time Quantitative Reverse Transcription PCR (RT-qPCR) method revealed that two of these miRNAs were significantly downregulated (i.e., the expression of miR-30b and miR-133b differed between stage B and A and between stage C and A, respectively) [18]. In 2015, Li and co-workers [20] investigated a larger number of miRNAs (n. 277) by RT-qPCR in dogs belonging to various breeds and affected by MMVD at different ACVIM stages (A = 6; B1+B2 = 6; C+D = 6). Stage B1+B2 and C+D showed 11 dysregulated miRNAs (four upregulated and seven downregulated) as compared to stage A. Moreover, the expression of six out of 11 miRNAs was significantly different between stages B1+B2 and C+D [20]. In 2017, Yang and colleagues [22] for the first time evaluated miRNAs from exosomes (i.e., extracellular vesicles of endosomal origin, enriched in cardioprotective miRNAs, which behave as key mediators of intercellular communication) using RT-qPCR in 13 healthy dogs and 14 affected by MMVD at various stages (asymptomatic = 7; symptomatic = 7). They demonstrated that exosomal miR-181 and miR-495 significantly increased in dogs with symptomatic MMVD compared to the others, whereas miR-9 and miR-599 increased and decreased, respectively, in all MMVD-affected dogs compared to healthy dogs. This led to the hypothesis that exosomal miRNA expression-level changes may be more specific to the disease status than total plasma miRNAs [22]. However, processing of exosomal miRNAs is not standardized yet, so a better consistency in the sample processing is still needed. For this reason, quantification of plasma miRNAs is preferable, at present [23]. In 2018, the expression of circulating miRNAs via genome-wide sequencing through the enrolment of 17 dogs of various breeds (8 dogs with congestive heart failure secondary to MMVD and 9 age-matched healthy dogs) was profiled by sequencing and qPCR validation [19]. The analysis showed the up-regulation of four (miR-133, miR-1, miR-let-7e, and miR-125) and the downregulation of four miRNAs (miR-30c, miR-128, miR-142, and miR-423) in MMVD-affected dogs [19]. In 2021, Ro et al. [21] used RT-qPCR to evaluate levels of 11 miRNAs in 82 dogs of various breeds, including 10 healthy dogs and 72 dogs with different heart diseases (including 35 affected by MMVD: stage B1 = 1; B2 = 7; C = 22; D = 5). In dogs affected by MMVD, upregulation of cfa-miR-130b was documented [21]. In 2022, Bagardi et al. [7] measured the concentration of miR-1-3p, miR-30b-5p and miR-128-3p in 33 CKCSs of ACVIM stage B1 divided into three age-related groups, and 11 ACVIM stage A CKCSs. The results of this study showed that only miR-30b-5p was significantly higher in ACVIM B1 dogs than in ACVIM A subjects. According to the age of dogs, the amount of miR-30b-5p was significantly higher in group B1<3y, B1 3-7y, and B1>7y than in group A. The areas under the receiver operating curves were fair in discriminating between group B1<3y and group A, and between B1 3-7y and A, and good in discriminating between group B1>7y and A [7]. This study demonstrates that the identification of dogs with early asymptomatic MMVD through the evaluation of miR-30b-5p represents an intriguing possibility that certainly merits further research.

Although these results could allow speculating that specific miRNAs expressions may accurately characterize different stages of MMVD, further research is necessary to confirm such a hypothesis, especially considering the remarkable limitations of the abovementioned studies (i.e., a limited sample size of each research, potential confounding effects associated with the heterogeneous breeds, age categories, MMVD stages as well as types of evaluated miRNAs, and lack of follow-up). Thus, the need arises to create studies that are statistically and methodologically more robust. To achieve this, the present work tries to assess the predictive value of one of the most promising of the miRNAs investigated to date [7,18], namely miR-30b-5p, by correlating its expression to the clinical and echocardiographic parameters of CKCSs affected by MMVD at different ACVIM stages. The choice to focus the attention exclusively on this miRNA has been dictated by the crossed findings reported in literature for the Dachshund by Hulanicka et al. [18] and for the CKCS by Bagardi et al. [7]. Both works, in fact, reported a significant downregulation of miR-30b-5p in MMVD-affected dogs belonging to ACVIM classes B (asymptomatic in general, with and without cardiac remodelling) and B1, respectively, when compared to healthy subjects.

The present study aims at (1) describing the clinical and echocardiographic parameters of healthy CKCSs (ACVIM A) and CKCSs affected by MMVD at different ACVIM stages (B1, B2 and C), comparing the obtained parameters at T0 (first examination) and T1 (end of the follow-up); and (2) assessing the association between the values of clinical and echocardiographic parameters at T1 and the expression profile of miR-30b-5p measured at T0.

Materials and methods

Clinical and echocardiographic examinations

This is a preliminary prospective pilot study. Owned CKCSs visited at the Cardiology Unit of the Department of Veterinary Medicine, University of Milan, between November 2019 and September 2021 were considered for enrollment. Informed consent was signed by the owners, according to the ethical committee statement of the University of Milan, number 2/2016.

Every dog was evaluated on two different occasions, a first examination (T0), and a follow-up clinic visit (T1). For this study, the time between T0 and T1 had not to be standardized but could vary among the recruited patients between 6 and 24 months.

At T0, animal history was collected and a complete physical examination was performed. The cardiovascular system was evaluated through cardiac auscultation performed by trained cardiologists (MB and PGB), who established the presence/absence of a cardiac murmur and its timing and intensity. Cardiac murmurs were classified as follows: 0 = absent; 1 = I-II/VI left apical systolic or soft; 2 = III-IV/VI bilateral systolic or moderate and loud, respectively; 3 = V-VI/VI bilateral systolic or palpable [24]. Blood pressure was indirectly measured through an oscillometric method according to the ACVIM consensus statement [25].

A standard transthoracic echocardiographic examination was performed by a single operator (MB) to detect MMVD using MyLabOmega Vet ultrasound machine (Esaote, Genova, Italy), equipped with multi-frequency phased array probes (1–5 and 2–9 MHz). Probes were chosen according to the weight of the subject. Echocardiographic examinations were performed with continuous electrocardiographic monitoring and video clips were obtained to make offline measurements [26]. Dogs were gently restrained during the examinations and were classified according to the ACVIM classification scheme [27].

The echocardiographic parameters considered pivotal to the MMVD diagnosis and monitoring over time and included in the study are described in Table 1.

[Figure omitted. See PDF.]

At the end of the echocardiographic examination, the Mitral INsufficiency Echocardiographic (MINE) score was calculated by assigning a score from 1 to 4 to the left ventricular end-diastolic diameter normalized for body weight (LVIDdN), left atrial-to-aortic root ratio (LA/Ao), E wave peak velocity (E) and fractional shortening (FS%) values. Every patient was classified on the base of the total reached score, as follows: mild (total score = 4–5), moderate (total score = 6–7), severe (total score = 8–12) and late stage (total score = 13–14) [32].

Peripheral venous blood sampling was performed at the end of the echocardiographic examination. Dogs were fasted for at least 12 hours. Blood was collected from the cephalic or saphenous vein in two 2.5-mL EDTA tubes and one serum collection tube. Blood from one EDTA tube was used for a complete blood count, while blood from the serum collection tube was used for biochemical analysis. The leftover was used for plasma miRNAs extraction. According to the study design, at T1, every procedure described for T0 was repeated, except for plasma miR-30b-5p extraction. Every variation of animal history, clinical signs, and echocardiographic parameters was carefully noted to monitor MMVD progression since T0.

CKCSs with other systemic diseases such as systemic hypertension, primary pulmonary hypertension, neoplasia, and other cardiac abnormalities such as dilated cardiomyopathy, congenital cardiac abnormalities, endocarditis, and severe arrhythmia were excluded from the study.

MiR-30b-5p extraction and RT-qPCR quantification

SmallRNA extraction was performed as previously reported for samples collected at T0 [7]. Briefly, within two hours from sampling, the leftover blood was centrifuged at 800 g for 15 minutes. Plasma was then stored at -80°C until RNA isolation.

An aliquot of 150 μL of plasma per sample was thawed on ice and centrifuged at 3000 × g for 5 min at 4 °C and smallRNA was extracted using miRNeasy Serum/Plasma Kits (Qiagen, Cat. No. 217184, Milano, Italy) following the manufacturer’s instructions. The sample was incubated at room temperature for 5 min and 3.75μL (25 fmol final concentration) of the synthetic spike-in control Caenorhabditis elegans miRNA cel-miR-39 (Qiagen, Cat. No. 219610) was spiked into samples. RNA extraction was then carried out according to the manufacturer’s instructions. RNA concentration was assayed using a NanoDrop ND-1000 spectrophotometer. Reverse transcription was carried out on 10 ng of total RNA using a TaqMan Advanced miRNA cDNA Synthesis Kit (Cat. No. A28007, Applied Biosystems) following the manufacturer’s instructions.

RT-qPCR was performed following the MIQE guidelines [33]. The small RNA TaqMan assays (ThermoFisher Scientific) included cel-miR-39-3p (assay ID 478293_mir) and miR-30b-5p (assay ID 478007_mir). The reference miRNA was miR-16-5p (assay ID rno481312_mir). Quantitation was performed in 15 μl in a CFX Connect Real-Time PCR Detection System (Bio-Rad) as previously reported [33]. No-RT controls and no-template controls were included. MiRNA expressions are reported in terms of fold change normalized to miR-16 using the formula 2^−ΔΔCq on Bio-Rad CFX Maestro Software.

Statistical analysis

Statistical analysis was performed using commercially available statistics software (SPSS^® 27.0, IBM^®, SPSS, USA). Descriptive statistics were generated. The distribution of data for continuous variables was assessed for normality using the Shapiro-Wilk test. Normally distributed parameters are expressed as mean ± standard deviation, while non-normally distributed parameters are reported as the median and interquartile range (25^th and 75^th percentile). Obtained values were defined as statistically significant at a P value < 0.05, and all the significance values were adjusted according to the Bonferroni post hoc correction.

The Inverse Probability Weighting (IPW) method was performed to face confounding biases using “R” software [34] to estimate the association of the exposure variable (miR-30b-5p) with the outcome variables (clinical and echocardiographic variables) considering the time. The regressors were miR-30b-5p, time, and their interaction. A gaussian model was used for continuous responses (LA/Ao, EDVI, ESVI, LVIDdN, FS% and E) and an ordinal model was used for ordinal responses (presence/intensity of a heart murmur, mitral regurgitant jet size, MINE score). For each outcome model and each response variable, the marginal effect of the exposure variable (miR-30b-5p) and its standard error were calculated. The marginal effect of miR-30b-5p for continuous responses was interpreted as the increase/decrease (depending on the positive or negative marginal effect value) of the response variable, corresponding to an increase of miR-30b-5p in one unit. The marginal effect of miR-30b-5p for ordinal responses was interpreted as the increase/decrease of the probability of the lowest-value category when miR-30p-5p increases by one unit [35].

Results

Demographics and characteristics of the study population

The study population was composed of 35 CKCS, of which 23 (65.7%) were females and 12 (34.3%) males. At T0, the median age was 4.29 years (2.07–7.66), and the body weight was 8.99 ± 2.08 kg. Clinical and echocardiographic characteristics and plasma miR-30b-5p expression of the population are summarized in Table 2.

[Figure omitted. See PDF.]

Data obtained at T0.

All the enrolled dogs had no abnormalities on complete blood count and biochemical analysis. All dogs were re-evaluated at a follow-up clinic visit (T1), although the time between T0 and T1 varied among the recruited patients as follows: between 6 and 12 months (14 dogs), between 12 and 18 months (19 dogs), and between 18 and 24 months (2 dogs). The ACVIM stage of some dogs evolved over time. Specifically, 5 of the 7 dogs that had been classified as ACVIM class A at T0 were classified as a class B1 at T1, all the B1 subjects remained in the same ACVIM class, one ACVIM class B2 dog at T0 was classified as a class C at T1, and one ACVIM class C dog at T0 was classified as a class D at T1. Similarly, also some demographic, clinical and echocardiographic variables changed from T0 to T1, as summarized in Table 3.

[Figure omitted. See PDF.]

Data obtained at T1.

IPW analysis for continuous variables

According to the IPW analysis, the obtained results were statistically significant for 4 of the selected echocardiographic parameters, namely LA/Ao, ESVI, EDVI, and LVIDdN. Specifically, the marginal effect of miR-30b-5p was negative for each of these parameters, which means that any increase in one unit of plasma miR-30b-5p expression corresponded to a decrease of the echocardiographic variable in a number denominated as “estimate”. In particular, the marginal effect of miR-30b-5p was as follows: LVIDdN (estimate: -0.0004, P = 0.0003), ESVI (estimate: -0.0124, P = 0.005), EDVI (estimate: -0.0503, P = 0.0005), LA/Ao (estimate: -0.0004, P = 0.048) (S1 Table).

Fig 1 shows the marginal effect of plasma miR-30b-5p expression on each echocardiographic variable at T0 and T1. For a better understanding of the figure, it should be considered that when T0 and T1 lines are less parallel, the tendency to an interaction effect between miR-30b-5p and the echocardiographic parameter is higher (i.e., ESVI).

[Figure omitted. See PDF.]

Black line: T0. Red line: T1. A) Marginal effect of miR-30b-5p on LVIDdN. B) Marginal effect of miR-30b-5p on ESVI: The interaction between the two parameters is high since T0 and T1 lines are poorly parallel. C) Marginal effect of miR-30b-5p on EDVI. D) Marginal effect of miR-30b-5p on LA/Ao.

Through the IPW analysis we also studied the marginal effect of the elapsed time between T0 and T1 on the echocardiographic variables. Interestingly, none of the results reached the statistical significance: LVIDdN (P = 0.194), ESVI (P = 0.534), EDVI (P = 0.121), LA/Ao (P = 0.602), FS% (P = 0.074), E (P = 0.785).

IPW analysis for ordinal variables

As previously explained, the IPW analysis tried to determine whether plasma miR-30b-5p expression has a marginal effect even on ordinal variables of our choice, specifically the presence/intensity of a heart murmur, the mitral regurgitant jet size, and the MINE score. None of the abovementioned parameters showed a statistically significant result when analyzed along the plasma miR-30b-5p expression; however, the marginal effect of the miRNA on these variables is always positive (S1 Table). Since these variables have been expressed in different classes to determine their severity, the positive marginal effect of miR-30b-5p examined in the IPW analysis means that any increase in one unit of plasma miR-30b-5p expression corresponds to an increase of the probability of the class with the lower value. Specifically, these classes are: 0 = absent, for the heart murmur; 0 = absent, for the mitral regurgitation; mild, for the MINE score. Fig 2 shows the marginal effect of plasma miR-30b-5p expression on each ordinal variable at T0 and T1.

[Figure omitted. See PDF.]

Black line: T0. Red line: T1. On the y-axis, the probability of the class with the lower value is reported. A) Marginal effect of miR-30b-5p on the presence/intensity of a heart murmur. B) Marginal effect of miR-30b-5p on the mitral regurgitant jet size. C) Marginal effect of miR-30b-5p on the MINE score.

Discussion

At present, the pathophysiological aspect of the MMVD is not fully elucidated, especially regarding the molecular background of the disease. The CKCS represent a perfect model for studies aimed at filling such a knowledge gap, considering the hereditability of the disease in this breed, the high prevalence, and the fact that in CKCS MMVD can develop at a relatively young age and show a fast progression [3–6,36–38]. Thus, CKCSs have been enrolled in the present study. In dogs, the diagnosis of MMVD is typically based on a combination of two-dimensional and Doppler echocardiographic findings [2]. However, this approach necessarily required specialized equipment and well-trained cardiologists, which may be not always available in all veterinary institutions. Moreover, since echocardiograms are often acquired only after the recognition of a heart murmur on physical examination, such an approach may lead to tardive diagnosis with potentially relevant clinical consequences, especially considering that in some dogs the disease progression may be particularly rapid. Investigation of reliable early biomarkers of MMVD, such miRNAs, could offer additional advantages in the recognition of the disease, with the aim of reaching prompt diagnosis and, consequently, trying to improve therapeutic timing and life expectancy.

A study published by Bagardi et al. [7] conducted on the CKCS recently provided interesting results for miR-30b-5p in ACVIM class B1 CKCSs, demonstrating that miR-30b-5p can discriminate between B1 subjects and healthy dogs, and B1 CKCSs younger than 3 years of age and healthy dogs [7]. Since a heart murmur is uncommon to detect in young patients and echocardiographic changes can be misclassified as normal findings by inexperienced sonographers, miR-30b-5p can support the early diagnosis of MMVD-affected dogs [34]. The results reported in the present pilot study may further strengthen the potential clinical utility of miR-30b-5p measurement in CKCSs affected by MMVD. We achieved this result by investigating, for the first time in veterinary medicine, the relationship between a plasma miRNA and the variation over time of several clinical and echocardiographic parameters in healthy and MMVD-affected CKCSs belonging to different ACVIM classes. Specifically, our results highlight a significant and negative marginal effect of miR-30b-5p on some important echocardiographic variables, since high levels of plasma miR-30b-5p corresponded to lower LVIDdN, ESVI, EDVI, and LA/Ao. As these parameters are pivotal to the MMVD diagnosis and monitoring because their increase over time reflects the progression of the disease [30,39–41], we could speculate that higher miR-30b-5p expression is related to an echocardiographically more stable form of MMVD over time in affected CKCSs.

These results are partially supported by the IPW analysis on the miRNA expression and the intensity of a heart murmur, the mitral regurgitation severity, and the MINE score. Indeed, although not statistically significant, these results suggest a positive marginal effect of miR-30b-5p on these three variables. High level of circulating miR-30b-5p corresponded to a higher probability of the class with the lower value (i.e., absence of a heart murmur, absence of a mitral regurgitation, and a mild MINE score). It is possible to hypothesize that the limited sample size of our study and the division of enrolled dogs into several classes of severity for each of these variables prevented us from obtaining statistically significant results. For this reason, studies enrolling a larger number of affected dogs are needed for validating and furtherly expand these results. Furthermore, it is fair to assume that the non-statistically significant results for the MINE score could be due to the fact that two of the parameters on which the score is calculated (i.e., E and FS%) showed no statistical significance on the IPW analysis.

Lastly, our results show a non-significant marginal effect of the elapsed time between T0 and T1 on all the examined variables. This is a somewhat surprising result, since MMVD is a disease that progresses over time. Our main hypothesis is that the lack of homogeneity of the follow-up period and its relatively limited extension (it ranged from 6 to 24 months) may have influenced the IPW analysis providing a hardly realistic result.

This pilot study presents some limitations, such as the small sample size and the limited number of subjects for each ACVIM class. Furthermore, the miR-30b-5p was evaluated only at T0 and not at T1. However, to fill such a gap, we are continuing to collect blood samples and data during the breed screenings, with the goal of setting up a larger prospective study that will include the evaluation of different miRNAs at established follow-up periods and, thus, trying to consolidate the preliminary findings from the present report. Another limitation was the variability of the follow-up, that did not always allow to have clinical and echocardiographic data associated to a fixed timing for each patient. The standardization of T1, and possibly subsequent timings, compared to T0 could clarify the variability of the biomarker due to the progression of the MMVD in this breed.

Conclusions

In conclusion, this pilot study allows to identify a relationship between high levels of miR-30b-5p and forms of MMVD that appear echocardiographically more stable over time in the CKCS.

Supporting information

S1 Table. Estimates of miR-30b-5p marginal effects for continuous and ordinal variables.

E = E wave velocity, EDVI = end-diastolic volume index, ESVI = end-systolic volume index, FS = fractional shortening, LA/Ao = left atrial-to-aortic root ratio, LVIDdN = left ventricular end-diastolic diameter normalized for body weight, Murmur = systolic heart murmur intensity. 95% Lower limit and 95% Upper limit: Lower and upper limits of 95% confidence interval for the mean marginal effect.

https://doi.org/10.1371/journal.pone.0274724.s001

(DOCX)

Acknowledgments

The Authors are thankful to the dog owners and breeders who trustingly and wholeheartedly participated to this research work.

Citation: Ghilardi S, Lecchi C, Bagardi M, Romito G, Colombo FM, Polli M, et al. (2022) Prospective pilot study on the predictive significance of plasma miR-30b-5p through the study of echocardiographic modifications in Cavalier King Charles Spaniels affected by different stages of myxomatous mitral valve disease: The PRIME study. PLoS ONE 17(12): e0274724. https://doi.org/10.1371/journal.pone.0274724

About the Authors:

Sara Ghilardi

Roles: Data curation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

Cristina Lecchi

Roles: Conceptualization, Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

Mara Bagardi

Roles: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

E-mail: [email protected], [email protected]

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

ORICD: https://orcid.org/0000-0002-5247-3796

Giovanni Romito

Roles: Validation, Writing – review & editing

Affiliation: Department of Veterinary Medical Sciences, Alma Mater Studiorum–University of Bologna, Bologna, Italy

Fabio M. Colombo

Roles: Formal analysis, Software

Affiliation: Department of Veterinary Sciences for Health, Animal Production and Food Security–Lodi (Italy), University of Milan, Milan, Italy

Michele Polli

Roles: Writing – review & editing

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

ORICD: https://orcid.org/0000-0003-0129-0161

Celeste Franco

Roles: Data curation

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

Paola G. Brambilla

Roles: Conceptualization, Project administration, Supervision, Validation, Writing – review & editing

Affiliation: Department of Veterinary Medicine and Animal Science, University of Milan, Lodi, Italy

ORICD: https://orcid.org/0000-0003-4255-555X

References

1. Borgarelli M, Buchanan JW. Historical review, epidemiology and natural history of degenerative mitral valve disease. J Vet Cardiol 2012;14:93–101. pmid:22386588

2. Ljungvall I, Häggström J. Adult-Onset Valvular Heart Disease. In: Ettinger SJ, Feldman EC, Côté E, editors. Textbook of Veterinary Internal Medicine. Diseases of the Dog and the Cat. Elsevier Saunders; Philadelphia, PA, USA: 2017. pp. 3033–57.

3. Lewis T, Swift S, Woolliams JA, Blott S. Heritability of premature mitral valve disease in Cavalier King Charles spaniels. Vet J 2011;188:73–6. pmid:20347358

4. Madsen MB, Olsen LH, Häggström J, Höglund K, Ljungvall I, Falk T, 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–7. pmid:21846748

5. Meurs KM, Friedenberg SG, Williams B, Keene BW, Atkins CE, Adin D, et al. Evaluation of genes associated with human myxomatous mitral valve disease in dogs with familial myxomatous mitral valve degeneration. Vet J 2018;232:16–9. pmid:29428085

6. O’Brien MJ, Beijerink NJ, Wade CM. Genetics of canine myxomatous mitral valve disease. Anim Genet. 2021 Aug;52(4):409–421. pmid:34028063

7. Bagardi M, Ghilardi S, Zamarian V, Ceciliani F, Brambilla PG, Lecchi C. Circulating MiR-30b-5p is upregulated in Cavalier King Charles Spaniels affected by early myxomatous mitral valve disease. PloS one 2022;17:e0266208. pmid:35816500

8. Wu Y, Li Q, Zhang R, Dai X, Chen W, Zing D. Circulating microRNAs: biomarkers of disease. Clin Chim Acta 2021 May;516:46–54. pmid:33485903

9. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell 2009;136:215–33. pmid:19167326

10. Zhang Y, Wang Z, Gemeinhart RA. Progress in microRNA delivery. J Control Release 2013;172:962–74. pmid:24075926

11. Gholaminejad A, Zare N, Dana N, Shafie D, Mani A, Javanmard SH. A meta-analysis of microRNA expression profiling studies in heart failure. Heart Fail Rev 2021;26:997–1021. pmid:33443726

12. Kalayinia S, Arjmand F, Maleki M, Malakootian M, Singh CP. MicroRNAs: roles in cardiovascular development and disease. Cardiovascular pathology: the official journal of the Society for Cardiov Pathol 2021;50,107296. pmid:33022373

13. Miretti S, Lecchi C, Ceciliani F, Baratta M. MicroRNAs as Biomarkers for Animal Health and Welfare in Livestock. Front Vet Sci 2020;18;7:578193. pmid:33392281

14. Nappi F, Iervolino A, Avtaar Singh SS, Chello M. MicroRNAs in Valvular Heart Diseases: Biological Regulators, Prognostic Markers and Therapeutical Targets. Int J Mol Sci 2021;22:12132. pmid:34830016

15. Lecchi C, Zamarian V, Gini C, Avanzini C, Polloni A, Rota Nodari S, et al. Salivary microRNAs are potential biomarkers for the accurate and precise identification of inflammatory response after tail docking and castration in piglets. J Anim Sci 2020;1;98:skaa153. pmid:32374403

16. Wang GK, Zhu JQ, Zhang JT, Li Q, Li Y, He J, et al. Circulating microRNA: a novel potential biomarker for early diagnosis of acute myocardial infarction in humans. Eur Heart J 2010;31:659–66. pmid:20159880

17. Wojciechowska A, Braniewska A, Kozar-Kamińska K. MicroRNA in cardiovascular biology and disease. Adv Clin Exp Med 2017;26:865–74. pmid:29068585

18. Hulanicka M, Garncarz M, Parzeniecka-Jaworska M, Jank M. Plasma miRNAs as potential biomarkers of chronic degenerative valvular disease in Dachshunds. BMC Vet Res 2014;26;10:205. pmid:25253075

19. Jung SW, Bohan A. Genome-wide sequencing and quantification of circulating microRNAs for dogs with congestive heart failure secondary to myxomatous mitral valve degeneration. Am J Vet Res 2018;79:163–9. pmid:29359980

20. Li Q, Freeman LM, Rush JE, Laflamme DP. Expression Profiling of Circulating MicroRNAs in Canine Myxomatous Mitral Valve Disease. Int J Mol Sci 2015;19;16:14098–108. pmid:26101868

21. Ro WB, Kang MH, Song DW, Kim HS, Lee GW, Park HM. Identification and Characterization of Circulating MicroRNAs as Novel Biomarkers in Dogs With Heart Diseases. Front Vet Sci 2021;8:729929. pmid:34708100

22. Yang VK, Loughran KA, Meola DM, Juhr CM, Thane KE, Davis AM, et al. Circulating exosome microRNA associated with heart failure secondary to myxomatous mitral valve disease in a naturally occurring canine model. J Extracell Vesicles. 2017;6:1350088. pmid:28804599

23. Reis-Ferreira A, Neto-Mendes J, Bras-Silva C, Lobo L, Fontes-Sousa AP. Emerging roles of microRNAs in Veterinary Cardiology. Vet Sci 2022 Sep 28;9:533. pmid:36288146

24. Rishniw M. Murmur grading in humans and animals: past and present. J Vet Cardiol. 2018;20:223–33. pmid:30017853

25. Acierno MJ et al. ACVIM consensus statement: Guidelines for the identification, evaluation, and management of systemic hypertension in dogs and cats. J Vet Intern Med. 2018;32:1803–22. pmid:30353952

26. Thomas WP, Gaber CE, Jacobs GJ, Kaplan PM, Lombard CW, Moise NS, et al. Recommendations for standards in transthoracic two-dimensional echocardiography in the dog and cat. J Vet Intern Med. 1993;7:247–52.

27. Keene BW, Atkins CE, Bonagura JD, Fox PR, Häggström J, Fuentes VL, et al. ACVIM consensus guidelines for the diagnosis and treatment of myxomatous mitral valve disease in dogs. J Vet Intern Med. 2019;33:1127–40. pmid:30974015

28. Cornell CC, Kittleson MD, Della Torre P, Häggström J, Lombard CW, Pedersen HD, et al. Allometric scaling of M-mode cardiac measurements in normal adult dogs. J Vet Intern Med. 2004;18:311–21. pmid:15188817

29. Marchesotti F, Vezzosi T, Tognetti R, Marchetti F, Patata V, Contieri B, et al. Left atrial anteroposterior diameter in dogs: reference interval, allometric scaling, and agreement with the left atrial-to-aortic root ratio. J Vet Med Sci. 2019;81:1655–62. pmid:31611490

30. Chetboul V, Tissier R. Echocardiographic assessment of canine degenerative mitral valve disease. J Vet Cardiol. 2012;14:127–48. pmid:22366573

31. Lombard CW. Normal values of the canine M-mode echocardiogram. Am J Vet Res. 1984;45:2015–18. pmid:6497098

32. Vezzosi T, Grosso G, Tognetti R, Meucci V, Patata V, Marchesotti F, et al. The Mitral INsufficiency Echocardiographic score: A severity classification of myxomatous mitral valve disease in dogs. J Vet Intern Med. 2021;35:1238–44. pmid:33951235

33. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4): 611–622. pmid:19246619

34. Bagardi M, Ghilardi S, Locatelli C, Bionda A, Polli M, Bussadori CM, et al. Influence of morphometry on echocardiographic measurements in Cavalier King Charles Spaniels: an inverse probability weighting analysis. Vet Sci. 2021;8:205. pmid:34679035

35. Leeper, TJ. Interpreting regression results using average marginal effects with R’s margins. 2017. https://www.semanticscholar.org/paper/Interpreting-Regression-Results-using-Average-with-Leeper/9615c76bd5d81f7ebbbdac9714619863dc3a2337.

36. Darke PG. Valvular incompetence in Cavalier King Charles Spaniels. Vet Rec. 1987;120:365–6. pmid:3590588

37. Häggström J, Hansson K, Kvart C, Swenson L. Chronic valvular disease in the Cavalier King Charles Spaniel in Sweden. Vet Rec. 1992;131:549–53. pmid:1481344

38. Pedersen HD, Lorentzen KA, Kristensen BO. Echocardiographic mitral valve prolapse in Cavalier King Charles Spaniels: epidemiology and prognostic significance for regurgitation. Vet Rec. 1999;144:315–20. pmid:10212505

39. Borgarelli M, Savarino P, Crosara S, Santilli RA, Chiavegato D, Poggi M, et al. Survival characteristics and prognostic variables of dogs with mitral regurgitation attributable to myxomatous valve disease. J Vet Intern Med. 2008;22:120–128. pmid:18289298

40. Häggström J, Boswood A, O’Grady M, Jöns O, Smith S, Swift S, et al. Effect of pimobendan or benazepril hydrochloride on survival times in dogs with congestive heart failure caused by naturally occurring myxomatous mitral valve disease: the QUEST study. J Vet Intern Med. 2008; 22:1124–35. pmid:18638016

41. Borgarelli M, Crosara S, Lamb K, Savarino P, La Rosa G, Tarducci A, et al. 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. pmid:22211523