- 1Institute of Molecular Biology and Genetics, Almazov National Medical Research Centre, Saint Petersburg, Russia
- 2Department of Bioengineering Systems, Saint Petersburg Electrotechnical University “LETI”, Saint Petersburg, Russia
- 3Department of Women's and Children's Health, Center for Molecular Medicine, Karolinska Institute, Stockholm, Sweden
Introduction: Pericardial fluid is enriched with biologically active molecules of cardiovascular origin including microRNAs. Investigation of the disease-specific extracellular microRNAs could shed light on the molecular processes underlying disease development. Arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart disease characterized by life-threatening arrhythmias and progressive heart failure development. The current data about the association between microRNAs and ARVC development are limited.
Methods and Results: We performed small RNA sequence analysis of microRNAs of pericardial fluid samples obtained during transcutaneous epicardial access for ventricular tachycardia (VT) ablation of six patients with definite ARVC and three post-infarction VT patients. Disease-associated microRNAs of pericardial fluid were identified. Five microRNAs (hsa-miR-1-3p, hsa-miR-21-5p, hsa-miR-122-5p, hsa-miR-206, and hsa-miR-3679-5p) were found to be differentially expressed between patients with ARVC and patients with post-infarction VT. Enrichment analysis of differentially expressed microRNAs revealed their close linkage to cardiac diseases.
Conclusion: Our data extend the knowledge of pericardial fluid microRNA composition and highlight five pericardial fluid microRNAs potentially linked to ARVC pathogenesis. Further studies are required to confirm the use of pericardial fluid RNA sequencing in differential diagnosis of ARVC.
Introduction
Pericardial fluid is a plasma ultrafiltrate found between the visceral and parietal pericardium, a double-walled sac surrounding the heart and roots of great vessels (1). Pericardial fluid is formed by the diffusion from both pericardial and epicardial vessels, as well as trans myocardial diffusion and acts as a lubricant, isolating, and protective layer for the beating heart (2). Due to a low clearance rate (3), pericardial fluid is enriched with heart-derived biologically active molecules, including cytokines, hormones, and nucleic acids (4–9). Thus, the composition of pericardial fluid potentially could influence the heart physiology and reflect certain heart conditions.
Among nucleic acids circulating in the pericardial fluid, extracellular microRNAs are of special interest. These short regulatory RNAs are extremely stable in the extracellular space due to the formation of RNA–protein complexes or packaging into vesicles, which protect them from nuclease digestion (10–14). Currently, data about pericardial microRNAs in patients with different structural heart diseases are lacking. A few studies reported on the differential expression of pericardial microRNAs in some common cardiovascular conditions. For instance, miR-423-5p was found to be enriched in pericardial fluid compared to serum; moreover, its levels differed in stable and unstable angina pectoris and aortic stenosis (15). Kuosmanen et al. (9) profiled microRNAs from the pericardial fluid of heart failure patients undergoing open-heart surgery and found no associations between microRNA profile and the clinical phenotype. In another study, exosomes from the pericardial fluid of patients were reported to contain myocardial-derived microRNAs promoting angiogenesis in vitro and in vivo (16). Taken together, these data suggest that microRNAs of pericardial fluid are biologically active molecules and could participate in cell-to-cell crosstalk.
Here, we focused on arrhythmogenic right ventricular cardiomyopathy (ARVC), an inherited disease difficult to diagnose and prognosticate (17, 18). At the early stages, ARVC is often asymptomatic or is manifested by heart rhythm abnormalities. During the ARVC progression, heart failure develops as a result of substitution of myocardium with fibrous and fat tissues (19). At the molecular level, ARVC is accompanied by significant changes in the signaling pathway activity (20–22). MicroRNAs as transcriptional regulators were reported to be involved in this pathological signaling network. Expression levels of miR-21 and miR-135b were found to be upregulated and downregulated correspondingly in the myocardium of ARVC patients (23). Transcriptome analysis of the HL-1 cells with PKP2 knockdown representing ARVC in vitro model revealed the downregulation of miR-184 (24). MicroRNA expression screening in plasma samples of ARVC patients and patients with idiopathic ventricular tachycardia (VT) revealed decreased expression of miR-320a associated with ARVC (25). Study performed on cardiac stromal cells disclosed three microRNAs (hsa-miR-520c-3p, hsa-miR-29b-3p, and hsa-miR-1183) differentially expressed between ARVC and control condition (26). Reduced blood plasma level of miR-320a-3p and elevated plasma levels of miR-144-3p, miR-145-5p, miR-185-5p, and miR-494-3p were reported to be associated with ARVC (27). A study involving transgenic mice carrying human DSG Q558* gene revealed different patterns of miRNA expression between the right and left ventricles: miR-217-5p and miR-708-5p were found to be upregulated and miR-499-5p was found to be downregulated specifically in the right ventricle (28). A recent study analyzing the microRNA expression profiles in blood and right ventricle tissue samples revealed that the expression of six microRNAs (miR-122-5p, miR-133a-3p, miR-133b, miR-142-3p, miR-182-5p, and miR-183-5p) was able to discriminate ARVC samples from healthy ones or other cardiomyopathy samples (29).
Although serum and myocardial microRNA expression evaluation has been reported, no common microRNA expression signature for ARVC is known. We suggested that microRNA composition of pericardial fluid in patients with ARVC might be specific, reflecting myocardial ARVC-related structural and molecular changes. We performed sequencing of microRNAs circulating in pericardial fluid of ARVC patients and patients with post-infarction VT (control group). We described the microRNA composition of pericardial fluid and performed an analysis of differentially expressed microRNAs.
Materials and Methods
Patient Characteristics
Patients referred for epicardial VT mapping and ablation between January 2019 and November 2020 were prospectively screened for inclusion into the study. Inclusion criteria were the following: indication to VT ablation; presumably epicardial VT exit site; a definite ARVC diagnosis, or the presence of a proven post-myocardial infarction scar; signed informed consent for the study. General exclusion criteria were the following: previous cardiothoracic surgery that potentially prevented pericardial manipulations; ongoing electrical storm (multiple defibrillation shocks in a short period of time); unsuccessful epicardial access; a previous ablation procedure <3 months ago; previous epicardial ablation; inadvertent right ventricle puncture during epicardial access, and/or visible blood in the pericardial fluid sample. Specific exclusion criteria for ARVC patients were other from “definite ARVC” diagnosis according to the established criteria; predominantly left ventricle disease and/or severe left ventricle systolic dysfunction (<40%); stenotic coronary artery disease; the history of angina and/or myocardial infarction. Specific exclusion criteria for ischemic group patients were severe right ventricle systolic dysfunction; indefinite myocardial scar; a coronary artery lesion requiring intervention. Patient clinical characteristics are summarized in Table 1.
Pericardial Fluid Collection
Epicardial access was performed under general anesthesia via a subxiphoid transcutaneous puncture under fluoroscopic guidance, as described in detail earlier (30). Special attention was paid to enter the pericardial space without any damage to the right ventricle and to a minimal use of contrast media. Once a sheath was introduced into the pericardial space, pericardial fluid was aspirated into an empty sterile syringe. The fluid was visually assessed for the presence of blood; blood-contaminated samples were discarded. A blood sample was collected from a femoral vein sheath immediately after the pericardial access. Collected pericardial fluid was centrifuged at 3,000 g at 4°C for 15 min, and then supernatant was collected, aliquoted, and stored at −80°C.
Genetic Testing
Genomic DNA was extracted from blood using FlexiGene DNA Kit (Qiagen). Target sequencing of 108 cardiomyopathy- and arrhythmia-associated genes was performed using Haloplex target enrichment (Agilent) with subsequent sequencing on MiSeq instrument (Illumina) as previously described (31). If disease-related genetic variants had not been identified, exome sequencing was performed as previously described (32). For the genetic variant verification, Sanger sequencing using a BigDye Terminator v3.1 kit and a 3,500 Genetic Analyzer (Applied Biosystems) was performed.
RNA Extraction
Before RNA extraction, samples were additionally centrifuged at 3,000 g at 4°C for 15 min and the obtained supernatant was used for RNA extraction. Small RNAs were extracted using SPLIT RNA Extraction Kit (Lexogen) according to manufacturer's recommendations.
Small RNA Library Preparation and Sequencing
Small RNA libraries were generated using Small RNA-Seq Kit (Lexogen) according to manufacturer's recommendations. The number of amplification cycles was 20 for all samples. Libraries were quantified using capillary gel electrophoresis using Bioanalyzer 2,100 (Agilent), then pooled in equimolar ratios based on 143-bp peak area, and purified in 6% PAAG gel using Gel Extraction Module (Lexogen). Sequencing was performed using MiSeq Reagent Kit v3 2x75bp and MiSeq equipment (Illumina) according to manufacturer's recommendations.
Bioinformatic Analyses
Obtained paired-end reads were merged using FLASH tool (33), length filtered 15–31 bp using Geneious Prime 2020.0.5 (https://www.geneious.com), and then aligned to the mature microRNA database (miRbase, http://www.mirbase.org/) using Novoalign implemented in mirPRo tool (34). Counts were normalized, and differential expression was calculated using R Studio version 1.2.5019 (35) with R version 3.0.1 (36) DESeq2 package (37). Hierarchical clustering and data visualization were performed in Phantasus version 1.5.1 (https://artyomovlab.wustl.edu/phantasus/). Tissue-specific expression profile of microRNAs was determined using human microRNA tissue atlas (38). MicroRNA set enrichment analyses were performed using TAM 2.0 tool (http://www.lirmed.com/tam2/) (39). The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus (40) and are accessible through GEO Series accession number GSE164490.
Quantitative PCR
Levels of selected microRNA were evaluated by qPCR. To remove heparin traces, RNA was treated with heparinase (Sigma) as was described before (41). For reverse transcription and qPCR microRNA-specific TaqMan Assays, TaqMan MicroRNA Reverse Transcription Kit and TaqMan Universal Master Mix II no UNG (all Thermo Fisher Scientific) were used according to manufacturer's recommendations. For miR-3679-5p measurement, a reverse transcription stem-loop primer and a primer pair for amplification were designed using sRNAPrimerDB online service (42). In this case, reverse transcription was performed using TaqMan MicroRNA Reverse Transcription Kit (Thermo Fisher Scientific), and real-time PCR was performed using qPCRmix-HS SYBR (Evrogen). TaqMan Assays and designed primer sequences are indicated in Supplementary Table 1.
Correlation Analysis
Pearson's correlation coefficient (r) and p-value were calculated using GraphPad Prism v.5.00 to explore the association between small RNA sequencing data and real-time PCR data. Linear regression was used to plot the line of best fit shown in each graph.
Results
Patient Characteristics
Among 40 patients undergoing epicardial ablation during the study period, nine subjects were eligible according to the inclusion and exclusion criteria: six patients with definite ARVC according to 2010 ARVC Task Force Criteria (19) and three control patients with coronary artery disease and the history of previous myocardial infarction, with no suspicion toward inherited channelopathy syndromes (Table 1). Four out of six ARVC patients (patients 4, 5, 7, 9) carried genetic variants in genes (PKP2, DSG2, and JUP) coding for desmosomal proteins—plakophilin-2, desmoglein-2, and junctional plakoglobin. One ARVC patient (patient 8) carried a genetic variant in FLNC gene coding for actin-binding filamin C protein. Target sequencing and subsequent whole-exome sequencing did not reveal disease-linked genetic variants in patient 6.
Presence of MicroRNA in Pericardial Fluid
Sequencing of microRNA revealed its presence in all pericardial fluid samples. In each sample, 145–411 microRNAs were detected, with the average number of 269 microRNAs per sample. Here, 105 microRNAs were shared between all samples (Table 2).
Description of Pericardial MicroRNAs
We analyzed microRNAs found in pericardial fluid according to the following criteria: (1) affiliation to a particular microRNA family; (2) cell- and tissue-specific expression. We identified 19 microRNA families presented in pericardial fluid by two or more microRNAs (Table 3). The most abundant in pericardial fluid microRNA family was the ubiquitous let-7 family with 11 detected microRNAs. Since the pericardial fluid could be enriched by cardiac microRNAs, we evaluated the presence of known cardiomyocyte-specific microRNAs: hsa-miR-1-3p, -133a-3p, -208a-3p, -208b-3p, -486-5p, and -486-3p. Hsa-miR-486-5p was detected in all samples, hsa-miR-1-3p was presented in six samples, hsa-miR-133a-3p in four samples, hsa-miR-486-3p was detected only in two samples, and hsa-miR-208a-3p and hsa-miR-208b-3p were not detected. Besides cardiomyocytes, the heart also consists of fibroblasts and endothelial cells. Thus, we identified in all pericardial fluid samples the microRNAs from miR-29 and miR-30 families, also having a high expression level in fibroblasts. Moreover, hsa-miR-21-5p, known to be expressed in cardiac fibroblasts (43), was presented in all studied samples. Two microRNAs (hsa-miR-93-5p; hsa-miR-106b-3p) from endothelial-specific family miR-17 were also detected in all pericardial fluid samples. Blood cell-derived microRNAs are also likely to contribute to pericardial fluid microRNA profile. Erythrocyte-specific hsa-miR-144-3p and hsa-miR-451a were detected in five samples and in all samples correspondingly.
Clustering and Differential Expression Analyses
Analysis of sample similarity revealed the high similarity of microRNA profiles of control group (post-infarction VT) and ARVC samples (Figure 1A). Principal component analysis (PCA) was performed in order to determine whether microRNA expression pattern is able to separate control and ARVC pericardial fluid samples. PCA showed no segregation of control or ARVC samples over the first two principal components (Figure 1B). Then, differential expression analysis was performed to reveal microRNA expression levels that significantly differ between patients with ARVC and the control group. Only microRNAs with at least one non-zero count between samples were taken into analyses. Although none of the microRNAs passed the multiple testing correction, we used non-corrected p-values taking into account the small group size and pilot nature of the study (Figure 1C). We found five differentially expressed microRNAs: two were downregulated in the ARVC group compared to the control group (hsa-miR-3679-5p and hsa-miR-21-5p), and three were upregulated in the ARVC group compared to the control group (hsa-miR-122-5p, hsa-miR-206, and hsa-miR-1-3p).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/647812/fcvm-08-647812-HTML/image_m/fcvm-08-647812-g001.jpg)
Figure 1. Results of samples clustering and microRNA differential expression analysis. (A) Similarity matrix of pericardial fluid samples calculated based on their microRNA expression patterns. (B) Principal component analysis of pericardial fluid samples. (C) Expression heatmap and result of statistical analyses of five differentially expressed microRNAs between control and arrhythmogenic right ventricular cardiomyopathy (ARVC) pericardial fluid samples.
Enrichment Analysis of Differentially Expressed MicroRNAs
To provide the functional annotation of microRNAs differentially expressed between ARVC and control groups and to prove these microRNAs are associated with disease development, we performed microRNA set enrichment analysis. Comparison of differentially expressed microRNAs found in pericardial fluid with databases of disease-associated microRNAs showed overrepresentation of cardiovascular pathology terms—chronic atrial fibrillation, coronary heart disease, arrhythmia, heart diseases, and hypertension. Also, various muscular pathology terms were found including muscular dystrophy, musculoskeletal abnormalities, and distal myopathy (Figure 2A). MicroRNA set mapping against databases containing biological processes resulted in diverse terms including cell cycle, heart and muscle development, inflammation, hormone-mediated signaling pathway, T-helper 17 cell differentiation, muscle development, skeletal muscle cell differentiation, cell death, cell proliferation, and cardiogenesis (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/647812/fcvm-08-647812-HTML/image_m/fcvm-08-647812-g002.jpg)
Figure 2. Enrichment analyses of microRNAs differentially expressed between control and arrhythmogenic right ventricular cardiomyopathy (ARVC) pericardial fluid samples. (A) MicroRNA set mapping against database containing disease-associated microRNA sets. (B) MicroRNA set mapping against database containing biological processes-associated microRNA sets. The number of microRNAs overlapping between the datasets is indicated in round brackets.
Quantitative PCR Validation of Differentially Expressed MicroRNAs
To validate microRNA expression data obtained by small RNA sequencing technique, we measured levels of differentially expressed microRNAs using qPCR and performed correlation analyses. Expression levels of three differentially expressed microRNAs (hsa-miR-1-3p, hsa-miR-21-5p, and hsa-miR-122-5p) measured by qPCR strongly correlated with sequencing data, whether expression levels of two other microRNAs (hsa-miR-206 and hsa-miR-3679-5p) demonstrated inconsistency between qPCR and small RNA sequencing (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/647812/fcvm-08-647812-HTML/image_m/fcvm-08-647812-g003.jpg)
Figure 3. Correlation analyses of microRNA expression between qPCR and small RNA sequencing data. (A) hsa-miR-1-3p. (B) hsa-miR-21-5p. (C) hsa-miR-122-5p. (D) hsa-miR-206. (E) hsa-miR-3679-5p. RQ, relative quantity. Pearson's correlation coefficient (r) and p-value are indicated in each plot.
Discussion
In the present study, we performed profiling of microRNAs in pericardial fluid samples obtained from patients with ARVC or post-infarction VT patients. In contrast to previous studies that used open-heart surgery for pericardial fluid collection (9, 15, 16), we obtained pericardial fluid samples during transcutaneous epicardial access. This approach allowed us to obtain high-quality pericardial fluid samples from rare patient groups.
To date, qPCR is a gold standard approach to detect microRNAs including a circulating pool of microRNAs (44, 45). However, rapidly developing technologies of next-generation sequencing (NGS) allow to perform accurate quantitative and qualitative assessment of nucleic acids including small RNAs extracted from solid or liquid tissue samples (46–48). In the current study, we performed sequencing of microRNAs extracted from pericardial fluid samples, described the pericardial fluid microRNA composition, and performed the differential expression analysis with subsequent microRNA set enrichment analysis.
In line with previously published reports (9, 16), we detected microRNAs in all investigated samples and revealed a similar spectrum of pericardial fluid microRNAs, indicating the validity of NGS-based approaches for microRNA detection. Pericardial fluid is formed by the diffusion from both pericardial and epicardial vessels, as well as trans-myocardial diffusion. Consequently, pericardial fluid microRNA repertoire is at least partially formed by secretion of cardiac cells—cardiomyocytes, endothelial cells, and cardiac fibroblasts (9). Since most microRNAs are expressed in a broad spectrum of cell types and tissues, it is hardly possible to precisely determine their origin. However, we detected several microRNAs known to be expressed predominantly in fibroblasts (including cardiac fibroblasts), endothelial cells, or erythrocytes. As was previously reported (9, 15), cardiac-specific microRNAs were not present in all pericardial fluid samples and their average expression levels were low. This observation indicated the moderate release of these microRNAs from cardiomyocytes and likely the absence of acute myocardial injury in the studied patient group, as opposed to early stages of myocardial infarction, which is accompanied by the elevation of cardiac microRNA levels in serum (49).
It has been suggested that microRNA composition of pericardial fluid could reflect cellular and molecular events underlying cardiac pathologies (9). Consistent with a previous report (9), high levels of five microRNAs associated with cardiac disease (let-7b-5p, hsa-miR-16–5p, hsa-miR-21–5p, hsa-miR-125b-5p, and hsa-miR-451a) were found in pericardial fluid samples (Table 2). Some of these microRNAs could be potentially relevant for ARVC pathogenesis even in the absence of severe cardiac pathological remodeling.
Expression analysis revealed five microRNAs differentially expressed between ARVC and control groups. Among differentially expressed microRNAs, hsa-miR-1-3p and hsa-miR-21-5p were reported to be highly expressed by cardiomyocytes and cardiac fibroblast correspondingly (50, 51). These two microRNAs are well-known to contribute to various cardiovascular diseases including ischemic heart injury, atrial fibrillation, and cardiomyopathies of different origins (52). Hsa-miR-206 was reported to be highly expressed in skeletal muscle but can also be present in the myocardium. Heart-specific overexpression of hsa-miR-206 in transgenic mice led to Cx43 downregulation and subsequently contributed to abnormal heart rate and PR interval and shortened life span. At the same time, hsa-miR-122-5p and hsa-miR-3679-5p do not reveal any specific heart-expression profile; while the first one is present in the liver in high amounts (53) and in blood cells at lower levels (53), the second one does not have any tissue-specific expression pattern (38). Hsa-miR-122-5p is essential for embryonic liver development and also was reported to regulate multiple physiological and pathological processes in the adult liver (53, 54). Intriguingly, hsa-miR-122-5p levels in heart tissue and blood samples were shown to discriminate arrhythmogenic cardiomyopathy patients from unaffected family members and patients with other cardiomyopathies (29). In contrast, little is known about hsa-miR-3679-5p, the novel player in cardiovascular biology. Originally, it was discovered in peripheral blood (55). Later, two histone demethylases, KDM7A and KDM6A (UTX), were identified as hsa-miR-3679-5p direct targets in monocytes. Downregulation of these genes by hsa-miR-3679-5p led to the reduction of adhesion molecules and regulation of monocyte adhesion to endothelial cells, which could be linked to an inflammatory response (56). Future studies are needed to prove whether these processes are relevant to ARVC. Surprisingly, a link between differentially expressed microRNAs and fibrosis-related genes encoding for proteins responsible for extracellular matrix deposition was found. Four out of five differentially expressed microRNAs (hsa-miR-1-3p, hsa-miR-21-5p, hsa-miR-206, and hsa-miR-122-5p) were reported to regulate directly or indirectly matrix metalloproteinase 2 gene (43, 57–59) and genes coding for collagen isoforms (60–63). These microRNAs were also reported to regulate vimentin expression—a protein being a strong marker of mesenchymal cell- and fibroblast-specific intermediate filament (64–67).
In line with these data, the enrichment analysis found a strong association between differentially expressed microRNA set and cardiovascular diseases including persistent atrial fibrillation, coronary artery disease, unspecified heart disease, and arrhythmias. A number of associations with skeletal muscle pathologies were found likely due to the overlap of microRNA expression profiles between cardiac and skeletal muscle tissues. Analyses of associations with various biological processes found the associations with muscle development and differentiation and with basic biological processes like cell cycle, cell proliferation, and cell death. Moreover, an association with inflammation, which is frequently concomitant with the heart pathology, was observed.
Despite the fact that small RNA sequencing could be used for accurate microRNA quantification, the variation of detection levels between different methods and platforms was reported (68, 69). The combination of sequencing data with subsequent qPCR analyses of selected targets allows to take advantage of both techniques and validate the results using an independent approach (69). We performed correlation analyses of five differentially expressed microRNA levels measured by small RNA sequencing and qPCR. Surprisingly, only a partial correlation between small RNA sequencing and qPCR results was observed, and a similar fact was earlier reported in other studies (68, 69). In the present study, only three out of five microRNAs differentially expressed between ARVC and control group showed a strong correlation between the two techniques used for quantification.
Our study has several limitations. First, similarly to the previous studies, we were not able to profile pericardial fluid microRNAs of healthy subjects due to the invasive technique of sample collection. Consequently, as a result, there was an age difference between the experimental groups that could introduce additional bias to the microRNA expression. In addition, the study included a very limited number of patients meeting the inclusion criteria. At last, a low number of differentially expressed genes were identified, which restrict the power of enrichment analyses.
In conclusion, we performed microRNA profiling of pericardial fluid obtained from patients with recurrent VT due to ARVC or previous myocardial infarction using small RNA sequencing technique. We described the pericardial fluid microRNA composition and revealed five differentially expressed microRNAs. Once confirmed in future studies with a larger number of patients, these microRNAs might be used in differential diagnosis of structural heart diseases in patients undergoing invasive procedures involving epicardial access.
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 in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by Local ethical committee of Almazov National Medical Research Centre. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
AKh collected samples, performed experiments, analyzed the data, and wrote the paper. DP performed experiments and analyzed the data. YF collected samples and performed experiments. AKn collected samples and performed experiments. KS collected samples and co-wrote the paper. DL performed epicardial ablation, collected samples, and co-wrote the paper. EM performed epicardial ablation, collected samples, and co-wrote the paper. AKo supervised the research, acquired funding, and co-wrote the paper. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the Ministry of Science and Higher Education of the Russian Federation (Agreement No. 075-15-2020-901).
Conflict of Interest
The 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
The authors are grateful to Anton Fedorov and Kirill Kondratov for the discussion and valuable comments and to Natalia Smolina for manuscript proofreading.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.647812/full#supplementary-material
References
1. Shabetai R, Mangiardi L, Bhargava V, Ross JJ, Higgins CB. The pericardium and cardiac function. Prog Cardiovasc Dis. (1979) 22:107–34. doi: 10.1016/0033-0620(79)90017-3
2. Vogiatzidis K, Zarogiannis SG, Aidonidis I, Solenov EI, Molyvdas P-A, Gourgoulianis KI, et al. Physiology of pericardial fluid production and drainage. Front Physiol. (2015) 6:62. doi: 10.3389/fphys.2015.00062
3. Boulanger B, Yuan Z, Flessner M, Hay J, Johnston M. Pericardial fluid absorption into lymphatic vessels in sheep. Microvasc Res. (1999) 57:174–86. doi: 10.1006/mvre.1998.2127
4. Liou J-Y, Kuan P-L, Shyu K-G, Lu M-J, Chao H-H, Wang B-W. Pericardial fluid and serum levels of vascular endothelial growth factor and endostatin in patients with or without coronary artery disease. J Formos Med Assoc. (2006) 105:377–83. doi: 10.1016/S0929-6646(09)60133-9
5. Santiago J-J, McNaughton LJ, Koleini N, Ma X, Bestvater B, Nickel BE, et al. High molecular weight fibroblast growth factor-2 in the human heart is a potential target for prevention of cardiac remodeling. PLoS ONE. (2014) 9:e97281. doi: 10.1371/journal.pone.0097281
6. Nishikimi T, Asakawa H, Iida H, Matsushita Y, Shibasaki I, Tadokoro K, et al. Different secretion patterns of two molecular forms of cardiac adrenomedullin in pressure- and volume-overloaded human heart failure. J Card Fail. (2004) 10:321–7. doi: 10.1016/j.cardfail.2003.12.002
7. Kuh JH, Kim KS, Kim SH, Cho KW, Seul KH, Koh GY. Presence of immunoreactive atrial natriuretic peptides in pericardial fluid of human subjects with congenital heart diseases. Life Sci. (1990) 46:1977–83. doi: 10.1016/0024-3205(90)90514-R
8. Tanaka T, Hasegawa K, Fujita M, Tamaki SI, Yamazato A, Kihara Y, et al. Marked elevation of brain natriuretic peptide levels in pericardial fluid is closely associated with left ventricular dysfunction. J Am Coll Cardiol. (1998) 31:399–403. doi: 10.1016/S0735-1097(97)00482-8
9. Kuosmanen SM, Hartikainen J, Hippeläinen M, Kokki H, Levonen AL, Tavi P. MicroRNA profiling of pericardial fluid samples from patients with heart failure. PLoS ONE. (2015) 10:1–23. doi: 10.1371/journal.pone.0119646
10. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci. (2011) 108:5003–8. doi: 10.1073/pnas.1019055108
11. Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. (2011) 39:7223–33. doi: 10.1093/nar/gkr254
12. Wagner J, Riwanto M, Besler C, Knau A, Fichtlscherer S, Röxe T, et al. Characterization of levels and cellular transfer of circulating lipoprotein-bound microRNAs. Arterioscler Thromb Vasc Biol. (2013) 33:1392–400. doi: 10.1161/ATVBAHA.112.300741
13. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. (2011) 13:423–33. doi: 10.1038/ncb2210
14. Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, et al. Detection of microRNA expression in human peripheral blood microvesicles. PLoS ONE. (2008) 3:e3694. doi: 10.1371/journal.pone.0003694
15. Miyamoto S, Usami S, Kuwabara Y, Horie T, Baba O, Hakuno D, et al. Expression patterns of miRNA-423-5p in the serum and pericardial fluid in patients undergoing cardiac surgery. PLoS ONE. (2015) 10:1–12. doi: 10.1371/journal.pone.0142904
16. Beltrami C, Besnier M, Shantikumar S, Shearn AIUU, Rajakaruna C, Laftah A, et al. Human pericardial fluid contains exosomes enriched with cardiovascular-expressed microRNAs and promotes therapeutic angiogenesis. Mol Ther. (2017) 25:679–93. doi: 10.1016/j.ymthe.2016.12.022
17. Corrado D, Basso C, Thiene G. Arrhythmogenic right ventricular cardiomyopathy: an update. Heart. (2009) 95:766–73. doi: 10.1136/hrt.2008.149823
18. Corrado D, Basso C, Judge DP. Arrhythmogenic cardiomyopathy. Circ Res. (2017) 121:784–802. doi: 10.1161/CIRCRESAHA.117.309345
19. Marcus FI, McKenna WJ, Sherrill D, Basso C, Bauce B, Bluemke DA, et al. Diagnosis of arrhythmogenic right ventricular cardiomyopathy/dysplasia: proposed modification of the task force criteria. Eur Heart J. (2010) 31:806–14. doi: 10.1093/eurheartj/ehq025
20. Dubash AD, Kam CY, Aguado BA, Patel DM, Delmar M, Shea LD, et al. Plakophilin-2 loss promotes TGF-β1/p38 MAPK-dependent fibrotic gene expression in cardiomyocytes. J Cell Biol. (2016) 212:425–38. doi: 10.1083/jcb.201507018
21. Chen SN, Gurha P, Lombardi R, Ruggiero A, Willerson JT, Marian AJ. The hippo pathway is activated and is a causal mechanism for adipogenesis in arrhythmogenic cardiomyopathy. Circ Res. (2014) 114:454–68. doi: 10.1161/CIRCRESAHA.114.302810
22. Kim C, Wong J, Wen J, Wang S, Wang C, Spiering S, et al. Studying arrhythmogenic right ventricular dysplasia with patient-specific iPSCs. Nature. (2013) 494:105–10. doi: 10.1038/nature11799
23. Zhang H, Liu S, Dong T, Yang J, Xie Y, Wu Y, et al. Profiling of differentially expressed microRNAs in arrhythmogenic right ventricular cardiomyopathy. Sci Rep. (2016) 6:28101. doi: 10.1038/srep28101
24. Gurha P, Chen X, Lombardi R, Willerson JT, Marian AJ. Knockdown of plakophilin 2 downregulates miR-184 through CpG hypermethylation and suppression of the E2F1 pathway and leads to enhanced adipogenesis in vitro. Circ Res. (2016) 119:731–50. doi: 10.1161/CIRCRESAHA.116.308422
25. Sommariva E, D'Alessandra Y, Farina FM, Casella M, Cattaneo F, Catto V, et al. MiR-320a as a potential novel circulating biomarker of arrhythmogenic cardiomyopathy. Sci Rep. (2017) 7:4802. doi: 10.1038/s41598-017-05001-z
26. Rainer J, Meraviglia V, Blankenburg H, Piubelli C, Pramstaller PP, Paolin A, et al. The arrhythmogenic cardiomyopathy-specific coding and non-coding transcriptome in human cardiac stromal cells. BMC Genomics. (2018) 19:1–12. doi: 10.1186/s12864-018-4876-6
27. Yamada S, Hsiao Y-W, Chang S-L, Lin Y-J, Lo L-W, Chung F-P, et al. Circulating microRNAs in arrhythmogenic right ventricular cardiomyopathy with ventricular arrhythmia. EP Eur. (2018) 20:f37–45. doi: 10.1093/europace/eux289
28. Calore M, Lorenzon A, Vitiello L, Poloni G, Khan MAF, Beffagna G, et al. A novel murine model for arrhythmogenic cardiomyopathy points to a pathogenic role of Wnt signalling and miRNA dysregulation. Cardiovasc Res. (2019) 115:739–51. doi: 10.1093/cvr/cvy253
29. Marinas MB, Celeghin R, Cason M, Bariani R, Frigo AC, Jager J, et al. A microRNA expression profile as non-invasive biomarker in a large arrhythmogenic cardiomyopathy cohort. Int J Mol Sci. (2020) 21:1536. doi: 10.3390/ijms21041536
30. Simonova KA, Mikhaylov EN, Tatarskiy RB, Kamenev AV, Panin DV, Orshanskaya VS, et al. Epicardial voltage mapping in patients with postinfarction ventricular tachycardia: a pilot study. J Arrhythmol. (2020) 27:22–7. doi: 10.35336/VA-2020-E-22-27
31. Kostareva A, Kiselev A, Gudkova A, Frishman G, Ruepp A, Frishman D, et al. Genetic spectrum of idiopathic restrictive cardiomyopathy uncovered by next-generation sequencing. PLoS One. (2016) 11:e0163362. doi: 10.1371/journal.pone.0163362
32. Kiselev A, Vaz R, Knyazeva A, Khudiakov A, Tarnovskaya S, Liu J, et al. De novo mutations in FLNC leading to early-onset restrictive cardiomyopathy and congenital myopathy. Hum Mutat. (2018) 39:1161–72. doi: 10.1002/humu.23559
33. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. doi: 10.1093/bioinformatics/btr507
34. Shi J, Dong M, Li L, Liu L, Luz-Madrigal A, Tsonis PA, et al. mirPRo–a novel standalone program for differential expression and variation analysis of miRNAs. Sci Rep. (2015) 5:14617. doi: 10.1038/srep14617
36. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Found Stat Comput (2013). Available online at: http://www.r-project.org/
37. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
38. Ludwig N, Leidinger P, Becker K, Backes C, Fehlmann T, Pallasch C, et al. Distribution of miRNA expression across human tissues. Nucleic Acids Res. (2016) 44:3865–77. doi: 10.1093/nar/gkw116
39. Li J, Han X, Wan Y, Zhang S, Zhao Y, Fan R, et al. TAM 2.0: tool for microRNA set analysis. Nucleic Acids Res. (2018) 46:W180–5. doi: 10.1093/nar/gky509
40. Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. (2002) 30:207–10. doi: 10.1093/nar/30.1.207
41. Kondratov K, Kurapeev D, Popov M, Sidorova M, Minasian S, Galagudza M, et al. Heparinase treatment of heparin-contaminated plasma from coronary artery bypass grafting patients enables reliable quantification of microRNAs. Biomol Detect Quantif. (2016) 8:9–14. doi: 10.1016/j.bdq.2016.03.001
42. Xie S, Zhu Q, Qu W, Xu Z, Liu X, Li X, et al. sRNAPrimerDB: comprehensive primer design and search web service for small non-coding RNAs. Bioinformatics. (2019) 35:1566–72. doi: 10.1093/bioinformatics/bty852
43. Roy S, Khanna S, Hussain S-RA, Biswas S, Azad A, Rink C, Gnyawali S, Shilo S, Nuovo GJ, Sen CK. MicroRNA expression in response to murine myocardial infarction: miR-21 regulates fibroblast metalloprotease-2 via phosphatase and tensin homologue. Cardiovasc Res. (2009) 82:21–29. doi: 10.1093/cvr/cvp015
44. Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. (2005) 33:e179. doi: 10.1093/nar/gni178
45. Jiang J, Lee EJ, Gusev Y, Schmittgen TD. Real-time expression profiling of microRNA precursors in human cancer cell lines. Nucleic Acids Res. (2005) 33:5394–403. doi: 10.1093/nar/gki863
46. Wyman SK, Parkin RK, Mitchell PS, Fritz BR, O'Briant K, Godwin AK, et al. Repertoire of microRNAs in epithelial ovarian cancer as determined by next generation sequencing of small RNA cDNA libraries. PLoS ONE. (2009) 4:e5311. doi: 10.1371/journal.pone.0005311
47. Willenbrock H, Salomon J, Søkilde R, Barken KB, Hansen TN, Nielsen FC, et al. Quantitative miRNA expression analysis: comparing microarrays with next-generation sequencing. RNA. (2009) 15:2028–34. doi: 10.1261/rna.1699809
48. Wang K, Li H, Yuan Y, Etheridge A, Zhou Y, Huang D, et al. The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PLoS ONE. (2012) 7:e51009. doi: 10.1371/journal.pone.0051009
49. Kuwabara Y, Ono K, Horie T, Nishi H, Nagao K, Kinoshita M, et al. Increased microRNA-1 and microRNA-133a levels in serum of patients with cardiovascular disease indicate myocardial damage. Circ Cardiovasc Genet. (2011) 4:446–54. doi: 10.1161/CIRCGENETICS.110.958975
50. Zhao Y, Ransom JF, Li A, Vedantham V, von Drehle M, Muth AN, et al. Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking miRNA-1-2. Cell. (2007) 129:303–17. doi: 10.1016/j.cell.2007.03.030
51. van Rooij E, Sutherland LB, Liu N, Williams AH, McAnally J, Gerard RD, et al. A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc Natl Acad Sci. (2006) 103:18255–60. doi: 10.1073/pnas.0608791103
52. Kura B, Kalocayova B, Devaux Y, Bartekova M. Potential Clinical Implications of miR-1 and miR-21 in Heart Disease and Cardioprotection. Int J Mol Sci. (2020) 21: doi: 10.3390/ijms21030700
53. Jopling C. Liver-specific microRNA-122: Biogenesis and function. RNA Biol. (2012) 9:137–42. doi: 10.4161/rna.18827
54. Bandiera S, Pfeffer S, Baumert TF, Zeisel MB. miR-122–a key factor and therapeutic target in liver disease. J Hepatol. (2015) 62:448–57. doi: 10.1016/j.jhep.2014.10.004
55. Vaz C, Ahmad HM, Sharma P, Gupta R, Kumar L, Kulshreshtha R, et al. Analysis of microRNA transcriptome by deep sequencing of small RNA libraries of peripheral blood. BMC Genomics. (2010) 11:288. doi: 10.1186/1471-2164-11-288
56. Higashijima Y, Matsui Y, Shimamura T, Nakaki R, Nagai N, Tsutsumi S, et al. Coordinated demethylation of H3K9 and H3K27 is required for rapid inflammatory responses of endothelial cells. EMBO J. (2020) 39:e103949. doi: 10.15252/embj.2019103949
57. Lv Y, Gao S, Zhang Y, Wang L, Chen X, Wang Y. miRNA and target gene expression in menstrual endometria and early pregnancy decidua. Eur J Obstet Gynecol Reprod Biol. (2016) 197:27–30. doi: 10.1016/j.ejogrb.2015.11.003
58. Yang Y, Li Q, Guo L. MicroRNA-122 acts as tumor suppressor by targeting TRIM29 and blocking the activity of PI3K/AKT signaling in nasopharyngeal carcinoma in vitro. Mol Med Rep. (2018) 17:8244–52. doi: 10.3892/mmr.2018.8894
59. Liu H, Cao Y-D, Ye W-X, Sun Y-Y. Effect of microRNA-206 on cytoskeleton remodelling by downregulating Cdc42 in MDA-MB-231 cells. Tumori. (2010) 96:751–5. doi: 10.1177/030089161009600518
60. Beaumont J, Lopez B, Hermida N, Schroen B, San Jose G, Heymans S, et al. microRNA-122 down-regulation may play a role in severe myocardial fibrosis in human aortic stenosis through TGF-beta1 up-regulation. Clin Sci. (2014) 126:497–506. doi: 10.1042/CS20130538
61. Wei L, Zhang Y, Qi X, Sun X, Li Y, Xu Y. Ubiquitin-proteasomes are the dominant mediators of the regulatory effect of microRNA-1 on cardiac remodeling after myocardial infarction. Int J Mol Med. (2019) 44:1899–907. doi: 10.3892/ijmm.2019.4330
62. Fry CS, Kirby TJ, Kosmac K, McCarthy JJ, Peterson CA. Myogenic progenitor cells control extracellular matrix production by fibroblasts during skeletal muscle hypertrophy. Cell Stem Cell. (2017) 20:56–69. doi: 10.1016/j.stem.2016.09.010
63. Liu S, Li W, Xu M, Huang H, Wang J, Chen X. Micro-RNA 21Targets dual specific phosphatase 8 to promote collagen synthesis in high glucose-treated primary cardiac fibroblasts. Can J Cardiol. (2014) 30:1689–99. doi: 10.1016/j.cjca.2014.07.747
64. Park YR, Seo SY, Kim SL, Zhu SM, Chun S, Oh J-M, et al. MiRNA-206 suppresses PGE2-induced colorectal cancer cell proliferation, migration, and invasion by targetting TM4SF1. Biosci Rep. (2018) 38: BSR20180664. doi: 10.1042/BSR20180664
65. Chou J, Wang B, Zheng T, Li X, Zheng L, Hu J, et al. MALAT1 induced migration and invasion of human breast cancer cells by competitively binding miR-1 with cdc42. Biochem Biophys Res Commun. (2016) 472:262–9. doi: 10.1016/j.bbrc.2016.02.102
66. Luo F, Ji J, Liu Y, Xu Y, Zheng G, Jing J, et al. MicroRNA-21, up-regulated by arsenite, directs the epithelial-mesenchymal transition and enhances the invasive potential of transformed human bronchial epithelial cells by targeting PDCD4. Toxicol Lett. (2015) 232:301–9. doi: 10.1016/j.toxlet.2014.11.001
67. Jin Y, Wang J, Han J, Luo D, Sun Z. MiR-122 inhibits epithelial-mesenchymal transition in hepatocellular carcinoma by targeting Snail1 and Snail2 and suppressing WNT/beta-cadherin signaling pathway. Exp Cell Res. (2017) 360:210–7. doi: 10.1016/j.yexcr.2017.09.010
68. Mestdagh P, Hartmann N, Baeriswyl L, Andreasen D, Bernard N, Chen C, et al. Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. (2014) 11:809–15. doi: 10.1038/nmeth.3014
Keywords: microRNA, small RNA sequencing, pericardial fluid, arrhythmogenic (right ventricular) cardiomyopathy, expression analyses
Citation: Khudiakov AA, Panshin DD, Fomicheva YV, Knyazeva AA, Simonova KA, Lebedev DS, Mikhaylov EN and Kostareva AA (2021) Different Expressions of Pericardial Fluid MicroRNAs in Patients With Arrhythmogenic Right Ventricular Cardiomyopathy and Ischemic Heart Disease Undergoing Ventricular Tachycardia Ablation. Front. Cardiovasc. Med. 8:647812. doi: 10.3389/fcvm.2021.647812
Received: 30 December 2020; Accepted: 15 February 2021;
Published: 19 March 2021.
Edited by:
Kun Wang, Qingdao University, ChinaCopyright © 2021 Khudiakov, Panshin, Fomicheva, Knyazeva, Simonova, Lebedev, Mikhaylov and Kostareva. 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: Aleksandr A. Khudiakov, khudyakov_aa@almazovcentre.ru
†These authors have contributed equally to this work