Next Article in Journal
Recent Insights into the Rhythmogenic Core of the Locomotor CPG
Next Article in Special Issue
Ivabradine Induces Cardiac Protection against Myocardial Infarction by Preventing Cyclophilin-A Secretion in Pigs under Coronary Ischemia/Reperfusion
Previous Article in Journal
Expression Signature of lncRNAs and mRNAs in Sevoflurane-Induced Mouse Brain Injury: Implication of Involvement of Wide Molecular Networks and Pathways
Previous Article in Special Issue
CD40/CD40L and Related Signaling Pathways in Cardiovascular Health and Disease—The Pros and Cons for Cardioprotection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Identified Circular Transcript of RCAN2, circ-RCAN2, Shows Deviated Expression Pattern in Pig Reperfused Infarcted Myocardium and Hypoxic Porcine Cardiac Progenitor Cells In Vitro

1
Department of Internal Medicine II, Division of Cardiology, Medical University of Vienna, 1090 Vienna, Austria
2
Institute of Information Systems Engineering, Research Unit of Information and Software Engineering, Vienna University of Technology, 1040 Vienna, Austria
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(3), 1390; https://doi.org/10.3390/ijms22031390
Submission received: 23 December 2020 / Revised: 27 January 2021 / Accepted: 28 January 2021 / Published: 30 January 2021

Abstract

:
Circular RNAs (circRNAs) are crucial in gene regulatory networks and disease development, yet circRNA expression in myocardial infarction (MI) is poorly understood. Here, we harvested myocardium samples from domestic pigs 3 days after closed-chest reperfused MI or sham surgery. Cardiac circRNAs were identified by RNA-sequencing of rRNA-depleted RNA from infarcted and healthy myocardium tissue samples. Bioinformatics analysis was performed using the CIRIfull and KNIFE algorithms, and circRNAs identified with both algorithms were subjected to differential expression (DE) analysis and validation by qPCR. Circ-RCAN2 and circ-C12orf29 expressions were significantly downregulated in infarcted tissue compared to healthy pig heart. Sanger sequencing was performed to identify the backsplice junctions of circular transcripts. Finally, we compared the expressions of circ-C12orf29 and circ-RCAN2 between porcine cardiac progenitor cells (pCPCs) that were incubated in a hypoxia chamber for different time periods versus normoxic pCPCs. Circ-C12orf29 did not show significant DE in vitro, whereas circ-RCAN2 exhibited significant ischemia-time-dependent upregulation in hypoxic pCPCs. Overall, our results revealed novel cardiac circRNAs with DE patterns in pCPCs, and in infarcted and healthy myocardium. Circ-RCAN2 exhibited differential regulation by myocardial infarction in vivo and by hypoxia in vitro. These results will improve our understanding of circRNA regulation during acute MI.

1. Introduction

Circular RNAs (circRNAs) are long non-coding RNAs (lncRNAs) characterized by a closed-loop structure that lacks 3′ poly(A) tails and 5′ cap structures [1]. Since their initial discovery in 1976 in viroids (i.e., circular single-stranded RNA pathogens of higher plants) [2], circRNAs have long been thought to be splicing by-products or generated through splicing errors [3,4]. The development of high-throughput sequencing has enabled increasing circRNA research. Over the past decade, next-generation sequencing (NGS) has led to the identification of thousands of circRNAs in humans, nematodes, mice, and pigs [5,6,7,8], and their functions and roles in diseases have been investigated in neurologic disorders, cancer, and cardiovascular diseases (CVD) [9,10,11]. CircRNAs reportedly act as microRNA (miRNA) sponges by inhibiting miRNA-mRNA binding [12,13], and can also interact with RNA-binding proteins [14,15]. Despite their classification as non-coding molecules, some circRNAs can be translated in a cap-independent manner [16,17]. However, the exact mechanism of circRNA translation is not yet fully understood. Furthermore, circRNAs show tissue- and development-specific expression patterns [18,19]. Investigating circRNA expression patterns is essential for improving our understanding of their functions in diseases. Recent studies indicate imperative roles of circRNAs in cellular processes and disease development [20,21]. A 2019 study demonstrated circFOXP1 involvement in the transcriptional preservation of stem cell identity, and showed that silencing circFOXP1 reduces human mesenchymal stem cell (huMSC) growth and proliferation [22].
In CVD, circFndc3b modulates cardiac function after myocardial infarction (MI) via FUS/VEGF signaling [15]. Additionally, the circRNAs ZFAS1 and CDR1as may act as novel whole blood biomarkers that can distinguish between acute MI (AMI) and control patients [23]. Moreover, a 2019 study identified circNfix as a super-enhancer-regulated circRNA involved in cardiac regeneration after MI [8]. However, this area has not yet been sufficiently studied and there remains a need to identify altered circRNA expression patterns in AMI and to elaborate their molecular functions, which may lead to the development of new diagnostic and therapeutic treatment methods. In the present study, we performed RNA sequencing analysis (RNA-seq) and used two state-of-the-art algorithms to identify circRNA expressions in infarcted and control pig hearts. We then performed qPCR with divergent primers to validate our RNA-seq data in tissue samples. Finally, we further analyzed the expression patterns of the two most promising circRNAs, circ-RCAN2 and circ-C12orf29, in porcine cardiac progenitor cells (pCPCs) under hypoxic and normoxic conditions. Our findings will advance our knowledge of circRNAs in relation to AMI, and improve our understanding of their roles in CVD.

2. Results

2.1. Bioinformatics Analysis of Cardiac-Related circRNAs in Healthy and Infarcted Pig Hearts

We analyzed the expression of cardiac circRNAs in infarcted and healthy myocardium tissue of domestic pigs, three days after closed-chest reperfused MI (or sham operation). To identify cardiac circRNAs in rRNA-depleted RNA of infarcted and healthy porcine myocardium, we used the CIRIfull [24] and KNIFE [25] algorithms. KNIFE is pre-packaged with annotated junction indices for human, mouse, rat, and drosophila, but not for Sus scrofa. Therefore, we created a new junction index using the Ensemble reference genome Sus scrofa 11.1 and annotations (downloaded from iGenome, https://support.illumina.com/sequencing/sequencing_software/igenome.html). Transcriptomic sequences were extracted with gffread, and ribosomal sequences were acquired from NCBI. From these input data, we prepared Bowtie indices, as well as linear and scrambled junction indices, following the instructions at https://github.com/lindaszabo/KNIFE/tree/master/createJunctionIndex. CIRIfull and KNIFE were used independently from each other to ensure correct backsplice junction (BSJ) identification in circular transcripts and to eliminate false positives [26]. Figure 1 shows the study design and bioinformatics pipeline.
CIRIfull identified 5487 fully reconstructed circRNAs, while KNIFE identified the BSJs of 8668 circRNAs. We additionally analyzed the types of circRNAs expressed in porcine hearts, and found that most of the candidates were exonic.
Using the fully reconstructed circRNAs detected by CIRIfull, we performed differential expression (DE) analysis with DESeq2. This revealed significant DE of eight circRNAs in infarcted pig hearts compared to healthy controls. In Figure 2A, a volcano plot shows the p values and log2 fold changes. In Figure 2B, a heat-map displays the read counts of our RNA-seq data. Figure 2C illustrates the circRNA distribution on the chromosomes, indicating that the most circRNAs in pig hearts were located on chromosome 1.

2.2. Validation of Our RNA-Seq Data

Among the total identified circRNAs, we selected the circRNAs that were significantly dysregulated between the two study groups for validation in porcine heart tissue using quantitative PCR (qPCR). For these eight circRNAs, we designed divergent primers spanning the BSJ (Figure 3A). The qPCR data validated four of the circRNAs that showed significantly deviating expression patterns between AMI and healthy pig myocardium in our RNA-seq data: circ-ZNF644 (p = 2.25 × 10−5), circ-SLCO5A1 (p = 0.021), circ-RCAN2 (p = 0.027), and circ-C12orf29 (p = 0.024) (Figure 3B). Of these four circRNAs, only two showed significant downregulation in reperfused AMI compared to control tissue based on qPCR: circ-C12orf29 (p = 0.004) and circ-RCAN2 (p = 0.044) (Figure 3B). After qPCR, Sanger sequencing was performed, confirming the BSJ of the circRNAs (Figure 3C).

2.3. CircRNA Expression Pattern in Hypoxic and Normoxic pCPCs In Vitro

Based on our qPCR results, we selected circ-C12orf29 and circ-RCAN2 for further investigation in pCPCs under normoxic and hypoxic conditions in vitro. The genomic locus of circ-C12orf29 in pigs is on chromosome 5. To our knowledge, no further function has been described for circC12orf29. Our bioinformatics analysis revealed that the genomic locus of circ-RCAN2 in pigs lies on chromosome 7 on the regulator of the calcineurin 2 gene (RCAN2), and that the BSJ is located at 7:41003826-41008498.
We isolated pCPCs as previously described [27], and performed immunofluorescence staining. Our results showed positive expression of connexin-43 (Cx43) (Figure 4a, green), alpha smooth muscle actin (αSMA) (Figure 4b, green), islet-1 (Isl-1) (Figure 4c, green), pro-BNP (Figure 4d, green), and stem cells antigen-1 (Sca-1) (Figure 4e, green), confirming successful isolation of Sca-1+, and Isl1+ pCPCs.
In our in vitro experiment, pCPCs were cultured either under normoxic conditions, or in a hypoxic chamber for 60 min, 90 min, 2 h, 4 h, or 8 h without further normoxic culturing of the cells, thus without mimicking reperfusion in vivo. Then, we measured the expressions of circ-C12orf29 and circ-RCAN2. Circ-C12orf29 did not exhibit any significant up- or downregulation in vitro. However, circ-RCAN2 showed a deviated expression pattern in pCPCs, with upregulation in hypoxic pCPCs compared to in normoxic cells (Figure 5). In contrast with the circ-RCAN2 downregulation observed in reperfused AMI tissue at 3 days after infarction, constant hypoxia of pCPCs led to significant circ-RCAN2 upregulation, which was dependent on ischemia time, compared to control cells (Figure 5).

3. Discussion

In our present study, we revealed novel cardiac circRNAs in a pig model using deep RNA-sequencing paired with state-of-the-art bioinformatics algorithms. We also created annotated junction indices for Sus scrofa for the KNIFE algorithm, as KNIFE only had annotated junction indices for mouse, rat, human, and drosophila. Our expression pattern analysis revealed decreased circ-RCAN2 expression in porcine reperfused ischemic myocardium compared to healthy porcine heart tissue. On the other hand, within a constant hypoxic milieu, pCPCs showed a time-dependent increase in circ-RCAN2 expression. These results suggest that ischemia time and restoration of normoxic conditions play a role in circRNA regulation.
The literature includes only scarce information regarding circRNA expression patters in infarcted myocardium; therefore, here we focused on the detection of novel circRNAs linked to MI in pig hearts. We found that the genomic locus of circ-RCAN2 in pigs lies on chromosome 7 on RCAN2, and the BSJ is located at 7:41003826-41008498. Furthermore, we demonstrated that circ-RCAN2 was significantly downregulated in a porcine model of reperfused MI compared to healthy hearts. To our knowledge, this is the first description of a circular transcript of RCAN2.
The RCAN gene family members RCAN1, RCAN2, and RCAN3—also known as adapt78 [28], DSCR1 [29], and MCIP1 [30]—are inhibitors of the Ca2+-activated protein phosphatase calcineurin [29,31,32]. Calcineurin is crucial for cardiac remodeling [33], and is present at notably high levels in excitatory cells, such as cardiomyocytes, that exhibit frequent and rapid changes in cytoplasmic calcium, indicating that such cells require high calcineurin activity [34]. Upregulated expression levels of RCAN homologs have been reported in various diseases, including cardiac hypertrophy and diabetes [35].
One member of the RCAN gene family, RCAN1, has been investigated for its role in ischemia/reperfusion (I/R) injury [36]. In an in vitro study, Parra et al. found that RCAN1-depleted cardiomyocytes were more susceptible to I/R damage [33]. However, our present bioinformatics analysis did not reveal a circular transcript of RCAN1. One possible explanation is that no circular transcript of RCAN1 exists. In this study, we detected circ-RCAN2. Little is presently known regarding the possible function of RCAN2 in the heart and, to our knowledge, no circular transcript of RCAN2 has previously been identified. We found that circ-RCAN2 was significantly downregulated in the reperfused infarcted tissue after 3 days. In contrast, in pCPCs incubated under constant hypoxia (without mimicking reperfusion) for different time periods, we observed time-dependent significant upregulation of circ-RCAN2 compared to in normoxic cells. This may be related to the fact that the pigs underwent 90-min occlusion of the left anterior descending coronary artery (LAD), followed by reperfusion, and tissue was harvested 3 days later. In contrast, in the in vitro experiment, we applied constant hypoxic conditions to our pCPCs for different time periods, without subsequent normoxia. This difference might explain the biphasic expression pattern of circ-RCAN2, and emphasizes the importance of reperfusion and ischemia time in AMI, and the differences between in vivo and in vitro conditions.
Our study has several limitations. One is the low number of animals. However, this was a pilot study, and the results are the first documentation of the existence of circ-RCAN2 and its dysregulation in infarcted cardiac tissue in a translational large animal experiment. Additionally, in our experiments, we harvested the transmural myocardial area for analysis at 3 days post-infarction, since the local inflammation and edema peaks at 3 days post-infarction onset. Despite the well-known heterogeneity of transmural myocardial ischemia, we collected transmural tissue samples because it is not possible to macroscopically separate subendocardial tissue from the mid-cardiac and subepicardial tissues, and a transmural section of the infarcted heart mirrors the complexity of the global pathophysiology processes and similarities to human conditions. In accordance with our previous findings regarding the intrinsic remote conditioning of ischemia-nonaffected regions after ischemia/reperfusion [37], we used non-infarcted animals as controls.
Although thousands of circRNAs have been described in the mammalian heart [7], it remains largely unknown which circRNAs may be linked to AMI or hypoxia. CircRNAs can possess gene regulatory functions, reportedly acting as miRNA sponges [12,13] or interacting with RNA-binding proteins (RBPs) [38,39]. Additionally, they show the potential to be used as diagnostic or therapeutic biomarkers in diseases, such as cancer and CVD [40,41,42]. Here, we investigated circRNA expression in porcine heart tissue to unravel novel circRNAs associated with AMI. To our knowledge, our results constitute the first reported detection of a circular transcript of RCAN2, and we further show that circ-RCAN2 is important in infarcted tissue and hypoxic pCPCs in vitro.

4. Materials and Methods

4.1. Animal Study

This study included six female domestic pigs (Sus scrofa domesticus), which were 6 months old and had weights ranging between 25–30 kg at the time of the experiment. These pigs were randomized to either the AMI (n = 3) or sham-operated control (n = 3) group. After overnight fasting, the animals were anesthetized via intramuscular injection using 12 mg/kg ketamine hydrochloride, 1 mg/kg xylazine, and 0.04 mg/kg atropine. Then, coronary angiography was performed using a 6F guiding catheter (Terumo Medical Corporation, Somerset, New Jersey, USA). For pigs in the AMI group, a balloon catheter (Maverick XL Monorail Balloon Catheter, Boston Scientific) was inserted into the LAD after the first diagonal branch. Then, reperfused AMI was induced by balloon occlusion of the mid LAD at 5 atm for 90 min, followed by balloon deflation and reperfusion. The occlusion was confirmed by control angiography. Three days after AMI or sham procedure, all animals were euthanized under general anesthesia using 10 mL intravenous saturated KCl (10%).

4.2. Tissue Sample Collection

During autopsy, we collected tissue samples (150–250 mg) from the transmural infarcted area of the pig hearts (distal anterior tissue, distal to the LAD occlusion site), and from the same location of the left ventricle from healthy pig hearts. In the early phase of infarction, the infarcted tissue was characterized by edema and inflammation, peaking at day 2–3 post onset. For that reason, we chose day 3 for tissue collection. Tissue samples were stored at −80 °C in RNAlater (Thermo Fisher Scientific, Waltham, MA, USA).

4.3. RNA Extraction from Porcine Heart Tissue

Tissue samples were thawed at room temperature, and then 25–30 mg of tissue was cut and placed in Precellys CK28 tubes containing Qiazol Lysis buffer (Qiagen, Hilden, Germany). These tissue samples were disrupted and homogenized at room temperature using the 5000:3×20 - 20 s wait program on the Precellys 24 (Bertin, Rockville, MD, USA). Next, the supernatant was transferred into new 2-mL tubes (Eppendorf, Hamburg, Germany) and 200 µL chloroform was added. This mixture was incubated at room temperature for 2 min, and then the tubes were centrifuged at 11.600× g at 4 °C for 15 min. The upper aqueous phase was collected and transferred to 2-mL tubes (Eppendorf), and total RNA was isolated using the miRNeasy Mini kit (Qiagen, Hilden, Germany), including a DNase I digestion, following the manufacturer’s protocol. Total RNA quality and quantity were assessed using the Nanodrop One (Thermo Fisher Scientific).

4.4. RNA Sequencing of circRNAs

Isolated total RNA was sequenced using the Illumina NextSeq550 sequencing approach. We performed mRNA enrichment using the NEBNext rRNA depletion kit (New England Biolabs, Ipswich, MA, USA). No poly(A) enrichment was performed, to enable accurate circRNA sequencing. The sequencing samples had a mean depth of 140 million paired-end (150 bp) reads. Library preparation and sequencing were performed at the VBCF NGS Unit (www.viennabiocenter.org/facilities).

4.5. Bioinformatics

The raw sequencing reads (BAM files) were sorted using samtools, and then converted to FASTQ files using bedtools. Since we had paired-end data, this resulted in two files per sample (with ending R1 and R2). We analyzed our data using two algorithms for circular RNA identification: CIRI-full and KNIFE. Further analyses were performed using the results of CIRI-full, but we checked whether the circular RNA candidates were also identified using KNIFE. Circular RNA candidates identified by both algorithms are far less likely to be false positives, as demonstrated by Hansen et al. [26]. When using the KNIFE algorithm, we trimmed the reads with Trim Galore (default values). On the other hand, CIRI-full requires reads of equal length, and thus we had to use the untrimmed reads. For both algorithms, quality control was performed using FastQC. As further input for the two algorithms, we used reference genomes and annotation files from iGenomes (Ensembl Sus Scrofa 11.1). CIRI-full requires prior alignment using the Burrows-Wheeler Alignment Tool (bwa mem), which we performed with option T = 19 (i.e., alignment with a score < 19 is skipped in the output). For KNIFE, we used the recommended values of 13 for junction overlap (because the reads were paired-end), and 100 for the “ntrim” parameter (because the original read length was 150, and the recommended value is 2/3 of this).

4.6. Primers

CircRNA primers were designed in accordance with the results of our bioinformatics analyses. All circRNA primers were designed to amplify the backsplice junction. The primers were purified by HPLC, and were purchased from Microsynth (Balgach, Switzerland). The primer sequences are listed in the Supplementary Material 1, Table S1. We used β-actin and hypoxanthine guanine phosphoribosyltransferase (HPRT) as housekeeping genes in porcine heart tissue, and HPRT as a housekeeping gene in pCPCs.

4.7. cDNA Synthesis

We performed cDNA synthesis using the QuantiTect Reverse Transcription kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol.

4.8. Real-Time Quantitative PCR Analysis

Experiments were each performed twice. For tissue sample analysis, the circRNA expression levels were normalized to the geometric mean of two housekeeping genes: β-actin and HPRT. DE analysis was performed using the ΔΔCt method. We performed qPCR analyses of pig heart tissue in duplicate with a standard curve. For the pCPC analysis, the circRNA expression levels were normalized to HPRT. We performed qPCR analyses of pCPCs in triplicate, and the primers were diluted with RNase-free H2O to a final concentration of 10 µM. No-template controls (run in duplicate) served as control wells. The PCR plates were sealed using transparent PCR foil, and qPCR was run using Quantstudio 5 (Thermo Fisher Scientific, Waltham, MA, USA).

4.9. Sanger Sequencing

To confirm our sequencing results and the head-to-tail splicing of the selected circRNAs, Sanger sequencing was performed by Microsynth (Balgach, Switzerland) using our designed circRNA primers.

4.10. Porcine Cardiac Progenitor Cells

Porcine cardiac progenitor cells were isolated from pig hearts, and qPCR was performed to assess the relative gene expressions of ISL-1, SCA-1, c-Kit, and BNP, as previously described [27]. To further evaluate the expression of progenitor cell markers, such as islet-1 (Isl-1) and stem cell antigen (Sca-1), we also performed immunofluorescence staining. The cells were cultured in a 1:1 ratio of M199 and Dulbecco’s modified Eagle’s medium (DMEM) (Sigma Aldrich, St. Louis, MO, USA), containing 10% fetal bovine serum (FBS, Biochrom Ltd., Cambridge, UK) and 1% pen-strep. Cells were grown to 90% confluence, and then seeded in 96-well plates containing fresh medium and cultivated overnight in a CO2 incubator at 37 °C. The next day, the medium was removed, and the cells were washed with 1× D-PBS (Sigma Aldrich) and then fixed with 4% paraformaldehyde (Carl Roth GmbH, Karlsruhe, Germany) in 1× D-PBS for 15 min. Then, the paraformaldehyde was removed, and 1× PBS (Sigma Aldrich) was added to each well. The cells were incubated at room temperature for 5 min, and then PBS was removed and each well was washed twice with 1× D-PBS. The primary antibodies were diluted following the manufacturer’s protocol, with permeabilization buffer comprising 3% non-fat dry milk and 0.1% Triton X-100 (both from Sigma Aldrich). The primary antibodies were added (50 µL per well), followed by a 2-h incubation at room temperature, after which the primary antibodies were aspirated, and each well was washed with 1× D-PBS. Next, the secondary antibodies were diluted with permeabilization buffer following the manufacturer’s protocol, and added to the wells, followed by a 2-h incubation at room temperature protected from light. Finally, the cells were counterstained using Hoechst (1:5000 in D-PBS, Sigma Aldrich) and Phalloidin (1:40 in PBS, Thermo Fisher Scientific). Immunofluorescence staining was analyzed using an Olympus IX83 Inverted Microscope with the cellSens imaging software (Olympus, Tokyo, Japan). The primary and secondary antibodies, including their dilutions, are listed in the Supplementary materials, Table S2.

4.11. Hypoxic and Normoxic Conditions of pCPCs

Next, we investigated circRNA expression in pCPCs under normoxic and hypoxic conditions. We cultured pCPCs in a 1:1 ratio of M199 and DMEM (Sigma Aldrich) with 10% FBS and 1% pen-strep, in 6-well plates, with incubation at 37 °C. The media was changed every 2 days until the cells reached the 90% confluence, at which time the medium was changed to serum-free medium. After an overnight incubation at 37 °C, the cells were placed in a hypoxia incubator chamber (Stemcell Technologies, Vancouver, Canada) for different time periods (60 min, 90 min, 2 h, 4 h, and 8 h). Control cells were not put in the hypoxia incubator chamber (normoxic conditions). To harvest the cells, the medium was removed, the wells were washed with 1× D-PBS, and 350 µL RLT buffer containing β-mercaptoethanol (1:100, Sigma Aldrich) was added to each well. Finally, the cells were transferred into 2-mL tubes (Eppendorf) and stored at −80 °C until further use.

4.12. RNA Extraction from pCPCs

Cells in RLT buffer were thawed at room temperature until the solution was completely liquid, and then transferred into QIAshredder tubes (Qiagen). These tubes were then centrifuged at 11.600× g for 2 min at 4 °C, and the flow-through was transferred into 2-mL collection tubes (Qiagen). RNA was isolated using the QIAcube and the miRNeasy Mini kit (Qiagen). RNA quantity and quality was assessed using the Nanodrop One (Thermo Fisher Scientific).

4.13. Statistics

We performed differential expression analysis of the RNA-seq data (cases vs. controls) using the package DESeq2 [43], which applies a model based on the negative binomial distribution. Prior to the analysis, we filtered transcripts with low expression by requiring counts to be over 3 in more than two of the samples. Results with a p value of <0.05 were considered significant. We performed differential expression analyses of the qPCR results for cardiac tissue (cases vs. control) and pCPCs (hypoxia of 60 min, 90 min, 2 h, 4 h, or 8 h vs. normoxia) using the ΔΔCt method, with unequal variance t-tests applied to the log fold changes.

4.14. Ethics Approval

Animal handling procedures were conducted in accordance with the “Guide for The Care and Use of Laboratory Animals“ [44], and were approved by the ethical review committee of the University of Kaposvar (SOI/31/26-11/2014).

5. Conclusions

In conclusion, here we identified and validated novel circRNAs, and examined their variable expression patterns in infarcted and healthy pig hearts. We confirmed the presence of circ-RCAN2 in pCPCs in vitro, and revealed important differences in deregulation under hypoxic conditions both in vivo and in vitro. Our findings elucidate an interesting new target for MI.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/22/3/1390/s1.

Author Contributions

Conceptualization, J.M.-T., P.E., J.W., and M.G.; methodology, J.M.-T., P.E., K.Z., and D.L.; software, P.E.; validation, J.M.-T., N.K., and A.S.; formal analysis, J.M.-T. and P.E.; investigation, J.M.-T.; resources, M.G.; data curation, J.M-T.; writing—original draft preparation, J.M.-T., P.E., and M.G.; writing—review and editing, J.W., N.K., A.S., K.Z., D.T., D.L., E.H., G.G., N.P., and M.G.; visualization, J.M.-T. and P.E.; supervision, J.W. and M.G.; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CResPace (grant number 732170), SCIENCE (grant number 643478), and REGENHEART (grant number 731532).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the study design; in the data collection, analyses, or interpretation; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

AMIAcute myocardial infarction
BSJBacksplice junction
circRNACircular RNA
CVDCardiovascular diseases
Cx43Connexin-43
DEDifferential expression
I/R InjuryIschemia/reperfusion injury
HPRTHypoxanthine-guanine phosphoribosyltransferase
Isl-1Islet-1
lncRNALong non-coding RNA
NGSNext-generation sequencing
pCPCspro-BNPPorcine cardiac progenitor cellsPro-brain natriuretic peptide
qPCRReal-time quantitative polymerase chain reaction
RNA-seqRNA sequencing
Sca-1Stem cells antigen-1
αSMAAlpha smooth muscle actin

References

  1. Barrett, S.P.; Salzman, J. Circular RNAs: Analysis, expression and potential functions. Development 2016, 143, 1838–1847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sanger, H.L.; Klotz, G.; Riesner, D.; Gross, H.J.; Kleinschmidt, A.K. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc. Natl. Acad. Sci. USA 1976, 73, 3852–3856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Capel, B.; Swain, A.; Nicolis, S.; Hacker, A.; Walter, M.; Koopman, P.; Goodfellow, P.; Lovell-Badge, R. Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 1993, 73, 1019–1030. [Google Scholar] [CrossRef]
  4. Cocquerelle, C.; Mascrez, B.; Hétuin, D.; Bailleul, B. Mis-splicing yields circular RNA molecules. FASEB J. 1993, 7, 155–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Liang, G.; Yang, Y.; Niu, G.; Tang, Z.; Li, K. Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages. DNA Res. 2017, 24, 523–535. [Google Scholar] [CrossRef] [Green Version]
  6. Ivanov, A.; Memczak, S.; Wyler, E.; Torti, F.; Porath, H.T.; Orejuela, M.R.; Piechotta, M.; Levanon, E.Y.; Landthaler, M.; Dieterich, C.; et al. Analysis of intron sequences reveals hallmarks of circular rna biogenesis in animals. Cell Rep. 2015, 10, 170–177. [Google Scholar] [CrossRef] [Green Version]
  7. Tan, W.L.; Lim, B.T.; Anene-Nzelu, C.G.; Ackers-Johnson, M.; Dashi, A.; See, K.; Tiang, Z.; Lee, D.P.; Chua, W.W.; Luu, T.D.; et al. A landscape of circular RNA expression in the human heart. Cardiovasc. Res. 2016, 113, 298–309. [Google Scholar] [CrossRef]
  8. Huang, S.; Li, X.; Zheng, H.; Si, X.; Li, B.; Wei, G.; Li, C.; Chen, Y.; Chen, Y.; Liao, W.; et al. Loss of Super-enhancer-regulated circRNA Nfix induces cardiac regeneration after myocardial infarction in adult mice. Circulation 2019, 139, 2857–2876. [Google Scholar] [CrossRef]
  9. Song, T.; Xu, A.; Zhang, Z.; Gao, F.; Zhao, L.; Chen, X.; Gao, J.; Kong, X. CircRNA hsa_circRNA_101996 increases cervical cancer proliferation and invasion through activating TPX2 expression by restraining miR-8075. J. Cell. Physiol. 2019, 234, 14296–14305. [Google Scholar] [CrossRef]
  10. Wu, F.; Han, B.; Wu, S.; Yang, L.; Leng, S.; Li, M.; Liao, J.; Wang, G.; Ye, Q.; Zhang, Y.; et al. Circular RNA TLK1 aggravates neuronal injury and neurological deficits after ischemic stroke via miR-335-3p/TIPARP. J. Neurosci. 2019, 39, 7369–7393. [Google Scholar] [CrossRef]
  11. Zhou, L.-Y.; Zhai, M.; Huang, Y.; Xu, S.; An, T.; Wang, Y.-H.; Zhang, R.-C.; Liu, C.; Dong, Y.-H.; Wang, M.; et al. The circular RNA ACR attenuates myocardial ischemia/reperfusion injury by suppressing autophagy via modulation of the Pink1/ FAM65B pathway. Cell Death Differ. 2019, 26, 1299–1315. [Google Scholar] [CrossRef] [PubMed]
  12. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef] [PubMed]
  13. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef] [PubMed]
  14. Du, W.W.; Yang, W.; Chen, Y.; Wu, Z.-K.; Foster, F.S.; Yang, Z.; Li, X.; Yang, B.B. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur. Heart J. 2016, 38, 1402–1412. [Google Scholar] [CrossRef] [Green Version]
  15. Garikipati, V.N.S.; Verma, S.K.; Cheng, Z.; Liang, D.; Truongcao, M.M.; Cimini, M.; Yue, Y.; Huang, G.; Wang, C.; Benedict, C.; et al. Circular RNA CircFndc3b modulates cardiac repair after myocardial infarction via FUS/VEGF-A axis. Nat. Commun. 2019, 10, 1–14. [Google Scholar] [CrossRef]
  16. Legnini, I.; Di Timoteo, G.; Rossi, F.; Morlando, M.; Briganti, F.; Sthandier, O.; Fatica, A.; Santini, T.; Andronache, A.; Wade, M.; et al. Circ-ZNF609 Is a Circular RNA that can be translated and functions in myogenesis. Mol. Cell 2017, 66, 22–37. [Google Scholar] [CrossRef] [Green Version]
  17. Yang, Y.; Fan, X.; Mao, M.; Song, X.; Wu, P.; Zhang, Y.; Jin, Y.; Yang, Y.; Chen, L.-L.; Wang, Y.; et al. Extensive translation of circular RNAs driven by N6-methyladenosine. Cell Res. 2017, 27, 626–641. [Google Scholar] [CrossRef] [Green Version]
  18. Xia, S.; Feng, J.; Lei, L.; Hu, J.; Xia, L.; Wang, J.; Xiang, Y.; Liu, L.; Zhong, S.; Han, L.; et al. Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes. Brief. Bioinform. 2016, 18, 984–992. [Google Scholar] [CrossRef] [Green Version]
  19. Salzman, J.; Chen, R.E.; Olsen, M.N.; Wang, P.L.; Brown, P.O. Cell-type specific features of circular RNA expression. PLoS Genet 2013, 9, e1003777. [Google Scholar] [CrossRef]
  20. Zhang, M.; Wang, Z.; Cheng, Q.; Wang, Z.; Lv, X.; Li, N. Circular RNA (circRNA) CDYL induces myocardial regeneration by ceRNA after myocardial infarction. Med Sci. Monit. 2020, 26, e923188. [Google Scholar] [CrossRef] [Green Version]
  21. Zhu, Y.; Pan, W.; Yang, T.; Meng, X.; Jiang, Z.; Tao, L.; Wang, L. Upregulation of circular RNA CircNFIB Attenuates cardiac fibrosis by sponging miR-433. Front. Genet. 2019, 10, 564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cherubini, A.; Barilani, M.; Rossi, R.L.; Jalal, M.M.K.; Rusconi, F.; Buono, G.; Ragni, E.; Cantarella, G.; Simpson, A.H.W.R.; Péault, B.; et al. FOXP1 circular RNA sustains mesenchymal stem cell identity via microRNA inhibition. Nucleic Acids Res. 2019, 47, 5325–5340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zhang, Y.; Sun, L.; Xuan, L.; Pan, Z.; Li, K.; Liu, S.; Huang, Y.; Zhao, X.; Huang, L.; Wang, Z.; et al. Reciprocal changes of circulating long non-coding RNAs ZFAS1 and CDR1AS predict acute myocardial infarction. Sci. Rep. 2016, 6, 1–11. [Google Scholar] [CrossRef] [PubMed]
  24. Zheng, Y.; Ji, P.; Chen, S.; Hou, L.; Zhao, F. Reconstruction of full-length circular RNAs enables isoform-level quantification. Genome Med. 2019, 11, 1–20. [Google Scholar] [CrossRef]
  25. Szabo, L.; Morey, R.; Palpant, N.J.; Wang, P.L.; Afari, N.; Jiang, C.; Parast, M.M.; Murry, C.E.; Laurent, L.C.; Salzman, J. Statistically based splicing detection reveals neural enrichment and tissue-specific induction of circular RNA during human fetal development. Genome Biol. 2015, 16, 1–26. [Google Scholar] [CrossRef] [Green Version]
  26. Hansen, T.B. Improved circRNA Identification by Combining Prediction Algorithms. Front. Cell Dev. Biol. 2018, 6, 20. [Google Scholar] [CrossRef]
  27. Zlabinger, K.; Spannbauer, A.; Traxler, D.; Gugerell, A.; Lukovic, D.; Winkler, J.; Mester-Tonczar, J.; Podesser, B.K.; Gyongyosi, M. MiR-21, MiR-29a, GATA4, and MEF2c expression changes in endothelin-1 and angiotensin II Cardiac hypertrophy stimulated Isl-1+Sca-1+c-kit+ porcine cardiac progenitor cells in vitro. Cells 2019, 8, 1416. [Google Scholar] [CrossRef] [Green Version]
  28. Crawford, D.R.; Leahy, K.P.; Abramova, N.; Lan, L.; Wang, Y.; Davies, K.J. Hamsteradapt78mRNA Is a down syndrome critical region homologue that is inducible by oxidative stress. Arch. Biochem. Biophys. 1997, 342, 6–12. [Google Scholar] [CrossRef]
  29. Kingsbury, T.J.; Cunningham, K.W. A conserved family of calcineurin regulators. Genes Dev. 2000, 14, 1595–1604. [Google Scholar]
  30. Yang, J.; Rothermel, B.; Vega, R.B.; Frey, N.; McKinsey, T.A.; Olson, E.N.; Bassel-Duby, R.; Williams, R.S. Independent signals control expression of the calcineurin inhibitory proteins MCIP1 and MCIP2 in striated muscles. Circ. Res. 2000, 87, E61–E68. [Google Scholar] [CrossRef] [Green Version]
  31. Strippoli, P.; Lenzi, L.; Petrini, M.; Carinci, P.; Zannotti, M. A New gene family including DSCR1 (down syndrome candidate region 1) and ZAKI-4: Characterization from Yeast to human and identification of DSCR1-like 2, a novel human member (DSCR1L2). Genomics 2000, 64, 252–263. [Google Scholar] [CrossRef] [PubMed]
  32. Rothermel, B.; Vega, R.B.; Yang, J.; Wu, H.; Bassel-Duby, R.; Williams, R.S. A protein encoded within the down syndrome critical region is enriched in striated muscles and inhibits calcineurin signaling. J. Biol. Chem. 2000, 275, 8719–8725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Parra, V.; Altamirano, F.; Hernández-Fuentes, C.P.; Tong, D.; Kyrychenko, V.; Rotter, D.; Pedrozo, Z.; Hill, J.A.; Eisner, V.; Lavandero, S.; et al. Down syndrome critical region 1 Gene, Rcan1, helps maintain a more fused mitochondrial network. Circ. Res. 2018, 122, e20–e33. [Google Scholar] [CrossRef] [PubMed]
  34. Parra, V.; Rothermel, B.A. Calcineurin signaling in the heart: The importance of time and place. J. Mol. Cell. Cardiol. 2017, 103, 121–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lee, S.-K.; Ahnn, J. Regulator of calcineurin (RCAN): Beyond down syndrome critical region. Mol Cells 2020, 43, 671–685. [Google Scholar] [CrossRef]
  36. Rotter, D.; Grinsfelder, D.B.; Parra, V.; Pedrozo, Z.; Singh, S.; Sachan, N.; Rothermel, B.A. Calcineurin and its regulator, RCAN1, confer time-of-day changes in susceptibility of the heart to ischemia/reperfusion. J. Mol. Cell. Cardiol. 2014, 74, 103–111. [Google Scholar] [CrossRef] [Green Version]
  37. Pavo, N.; Lukovic, D.; Zlabinger, K.; Zimba, A.; Lorant, D.; Goliasch, G.; Winkler, J.; Pils, D.; Auer, K.; Ankersmit, H.J.; et al. Sequential activation of different pathway networks in ischemia-affected and non-affected myocardium, inducing intrinsic remote conditioning to prevent left ventricular remodeling. Sci. Rep. 2017, 7, srep43958. [Google Scholar] [CrossRef] [Green Version]
  38. Conn, S.J.; Pillman, K.A.; Toubia, J.; Conn, V.M.; Salmanidis, M.; Phillips, C.A.; Roslan, S.; Schreiber, A.W.; Gregory, P.A.; Goodall, G.J. The RNA binding protein quaking regulates formation of circRNAs. Cell 2015, 160, 1125–1134. [Google Scholar] [CrossRef] [Green Version]
  39. Dong, W.; Dai, Z.-H.; Liu, F.-C.; Guo, X.-G.; Ge, C.-M.; Ding, J.; Liu, H.; Yang, F. The RNA-binding protein RBM3 promotes cell proliferation in hepatocellular carcinoma by regulating circular RNA SCD-circRNA 2 production. Ebiomedicine 2019, 45, 155–167. [Google Scholar] [CrossRef] [Green Version]
  40. Li, Y.; Zheng, Q.; Bao, C.; Li, S.; Guo, W.; Zhao, J.; Chen, D.; Gu, J.; He, X.; Huang, S. Circular RNA is enriched and stable in exosomes: A promising biomarker for cancer diagnosis. Cell Res. 2015, 25, 981–984. [Google Scholar] [CrossRef] [Green Version]
  41. Holdt, L.M.; Stahringer, A.; Sass, K.; Pichler, G.; Kulak, N.A.; Wilfert, W.; Kohlmaier, A.; Herbst, A.; Northoff, B.H.; Nicolaou, A.; et al. Circular non-coding RNA ANRIL modulates ribosomal RNA maturation and atherosclerosis in humans. Nat. Commun. 2016, 7, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Wang, K.; Long, B.; Liu, F.; Wang, J.-X.; Liu, C.-Y.; Zhao, B.; Zhou, L.-Y.; Sun, T.; Wang, M.; Cui-Yun, L.; et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur. Heart J. 2016, 37, 2602–2611. [Google Scholar] [CrossRef] [PubMed]
  43. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. National Research Council. Guide for the Care and Use of Laboratory Animals, 8th ed.; National Academies Press (US): Washington, DC, USA, 2011; ISBN 978-0-309-15400-0. [Google Scholar]
Figure 1. Identification of circRNAs in pig hearts. Study design and bioinformatics pipeline for the identification of circRNAs using reads of rRNA-depleted RNA from pig hearts at three days after reperfused acute myocardial infarction (AMI) (n = 3) or sham operation (healthy pig hearts) (n = 3).
Figure 1. Identification of circRNAs in pig hearts. Study design and bioinformatics pipeline for the identification of circRNAs using reads of rRNA-depleted RNA from pig hearts at three days after reperfused acute myocardial infarction (AMI) (n = 3) or sham operation (healthy pig hearts) (n = 3).
Ijms 22 01390 g001
Figure 2. Regulation of circRNAs in infarcted (AMI) and healthy (control) pig hearts. (a) Volcano plot shows identified circRNAs. Red dots indicate significantly dysregulated circRNAs. (b) Heat-map shows the eight circRNAs with deviating expression patterns between AMI and control tissue. (c) Distribution of circRNAs among chromosomes.
Figure 2. Regulation of circRNAs in infarcted (AMI) and healthy (control) pig hearts. (a) Volcano plot shows identified circRNAs. Red dots indicate significantly dysregulated circRNAs. (b) Heat-map shows the eight circRNAs with deviating expression patterns between AMI and control tissue. (c) Distribution of circRNAs among chromosomes.
Ijms 22 01390 g002
Figure 3. Validation of the RNA-seq data in porcine heart tissue. (a) Divergent primers were used to amplify the backsplice junction of the circRNAs. Arrows show the directions of the forward/reverse primers. Numbers inside the boxes indicate the exons. (b) Differential expression analysis of our four circRNAs between reperfused AMI and healthy heart using qPCR. * p < 0.05, ** p < 0.01. Error bars show standard deviation (SD). (c) Sanger sequencing chromatogram shows the backsplice junction.
Figure 3. Validation of the RNA-seq data in porcine heart tissue. (a) Divergent primers were used to amplify the backsplice junction of the circRNAs. Arrows show the directions of the forward/reverse primers. Numbers inside the boxes indicate the exons. (b) Differential expression analysis of our four circRNAs between reperfused AMI and healthy heart using qPCR. * p < 0.05, ** p < 0.01. Error bars show standard deviation (SD). (c) Sanger sequencing chromatogram shows the backsplice junction.
Ijms 22 01390 g003
Figure 4. Characterization of pCPCs using immunofluorescence staining. Results showed positive expression of Cx43 (a), αSMA (b), and pro-BNP (d), and of markers associated with cardiac progenitor cells, such as Isl-1 (c) and Sca-1 (e). In panels (d,e), the cytoplasm was counterstained with Phalloidin. Nuclei were stained with Hoechst.
Figure 4. Characterization of pCPCs using immunofluorescence staining. Results showed positive expression of Cx43 (a), αSMA (b), and pro-BNP (d), and of markers associated with cardiac progenitor cells, such as Isl-1 (c) and Sca-1 (e). In panels (d,e), the cytoplasm was counterstained with Phalloidin. Nuclei were stained with Hoechst.
Ijms 22 01390 g004
Figure 5. Differential expression analysis of circ-RCAN2 in hypoxic and normoxic pCPCs in vitro. * p < 0.05, ** p < 0.01, *** p < 0.001. Error bars show SD.
Figure 5. Differential expression analysis of circ-RCAN2 in hypoxic and normoxic pCPCs in vitro. * p < 0.05, ** p < 0.01, *** p < 0.001. Error bars show SD.
Ijms 22 01390 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mester-Tonczar, J.; Einzinger, P.; Winkler, J.; Kastner, N.; Spannbauer, A.; Zlabinger, K.; Traxler, D.; Lukovic, D.; Hasimbegovic, E.; Goliasch, G.; et al. Novel Identified Circular Transcript of RCAN2, circ-RCAN2, Shows Deviated Expression Pattern in Pig Reperfused Infarcted Myocardium and Hypoxic Porcine Cardiac Progenitor Cells In Vitro. Int. J. Mol. Sci. 2021, 22, 1390. https://doi.org/10.3390/ijms22031390

AMA Style

Mester-Tonczar J, Einzinger P, Winkler J, Kastner N, Spannbauer A, Zlabinger K, Traxler D, Lukovic D, Hasimbegovic E, Goliasch G, et al. Novel Identified Circular Transcript of RCAN2, circ-RCAN2, Shows Deviated Expression Pattern in Pig Reperfused Infarcted Myocardium and Hypoxic Porcine Cardiac Progenitor Cells In Vitro. International Journal of Molecular Sciences. 2021; 22(3):1390. https://doi.org/10.3390/ijms22031390

Chicago/Turabian Style

Mester-Tonczar, Julia, Patrick Einzinger, Johannes Winkler, Nina Kastner, Andreas Spannbauer, Katrin Zlabinger, Denise Traxler, Dominika Lukovic, Ena Hasimbegovic, Georg Goliasch, and et al. 2021. "Novel Identified Circular Transcript of RCAN2, circ-RCAN2, Shows Deviated Expression Pattern in Pig Reperfused Infarcted Myocardium and Hypoxic Porcine Cardiac Progenitor Cells In Vitro" International Journal of Molecular Sciences 22, no. 3: 1390. https://doi.org/10.3390/ijms22031390

APA Style

Mester-Tonczar, J., Einzinger, P., Winkler, J., Kastner, N., Spannbauer, A., Zlabinger, K., Traxler, D., Lukovic, D., Hasimbegovic, E., Goliasch, G., Pavo, N., & Gyöngyösi, M. (2021). Novel Identified Circular Transcript of RCAN2, circ-RCAN2, Shows Deviated Expression Pattern in Pig Reperfused Infarcted Myocardium and Hypoxic Porcine Cardiac Progenitor Cells In Vitro. International Journal of Molecular Sciences, 22(3), 1390. https://doi.org/10.3390/ijms22031390

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop