Parts of this article's content have been modified or rectified in:
Erratum: Expression Profile of Circulating MicroRNAs in Dogs With Cardiac Hypertrophy: A Pilot Study
Woong-Bin Ro
Min-Hee Kang
Doo-Won Song1- 1Department of Veterinary Internal Medicine, College of Veterinary Medicine, Konkuk University, Seoul, South Korea
- 2Department of Cancer Genome Research, Cancer Research Institute, Clinomics Inc., Seoul, South Korea
This study aimed to identify the expression profile of circulating microRNAs in dogs with eccentric or concentric cardiac hypertrophy. A total of 291 microRNAs in serum samples of five dogs with myxomatous mitral valve degeneration (MMVD) and five dogs with pulmonic stenosis (PS) were compared with those of five healthy dogs using microarray analysis. Results of microarray analysis revealed up-regulation of cfa-miR-130b [fold change (FC) = 2.13, p = 0.014), down-regulation of cfa-miR-375 (FC = 1.51, p = 0.014), cfa-miR-425 (FC = 2.56, p = 0.045), cfa-miR-30d (FC = 3.02, p = 0.047), cfa-miR-151 (FC = 1.89, p = 0.023), cfa-miR-19b (FC = 3.01, p = 0.008), and cfa-let-7g (FC = 2.53, p = 0.015) in MMVD group which showed eccentric cardiac hypertrophy, up-regulation of cfa-miR-346 (FC = 2.74, p = 0.032), down-regulation of cfa-miR-505 (FC = 1.56, p = 0.016) in PS group which showed concentric cardiac hypertrophy, and down-regulation of cfa-miR-30c (FC = 3.45, p = 0.013 in MMVD group; FC = 3.31, p = 0.014 in PS group) and cfa-let-7b (FC = 11.42, p = 0.049 in MMVD group; FC = 5.88, p = 0.01 in PS group) in both MMVD and PS groups. In addition, the unsupervised hierarchical clustering of differentially expressed microRNAs in each group resulted in complete separation of healthy dogs from dogs with heart diseases. Therefore, eleven microRNAs among 291 microRNAs were identified as differentially expressed circulating microRNAs related to MMVD or PS in dogs. This pilot study demonstrates that the microRNAs identified in this study could be possible candidates for novel biomarker or therapeutic target related to cardiac hypertrophy in dogs.
Introduction
Cardiac hypertrophy is an enlargement of the heart, which occurs in response to increased cardiac workload caused by various heart diseases (1). There are two types of cardiac hypertrophy such as eccentric and concentric cardiac hypertrophy. The eccentric cardiac hypertrophy is caused by volume overload diseases such as myxomatous mitral valve degeneration (MMVD), patent ductus arteriosus, and dilated cardiomyopathy in dogs, and results in increased myocyte length and dilated chambers with variable relative wall thickness (RWT) (2). The concentric cardiac hypertrophy is caused by pressure overload diseases such as pulmonic stenosis (PS) and aortic stenosis in dogs, and results in increased myocyte width and RWT with normal to decreased chamber volume (2).
MicroRNAs (miRNAs) are single-stranded, small non-coding RNAs containing 18–22 nucleotides that pair with specific target mRNAs to regulate and inhibit gene expression by degradation or translation inhibition (3). Previous studies have reported that miRNA play crucial roles in heart development and pathophysiology of cardiovascular diseases (4). In addition, circulating microRNAs are increasingly recognized as promising biomarkers for heart diseases because of their stability in peripheral blood (5). In humans, significant dysregulations of circulating microRNAs in cardiovascular diseases have been reported (5), and the role of microRNAs in cardiac hypertrophy has been identified (6–8). Moreover, the microRNAs related to cardiac hypertrophy are increasingly considered as potential therapeutic targets for various heart diseases (9).
Previous studies on circulating microRNAs in dogs with naturally occurring heart diseases are limited in number, and the studies have been conducted only in dogs with eccentric cardiac hypertrophy, such as MMVD and dilated cardiomyopathy (10–14). Therefore, the expression of microRNAs specific for concentric cardiac hypertrophy such as PS has not been studied yet in dogs. In this study, we aimed to demonstrate expression profiles of circulating microRNAs in dogs with eccentric or concentric cardiac hypertrophy, and investigate whether there is a difference in microRNA expression according to the type of cardiac hypertrophy. Dogs with MMVD or PS were included in this study, which are the two representative diseases of eccentric or concentric cardiac hypertrophy in dogs, respectively. We hypothesized that circulating microRNAs would be differentially expressed according to the type of cardiac hypertrophy or type of heart disease.
Materials and Methods
Serum Samples
Stored serum samples of five dogs with MMVD and with PS each were retrospectively retrieved and used in this study.
The inclusion criteria for MMVD group were dogs with MMVD of American College of Veterinary Internal Medicine (ACVIM) stage C or D, which were diagnosed and staged according to the guidelines of the ACVIM Consensus Statement for canine mitral valve disease, as previously described (15). According to the ACVIM guideline (15), stage C referred to dogs with past or current clinical signs by MMVD, and stage D indicated end-stage MMVD dogs, in which the clinical signs of heart failure were refractory to standard treatment used in stage C.
The inclusion criteria for PS group were dogs with PS of moderate or severe stenosis, which were diagnosed and staged according to the echocardiographic criteria for PS in dogs (16). According to the criteria (16), moderate stenosis referred to right ventricular outflow tract (RVOT) peak pressure gradient from 50 to 80 mmHg, and severe stenosis referred to RVOT peak pressure gradient above 80 mmHg.
Although the stages of MMVD and PS cannot be compared directly, progressed stages of the two diseases (stage C or D in MMVD, moderate or severe stenosis in PS) were included in this study in order to match the severity of the two diseases, as well as to ensure the eccentric and concentric cardiac hypertrophy in MMVD and PS groups, respectively.
The exclusion criteria for both MMVD and PS groups were dogs with any diseases other than MMVD or PS that could possibly affect circulating microRNA expression or cardiac hypertrophy, such as concurrent heart diseases, neoplasia, hypothyroidism, hyperadrenocorticism, diabetes mellitus, and systemic hypertension.
Five stored serum samples from healthy beagle dogs were used and they were retrieved from the remaining of a previous study (17). The previous study was approved by KBNP Institutional Animal Care and Use Committee (R0006046). The serum samples were not affected by the previous study because the blood was collected before initiation of the previous study, and were stored for 3 months before analysis. All healthy dogs in the previous study were examined to be healthy by physical examination, and there were no remarkable findings in complete blood count, serum chemistry, and urinalysis.
Blood was collected from jugular vein into 5-ml serum separating tubes (BD Vacutainer® SST™ Tube, USA). The tube was gently inverted about 5 times and was allowed to stand for 20–30 min at 4°C before centrifugation at 3,000 rpm for 15 min. The serum was aliquoted into cryovials (SPL Life Sciences, Korea) and stored at−70°C until use.
A total of 15 clinical data including breed, age, sex, heart rate, blood pressure, and results of physical examination were retrieved and evaluated from medical records from December 2014 to September 2018.
Echocardiographic Assessments
Echocardiographic data were acquired from the medical records. The EPIQ 7 ultrasound system (Philips Medical Systems, Andover, MA, USA) was utilized. A short-axis images of left ventricle (LV) were obtained at the left atrium and papillary muscle levels from the right-parasternal view. The values of LV measurements were normalized to the body weight according to the previous study (18). The fractional shortening (FS), ejection fraction (EF), end-diastolic volume index (EDVI), end-systolic volume index (ESVI), ratio of left atrial to aortic root diameter (LA/AO), normalized value of end-diastolic LV internal dimension (LVIDDN), normalized value of end-systolic LV internal dimension (LVIDSN), normalized value of end-diastolic LV free wall thickness (LVPWDN), normalized value of end-systolic LV free wall thickness (LVPWSN), normalized value of end-diastolic interventricular septal thickness (IVSDN), normalized value of end-systolic interventricular septal thickness (IVSSN), and RWT were analyzed for evaluation of cardiac function, hypertrophy, and remodeling of dogs included in this study. The RWT was calculated by following formula: RWT = (end-diastolic interventricular septal thickness + end-diastolic LV free wall thickness)/end-diastolic LV internal dimension (19).
RNA Preparation
RNA was extracted from serum using the miRNeasy Serum/Plasma Kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol. RNA integrity and quantity were evaluated with nanodrop 1000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), Quant-IT microRNA assay kit (Thermo Fisher Scientific) by QuantusTM Fluorometer (Promega, Madison, WI, USA), and Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, USA).
Gene Microarray Hybridization
The RNA was labeled using the FlashTagTM Biotin RNA Labeling Kit (Affymetrix, Santa Clara, CA, USA) and hybridized to GeneChip miRNA 4.0 microarrays (Affymetrix).
A hybridization mixture was added to samples with controls including control oligo B2, 20 X hybridization controls (bioB, bioC, bioD, cre), 27.5% Formamide, DMSO, 2 X hybridization buffer and water. 110.5 μl of this mixture was injected into GeneChip miRNA 4.0 arrays containing 291 canine probes annotated in miRBase Release 20, and placed in the Affymetrix GeneChip Hybridization Oven 640 at 48 °C and 60 rpm for 16 h overnight. The information on 291 microRNAs investigated in this study are shown in Supplementary Table 1. To amplify singal intensities, stain cocktails (stain cocktail 1 and 2) were added. Arrays were stained and washed using the Affymetrix GeneChip Fluidics Station 450 according to the FS450_0002 fluidics protocol.
Microarray Scanning and Data Processing
All arrays were scanned using the Affymetrix GeneChip Scanner 3000 (Affymetrix) and raw analysis performed with Transcriptome Analysis ConsoleTM (TAC) software (Affymetrix). The raw data images generated from the scanner were processed into CEL files, that contained measured intensities for each probe on the array. The fold change of > 1.5 (up or down regulation) and p-value of < 0.05 was used as a cut-off to screen wide range of differentially expressed microRNAs.
Statistical Analysis
The clinical data of dogs included in this study were presented as the mean ± standard deviation. The different expressions of microRNAs were compared between healthy vs. MMVD groups and healthy vs. PS groups, using Mann-Whitney test. The values of clinical and echocardiographic data for three groups were compared using Kruskal-Wallis test and post-hoc test. Statistical analysis was performed using the SPSS 25.0 software (SPSS Inc., Chicago, IL, USA). A p-value was considered statistically significant at the less of 0.05. Unsupervised hierarchical clustering of differentially expressed microRNAs was performed by Multi Experiment Viewer (MeV) software version 4.9.0, using average link hierarchical clustering and Euclidian distance measure.
Results
The clinical characteristics of dogs included in this study are shown in Table 1. The breed of dogs included in this study were five Beagle in healthy group, four Maltese and one Shih Tzu in MMVD group, and two Pomeranian, one French Bulldog, one Mixed-breed, and one Poodle in PS group. The mean age of dogs with MMVD was significantly higher than that of other groups (p = 0.007), and the mean diastolic blood pressure of dogs with PS was significantly higher when compared with healthy dogs (p = 0.039). In echocardiographic analysis, the FS and EF were significantly increased in both MMVD and PS group compared with healthy group (p = 0.006 and p = 0.006, respectively). The EDVI of PS group was significantly lower than that of healthy group (p = 0.039), and the ESVI was significantly decreased in both MMVD and PS group compared with healthy group (p = 0.004). The LA/AO of MMVD group was significantly higher than that of other groups (p = 0.002). The LVIDDN was significantly different between the three groups (p = 0.002), the highest in MMVD group and the lowest in PS group. The LVIDSN of PS group was significantly lower than that of other groups (p = 0.007). The IVSDN and IVSSN of PS group were significantly higher than those of healthy group (p = 0.030 and p = 0.041, respectively), and the RWT of dogs with PS was significantly higher than that of other groups (p = 0.007). Of the five dogs with MMVD, two dogs were stage C and three dogs were stage D based on the ACVIM staging system. All dogs with PS presented severe stenosis in valvular form. The representative echocardiographic images of dogs included in this study are shown in Figure 1.
TABLE 1
Table 1. Clinical characteristics of dogs included in this study.
FIGURE 1
Figure 1. Representative echocardiographic images of dogs included in (A) healthy, (B) MMVD, and (C) PS groups. Upper images are right parasternal 4-chamber view, and lower images are right parasternal short-axis view. All images were taken during the diastole. (B) Note the enlargement of chambers (LV and LA), eccentric hypertrophy of ventricular walls (white asterisk), and thickened mitral valve (arrowhead) in a dog with MMVD. (C) Also, note the concentric hypertrophy of ventricular walls (yellow asterisk) in a dog with PS. LA, left atrium; LV, left ventricle; MMVD, myxomatous mitral valve degeneration; PS, pulmonic stenosis; RA, right atrium; RV, right ventricle.
The circulating microRNA expression levels of dogs with MMVD or PS, compared with those of the healthy dogs are illustrated in the volcano plots in Figure 2. Using two independent parameters, a fold change cut-off of > 1.5 (up or down regulation) and p-value of <0.05, one microRNA (cfa-miR-130b) was significantly up-regulated and eight microRNAs (cfa-miR-375, cfa-miR-425, cfa-miR-30d, cfa-miR-30c, cfa-miR-151, cfa-let-7b, cfa-miR-19b, cfa-let-7g) were significantly down-regulated in dogs with MMVD when compared to healthy dogs (Figure 2A, Supplementary Table 2). In addition, one microRNA (cfa-miR-346) was significantly up-regulated and three microRNAs (cfa-miR-30c, cfa-let-7b, cfa-miR-505) were significantly down-regulated in dogs with PS when compared to healthy dogs (Figure 2B, Supplementary Table 2). The unsupervised hierarchical clustering of differentially expressed microRNAs resulted in complete separation of healthy dogs from dogs with heart diseases, exhibiting distinct microRNA expression patterns according to the sample types (Figure 3). Therefore, eleven microRNAs among 291 microRNAs were identified as differentially expressed circulating microRNAs related to MMVD or PS in dogs.
FIGURE 2
Figure 2. Volcano plots of circulating miRNA in dogs with (A) MMVD or (B) PS, compared with healthy dogs. The volcano plot shows the fold change and p-value between the dogs with heart diseases and healthy dogs for 291 miRNAs. The vertical lines correspond to a 1.5-fold differences, and the horizontal line displayed a p-value of 0.05. Therefore, the blue and red points in the plot represent the miRNAs with statistically significant differential expression in dogs with MMVD or PS (A or B, respectively) compared with healthy dogs. miRNA, microRNA; MMVD, myxomatous mitral valve degeneration; PS, pulmonic stenosis.
FIGURE 3
Figure 3. Heat map of significantly altered miRNA microarray expression data from serum samples of dogs with (A) MMVD and (B) PS compared with healthy dogs. Sample species are shown at the top and the miRNA species are shown on the right. Unsupervised hierarchical clustering of expression levels of nine and four (A and B, respectively) differentially expressed microRNAs using the Euclidian distance measure. The hierarchical clustering completely separated healthy dogs from dogs with heart diseases, showing distinct microRNA patterns among sample types. Red and green colors indicate relatively high and low expressions of microRNA, respectively. miRNA, microRNA; MMVD, myxomatous mitral valve degeneration; PS, pulmonic stenosis.
Discussion
Myxomatous mitral valve degeneration is the most common heart disease in dogs, and considered as excellent disease model for studying mechanisms that cause eccentric cardiac hypertrophy (24, 25). Pulmonic stenosis is the most common among heart diseases that cause concentric cardiac hypertrophy in dogs, and is characterized by right ventricular (RV) concentric hypertrophy, often RV dilation due to tricuspid regurgitation, and LV pseudo-hypertrophy secondary to decreased LV filling (26–28).
In this study, the distinct characteristics of eccentric and concentric cardiac hypertrophy of dogs with MMVD and PS, respectively, were confirmed by the significant differences in echocardiographic indices between the three groups (Table 1). In addition, similarly matched severity of MMVD and PS in this study ensured that substantial and similar degree of eccentric and concentric cardiac hypertrophy occurred in MMVD and PS groups, respectively. The remarkable increase of LA/AO and LVIDDN in MMVD group indicated prominent cardiac remodeling by volume overload. The RV concentric hypertrophy in PS group was confirmed by septal hypertrophy (increased IVSDN and IVSSN), which is known to occur secondary to increased RV afterload in dogs with PS (29). This septal hypertrophy in PS group may have been also partially contributed by LV pseudo-hypertrophy, which is characterized by decreased LV filling (decreased EDVI, ESVI, LA/AO, LVIDDN, LVIDSN) and consequent pseudo-increased LV wall thickness (IVSDN, IVSSN and RWT). The increase of mean systolic and diastolic blood pressure in PS group may have been due to increased systemic vascular resistance, as previously described in dogs with pulmonary artery occlusion (30) (Table 1). However, there was no statistical significance in systolic blood pressure, which is thought to be due to small sample size.
Of the eleven significantly dysregulated microRNAs in this study, miR-30c and let-7b were commonly dysregulated in both MMVD and PS groups. Therefore, these two microRNAs are thought to be involved in common pathway related to both eccentric and concentric cardiac hypertrophy in dogs.
In a previous study (31), miR-30c was down-regulated in pathological concentric cardiac hypertrophy, and the direct target gene of miR-30c was mRNA of connective tissue growth factor (CTGF), a key pro-fibrotic protein. In addition, down-regulation of circulating miR-30c was also observed in both eccentric and concentric cardiac hypertrophy, in dogs with MMVD and in heart failure patients with LV concentric hypertrophy, respectively (11, 32). Taken together, these results indicate that the down-regulation of miR-30c is associated with both eccentric and concentric cardiac hypertrophy, and the decrease of miR-30c may cause accumulation of CTGF, contributing to disease progression by pro-fibrotic effect. These results are consistent with the results of the present study, which showed down-regulation of miR-30c in both eccentric and concentric cardiac hypertrophy groups.
According to previous studies (33, 34), let-7b is known to be up-regulated by thioredoxin1 (Trx1), which is increased in response to various stress to heart including ischemia, heart failure, and pressure overload. These heart conditions such as ischemia and heart failure are common pathologic states that can be resulted from both eccentric and concentric cardiac hypertrophy, which explains the dysregulation of let-7b in both MMVD and PS groups in this study. In addition, it is reported that the let-7b up-regulated by Trx1 inhibits angiotensin-II induced cardiac hypertrophy by regulation of target, cyclin D2 (33). This anti-hypertrophic property of let-7b implicates its potential as promising therapeutic target in dogs with cardiac hypertrophy.
Meanwhile, let-7b was up-regulated in MMVD dogs in a previous study (12), whereas down-regulation was identified in this study. Similarly, discrepant expression of miR-30c was reported in a previous study in rats, in which miR-30c was down-regulated in heart failure but also showed up-regulation in response to therapeutic treatment (35). This opposite direction of regulation related to therapeutic treatment was also observed in previous studies on let-7b, which showed up-regulation in heart failure and down-regulation in response to treatment (35, 36). Based on these results, miR-30c and let-7b are thought to be associated with not only cardiac hypertrophy but also recovery of cardiac condition by treatment, and the expression of those microRNAs in heart diseases may be different depending on the stage of the heart disease or treatment status. For example, the let-7b may be up-regulated when the anti-hypertrophic compensatory mechanism is still working, however down-regulation of let-7b would be observed in non-compensated stage of heart disease. However, the reason of down-regulation of let-7b in response to treatment is required to be clarified in further studies. Meanwhile, the previous study in MMVD dogs reported that the expression level of let-7b showed tendency to be positively correlated to severity of MMVD (12). Although several mechanisms of gene regulation by miR-30c and let-7b in heart diseases have been identified in previous studies (31, 37–39), further studies on specific roles of these microRNAs in dogs with various cardiac states are required.
In this study, seven microRNAs (miR-130b, miR-375, miR-425, miR-30d, miR-151, miR-19b, let-7g) were significantly dysregulated in MMVD group. According to the previous studies (40, 41), miR-130b was reported to inhibit CYLD, which mediates lesion formation of vessels and cardiomyocytes, leading to cardiac dysfunction. The recovery of heart condition by alteration of microRNA was also reported in previous studies on miR-375 (42, 43), in which the inhibition of miR-375 attenuated inflammation and LV dysfunction, and also enhanced myocardial repair and function after myocardial infarction. The miR-425 has been reported to have anti-fibrotic effect on heart by inhibiting TGFβ1 (44). The expression of miR-30d, which is selectively enriched in cardiomyocytes, is induced by hypoxic stress and inhibits MAP4K4 to reduce cardiac apoptosis (45). The miR-151 targets PLM and is known to prevent ischemia-induced arrhythmias (46). The miR-19b also has cardio-protective property to regenerate cardiomyocyte in response to myocardial infarction (47). Lastly, let-7g inhibits cell adhesion and inflammation, while enhancing angiogenesis by targeting THBS1, which is involved in TGF-β pathway (48).
All these microRNAs dysregulated in MMVD group were reported to be dysregulated in previous studies on heart diseases in humans and rats (32, 35, 36, 42–44, 49–51), and four microRNAs among them (miR-130b, miR-425, miR-19b, let-7g) were identified to be dysregulated in previous studies in dogs with MMVD (11, 13). In the previous studies in dogs (11, 13), these four microRNAs showed down-regulation, which is consistent with the results of the present study except for miR-130b.
Contrary to the findings from previous studies in dogs (11, 13), up-regulation of miR-130b was observed in MMVD group in this study, and also in previous studies in rats and mice with heart diseases (49, 51). This discrepant expressions of miR-130b may be explained by the same reason as miR-30c and let-7b mentioned above. In the previous study in rats with induced cardiac hypertrophy (36), miR-130b was initially down-regulated in cardiac hypertrophy, but showed up-regulation in reverse-remodeling process in which the heart condition was restored. Therefore, similar to miR-30c and let-7b, miR-130b may be associated with protection of heart from cardiac hypertrophy and also recovery of cardiac condition by reverse-remodeling. In a previous human study (52), the miR-130b was down-regulated in the plasma of obese heart failure patients, and the biological pathways of predicted target genes of miR-130b were associated with prevention of cardiac hypertrophy, cardiomyocyte differentiation, and cardiac remodeling. This result corresponds to the results of the present and previous studies (11, 36, 49, 51). As mentioned above, miR-130b was reported to be protective for cardiac lesion and heart dysfunction (40, 41). Therefore, the up-regulation of miR-130b in MMVD group in this study may indicate that the miR-130b related protective pathway was activated in order to attenuate the progression of disease. However, this is difficult to be confirmed from the results of this study because serial evaluation of heart condition was not conducted. In addition, pathway causing miR-130b to be dysregulated is not identified clearly in dogs with heart diseases, therefore these speculations are needed to be verified in future studies.
Meanwhile, six microRNAs that were dysregulated only in MMVD group in this study (miR-130b, miR-375, miR-30d, miR-151, miR-19b, let-7g) were dysregulated also in concentric cardiac hypertrophy in previous murine studies (35, 36). In addition, as mentioned above, miR-425 was associated with cardiac fibrosis, which is a common pathologic change that can be resulted from both eccentric and concentric cardiac hypertrophy (44, 53). Therefore, although these seven microRNAs were significantly dysregulated only in MMVD group in this study, they are considered to be microRNAs that can be observed in both eccentric and concentric cardiac hypertrophy. The possible reasons why these microRNAs were not identified in the PS group in this study may include small number of samples and variable expression of microRNAs depending on the various factors, such as disease states as well as disease type.
The miR-346 and miR-505, which showed dysregulations only in the PS group, have not been reported in dogs until now. The target of miR-346 is Bax, which is a regulator of apoptosis, thus miR-346 inhibits apoptosis in myocardial ischemia-reperfusion injury (54). The up-regulation of miR-346 was observed in dogs with PS in this study, which is consistent with the result of previous study in which miR-346 was significantly up-regulated in ischemia reperfusion injured mouse hearts (55). Meanwhile, a previous study reported that miR-346 was up-regulated in ischemic myocardium of rats with myocardial infarction, and silencing of miR-346 significantly inhibited inflammatory response and apoptosis of cardiomyocytes (56). Since hearts with concentric hypertrophy are known to be susceptible to ischemic injury (57), miR-346 may have been predominantly expressed in dogs with PS in this study. The miR-505 was down-regulated in dogs with PS in this study. This microRNA interferes vascular regeneration by targeting FGF18 (58). In a previous study in neonatal mouse heart (59), miR-505-5p was down-regulated in heart tissue of 7-day-old mice compared to heart tissue of 1- and 6-day-old mice, and the target genes of miR-505-5p was associated with myocardial regenerative process. On the other hand, the up-regulation of miR-505-5p was observed in human patients with degenerative aortic stenosis. Since miR-346 and miR-505 were dysregulated only in PS group in this study and not identified in previous studies in dogs with MMVD, these microRNAs can be possibly associated with specific changes of PS or concentric cardiac hypertrophy.
It should be noted that the target genes or proteins, and related pathways of microRNAs described in this study were previously confirmed in other species, rather than dogs. Although most microRNAs and their targets are conserved among animals with different species, the conserved microRNAs do not always show the same expression patterns or levels in different species (60). In addition, multiple microRNAs are related to multiple target genes, thus identification of disease-specific target gene that is actually related to pathologic pathways is challenging (61). Therefore, further identification and validation of target genes in dogs are warranted.
The relationship of circulating microRNA with cardiac tissue microRNA remains largely unknown. However, recent study in humans showed strong positive correlation of circulating microRNA with pericardial microRNAs (62). In addition, several previous studies revealed that cardiac-enriched microRNAs were released to circulation in response to myocardial injury (63, 64). These results have raised potential of circulating microRNAs as novel diagnostic and prognostic biomarker. Nevertheless, there are still many challenges to use microRNAs as useful biomarker or therapeutic targets.
Despite of promising characteristics of circulating microRNAs, identifying the source of circulating microRNAs is very difficult, because many microRNAs have ubiquitous characteristics, and are related to variety of different clinical states (65). Therefore, it is unsure if circulating microRNAs identified in this study are representatives for cardiac microRNA. In addition, low expression of circulating microRNAs makes it more challenging to detect and accurately quantify them (66). Therefore, still large gap remains before clinical application of circulating microRNAs, which requires accumulation of further studies to identify and validate specific characteristics of circulating microRNAs.
This study has several limitations. First, the breed, age, and sex of healthy dogs were not matched with those of dogs with heart diseases, not satisfying the homogeneity of sample groups. In previous human studies (67–70), racial or ethnic difference, age, and sex were reported to be related to microRNA expression. However, association of microRNA expression with breed, age, and sex is not identified in dogs yet, which needs to be clarified in the future. Second, primary cardiac hypertrophy due to genetic cause could not be evaluated, because dogs with congenital cardiomyopathy were not included in this study. Third, the cardiac hypertrophy may not be the only factor associated with expression of circulating microRNAs in this study. Since the cardiac hypertrophy of dogs included in this study was one of various secondary changes to primary heart diseases, the changes in microRNA expression can be associated with not only cardiac hypertrophy but also other conditions such as disease type, myocardial infarction, cardiac fibrosis, and heart failure, as previously reported (71–73). In addition, the different location of cardiac hypertrophy between the two groups (LV hypertrophy in MMVD and RV hypertrophy in PS) should also be considered as another factor that can be associated with different microRNA expression. The factors associated with microRNAs are expected to be clarified by further studies on target genes and specific pathways of each candidate microRNA identified in this study. Lastly, due to small sample size and variable expression of microRNAs in different disease states, there may be unidentified microRNAs in either MMVD or PS groups in this study. Although the design of this study and analysis on the results were conducted according to pre-specified protocol, the results of this study should be verified by more specific study design. Further large-scale follow-up studies are required to validate and evaluate the microRNAs identified in this study as possible biomarkers and therapeutic targets in dogs with cardiac hypertrophy.
This study identified and demonstrated dysregulations of circulating microRNAs in dogs with two different types of cardiac hypertrophy. In addition, the expression profile of circulating microRNAs in dogs with PS was investigated for the first time in this study. This study provides a useful preliminary result for further validation and evaluation of microRNAs as biomarkers or therapeutic targets in dogs with eccentric or concentric cardiac hypertrophy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below:ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number “E-MTAB-10029” (https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-10029).
Ethics Statement
Ethical review and approval was not required for the animal study because stored serum samples were retrospectively used in this study. Informed consent was obtained from the owner for sample collection of client-owned dog. The serum samples of healthy dogs were retrieved from the remaining of a previous study (Ro W-B, Song D-W, Kim KH, Jeong SH, Kang MH. Pharmacokinetics and Pharmacodynamics of Pimobendan-Pentoxifylline Liquid Mixture After Oral Administration in Dogs. J Vet Clin (2019) 36(3):159-65.), which was approved by KBNP Institutional Animal Care and Use Committee (R0006046).
Author Contributions
W-BR, M-HK, S-HL, and H-MP conceived and designed the study. W-BR and D-WS curated the data and carried out the research. W-BR, M-HK, D-WS and H-MP analyzed the data. W-BR wrote the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of Interest
S-HL was employed by the company Clinomics Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This manuscript included a portion of a dissertation submitted by Woong-Bin Ro to the Department of Veterinary Medicine and the Graduate School of Konkuk University in partial fulfillment of the requirements for the degree of Doctor of Philosophy Veterinary Medicine.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.652224/full#supplementary-material
References
1. Barry SP, Davidson SM, Townsend PA. Molecular regulation of cardiac hypertrophy. Int J Bio Cell Biol. (2008) 40:2023–39. doi: 10.1016/j.biocel.2008.02.020
2. Bernardo BC, Weeks KL, Pretorius L, McMullen JR. Molecular distinction between physiological and pathological cardiac hypertrophy: experimental findings and therapeutic strategies. Pharmacol Ther. (2010) 128:191–227. doi: 10.1016/j.pharmthera.2010.04.005
3. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. (2004) 116:281–97. doi: 10.1016/S0092-8674(04)00045-5
4. Vegter EL, van der Meer P, de Windt LJ, Pinto YM, Voors AA. MicroRNAs in heart failure: from biomarker to target for therapy. Eur J Heart Fail. (2016) 18:457–68. doi: 10.1002/ejhf.495
5. Xu J, Zhao J, Evan G, Xiao C, Cheng Y, Xiao J. Circulating microRNAs: novel biomarkers for cardiovascular diseases. J Mol Med. (2012) 90:865–75. doi: 10.1007/s00109-011-0840-5
6. Wang J, Liew OW, Richards AM, Chen Y-T. Overview of microRNAs in cardiac hypertrophy, fibrosis, and apoptosis. Int J Mol Sci. (2016) 17:749. doi: 10.3390/ijms17050749
7. Zhou J, Zhang X, Liu Y, Gu S, Zhang X, An X, et al. microRNA-340-5p functions downstream of cardiotrophin-1 to regulate cardiac eccentric hypertrophy and heart failure via target gene dystrophin. Int Heart J. (2015) 56:454–8. doi: 10.1536/ihj.14-386
8. Li Q, Xie J, Wang B, Li R, Bai J, Ding L, et al. Overexpression of microRNA-99a attenuates cardiac hypertrophy. PLoS ONE. (2016) 11:e0148480. doi: 10.1371/journal.pone.0148480
9. Ooi JY, Bernardo BC, McMullen JR. The therapeutic potential of microRNAs regulated in settings of physiological cardiac hypertrophy. Future Med Chem. (2014) 6:205–22. doi: 10.4155/fmc.13.196
10. Hulanicka M, Garncarz M, Parzeniecka-Jaworska M, Jank M. Plasma microRNAs as potential biomarkers of chronic degenerative valvular disease in dachshunds. BMC Vet Res. (2014) 10:205. doi: 10.1186/s12917-014-0205-8
11. Jung S, 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. doi: 10.2460/ajvr.79.2.163
12. Li Q, Freeman L, Rush J, Laflamme D. Expression profiling of circulating microRNAs in canine myxomatous mitral valve disease. Int J Mol Sci. (2015) 16:14098–108. doi: 10.3390/ijms160614098
13. 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 Ext Ves. (2017) 6:1350088. doi: 10.1080/20013078.2017.1350088
14. Steudemann C, Bauersachs S, Weber K, Wess G. Detection and comparison of microRNA expression in the serum of doberman pinschers with dilated cardiomyopathy and healthy controls. BMC Vet Res. (2013) 9:12. doi: 10.1186/1746-6148-9-12
15. Atkins C, Bonagura J, Ettinger S, Fox P, Gordon S, Haggstrom J, et al. Guidelines for the diagnosis and treatment of canine chronic valvular heart disease. J Vet Intern Med. (2009) 23:1142–50. doi: 10.1111/j.1939-1676.2009.0392.x
16. Bussadori C, Amberger C, Le Bobinnec G, Lombard CW. Guidelines for the echocardiographic studies of suspected subaortic and pulmonic stenosis. J Vet Cardiol. (2000) 2:15–22. doi: 10.1016/S1760-2734(06)70007-8
17. Ro W-B, Song D-W, Kim K-H, Jeong S-H, Kang M-H. Pharmacokinetics and pharmacodynamics of pimobendan-pentoxifylline liquid mixture after oral administration in dogs. J Vet Clin. (2019) 36:159–65. doi: 10.17555/jvc.2019.06.36.3.159
18. Cornell CC, Kittleson MD, Torre PD, 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. doi: 10.1111/j.1939-1676.2004.tb02551.x
19. Foppa M, Duncan BB, Rohde LE. Echocardiography-based left ventricular mass estimation. How should we define hypertrophy? Card Ultr. (2005) 3:17. doi: 10.1186/1476-7120-3-17
20. Nelson RW, Couto CG. Small Animal Internal Medicine-E-Book. St. Louis, MO: Elsevier Health Sciences (2019).
21. Burkitt-Creedon JM, Davis H. Advanced Monitoring and Procedures for Small Animal Emergency and Critical Care. Ames, IA: John Wiley & Sons (2012).
22. Brown DJ, Rush JE, MacGregor J, Ross Jr JN, Brewer B, Rand WM. M-mode echocardiographic ratio indices in normal dogs, cats, and horses: a novel quantitative method. J Vet Intern Med. (2003) 17:653–62. doi: 10.1111/j.1939-1676.2003.tb02496.x
23. Serres F, Chetboul V, Tissier R, Poujol L, Gouni V, Carlos Sampedrano C, et al. Comparison of 3 ultrasound methods for quantifying left ventricular systolic function: correlation with disease severity and prognostic value in dogs with mitral valve disease. J Vet Intern Med. (2008) 22:566–77. doi: 10.1111/j.1939-1676.2008.0097.x
24. Matsuo T, Carabello BA, Nagatomo Y, Koide M, Hamawaki M, Zile MR, et al. Mechanisms of cardiac hypertrophy in canine volume overload. Am J Physiol Heart Circ Physiol. (1998) 275:H65–74. doi: 10.1152/ajpheart.1998.275.1.H65
25. Borgarelli M, Buchanan JW. Historical review, epidemiology and natural history of degenerative mitral valve disease. J Vet Cardiol. (2012) 14:93–101. doi: 10.1016/j.jvc.2012.01.011
26. Brambilla PG, Polli M, Pradelli D, Papa M, Rizzi R, Bagardi M, et al. Epidemiological study of congenital heart diseases in dogs: prevalence, popularity, and volatility throughout twenty years of clinical practice. PLoS ONE. (2020) 15:e0230160. doi: 10.1371/journal.pone.0230160
27. Ontiveros ES, Fousse SL, Crofton AE, Hodge TE, Gunther-Harrington CT, Visser LC, et al. Congenital cardiac outflow tract abnormalities in dogs: prevalence and pattern of inheritance from 2008 to 2017. Front Vet Sci. (2019) 6:52. doi: 10.3389/fvets.2019.00052
29. Tidholm A. Retrospective study of congenital heart defects in 151 dogs. J Small Anim Pract. (1997) 38:94–8. doi: 10.1111/j.1748-5827.1997.tb03326.x
30. Pinkerson AL, Kot PA. Systemic blood pressure response to changes in right ventricular function. Circ Res. (1964) 14:461–6. doi: 10.1161/01.RES.14.5.461
31. Duisters RF, Tijsen AJ, Schroen B, Leenders JJ, Lentink V, van der Made I, et al. miR-133 and miR-30 regulate connective tissue growth factor: implications for a role of microRNAs in myocardial matrix remodeling. Circ Res. (2009) 104:170–8. doi: 10.1161/CIRCRESAHA.108.182535
32. Watson CJ, Gupta SK, O'Connell E, Thum S, Glezeva N, Fendrich J, et al. MicroRNA signatures differentiate preserved from reduced ejection fraction heart failure. Eur J Heart Fail. (2015) 17:405–15. doi: 10.1002/ejhf.244
33. Yang Y, Ago T, Zhai P, Abdellatif M, Sadoshima J. Thioredoxin 1 negatively regulates angiotensin II–induced cardiac hypertrophy through upregulation of miR-98/let-7. Circ Res. (2011) 108:305–13. doi: 10.1161/CIRCRESAHA.110.228437
34. Ago T, Sadoshima J. Thioredoxin and ventricular remodeling. J Mol Cell Cardiol. (2006) 41:762–73. doi: 10.1016/j.yjmcc.2006.08.006
35. Dickinson BA, Semus HM, Montgomery RL, Stack C, Latimer PA, Lewton SM, et al. Plasma microRNAs serve as biomarkers of therapeutic efficacy and disease progression in hypertension-induced heart failure. Eur J Heart Fail. (2013) 15:650–9. doi: 10.1093/eurjhf/hft018
36. Wang J, Xu R, Lin F, Zhang S, Zhang G, Hu S, et al. MicroRNA: novel regulators involved in the remodeling and reverse remodeling of the heart. Cardiology. (2009) 113:81–8. doi: 10.1159/000172616
37. Mao L, Liu S, Hu L, Jia L, Wang H, Guo M. miR-30 family: a promising regulator in development and disease. Bio Med Res Int. (2018) 2018:9623412. doi: 10.1155/2018/9623412
38. Liu X, Li M, Peng Y, Hu X, Xu J, Zhu S, et al. miR-30c regulates proliferation, apoptosis and differentiation via the Shh signaling pathway in P19 cells. Exp Mol Med. (2016) 48:e248. doi: 10.1038/emm.2016.57
39. Bao M-H, Feng X, Zhang Y-W, Lou X-Y, Cheng Y, Zhou H-H. Let-7 in cardiovascular diseases, heart development and cardiovascular differentiation from stem cells. Int J Mol Sci. (2013) 14:23086–102. doi: 10.3390/ijms141123086
40. Xiao Z-Q, Yin T-K, Li Y-X, Zhang J-H, Gu J-J. miR-130b regulates the proliferation, invasion and apoptosis of glioma cells via targeting of CYLD. Oncol Rep. (2017) 38:167–74. doi: 10.3892/or.2017.5651
41. J Mathis B, Lai Y, Qu C, S Janicki J, Cui T. CYLD-mediated signaling and diseases. Curr Drug Targets. (2015) 16:284–94. doi: 10.2174/1389450115666141024152421
42. Garikipati VN, Verma SK, Jolardarashi D, Cheng Z, Ibetti J, Cimini M, et al. Therapeutic inhibition of miR-375 attenuates post-myocardial infarction inflammatory response and left ventricular dysfunction via PDK-1-AKT signalling axis. Cardiovasc Res. (2017) 113:938–49. doi: 10.1093/cvr/cvx052
43. Garikipati VNS, Krishnamurthy P, Verma SK, Khan M, Abramova T, Mackie AR, et al. Negative regulation of miR-375 by interleukin-10 enhances bone marrow-derived progenitor cell-mediated myocardial repair and function after myocardial infarction. Stem Cells. (2015) 33:3519–29. doi: 10.1002/stem.2121
44. Wang L, Liu J, Xu B, Liu Y-L, Liu Z. Reduced exosome miR-425 and miR-744 in the plasma represents the progression of fibrosis and heart failure. Kaohsiung J Med Sci. (2018) 34:626–33. doi: 10.1016/j.kjms.2018.05.008
45. Li J, Salvador AM, Li G, Valkov N, Ziegler O, Yeri A, et al. Mir-30d regulates cardiac remodeling by intracellular and paracrine signaling. Circ Res. (2021) 128:e1–23. doi: 10.1161/CIRCRESAHA.120.317244
46. Zhang Y, Wang R, Du W, Wang S, Yang L, Pan Z, et al. Downregulation of miR-151-5p contributes to increased susceptibility to arrhythmogenesis during myocardial infarction with estrogen deprivation. PLoS ONE. (2013) 8:e72985. doi: 10.1371/journal.pone.0072985
47. Gao F, Kataoka M, Liu N, Liang T, Huang Z-P, Gu F, et al. Therapeutic role of miR-19a/19b in cardiac regeneration and protection from myocardial infarction. Nat Comm. (2019) 10:1–15. doi: 10.1038/s41467-019-09530-1
48. Liao Y-C, Wang Y-S, Guo Y-C, Lin W-L, Chang M-H, Juo S-HH. Let-7g improves multiple endothelial functions through targeting transforming growth factor-beta and SIRT-1 signaling. J Am Coll Cardiol. (2014) 63:1685–94. doi: 10.1016/j.jacc.2013.09.069
49. Endo K, Naito Y, Ji X, Nakanishi M, Noguchi T, Goto Y, et al. MicroRNA 210 as a biomarker for congestive heart failure. Biol Pharm Bull. (2013) 36:48–54. doi: 10.1248/bpb.b12-00578
50. Zhu S, Cao L, Zhu J, Kong L, Jin J, Qian L, et al. Identification of maternal serum microRNAs as novel non-invasive biomarkers for prenatal detection of fetal congenital heart defects. Clin Chim Acta. (2013) 424:66–72. doi: 10.1016/j.cca.2013.05.010
51. Li L, Bounds KR, Chatterjee P, Gupta S. Micro RNA-130a, a potential antifibrotic target in cardiac fibrosis. J Am Heart Ass. (2017) 6:e006763. doi: 10.1161/JAHA.117.006763
52. Thomé JG, Mendoza MR, Cheuiche AV, La Porta VL, Silvello D, dos Santos KG, et al. Circulating microRNAs in obese and lean heart failure patients: a case–control study with computational target prediction analysis. Gene. (2015) 574:1–10. doi: 10.1016/j.gene.2015.07.068
53. Frey N, Katus HA, Olson EN, Hill JA. Hypertrophy of the heart: a new therapeutic target? Circulation. (2004) 109:1580–9. doi: 10.1161/01.CIR.0000120390.68287.BB
54. Lv X, Lu P, Hu Y, Xu T. miR-346 inhibited apoptosis against myocardial ischemia-reperfusion injury via targeting bax in rats. Drug Des Devel Ther. (2020) 14:895. doi: 10.2147/DDDT.S245193
55. Zhou L, Zang G, Zhang G, Wang H, Zhang X, Johnston N, et al. MicroRNA and mRNA signatures in ischemia reperfusion injury in heart transplantation. PLoS ONE. (2013) 8:e79805. doi: 10.1371/journal.pone.0079805
56. Yang B, Dong R, Zhao H. Inhibition of microRNA-346 inhibits myocardial inflammation and apoptosis after myocardial infarction via targeting NFIB. Eur Rev Med Pharmacol Sci. (2020) 24:11752–60.
57. Friehs I, del Nido PJ. Increased susceptibility of hypertrophied hearts to ischemic injury. Ann Thor Surg. (2003) 75:S678–84. doi: 10.1016/S0003-4975(02)04692-1
58. Huang YM, Li WW, Wu J, Han M, Li BH. The diagnostic value of circulating microRNAs in heart failure. Exp Ther Med. (2019) 17:1985–2003. doi: 10.3892/etm.2019.7177
59. Liu H, Zhu J, Liu Y, Fan Z, Zhu C, Qian L. Identification of the microRNA expression profile in the regenerative neonatal mouse heart by deep sequencing. Cell Biochem Biophys. (2014) 70:635–42. doi: 10.1007/s12013-014-9967-7
60. Ha M, Pang M, Agarwal V, Chen ZJ. Interspecies regulation of microRNAs and their targets. Bio Biop Acta Gene Regul Mech. (2008) 1779:735–42. doi: 10.1016/j.bbagrm.2008.03.004
61. Moqadam FA, Pieters R, Den Boer M. The hunting of targets: challenge in miRNA research. Leukemia. (2013) 27:16–23. doi: 10.1038/leu.2012.179
62. Cui C, Cui Q. The relationship of human tissue microRNAs with those from body fluids. Sci Rep. (2020) 10:1–7. doi: 10.1038/s41598-020-62534-6
63. De Rosa S, Eposito F, Carella C, Strangio A, Ammirati G, Sabatino J, et al. Transcoronary concentration gradients of circulating microRNAs in heart failure. Eur J Heart Fail. (2018) 20:1000–10. doi: 10.1002/ejhf.1119
64. Gidlöf O, Smith JG, Miyazu K, Gilje P, Spencer A, Blomquist S, et al. Circulating cardio-enriched microRNAs are associated with long-term prognosis following myocardial infarction. BMC Cardiovasc Disord. (2013) 13:1–9. doi: 10.1186/1471-2261-13-12
65. Small EM, Olson EN. Pervasive roles of microRNAs in cardiovascular biology. Nature. (2011) 469:336–42. doi: 10.1038/nature09783
66. Min PK, Chan SY. The biology of circulating micro RNA s in cardiovascular disease. Eur J Clin Invest. (2015) 45:860–74. doi: 10.1111/eci.12475
67. Dluzen DF, Hooten NN, Zhang Y, Kim Y, Glover FE, Tajuddin SM, et al. Racial differences in microRNA and gene expression in hypertensive women. Sci Rep. (2016) 6:1–14. doi: 10.1038/srep35815
68. Huang RS, Gamazon ER, Ziliak D, Wen Y, Im HK, Zhang W, et al. Population differences in microRNA expression and biological implications. RNA Biol. (2011) 8:692–701. doi: 10.4161/rna.8.4.16029
69. Huan T, Chen G, Liu C, Bhattacharya A, Rong J, Chen BH, et al. Age-associated micro RNA expression in human peripheral blood is associated with all-cause mortality and age-related traits. Aging Cell. (2018) 17:e12687. doi: 10.1111/acel.12687
70. Guo L, Zhang Q, Ma X, Wang J, Liang T. miRNA and mRNA expression analysis reveals potential sex-biased miRNA expression. Sci Rep. (2017) 7:1–10. doi: 10.1038/srep39812
71. Colpaert RM, Calore M. MicroRNAs in cardiac diseases. Cells. (2019) 8:737. doi: 10.3390/cells8070737
72. Wehbe N, Nasser SA, Pintus G, Badran A, Eid AH, Baydoun E. MicroRNAs in cardiac hypertrophy. Int J Mol Sci. (2019) 20:4714. doi: 10.3390/ijms20194714
Keywords: microRNA, dog, cardiac hypertrophy, myxomatous mitral valve degeneration, pulmonic stenosis, serum, novel biomarker, therapeutic target
Citation: Ro W-B, Kang M-H, Song D-W, Lee S-H and Park H-M (2021) Expression Profile of Circulating MicroRNAs in Dogs With Cardiac Hypertrophy: A Pilot Study. Front. Vet. Sci. 8:652224. doi: 10.3389/fvets.2021.652224
Received: 12 January 2021; Accepted: 15 March 2021;
Published: 09 April 2021.
Edited by:
Korakot Nganvongpanit, Chiang Mai University, ThailandReviewed by:
Sirilak Disatian Surachetpong, Chulalongkorn University, ThailandWanpitak Pongkan, Chiang Mai University, Thailand
Anudep Rungsipipat, Chulalongkorn University, Thailand
Copyright © 2021 Ro, Kang, Song, Lee and Park. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hee-Myung Park, parkhee@konkuk.ac.kr