Front. Cardiovasc. Med. Frontiers in Cardiovascular Medicine Front. Cardiovasc. Med. 2297-055X Frontiers Media S.A. 10.3389/fcvm.2021.647812 Cardiovascular Medicine Original Research Different Expressions of Pericardial Fluid MicroRNAs in Patients With Arrhythmogenic Right Ventricular Cardiomyopathy and Ischemic Heart Disease Undergoing Ventricular Tachycardia Ablation Khudiakov Aleksandr A. 1 * Panshin Daniil D. 1 Fomicheva Yulia V. 1 Knyazeva Anastasia A. 1 Simonova Ksenia A. 1 Lebedev Dmitry S. 1 2 Mikhaylov Evgeny N. 1 2 † Kostareva Anna A. 1 3 † 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

Edited by: Kun Wang, Qingdao University, China

Reviewed by: Zhihua Wang, Wuhan University, China; Luyu Zhou, Qingdao University, China

 *Correspondence: Aleksandr A. Khudiakov khudyakov_aa@almazovcentre.ru

This article was submitted to General Cardiovascular Medicine, a section of the journal Frontiers in Cardiovascular Medicine

†These authors have contributed equally to this work

 19 03 2021 2021 8 647812 30 12 2020 15 02 2021 Copyright © 2021 Khudiakov, Panshin, Fomicheva, Knyazeva, Simonova, Lebedev, Mikhaylov and Kostareva. 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.

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.

 microRNA small RNA sequencing pericardial fluid arrhythmogenic (right ventricular) cardiomyopathy expression analyses Ministry of Science and Higher Education of the Russian Federation10.13039/501100012190 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.

Patient demographics and clinical characteristics.

Patient Age Sex Disease-associated genetic variants Diagnosis Cardiac arrest ICD implanted LV EF, % RV dysfunction Heart failure, functional class Antiarrhythmic drugs* Number of VT induced Previous ablation, >3 months Previous ablation, <3 months Acute ablation results
1 73 M CAD, post-MI No No 23% No II BB 1 0 1 VT non-inducible
2 70 M CAD, post-MI No Yes 40% No II Amiodarone + BB 2 0 1 VT non-inducible
3 53 M CAD, post-MI No Yes 30% No III Sotalol 1 0 0 VT non-inducible
4 30 F PKP2 c.C235T, p.R79X, rs121434420 ARVC, definite No Yes 57% No I BB 0 0 1 VT non-inducible
5 53 M PKP2 c.2509delA, p.S837Vfs, rs727504432, DSG2 c.T3352A, p.S1118T ARVC, definite No Yes 64% Yes I Amiodarone + BB 1 0 0 VT non-inducible
6 59 F ARVC, definite No Yes 67% No I Amiodarone 0 0 0 VT non-inducible
7 38 M JUP c.G2105A, p.R702H, rs200690479 ARVC, definite No Yes 60% Yes I BB 0 0 0 VT non-inducible
8 20 M FLNC c.G3800A, p.R1267Q, rs768767784 ARVC, definite No Yes 53% No I BB VF 3 1 Clinical VT non-inducible
9 24 M DSG2 c.G671A, p.S224N ARVC, definite No No 50% No I BB 0 0 0 VT non-inducible

CAD, coronary artery disease; MI, myocardial infarction; ARVC, arrhythmogenic right ventricular cardiomyopathy; ICD, implantable cardioverter-defibrillator; VT, ventricular tachycardia; BB, beta-blocker; VF, ventricular fibrillation.

Antiarrhythmic drugs present during the ablation procedure.

 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).

List of microRNAs shared between all samples.

hsa-let-7a-3p hsa-miR-26a-5p hsa-miR-99b-5p hsa-miR-181a-5p hsa-miR-320c
hsa-let-7a-5p hsa-miR-26b-3p hsa-miR-99a-3p hsa-miR-181b-5p hsa-miR-335-3p
hsa-let-7b-5p hsa-miR-26b-5p hsa-miR-99a-5p hsa-miR-181c-5p hsa-miR-335-5p
hsa-let-7c-5p hsa-miR-27a-3p hsa-miR-100-5p hsa-miR-181d-5p hsa-miR-345-5p
hsa-let-7d-3p hsa-miR-27b-3p hsa-miR-101-3p hsa-miR-182-5p hsa-miR-361-3p
hsa-let-7d-5p hsa-miR-28-3p hsa-miR-106b-3p hsa-miR-183-5p hsa-miR-378a-3p
hsa-let-7e-5p hsa-miR-28-5p hsa-miR-122-5p hsa-miR-186-5p hsa-miR-421
hsa-let-7f-5p hsa-miR-29a-3p hsa-miR-125a-5p hsa-miR-191-5p hsa-miR-423-5p
hsa-let-7g-5p hsa-miR-29c-3p hsa-miR-125b-5p hsa-miR-192-5p hsa-miR-451a
hsa-let-7i-5p hsa-miR-30a-3p hsa-miR-130a-3p hsa-miR-193a-5p hsa-miR-484
hsa-miR-10a-3p hsa-miR-30a-5p hsa-miR-140-3p hsa-miR-195-5p hsa-miR-486-5p
hsa-miR-10a-5p hsa-miR-30b-5p hsa-miR-141-3p hsa-miR-199a-3p hsa-miR-497-5p
hsa-miR-10b-5p hsa-miR-30c-5p hsa-miR-143-3p hsa-miR-199b-3p hsa-miR-532-5p
hsa-miR-15b-5p hsa-miR-30d-5p hsa-miR-146a-5p hsa-miR-200a-3p hsa-miR-574-5p
hsa-miR-16-5p hsa-miR-30e-3p hsa-miR-146b-5p hsa-miR-200b-3p hsa-miR-652-3p
hsa-miR-19b-3p hsa-miR-30e-5p hsa-miR-148a-3p hsa-miR-203a-3p hsa-miR-744-5p
hsa-miR-21-5p hsa-miR-34a-5p hsa-miR-148b-3p hsa-miR-204-5p hsa-miR-769-5p
hsa-miR-22-3p hsa-miR-92a-3p hsa-miR-151a-3p hsa-miR-221-3p hsa-miR-888-5p
hsa-miR-23b-3p hsa-miR-93-5p hsa-miR-151a-5p hsa-miR-222-3p hsa-miR-1180-3p
hsa-miR-24-3p hsa-miR-95-3p hsa-miR-152-3p hsa-miR-320a-3p hsa-miR-1246
hsa-miR-25-3p hsa-miR-98-5p hsa-miR-181a-3p hsa-miR-320b hsa-miR-4286
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.

MicroRNA families with two or more microRNAs detected in all pericardial fluid samples.

MicroRNA family Number of microRNAs MicroRNAs
let-7 11 hsa-let-7a-3p; hsa-let-7a-5p; hsa-let-7b-5p; hsa-let-7c-5p; hsa-let-7d-3p; hsa-let-7d-5p; hsa-let-7e-5p; hsa-let-7f-5p; hsa-let-7g-5p; hsa-let-7i-5p; hsa-miR-98-5p
miR-10 9 hsa-miR-10a-3p; hsa-miR-10a-5p; hsa-miR-10b-5p; hsa-miR-99a-3p; hsa-miR-99a-5p; hsa-miR-99b-5p; hsa-miR-100-5p; hsa-miR-125a-5p; hsa-miR-125b-5p
miR-30 7 hsa-miR-30a-3p; hsa-miR-30a-5p; hsa-miR-30b-5p; hsa-miR-30c-5p; hsa-miR-30d-5p; hsa-miR-30e-3p; hsa-miR-30e-5p
miR-181 5 hsa-miR-181a-3p; hsa-miR-181a-5p; hsa-miR-181b-5p; hsa-miR-181c-5p; hsa-miR-181d-5p
miR-28 4 hsa-miR-28-3p; hsa-miR-28-5p; hsa-miR-151a-3p; hsa-miR-151a-5p
miR-8 3 hsa-miR-141-3p; hsa-miR-200a-3p; hsa-miR-200b-3p
miR-15 3 hsa-miR-15b-5p; hsa-miR-16-5p; hsa-miR-195-5p
miR-26 3 hsa-miR-26a-5p; hsa-miR-26b-3p; hsa-miR-26b-5p
miR-148 3 hsa-miR-148a-3p; hsa-miR-148b-3p; hsa-miR-152-3p
miR-320 3 hsa-miR-320a-3p; hsa-miR-320b; hsa-miR-320c
miR-17 2 hsa-miR-93-5p; hsa-miR-106b-3p
miR-25 2 hsa-miR-25-3p; hsa-miR-92a-3p
miR-27 2 hsa-miR-27a-3p; hsa-miR-27b-3p
miR-29 2 hsa-miR-29a-3p; hsa-miR-29c-3p
miR-95 2 hsa-miR-95-3p; hsa-miR-421
miR-146 2 hsa-miR-146a-5p; hsa-miR-146b-5p
miR-199 2 hsa-miR-199a-3p; hsa-miR-199b-3p
miR-221 2 hsa-miR-221-3p; hsa-miR-222-3p
miR-335 2 hsa-miR-335-3p; hsa-miR-335-5p
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).

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).

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).

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.

 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.

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 Shabetai R Mangiardi L Bhargava V Ross JJ Higgins CB. The pericardium and cardiac function. Prog Cardiovasc Dis. (1979) 22:107–34. 10.1016/0033-0620(79)90017-3 Vogiatzidis K Zarogiannis SG Aidonidis I Solenov EI Molyvdas P-A Gourgoulianis KI . Physiology of pericardial fluid production and drainage. Front Physiol. (2015) 6:62. 10.3389/fphys.2015.0006225852564 Boulanger B Yuan Z Flessner M Hay J Johnston M. Pericardial fluid absorption into lymphatic vessels in sheep. Microvasc Res. (1999) 57:174–86. 10.1006/mvre.1998.212710049665 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. 10.1016/S0929-6646(09)60133-916638647 Santiago J-J McNaughton LJ Koleini N Ma X Bestvater B Nickel BE . 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. 10.1371/journal.pone.009728124827991 Nishikimi T Asakawa H Iida H Matsushita Y Shibasaki I Tadokoro K . 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. 10.1016/j.cardfail.2003.12.00215309699 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. 10.1016/0024-3205(90)90514-R2141890 Tanaka T Hasegawa K Fujita M Tamaki SI Yamazato A Kihara Y . 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. 10.1016/S0735-1097(97)00482-89462585 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. 10.1371/journal.pone.011964625763857 Arroyo JD Chevillet JR Kroh EM Ruf IK Pritchard CC Gibson DF . Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci. (2011) 108:5003–8. 10.1073/pnas.101905510821383194 Turchinovich A Weiz L Langheinz A Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. (2011) 39:7223–33. 10.1093/nar/gkr254 Wagner J Riwanto M Besler C Knau A Fichtlscherer S Röxe T . Characterization of levels and cellular transfer of circulating lipoprotein-bound microRNAs. Arterioscler Thromb Vasc Biol. (2013) 33:1392–400. 10.1161/ATVBAHA.112.30074123559634 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. 10.1038/ncb221021423178 Hunter MP Ismail N Zhang X Aguda BD Lee EJ Yu L . Detection of microRNA expression in human peripheral blood microvesicles. PLoS ONE. (2008) 3:e3694. 10.1371/journal.pone.000369419002258 Miyamoto S Usami S Kuwabara Y Horie T Baba O Hakuno D . Expression patterns of miRNA-423-5p in the serum and pericardial fluid in patients undergoing cardiac surgery. PLoS ONE. (2015) 10:1–12. 10.1371/journal.pone.014290426562412 Beltrami C Besnier M Shantikumar S Shearn AIUU Rajakaruna C Laftah A . Human pericardial fluid contains exosomes enriched with cardiovascular-expressed microRNAs and promotes therapeutic angiogenesis. Mol Ther. (2017) 25:679–93. 10.1016/j.ymthe.2016.12.02228159509 Corrado D Basso C Thiene G. Arrhythmogenic right ventricular cardiomyopathy: an update. Heart. (2009) 95:766–73. 10.1136/hrt.2008.149823 Corrado D Basso C Judge DP. Arrhythmogenic cardiomyopathy. Circ Res. (2017) 121:784–802. 10.1161/CIRCRESAHA.117.309345 Marcus FI McKenna WJ Sherrill D Basso C Bauce B Bluemke DA . Diagnosis of arrhythmogenic right ventricular cardiomyopathy/dysplasia: proposed modification of the task force criteria. Eur Heart J. (2010) 31:806–14. 10.1093/eurheartj/ehq02520172912 Dubash AD Kam CY Aguado BA Patel DM Delmar M Shea LD . Plakophilin-2 loss promotes TGF-β1/p38 MAPK-dependent fibrotic gene expression in cardiomyocytes. J Cell Biol. (2016) 212:425–38. 10.1083/jcb.20150701826858265 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. 10.1161/CIRCRESAHA.114.302810 Kim C Wong J Wen J Wang S Wang C Spiering S . Studying arrhythmogenic right ventricular dysplasia with patient-specific iPSCs. Nature. (2013) 494:105–10. 10.1038/nature11799 Zhang H Liu S Dong T Yang J Xie Y Wu Y . Profiling of differentially expressed microRNAs in arrhythmogenic right ventricular cardiomyopathy. Sci Rep. (2016) 6:28101. 10.1038/srep28101 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. 10.1161/CIRCRESAHA.116.308422 Sommariva E D'Alessandra Y Farina FM Casella M Cattaneo F Catto V . MiR-320a as a potential novel circulating biomarker of arrhythmogenic cardiomyopathy. Sci Rep. (2017) 7:4802. 10.1038/s41598-017-05001-z28684747 Rainer J Meraviglia V Blankenburg H Piubelli C Pramstaller PP Paolin A . The arrhythmogenic cardiomyopathy-specific coding and non-coding transcriptome in human cardiac stromal cells. BMC Genomics. (2018) 19:1–12. 10.1186/s12864-018-4876-629940860 Yamada S Hsiao Y-W Chang S-L Lin Y-J Lo L-W Chung F-P . Circulating microRNAs in arrhythmogenic right ventricular cardiomyopathy with ventricular arrhythmia. EP Eur. (2018) 20:f37–45. 10.1093/europace/eux28929036525 Calore M Lorenzon A Vitiello L Poloni G Khan MAF Beffagna G . A novel murine model for arrhythmogenic cardiomyopathy points to a pathogenic role of Wnt signalling and miRNA dysregulation. Cardiovasc Res. (2019) 115:739–51. 10.1093/cvr/cvy25330304392 Marinas MB Celeghin R Cason M Bariani R Frigo AC Jager J . A microRNA expression profile as non-invasive biomarker in a large arrhythmogenic cardiomyopathy cohort. Int J Mol Sci. (2020) 21:1536. 10.3390/ijms2104153632102357 Simonova KA Mikhaylov EN Tatarskiy RB Kamenev AV Panin DV Orshanskaya VS . Epicardial voltage mapping in patients with postinfarction ventricular tachycardia: a pilot study. J Arrhythmol. (2020) 27:22–7. 10.35336/VA-2020-E-22-27 Kostareva A Kiselev A Gudkova A Frishman G Ruepp A Frishman D . Genetic spectrum of idiopathic restrictive cardiomyopathy uncovered by next-generation sequencing. PLoS One. (2016) 11:e0163362. 10.1371/journal.pone.016336227662471 Kiselev A Vaz R Knyazeva A Khudiakov A Tarnovskaya S Liu J . De novo mutations in FLNC leading to early-onset restrictive cardiomyopathy and congenital myopathy. Hum Mutat. (2018) 39:1161–72. 10.1002/humu.2355929858533 Magoc T Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. 10.1093/bioinformatics/btr50721903629 Shi J Dong M Li L Liu L Luz-Madrigal A Tsonis PA . mirPRo–a novel standalone program for differential expression and variation analysis of miRNAs. Sci Rep. (2015) 5:14617. 10.1038/srep1461726434581 RStudio Team. RStudio: Integrated Development for R. Boston, MA: RStudio, Inc. (2015). 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/ Love MI Huber W Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. 10.1186/s13059-014-0550-825516281 Ludwig N Leidinger P Becker K Backes C Fehlmann T Pallasch C . Distribution of miRNA expression across human tissues. Nucleic Acids Res. (2016) 44:3865–77. 10.1093/nar/gkw11626921406 Li J Han X Wan Y Zhang S Zhao Y Fan R . TAM 2.0: tool for microRNA set analysis. Nucleic Acids Res. (2018) 46:W180–5. 10.1093/nar/gky50929878154 Edgar R Domrachev M Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. (2002) 30:207–10. 10.1093/nar/30.1.20711752295 Kondratov K Kurapeev D Popov M Sidorova M Minasian S Galagudza M . Heparinase treatment of heparin-contaminated plasma from coronary artery bypass grafting patients enables reliable quantification of microRNAs. Biomol Detect Quantif . (2016) 8:9–14. 10.1016/j.bdq.2016.03.00127335806 Xie S Zhu Q Qu W Xu Z Liu X Li X . sRNAPrimerDB: comprehensive primer design and search web service for small non-coding RNAs. Bioinformatics. (2019) 35:1566–72. 10.1093/bioinformatics/bty85230295699 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. 10.1093/cvr/cvp01519147652 Chen C Ridzon DA Broomer AJ Zhou Z Lee DH Nguyen JT . Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. (2005) 33:e179. 10.1093/nar/gni17819374024 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. 10.1093/nar/gki86316192569 Wyman SK Parkin RK Mitchell PS Fritz BR O'Briant K Godwin AK . Repertoire of microRNAs in epithelial ovarian cancer as determined by next generation sequencing of small RNA cDNA libraries. PLoS ONE. (2009) 4:e5311. 10.1371/journal.pone.000531119390579 Willenbrock H Salomon J Søkilde R Barken KB Hansen TN Nielsen FC . Quantitative miRNA expression analysis: comparing microarrays with next-generation sequencing. RNA. (2009) 15:2028–34. 10.1261/rna.169980919745027 Wang K Li H Yuan Y Etheridge A Zhou Y Huang D . The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PLoS ONE. (2012) 7:e51009. 10.1371/journal.pone.005100923251414 Kuwabara Y Ono K Horie T Nishi H Nagao K Kinoshita M . Increased microRNA-1 and microRNA-133a levels in serum of patients with cardiovascular disease indicate myocardial damage. Circ Cardiovasc Genet. (2011) 4:446–54. 10.1161/CIRCGENETICS.110.95897521642241 Zhao Y Ransom JF Li A Vedantham V von Drehle M Muth AN . Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking miRNA-1-2. Cell. (2007) 129:303–17. 10.1016/j.cell.2007.03.03017397913 van Rooij E Sutherland LB Liu N Williams AH McAnally J Gerard RD . A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc Natl Acad Sci. (2006) 103:18255–60. 10.1073/pnas.060879110317108080 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: 10.3390/ijms2103070031973111 Jopling C. Liver-specific microRNA-122: Biogenesis and function. RNA Biol. (2012) 9:137–42. 10.4161/rna.1882722258222 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. 10.1016/j.jhep.2014.10.00425308172 Vaz C Ahmad HM Sharma P Gupta R Kumar L Kulshreshtha R . Analysis of microRNA transcriptome by deep sequencing of small RNA libraries of peripheral blood. BMC Genomics. (2010) 11:288. 10.1186/1471-2164-11-28820459673 Higashijima Y Matsui Y Shimamura T Nakaki R Nagai N Tsutsumi S . Coordinated demethylation of H3K9 and H3K27 is required for rapid inflammatory responses of endothelial cells. EMBO J. (2020) 39:e103949. 10.15252/embj.201910394932125007 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. 10.1016/j.ejogrb.2015.11.00326699100 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. 10.3892/mmr.2018.8894 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. 10.1177/03008916100960051821302623 Beaumont J Lopez B Hermida N Schroen B San Jose G Heymans S . 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. 10.1042/CS2013053824168656 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. 10.3892/ijmm.2019.4330 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. 10.1016/j.stem.2016.09.01027840022 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. 10.1016/j.cjca.2014.07.74725418215 Park YR Seo SY Kim SL Zhu SM Chun S Oh J-M . MiRNA-206 suppresses PGE2-induced colorectal cancer cell proliferation, migration, and invasion by targetting TM4SF1. Biosci Rep. (2018) 38: BSR20180664. 10.1042/BSR2018066430135139 Chou J Wang B Zheng T Li X Zheng L Hu J . 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. 10.1016/j.bbrc.2016.02.10231708098 Luo F Ji J Liu Y Xu Y Zheng G Jing J . 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. 10.1016/j.toxlet.2014.11.00125445583 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. 10.1016/j.yexcr.2017.09.01028890291 Mestdagh P Hartmann N Baeriswyl L Andreasen D Bernard N Chen C . Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. (2014) 11:809–15. 10.1038/nmeth.301424973947 Blondal T Brunetto MR Cavallone D Mikkelsen M Thorsen M Mang Y . Genome-wide comparison of next-generation sequencing and qPCR platforms for microRNA profiling in serum. Methods Mol Biol. (2017) 1580:21–44. 10.1007/978-1-4939-6866-4_328439824

Funding. This work was financially supported by the Ministry of Science and Higher Education of the Russian Federation (Agreement No. 075-15-2020-901).