- Open Access
Differential expression of microRNAs following cardiopulmonary bypass in children with congenital heart diseases
Journal of Translational Medicine volume 15, Article number: 117 (2017)
Children with congenital heart defects (CHDs) are at high risk for myocardial failure after operative procedures with cardiopulmonary bypass (CPB). Recent studies suggest that microRNAs (miRNA) are involved in the development of CHDs and myocardial failure. Therefore, the aim of this study was to determine alterations in the miRNA profile in heart tissue after cardiac surgery using CPB.
In total, 14 tissue samples from right atrium were collected from patients before and after connection of the CPB. SurePrint™ 8 × 60K Human v21 miRNA array and quantitative reverse transcription-polymerase chain reaction (RT-qPCR) were employed to determine the miRNA expression profile from three patients before and after connection of the CPB. Enrichment analyses of altered miRNA expression were predicted using bioinformatic tools.
According to miRNA array, a total of 90 miRNAs were significantly altered including 29 miRNAs with increased and 61 miRNAs with decreased expression after de-connection of CPB (n = 3) compared to before CPB (n = 3). Seven miRNAs had been validated using RT-qPCR in an independent cohort of 11 patients. Enrichment analyses applying the KEGG database displayed the highest correlation for signaling pathways, cellular community, cardiovascular disease and circulatory system.
Our result identified the overall changes of the miRNome in right atrium tissue of patients with CHDs after CPB. The differentially altered miRNAs lay a good foundation for further understanding of the molecular function of changed miRNAs in regulating CHDs and after CPB in particular.
Congenital heart diseases (CHDs) are the most common congenital organ malformations in newborns . Nearly 10 of 1000 newborns have congenital heart and vessel defects with a wide variation with regards to the severity of defects and the need for immediate therapeutic intervention [1,2,3]. Advances in diagnostic and therapeutic methods have improved the long-term survival and morbidities of these patients [3,4,5]. However, considerable mortality rates in neonates rather than older children in association with corrective cardiac surgery remain challenging [1,2,3,4,5]. Neonates needing surgical intervention in early life still have significantly higher mortality than infants and older children [6, 7]. Corrective cardiac surgery is one of the most major cause of mortality in neonates and infants . The morbidities and risk factors associated with cardiopulmonary bypass (CPB) are well documented in adults and children [9,10,11,12]. Data on the mechanisms of adaptation and remodeling in the myocardium of neonates and small children undergoing cardiac surgery are still scarce. Therefore, a better understanding of the molecular and genetic mechanisms underlying congenital heart diseases and the specific adapting mechanisms during stress and ischemia is fundamental. Recent studies in adults with different cardiac morbidities have suggested an important role for microRNAs (miRNAs) in the pathogenesis of numerous cardiovascular diseases [13, 14]. MiRNAs are a class of small non-coding RNAs that regulate gene expression post-transcriptionally via sequence-specific interaction with the 3′ UTR of target mRNAs, resulting in mRNA degradation and/or inhibition of translation . In signaling and transcriptional processes of cardiac biology, miRNAs can act as ‘molecular switches’ of post-transcriptional regulation of gene expression by regulating the expression of multiple proteins that function at different steps in cellular processes such as growth, differentiation, metabolism and apoptosis  as well as cardiac remodeling including increased angiogenesis and decreased fibrosis . These miRNAs-regulated pathways are contributing to the development of the cardiac remodeling process and their activation following a myocardial injury suggest a functional role for these signaling pathways in cardiac repair. MiRNAs known to be involved in myocardial development include miR-208a, miR-21, miR-145, miR-1, miR-133a and miR-29 [18,19,20]. Moreover, their aberrant expression has been linked to cardiac remodeling and hypertrophy [18,19,20]. However, the role of miRNAs in the myocardium under non-physiological conditions of CPB is not well understood. In this study, we aimed to characterize the miRNA profile in the atrial myocardium before and after cardiac surgery by means of CPB in patients (infants and children) with CHDs.
Study population and sample collection
Institutional Review Board approval (No. 156/14) was obtained before initiation of this study. Children aged from 5 days to 10.4 years who had cardiac surgery for CHDs were included in this pilot study after written informed consent was obtained. Selected demographics and preoperative characteristics of the children including diagnoses, surgical procedure and duration of hospitalization are shown in Table 1. Patients were excluded from the study when parents declined participation, complications occurred during surgery or scarring prevented secure extraction of a tissue sample. All patients received midazolam (0.1 mg/kg) for premedication. After preoxygenation with 100% oxygen, anesthesia was induced with thiopental and fentanyl in titrated doses. Atracurium (0.5 mg/kg) was given to facilitate nasotracheal intubation with a cuffed tube. Thereafter anesthesia was maintained with continuous infusions of fentanyl (5–10 µg/kg/h) and midazolam (0.1–0.2 mg/kg/h). The lungs were ventilated with oxygen/air to maintain normocapnia at pH of 7.35–7.45. Volume replacement was done to maintain physiological CVP and arterial pressure. Aortic and atrial cannulation were performed after systemic heparin administration with intravenous bolus. Before cannulation, the CPB was filled with priming solution which consisted of Sterofundin ISO, Gelafundin, Ringer’s lactate and tranexamic acid and then mixed with packed red blood cells (PRBCs) from donors. In every case the Heart–Lung Machine: JOSTRA HL 20 (Maquet, Rastatt, Germany) was used. Blood cardioplegia of 8–10 °C was given according to the Calafiore scheme . Atrial myocardial tissue samples from all patients were collected from the right auricular appendage or the right atrium before cannulation and after de-cannulation of the right atrial cannula. Tissue samples of interest were transferred to siliconized microcentrifuge tubes and immediately frozen in liquid nitrogen. These samples were stored in liquid nitrogen at −180 °C until processing.
Isolation of total RNA, including miRNAs, from myocardial tissue samples
Myocardial (right atrium) tissue samples were cut into small pieces, crushed and homogenized using the ball mill TissueLyser LT (Qiagen, Hilden, Germany) with an oscillation frequency of 50 Hz. Thereafter, total RNA including miRNA was extracted using the miRNeasy Mini kit on a QIAcube robot (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Concentration and purity of the samples were measured using the NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Massachusetts, USA). RNA integrity was assessed with an Agilent 2100 Bioanalyzer using RNA 6000 Nano kit (Agilent Technologies, California, USA). DNase I (Thermo Fisher Scientific, Massachusetts, USA) treatment was carried out according to the manufacturer’s instructions to remove any DNA contamination. Conventional PCR with exon spanning primers for Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was performed to exclude residual DNA in the samples .
MiRNA expression profiles of three samples from patients with CHD, before CPB (n = 3) and after CPB (n = 3) were hybridized on a SurePrint G3 miRNA Array v21.0 (Agilent Technologies, California, USA). The patients were chosen based on their age at the time of surgery to largely cover the range of ages included in our study. Briefly, 100 ng input of RNA from each sample was dephosphorylated by incubation with calf intestinal phosphatase and denatured using 100% DMSO. Samples were labeled with pCp-Cy3 by using T4 ligase. Each labeled RNA sample was hybridized onto an individual 8 × 60K format Agilent miRNA array slide v21. Arrays were washed and dried according to manufacturer’s recommendations and scanned at a resolution of 3 μm double-pass mode using an Agilent scanner. Data were acquired using Agilent AGW Feature Extraction software version 10.10.11 (Agilent Technologies, California, USA).
Reverse transcription qPCR of miRNA (array validation)
To validate the microarray results in the study, relative quantitative real-time PCR was performed on an ABI StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, USA) using SYBR Green I on nine differentially expressed miRNAs. Complementary DNA (cDNA) was generated by reverse transcription of 200 ng of total RNA using miScript RT II Kit (Qiagen, Hilden, Germany). Briefly, 200 ng of total RNA containing miRNAs was mixed with miScript 4 µL HiSpec Buffer, 2 µL nucleic mix, 2 µL miScript Reverse Transcriptase mix and RNase-free water to a final volume of 20 µL. Following the reverse transcription reaction, the cDNA was diluted 1:10 and then mixed with 10 µL QuantiTect SYBR Green PCR Master Mix, 2 µL miScript Universal Primer, 2 µL miScript Primer Assay for the selected nine miRNAs namely miR-328-5p, miR-4750-5p, miR-210-5p, miR-423-3p, miR-143-3p, miR-564, miR-770-5p, miR-874-5p, miR-93-5p and RNU6B as endogenous control and RNase-free water to a final volume of 20 µL. The primer sequences used in the study are shown in Additional file 1: Table S1. Reactions were assembled with the QIAgility automated pipetting system (Qiagen, Hilden, Germany). All PCR data were analyzed with SDS Relative Quantification Software version 2.3 (Applied Biosystems, Foster City, USA).
Target prediction and functional analysis
Enriched KEGG pathway analyses were performed using DIANA-miRPath v.3.0 software based on predicted targets by DIANA-microT-CDS . Targets of miRNAs with a score of more than 0.8 were selected. Only KEGG pathways with P value <0.05 and a false discovery rate (FDR) <0.05 were retained. The effect of miRNAs on target genes and networks has been evaluated using miRTargetLink software .
Reverse transcription qPCR of miRNA and mRNA (functional network validation)
We used miScript Primer Assays for 9 miRNAs (miR-744-5p, miR-648, miR-193b-3p, miR-212-3p, miR-143-3p, miR-93-5p, miR-222-3p, miR-423-3p and miR-766-3p) and QuantiTect Primer assays for 9 target genes (CDKN1A, MYC, PTEN, ESR1, ETS1, SOD2, MGMT, KRAS and HNF4A) (Qiagen, Hilden, Germany) to validate the different expression levels of the miRNA and their target genes, which are determined by miRTargetLink prediction software. In brief, 350 ng of total RNA including miRNA were converted into cDNA using the miScript II RT Kit. During the reverse transcription step, 5× miScript HiFlex Buffer was used to promote conversion of RNA into cDNA. The resulting cDNA was diluted to obtain a miRNA concentration of 1.5 ng/µL and a mRNA concentration of 5 ng/µL for. All reverse transcription PCR (RT-PCR) experiments were performed using the QIAgility™ automated PCR setup (Qiagen, Hilden, Germany). RT-qPCR analysis was done on a StepOnePlus™ Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). GAPDH and RNU6B were chosen as reference genes for mRNA and miRNA normalization, respectively. In addition, we included a no-template control (NTC) and no-reverse transcriptase control (NRT) in each run.
The freely available R software (R Development Core Team, 2010) was used to analyze the differences in miRNA expression in the atrial myocardial tissue samples from patients with CHD samples before and after CPB. After applying the Agilent Feature Extraction image analysis software on our hybridized microarray slides, we collected the computed total gene signals (TGS) for each miRNA, and performed quantile normalization and a log (base 2) transformation of the TGS values. Differential levels of miRNAs were analyzed by employing a paired two-tailed t-test for miRNAs that showed a significant change in the two groups. P-values below a threshold of 0.05 were considered statistically significant. We used the relative quantitative method of 2−ΔΔCq to measure the expression differences of specifically selected miRNAs in RT-qPCR . Two-tailed paired T-test was used to evaluate the fold change of miRNA levels.
A total of 14 patients who had cardiac surgery for CHDs were included in the present study. Preoperative characteristics of the children including diagnoses, surgical procedure and duration of hospitalization are shown in Table 1.
Differentially expressed miRNAs
In order to characterize miRNA expression profiles in the atrial myocardial tissue of patient after CPB, we performed microarray assays using SurePrint G3 miRNA arrays that contain 2549 human miRNAs annotated in miRBase version 21.0. Atrial myocardial tissue from patients with CHD before cardiac surgery by means of CPB and after CPB were collected (three patients in total). By applying a paired two-tailed t-test, microarray assays showed that miRNAs were expressed differentially in atrial myocardial tissue after surgical treatment with CPB. A total of 90 miRNAs were significantly altered in the two considered groups with a P value of <0.05. Out of the 90 altered miRNAs, 29 were up-regulated in the samples after CPB compared to samples before CPB, while 61 miRNA were down-regulated (Table 2). As shown in Table 2, the majority of altered miRNA (49 out of 90) fell into the range of 1.50–1.99 fold up- or down-regulation. In addition, 17 miRNAs including six up-regulated and another eleven down-regulated miRNAs displayed expression level with changes ≥2.0-fold and 24 miRNA including eight up-regulated and another sixteen down-regulated miRNAs displayed expression level with changes <1.5 fold between two groups. Using hierarchic clustering with the euclidian distance measure, we analyzed how the expression profiles of patients before and after CPB relate to each other. For this task, we used the 50 miRNAs with the highest variance of miRNA levels out of the 2.549 miRNAs. Figure 1 shows the heatmap of the hierarchic clustering from the differentially expressed miRNA by microarray. In general, the hierarchical clustering exhibited a clear separation of the examined groups based on miRNA expression profiles. The first cluster contains mostly patients before CPB and the second most of the patients after CPB.
Validation of differentially expressed miRNAs using RT-qPCR
To validate the data obtained from the miRNA microarray, RT-qPCR was performed to re-examine the expression level of nine miRNAs, namely miR-328-5p, miR-4750-5p, miR-210-5p, miR-423-3p, miR-143-3p, miR-564, miR-770-5p, miR-874-5p and miR-93-5p. These miRNAs were chosen based on their differential expression level in each patient group before and after CPB and/or on their known associations with heart diseases. In detail, we selected two miRNAs with highest (miR-4750-5p) and moderate (miR-328-5p) fold change among the up-regulated miRNAs and two miRNAs with highest (miR-770-5p) and moderate (miR-564) fold change among the down-regulated miRNAs. In addition, we selected six miRNAs (miR-328-5p, miR-210-5p, miR-423-3p, miR-143-3p, miR-874-5p and miR-93-5p) with low or moderate expression levels based on the array analysis and with known association with cardiac pathologies and its related process [26,27,28,29,30,31,32,33]. For this validation step, we analyzed a new cohort of 11 samples taken from patients with CHD undergoing corrective surgery by CPB. The RT-qPCR showed the same direction of expression changes as the microarray analysis for seven miRNAs (Fig. 2). The significance of the differences in the expression was confirmed for seven miRNAs namely miR-210-5p, miR-423-3p, miR-143-3p, miR-564, miR-770-5p, miR-874-5p and miR-93-5p) (P < 0.05) in the atrial myocardial samples from patients with CHDs after CPB compared to before CPB by RT-qPCR (Fig. 3). By including two out of the seven validated miRNAs namely miR-143-3p and miR-93-5p, a clear distinction between the groups based on the clustering dendrogram was, however, not possible (Additional file 2: Figure S1). We next analyzed the expression levels of the nine miRNAs before CPB and after CPB separately for neonate and infant patients. As for the atrial myocardial samples from neonates with CHDs, we found a significantly different expression for miR-874-5p with a P-value of 0.040 between prior CPB and after CPB. In addition, we found two miRNAs with border-line P-values (miR-423-3p, P-value = 0.060 and miR-93-5p, P-value = 0.070) (Fig. 4). As for atrial myocardial samples from infants with CHDs, we found significantly different expression of six miRNAs (miR-210-5p with a P-value of 0.014, miR-423-3p with a P-value of 0.034, miR-143-3p with a P-value of 0.018, miR-564 with a P-value of 0.024, miR-770-5p with a P-value of 0.025, and miR-874-5p with a P-value of 0.040) between prior and after CPB (Fig. 4). In addition, we analyzed the relative expression level of different miRNAs depending on the age of each patient before and after CPB. As shown in Additional file 3: Figure S2, the relative expression level (2−ΔΔCt) for certain miRNAs decreased with age. In addition, we found a remarkable difference between the miRNAs analyzed.
Comparative pathway analysis
We used DIANA-mirPath algorithm to gain insights into the biological pathways of the miRNAs that were altered in the atrial myocardial tissue from patients with CHD after CPB compared to before CPB. Based on the deregulated miRNAs by microarray analysis, we identified seventy KEGG pathways that were significantly enriched (P < 0.05, FDR corrected) for targets of the deregulated miRNAs, respectively (Additional file 1: Table S2). As shown in Table 3, the target genes of the altered miRNAs are mostly involved in “Signaling transduction”, “Cellular community”, “Cardiovascular disease”, “Circulatory system” and many others pathways. Using miRTargetLink, the resulting network is drawn schematically in Fig. 5. Besides the nine miRNAs, nine genes were targeted by one, two or more of the selected miRNAs in the “Strong” category. A strong interaction was observed between miRNAs and genes are highlighted in ‘‘Green’’ in the resulting network (Fig. 5).
Validation of miRNAs and mRNAs that were differentially expressed in atrial myocardial tissues
To confirm the differential expression of miRNAs and their targeted mRNAs before and after CPB, we analyzed eight of the down-regulated miRNAs, one of the up-regulated miRNAs and nine target mRNAs by RT-qPCR. With the exception of miR-222-3p, the remaining eight miRNAs showed the same direction of expression changes in the RT-PCR and in the microarray analysis. These findings are summarized in Additional file 4: Figure S3. Out of the eight miRNAs, miR-648 was up-regulated and the seven remaining miRNAs were down regulated (miR-744-5p, miR-193b-3p, miR-212-3p, miR-143-3p, miR-93-5p, miR-423-3p and miR-766-3p). As for the mRNA analysis, 6 out of 9 analyzed target mRNAs showed the expected expression changes. In detail, for five down-regulated miRNAs the according mRNA targets (ESR1, ETS1, MYC, PTEN and CDKN1A) were up-regulated. For the up-regulated miRNA miR-222-3p the according mRNA target (MGMT) was down-regulated. Due to low number of samples, we could not determine the significance of the expression differences. For the targets SOD2 and KRAS we did not find expression changes between prior CPB and after CPB. For the target HNF4A, we could not determine the fold changes.
In our study, we found evidence for an overall altered miRNA expression pattern in the atrial myocardial tissue of patients with CHD according to corrective cardiac surgery by using the cardiopulmonary bypass (CPB). MiRNA expression analysis using microarrays indicated 90 miRNAs with significantly differential expression, including 29 miRNAs with significantly up- and 61 miRNAs with significantly down-regulated expression after CPB. In the validation phase, by using a new cohort of samples, seven miRNAs were validated. These data show that miRNA expression levels in atrial myocardial tissue changes in the course of cardiac surgery using the CPB. MiRNAs specifically expressed or enriched in smooth, skeletal, and/or cardiac muscle that plays a physiological role in normal heart development and function are termed Myo-miRs . In particular, miR-208, a member of the miR-208 family, is differentially expressed during heart developmental and plays an essential role in normal cardiac conduction [34, 35]. Both over-expression and under-expression of miR-208a has been associated with cardiac arrhythmia , with apoptosis in ischemic cardiomyocytes  and generally with heart diseases . Down-regulation of miR-208a in atrial myocardial tissues of patients with CHD after CPB underscore the integral function of this miRNA in regulating cell morphology and contractility. The cardiomyogenesis regulator miR-143 is highly expressed during vascular smooth muscle cell differentiation and is involved in normal function and formation of the cardiac chamber via regulation of myocardial cell morphology . The down-regulation of miR-143 causes ventricular collapse by affecting cell morphology and contractility . Furthermore, miR-143 targets ETS Transcription Factor (Elk1), an activator of vascular smooth muscle cells proliferation .
Deregulation of a number of miRNAs in our study have been previously reported in right and left atrial appendages of patients with rheumatic mitral valve disease (RMVD), including miR-222-3p, miR-4484 and miR-940 in left atrial appendage , miR-5190 and miR-23c in right atrial appendages  and miR-143-3p in both, in the LAA and RAA of patients with RMVD . MiR-222 has a documented function in regulating cell proliferation and is involved in the vascular smooth muscle cells differentiation [38, 44]. In neonatal cardiomyocytes, miR-222 induces cellular hypertrophy and proliferation and inhibits apoptosis after ischemic injury . In addition, miR-222 increases cell proliferation and inhibits cardiomyogenic differentiation in right ventricular outflow tract (RVOT) myocardial tissues from infants with nonsyndromic tetralogy of fallot (TOF) . Our finding for miR-222 also suggests that this miRNA protects cardiac structure and functions after corrective surgery by CPB. Previous studies have identified that many miRNAs are expressed at low levels under normal condition and expressed strongly during pathological stress like miR-212, a cardiomyocyte-specific miRNA, which is strongly activated during heart failure. Its up-regulation expression leads to pathological cardiac hypertrophy . In agreement with our expression direction, Zhou et al.  reported down-regulation of miR-423-3p and miR-532-3p in ischemia–reperfusion injury heart grafts. MiR-423-5p was suggested as a diagnostic biomarker to distinguish heart failure patients from healthy controls. The plasma concentration patients with heart failure may reflect the severity of dilated cardiomyopathy . Similarly, reduced expression of miR-93 was observed in the left atrium compared to right atrial tissue of patients with sinus rhythm  indicating that miR-93 promotes angiogenesis in the ischemic tissue by coordinating the functional pathways of cell proliferation and apoptosis . The up-regulation of miR-328 expression suppresses its target gene SERCA2a (ATPase Sarcoplasmic/Endoplasmic Reticulum Ca2+ Transporting 2) in cardiac myocytes to indirectly activate the calcineurin/NFATc3 signaling pathway, leading to cardiac hypertrophy . In addition, miR-328 was up-regulated in the atria of patients with atrial fibrillation (AF) . Two miRNAs, miR-624 and miR-339, were deregulated in ischemic heart disease and coronary artery disease [53, 54]. MiR-744 was deregulated in patients with chronic congestive heart failure . Several miRNAs were detected in cardiac tissue at different stages of development and are highly expressed in the fetal heart, including miR-212, miR-210 and miR-423 . Moreover, Huang et al. found that miR-210 is involved in cellular hypoxia, regulation of angiogenesis and apoptosis . Similarly, miR-423-5p, miR-193b-3p, and miR-550a-5p showed an altered expression level in patients with heart failure [56, 57].
Bioinformatics analysis by DIANA-mirPath predicted several KEGG biological pathways that were significantly enriched for the differentially deregulated miRNAs in the atrial myocardial samples (Table 3). Enrichment analysis displayed the highest correlation for signaling pathways, providing further evidence that myocardial target proteins are involved in the signaling pathways. Pathways related to the cardiovascular system and vascular smooth muscle contraction and its related processes, including cardiomyocytes proliferation, differentiation and apoptosis have been identified [58,59,60,61,62,63]. Recent studies implicating the Hippo signaling pathway promotes cardiomyocyte proliferation by activating the insulin-like growth factor  and Wnt signaling pathways . Similarly, ErbB signaling pathway plays an important role in proper heart morphogenesis and also in all developmental stages of the heart . These results demonstrate that the deregulated miRNAs in atrial myocardial tissue play an important role in participating in many signaling pathways that control heart function. We employed bioinformatics tools to gain further insights into the impact of the deregulated miRNAs on target genes. As indicated in Fig. 5 and Additional file 4: Figure S3 specific genes are likely affected by the several deregulated miRNAs including the genes for O-6-methylguanine-DNA methyltransferase (MGMT) and superoxide dismutase 2 (SOD2). The expression of these genes is enhanced in myocardial remodeling, cardiac hypertrophy and/or failure indicating that severe stress effects plays a critical role in in the pathogenesis of myocardial remodeling and failure [67, 68]. Phosphatase and tensin homolog (PTEN) plays a role in promoting cardiomyocyte proliferation and regeneration and protect the heart from hypertrophy and heart failure under biomechanical stress [69,70,71]. The transcription factors v-Myc Avin myelocytomatosis viral oncogene homologue (c-Myc) also play a regulatory role in stimulating cardiac myocytes proliferation and differentiation during fetal development .
We would like to point out that this study has also a number of limitations. One important point is that the study population was small and heterogeneous; age at surgery ranged from 5 days of life to 10.4 years and we did not differ between male and female gender. Moreover, the examined CHDs reached from simple ASDs to complex CHDs like HLHS or TAPVC and complicated the interpretation of the data. In fact, we were able to find several miRNAs which were elevated in all patients despite the heterogeneous study population. Another limitation is that all samples were taken from the right atrium. Acquisition of the samples happened before cannulation and after de-cannulation. However, for example, it is not clear which impact the suture before excision of the sample on the miRNA expression has. A sampling of a myocardial specimen from the right or left ventricle before the connection to CPB is not possible. This can be done during the corrective surgery after the cardiac arrest and may not represent the physiological conditions before the connection to CPB and cardiac arrest. Thus, the tissue sampling from the atrial appendage may facilitate the examining of the myocardial conditions before and after surgery using the CPB without increasing risks. In addition, this side of the heart may represent nearly a uniform myocardial region in all patients with CHD. Other settings of myocardial sampling will not be approved by the ethic committee of our institution. The precise mechanisms underlying the expressions of such miRNA should be evaluated in a larger cohort of patients including subgroups of uniform diagnosis of congenital heart defects. Nevertheless, we conclude that the connection to CPB and hypothermic cardiac arrest may induce several physiological responses in the myocardium, which in part is reflected by changed miRNA expression levels. Furthermore, it is unclear which effect the cardiac surgery itself as well as the usage of the CPB has on the miRNA expression. Several non-physiological circumstances occur during this process; for example, cardioplegia induces cardiac arrest, hypothermia is used, the heart is mechanical unloaded and reperfusion occurs after CPB. In addition, the CPB means contact of the blood with a foreign body surface. Whether this circumstance alone already causes changes in the miRNA profile has to be cleared.
In summary, we report miRNAs in the atrial myocardial tissue with significantly altered expression levels in patients with CHD after cardiac surgery with CPB. These altered miRNAs include cardiac-specific miRNAs which have been described in various cardiac pathologies including congenital heart disease. The overlap between the miRNAs identified in our study and the miRNAs that are involved in various levels of cardiac development by either up- or down-regulation strongly supports the idea that these miRNAs play an essential role in the pathophysiology of congenital heart diseases. Although it is important to bear in mind, that the miRNA expression profiles were obtained from a small number of atrial myocardial tissue samples, the altered expression levels have been further confirmed RT-qPCR for nine deregulated in a cohort of 11 independent atrial myocardial tissue samples. The alteration of the expression of miRNAs may provide new insights into the underlying mechanisms of cardiac surgery with CPB and provide potential novel mechanism-based therapeutic strategies for CHD.
atrial septal defect
cyclin dependent kinase inhibitor 1A
congenital heart defect
v-Myc avin myelocytomastosis viral oncogene homologue
dextro-position of the great arteries
eTS transcription factor
estrogen receptor 1
eTS proto-oncogene 1
false discovery rate
glyceraldehyde 3-phosphate dehydrogenase
hypoplastic left heart syndrome
hepatocyte nuclear factor 4 alpha
Kyoto encyclopedia of genes and genomes
KRAS proto-oncogene, GTPase
left atrial appendages
miRNA enrichment analysis and annotation
packed red blood cells
phosphatase and tensin homolog
right atrial appendages
rheumatic mitral valve disease
quantitative reverse transcription-polymerase chain reaction
right ventricular outflow tract
ATPase sarcoplasmic/endoplasmic reticulum Ca2 + transporting 2
superoxide dismutase 2
total anomalous pulmonary venous connection
total gene signals
tetralogy of fallot
ventricular septal defect
van der Linde D, Konings EE, Slager MA, Witsenburg M, Helbing WA, Takkenberg JJ, Roos-Hesselink JW. Birth prevalence of congenital heart disease worldwide: a systematic review and meta-analysis. J Am Coll Cardiol. 2011;58(21):2241–7.
Schwedler G, Lindinger A, Lange PE, Sax U, Olchvary J, Peters B, Bauer U, Hense HW. Frequency and spectrum of congenital heart defects among live births in Germany: a study of the competence network for congenital heart defects. Clin Res Cardiol. 2011;100(12):1111–7.
Erikssen G, Liestol K, Seem E, Birkeland S, Saatvedt KJ, Hoel TN, Dohlen G, Skulstad H, Svennevig JL, Thaulow E, et al. Achievements in congenital heart defect surgery: a prospective, 40-year study of 7038 patients. Circulation. 2015;131(4):337–46 (discussion 346).
Raissadati A, Nieminen H, Jokinen E, Sairanen H. Progress in late results among pediatric cardiac surgery patients: a population-based 6-decade study with 98% follow-up. Circulation. 2015;131(4):347–53 (discussion 353).
Jones B, Muscara F, Lloyd O, McKinlay L, Justo R. Neurodevelopmental outcome following open heart surgery in infancy: 6-year follow-up. Cardiol Young. 2015;25(5):903–10.
Dollat C, Vergnat M, Laux D, Stos B, Baruteau A, Capderou A, Demontoux S, Hamann M, Mokhfi E, Van Aerschot I, et al. Critical congenital heart diseases in preterm neonates: is early cardiac surgery quite reasonable? Pediatr Cardiol. 2015;36(6):1279–86.
Kumar TK, Charpie JR, Ohye RG, Hirsch-Romano JC, Donohue JE, Yu S, Sood V, Wilkinson DA, Nelson K, Mitchell E, et al. Timing of neonatal cardiac surgery is not associated with perioperative outcomes. J Thorac Cardiovasc Surg. 2014;147(5):1573–9.
Ades AM, Dominguez TE, Nicolson SC, Gaynor JW, Spray TL, Wernovsky G, Tabbutt S. Morbidity and mortality after surgery for congenital cardiac disease in the infant born with low weight. Cardiol Young. 2010;20(1):8–17.
Cheng HH, Wypij D, Laussen PC, Bellinger DC, Stopp CD, Soul JS, Newburger JW, Kussman BD. Cerebral blood flow velocity and neurodevelopmental outcome in infants undergoing surgery for congenital heart disease. Ann Thorac Surg. 2014;98(1):125–32.
Kang N, Cole T, Tsang V, Elliott M, de Leval M. Risk stratification in paediatric open-heart surgery. Eur J Cardiothorac Surg. 2004;26(1):3–11.
Kussman BD, Wypij D, Laussen PC, Soul JS, Bellinger DC, DiNardo JA, Robertson R, Pigula FA, Jonas RA, Newburger JW. Relationship of intraoperative cerebral oxygen saturation to neurodevelopmental outcome and brain magnetic resonance imaging at 1 year of age in infants undergoing biventricular repair. Circulation. 2010;122(3):245–54.
Wypij D, Newburger JW, Rappaport LA, duPlessis AJ, Jonas RA, Wernovsky G, Lin M, Bellinger DC. The effect of duration of deep hypothermic circulatory arrest in infant heart surgery on late neurodevelopment: the Boston Circulatory Arrest Trial. J Thorac Cardiovasc Surg. 2003;126(5):1397–403.
Creemers EE, Tijsen AJ, Pinto YM. Circulating microRNAs: novel biomarkers and extracellular communicators in cardiovascular disease? Circ Res. 2012;110(3):483–95.
Mishra PK, Tyagi N, Kumar M, Tyagi SC. MicroRNAs as a therapeutic target for cardiovascular diseases. J Cell Mol Med. 2009;13(4):778–89.
Pasquinelli AE. MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship. Nat Rev Genet. 2012;13(4):271–82.
Liu N, Olson EN. MicroRNA regulatory networks in cardiovascular development. Dev Cell. 2010;18(4):510–25.
Agarwal U, George A, Bhutani S, Ghosh-Choudhary S, Maxwell JT, Brown ME, Mehta Y, Platt MO, Liang Y, Sahoo S, et al. Experimental, systems, and computational approaches to understanding the microRNA-mediated reparative potential of cardiac progenitor cell-derived exosomes From pediatric patients. Circ Res. 2017;120(4):701–12.
van Rooij E, Sutherland LB, Thatcher JE, DiMaio JM, Naseem RH, Marshall WS, Hill JA, Olson EN. Dysregulation of microRNAs after myocardial infarction reveals a role of miR-29 in cardiac fibrosis. Proc Natl Acad Sci USA. 2008;105(35):13027–32.
Thum T, Gross C, Fiedler J, Fischer T, Kissler S, Bussen M, Galuppo P, Just S, Rottbauer W, Frantz S, et al. MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts. Nature. 2008;456(7224):980–4.
Wang GK, Zhu JQ, Zhang JT, Li Q, Li Y, He J, Qin YW, Jing Q. Circulating microRNA: a novel potential biomarker for early diagnosis of acute myocardial infarction in humans. Eur Heart J. 2010;31(6):659–66.
Calafiore AM, Teodori G, Mezzetti A, Bosco G, Verna AM, Di Giammarco G, Lapenna D. Intermittent antegrade warm blood cardioplegia. Ann Thorac Surg. 1995;59(2):398–402.
Abu-Halima M, Ludwig N, Hart M, Leidinger P, Backes C, Keller A, Hammadeh M, Meese E. Altered micro-ribonucleic acid expression profiles of extracellular microvesicles in the seminal plasma of patients with oligoasthenozoospermia. Fertil Steril. 2016;106:1061–9.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, Dalamagas T, Hatzigeorgiou AG. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43(W1):W460–6.
Hamberg M, Backes C, Fehlmann T, Hart M, Meder B, Meese E, Keller A. MiRTargetLink–miRNAs, genes and interaction networks. Int J Mol Sci. 2016;17(4):564.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods. 2001;25(4):402–8.
Contu R, Latronico MV, Condorelli G. Circulating microRNAs as potential biomarkers of coronary artery disease: a promise to be fulfilled? Circ Res. 2010;107(5):573–4.
Dickinson BA, Semus HM, Montgomery RL, Stack C, Latimer PA, Lewton SM, Lynch JM, Hullinger TG, Seto AG, van Rooij E. Plasma microRNAs serve as biomarkers of therapeutic efficacy and disease progression in hypertension-induced heart failure. Eur J Heart Fail. 2013;15(6):650–9.
Endo K, Naito Y, Ji X, Nakanishi M, Noguchi T, Goto Y, Nonogi H, Ma X, Weng H, Hirokawa G, et al. MicroRNA 210 as a biomarker for congestive heart failure. Biol Pharm Bull. 2013;36(1):48–54.
Goren Y, Kushnir M, Zafrir B, Tabak S, Lewis BS, Amir O. Serum levels of microRNAs in patients with heart failure. Eur J Heart Fail. 2012;14(2):147–54.
Latronico MV, Catalucci D, Condorelli G. MicroRNA and cardiac pathologies. Physiol Genom. 2008;34(3):239–42.
Tijsen AJ, Creemers EE, Moerland PD, de Windt LJ, van der Wal AC, Kok WE, Pinto YM. MiR423-5p as a circulating biomarker for heart failure. Circ Res. 2010;106(6):1035–9.
Wang K, Liu F, Zhou LY, Ding SL, Long B, Liu CY, Sun T, Fan YY, Sun L, Li PF. miR-874 regulates myocardial necrosis by targeting caspase-8. Cell Death Dis. 2013;4:e709.
Li C, Li X, Gao X, Zhang R, Zhang Y, Liang H, Xu C, Du W, Zhang Y, Liu X, et al. MicroRNA-328 as a regulator of cardiac hypertrophy. Int J Cardiol. 2014;173(2):268–76.
van Rooij E, Quiat D, Johnson BA, Sutherland LB, Qi X, Richardson JA, Kelm RJ Jr, Olson EN. A family of microRNAs encoded by myosin genes governs myosin expression and muscle performance. Dev Cell. 2009;17(5):662–73.
Sayed D, Abdellatif M. MicroRNAs in development and disease. Physiol Rev. 2011;91(3):827–87.
Callis TE, Pandya K, Seok HY, Tang RH, Tatsuguchi M, Huang ZP, Chen JF, Deng Z, Gunn B, Shumate J, et al. MicroRNA-208a is a regulator of cardiac hypertrophy and conduction in mice. J Clin Investig. 2009;119(9):2772–86.
Tony H, Meng K, Wu B, Yu A, Zeng Q, Yu K, Zhong Y. MicroRNA-208a dysregulates apoptosis genes expression and promotes cardiomyocyte apoptosis during ischemia and its silencing improves cardiac function after myocardial infarction. Med Inflamm. 2015;2015:479123.
Winbanks CE, Ooi JY, Nguyen SS, McMullen JR, Bernardo BC. MicroRNAs differentially regulated in cardiac and skeletal muscle in health and disease: potential drug targets? Clin Exp Pharmacol Physiol. 2014;41(9):727–37.
Deacon DC, Nevis KR, Cashman TJ, Zhou Y, Zhao L, Washko D, Guner-Ataman B, Burns CG, Burns CE. The miR-143-adducin3 pathway is essential for cardiac chamber morphogenesis. Development. 2010;137(11):1887–96.
Cordes KR, Sheehy NT, White MP, Berry EC, Morton SU, Muth AN, Lee TH, Miano JM, Ivey KN, Srivastava D. miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature. 2009;460(7256):705–10.
Thum T, Galuppo P, Wolf C, Fiedler J, Kneitz S, van Laake LW, Doevendans PA, Mummery CL, Borlak J, Haverich A, et al. MicroRNAs in the human heart: a clue to fetal gene reprogramming in heart failure. Circulation. 2007;116(3):258–67.
Yu Z, Raabe T, Hecht NB. MicroRNA Mirn122a reduces expression of the posttranscriptionally regulated germ cell transition protein 2 (Tnp2) messenger RNA (mRNA) by mRNA cleavage. Biol Reprod. 2005;73(3):427–33.
Liu H, Qin H, Chen GX, Liang MY, Rong J, Yao JP, Wu ZK. Comparative expression profiles of microRNA in left and right atrial appendages from patients with rheumatic mitral valve disease exhibiting sinus rhythm or atrial fibrillation. J Transl Med. 2014;12:90.
Liu X, Cheng Y, Zhang S, Lin Y, Yang J, Zhang C. A necessary role of miR-221 and miR-222 in vascular smooth muscle cell proliferation and neointimal hyperplasia. Circ Res. 2009;104(4):476–87.
Liu X, Xiao J, Zhu H, Wei X, Platt C, Damilano F, Xiao C, Bezzerides V, Bostrom P, Che L, et al. miR-222 is necessary for exercise-induced cardiac growth and protects against pathological cardiac remodeling. Cell Metab. 2015;21(4):584–95.
Zhang J, Chang JJ, Xu F, Ma XJ, Wu Y, Li WC, Wang HJ, Huang GY, Ma D. MicroRNA deregulation in right ventricular outflow tract myocardium in nonsyndromic tetralogy of fallot. Can J Cardiol. 2013;29(12):1695–703.
Ucar A, Gupta SK, Fiedler J, Erikci E, Kardasinski M, Batkai S, Dangwal S, Kumarswamy R, Bang C, Holzmann A, et al. The miRNA-212/132 family regulates both cardiac hypertrophy and cardiomyocyte autophagy. Nat Commun. 1078;2012:3.
Zhou L, Zang G, Zhang G, Wang H, Zhang X, Johnston N, Min W, Luke P, Jevnikar A, Haig A, et al. MicroRNA and mRNA signatures in ischemia reperfusion injury in heart transplantation. PLoS ONE. 2013;8(11):e79805.
Fan KL, Zhang HF, Shen J, Zhang Q, Li XL. Circulating microRNAs levels in Chinese heart failure patients caused by dilated cardiomyopathy. Indian Heart J. 2013;65(1):12–6.
Slagsvold KH, Johnsen AB, Rognmo O, Hoydal M, Wisloff U, Wahba A. Comparison of left versus right atrial myocardium in patients with sinus rhythm or atrial fibrillation-an assessment of mitochondrial function and microRNA expression. Physiol Rep. 2014;2(8):e12124.
Hazarika S, Farber CR, Dokun AO, Pitsillides AN, Wang T, Lye RJ, Annex BH. MicroRNA-93 controls perfusion recovery after hindlimb ischemia by modulating expression of multiple genes in the cell cycle pathway. Circulation. 2013;127(17):1818–28.
McManus DD, Lin H, Tanriverdi K, Quercio M, Yin X, Larson MG, Ellinor PT, Levy D, Freedman JE, Benjamin EJ. Relations between circulating microRNAs and atrial fibrillation: data from the Framingham Offspring Study. Heart Rhythm. 2014;11(4):663–9.
Chen F, Zhao X, Peng J, Bo L, Fan B, Ma D. Integrated microRNA-mRNA analysis of coronary artery disease. Mol Biol Rep. 2014;41(8):5505–11.
Saddic LA, Chang TW, Sigurdsson MI, Heydarpour M, Raby BA, Shernan SK, Aranki SF, Body SC, Muehlschlegel JD. Integrated microRNA and mRNA responses to acute human left ventricular ischemia. Physiol Genom. 2015;47(10):455–62.
Cakmak HA, Coskunpinar E, Ikitimur B, Barman HA, Karadag B, Tiryakioglu NO, Kahraman K, Vural VA. The prognostic value of circulating microRNAs in heart failure: preliminary results from a genome-wide expression study. J Cardiovasc Med. 2015;16(6):431–7.
Huang X, Zuo J. Emerging roles of miR-210 and other non-coding RNAs in the hypoxic response. Acta Biochim Biophys Sin. 2014;46(3):220–32.
Wong LL, Armugam A, Sepramaniam S, Karolina DS, Lim KY, Lim JY, Chong JP, Ng JY, Chen YT, Chan MM, et al. Circulating microRNAs in heart failure with reduced and preserved left ventricular ejection fraction. Eur J Heart Fail. 2015;17(4):393–404.
Bruneau BG. Signaling and transcriptional networks in heart development and regeneration. Cold Spring Harbor Perspect Biol. 2013;5(3):a008292.
Galasso G, De Rosa R, Piscione F, Iaccarino G, Vosa C, Sorriento D, Piccolo R, Rapacciuolo A, Walsh K, Chiariello M. Myocardial expression of FOXO3a-Atrogin-1 pathway in human heart failure. Eur J Heart Fail. 2010;12(12):1290–6.
Holweg CT, Baan CC, Niesters HG, Vantrimpont PJ, Mulder PG, Maat AP, Weimar W, Balk AH. TGF-beta1 gene polymorphisms in patients with end-stage heart failure. J Heart Lung Transpl. 2001;20(9):979–84.
Tsai EJ, Kass DA. Cyclic GMP signaling in cardiovascular pathophysiology and therapeutics. Pharmacol Ther. 2009;122(3):216–38.
Xin M, Kim Y, Sutherland LB, Murakami M, Qi X, McAnally J, Porrello ER, Mahmoud AI, Tan W, Shelton JM, et al. Hippo pathway effector Yap promotes cardiac regeneration. Proc Natl Acad Sci USA. 2013;110(34):13839–44.
Communal C, Colucci WS. The control of cardiomyocyte apoptosis via the beta-adrenergic signaling pathways. Arch Mal Coeur Vaiss. 2005;98(3):236–41.
Heallen T, Zhang M, Wang J, Bonilla-Claudio M, Klysik E, Johnson RL, Martin JF. Hippo pathway inhibits Wnt signaling to restrain cardiomyocyte proliferation and heart size. Science. 2011;332(6028):458–61.
Xin M, Kim Y, Sutherland LB, Qi X, McAnally J, Schwartz RJ, Richardson JA, Bassel-Duby R, Olson EN. Regulation of insulin-like growth factor signaling by Yap governs cardiomyocyte proliferation and embryonic heart size. Sci Signal. 2011;4(196):ra70.
Sanchez-Soria P, Camenisch TD. ErbB signaling in cardiac development and disease. Semin Cell Dev Biol. 2010;21(9):929–35.
Fritz G, Kaina B. Stress factors affecting expression of O6-methylguanine-DNA methyltransferase mRNA in rat hepatoma cells. Biochem Biophys Acta. 1992;1171(1):35–40.
Siwik DA, Tzortzis JD, Pimental DR, Chang DL, Pagano PJ, Singh K, Sawyer DB, Colucci WS. Inhibition of copper-zinc superoxide dismutase induces cell growth, hypertrophic phenotype, and apoptosis in neonatal rat cardiac myocytes in vitro. Circ Res. 1999;85(2):147–53.
Chen J, Huang ZP, Seok HY, Ding J, Kataoka M, Zhang Z, Hu X, Wang G, Lin Z, Wang S, et al. mir-17-92 cluster is required for and sufficient to induce cardiomyocyte proliferation in postnatal and adult hearts. Circ Res. 2013;112(12):1557–66.
Hassink RJ, Pasumarthi KB, Nakajima H, Rubart M, Soonpaa MH, de la Riviere AB, Doevendans PA, Field LJ. Cardiomyocyte cell cycle activation improves cardiac function after myocardial infarction. Cardiovasc Res. 2008;78(1):18–25.
Sun F, Park KK, Belin S, Wang D, Lu T, Chen G, Zhang K, Yeung C, Feng G, Yankner BA, et al. Sustained axon regeneration induced by co-deletion of PTEN and SOCS3. Nature. 2011;480(7377):372–5.
Jackson T, Allard MF, Sreenan CM, Doss LK, Bishop SP, Swain JL. The c-myc proto-oncogene regulates cardiac development in transgenic mice. Mol Cell Biol. 1990;10(7):3709–16.
MA performed experimental work, particularly the miRNA array experiment, RT-qPCR validation and wrote the manuscript, MP, helped in the study design, diagnosed patients and helped in writing, NL, helped in RT-qPCR validation experiments and edited the manuscript. JM, helped in experimental work and samples collection. IM, isolated RNA, assessed the quality and quantity of RNA, GG, patients’ diagnosis and samples collection, JP, patients’ diagnosis and samples collection, HJS, patients’ diagnosis and samples collection, UG, patients’ diagnosis and samples collection, TP, patients’ diagnosis and samples collection, AK, performed bioinformatics analysis, EM, designed the study, coordinated the molecular biology experiment and article editing, HAK, designed the study, diagnosed patients and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
Consent for publication
Consent has been provided, the information letter and consent form are available in additional file 5.
Ethics approval and consent to participate
Institutional Review Board approval/Ethikvotum Ärztekammer des Saarlandes: Nr. 156/14.
This study was funded by the Fördergemeinschaft Kinderherzen in Bonn, the Competence Network for Congenital Heart Defects, which received funding from the Federal Ministry of Education and Research, Grant No. 01GI0601 (2014), and the German Centre for Cardiovascular Research (DZHK), Grant No. 81X2800112 (2015).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
RT-qPCR primer sequences used in the study. Table S2. MiRNA associated pathways in the atrial myocardial tissue of patients with CHD after CPB (n=3) compared to before CPB (n=3): KEGG Pathways with predicted interaction enrichments for deregulated miRNAs by microarray (P value < 0.05).
Unsupervised hierarchical clustering (Euclidian distance, complete linkage) of 3 patients with 6 samples including 3 samples before and 3 samples after CPB. The clustering was done based on the expression of the highest variance miRNAs as determined by RT-qPCR. MiRNAs with high expression are shown in blue and miRNAs with low expression in red. The blue and orange lines indicate the two main clusters of samples.
Influence of age on the expression of 9 differentially expressed miRNA as determined by RT-qPCR. The relative expression level of 2-ΔΔCt is indicated on the y-axis and the patient age at the time of surgery on the x-axis. The different miRNAs are color coded as indicated. RNAU6B was used as endogenous control for normalization and paired-two-tailed t-tests was used to evaluate differences in expression.
Functional target network validation of 9 differentially expressed miRNAs and 9 predicted target mRNAs. The comparison was done between atrial myocardial tissue of CHD patients before CBP (n=10) and after CPB (n=10). RNAU6B and GAPDH were used as an endogenous controls for normalization of miRNA and mRNA, respectively. Paired-two-tailed t-tests and ± standard deviation (STDV) were used to evaluate differences in expression.
1) Information letter to patients and parents, 2) consent form for anonymous publication, participation & storage of medical data and 3) tissue sample form.
About this article
Cite this article
Abu-Halima, M., Poryo, M., Ludwig, N. et al. Differential expression of microRNAs following cardiopulmonary bypass in children with congenital heart diseases. J Transl Med 15, 117 (2017). https://doi.org/10.1186/s12967-017-1213-9
- Congenital heart disease
- Atrial myocardium
- Cardiopulmonary bypass