1. Introduction
Cardiac development is a complex developmental process that initiates soon after gastrulation with the configuration of symmetrical pools of cardiomyogenic precursors that subsequently fused in the embryonic midline, leading to the formation of the cardiac straight tube [
1]. This early tube is composed of an external myocardial layer and an internal endocardial lining. Soon after, a rightward looping develops and embryonic atrial and ventricular chambers become progressively configured [
2]. At this stage, the embryonic heart becomes externally covered by a new epithelial lining, the embryonic epicardium (EE). The EE originates from the proepicardium (PE), a transitory structure with a cauliflower appearance located at the junction of the cardiac and hepatic anlagen within the
septum transversum [
3]. Cells emanating from the PE bridged to the naked embryonic myocardium providing an external lining [
4]. Soon after, a subset of these embryonic epicardial cells, undergoes an epithelial-to-mesenchymal transition (EMT), colonizing the subepicardial space and migrating into the embryonic myocardium [
5,
6]. In chicken, these epicardial derived cells (
i.e. EPDCs) have been consistently reported to contribute to distinct coronary vascular components, including vascular endothelium, smooth muscle and adventitial layers, as well as cardiac fibroblasts [
7,
8]. However, in mice, several lineage tracing experiments have provided evidence of a modest contribution of EPDCs to cardiac endothelium, while their contribution to coronary smooth muscle and adventitial vasculature and fibroskeleton is undisputed [
9].
In recent years, our understanding of the molecular mechanisms that govern proepicardial and epicardial formation has greatly emerged [
3]. Evidence highlighting the pivotal role of Bmp and Fgf during proepicardial and myocardial specification was reported by Kruithof
et al. [
10] demonstrating that Fgf enhances proepicardial specification while Bmp promotes cardiomyogenic differentiation. More recently, the role of distinct transcription factors has emerged as key players in epicardial development.
Wt1 is essential for epicardial EMT and maturation [
11,
12,
13,
14,
15,
16]. Similarly,
Tbx18 has also been reported to be expressed in the embryonic epicardium but its role in epicardial formation is controversial. Greulich
et al. [
17] reported that
Tbx18 is dispensable for epicardial formation, EMT and subsequent differentiation into smooth muscle cells and fibroblasts while Wu
et al. [
18] demonstrated impaired coronary plexus formation and identified impaired expression of several signaling pathways related to vascular development, such as Hedgehog, Vegf, Angiopoetin and Wnt signaling. Further support for the role of
Tbx18 in epicardial EMT was reported by Takeichi
et al. [
19] using murine primary epicardial cells
. Epicardin/Tcf21 is also discretely expressed in the PE and EE [
20,
21] and regulates the specification and maturation of the proepicardial cells [
22] in
Xenopus, while in mice it has been reported to be essential for epicardial derived fibroblast formation [
23] and the inhibition of smooth muscle differentiation [
24]. Other cardiac transcription factors with pivotal roles in cardiogenesis, including
Gata4,
Nkx2.5,
Isl1 and
Pitx2 have been reported during PE/EE formation [
25,
26], yet their functional contribution is still uncertain, except
Gata4 that is essential for PE formation [
27].
For decades, the adult epicardium was considered merely an external cardiac lining with no or limited physiological implications. However, it has become clearly demonstrated that upon injury, the epicardium is reactivated, providing essential cues to regenerate the damaged heart in experimental models such as zebrafish, medaka and mice [
28,
29,
30,
31,
32,
33]. While the direct epicardial cellular contribution is limited [
34], adult thymosin β4 primed epicardial cells can be converted into fully functionally integrated cardiomyocytes after myocardial infarction [
35,
36,
37], supporting the therapeutic potential of the epicardium to heal the damaged heart [
38,
39,
40,
41]. Furthermore, EPDCs play determinant roles in distinct cardiac pathological conditions, contributing to an increased fibrous response after myocardial infarction [
42] as well as providing cellular substrates for atrial fat deposition [
43] that correlates with an enhanced prevalence of arrhythmogenic diseases such atrial fibrillation [
44,
45,
46].
Over the last decade, we have witnessed the rise of a novel layer of complexity governing gene regulation with the discovery of non-coding RNAs [
47,
48]. Non-coding RNAs are broadly classified according to their length into small non-coding RNAs (<200 nt) and long non-coding RNAs (>200 nt) [
47,
48]. microRNAs constitute the most abundantly expressed and widely studied types of RNAs among small non-coding RNAs. microRNAs are 20-24 nucleotides in length molecules that play essential roles controlling post-transcriptional regulation of coding RNAs by base pair complementary binding within the 3´UTRs and promoting mRNA degradation and/or protein translation blockage [
47]. On the other hand, long non-coding RNAs displayed a similar biogenesis as coding RNAs but despite their length, no protein coding potential is identified [
48]. Several lncRNAs have been reported to exert pivotal roles during early embryonic development, such as
Fendrr and
Braveheart [
49,
50], but none of them have been identified during epicardial formation. Furthermore, some lncRNAs have been identified to be distinctly expressed in cardiac pathological conditions such as myocardial ischemia, heart failure or arrhythmogenic diseases, respectively [
51,
52,
53,
54].
Recently, we have identified three distinct lncRNAs, adjacent to the
Wt1,
Bmp4 and
Fgf8 chicken gene loci, with enhanced expression in the PE as compared to the embryonic myocardium in the developing embryonic chicken heart. These lncRNAs are distinctly regulated by Bmp, Fgf and thymosin β4, providing support of their plausible implication in epicardial cell lineage specification [
55]. Within this study, we therefore provided thorough characterization of their tissue and subcellular distribution, their transcriptional regulation, their conservation in mice and their functional characterization. We provide herein evidence that these lncRNAs are widely distributed and transcriptional regulated by cardiac enriched transcription factors. Furthermore, mouse
Gm14014, but not its chicken homologue
Wt1_76127, is required for epicardial cell migration.
3. Discussion
LncRNAs represent a large family of non-coding RNAs with a wide cellular and tissue distribution. We have previously characterized three newly identified lncRNAs in the chicken genome that display enhanced expression in the developing PE as compared to the embryonic ventricle [
55]. We now further elaborated into distinct characteristics of these lncRNAs, including a more detailed analysis of their tissue distribution, transcriptional regulation, functional role and evolutionary conservation in mice. Our data demonstrated that these lncRNAs displayed a broad tissue distribution, in line with other reports demonstrating a global expression of lncRNAs in different tissues [
87,
88,
89]. In addition, we also demonstrated that these lncRNAs are preferentially nuclearly located, supporting in first instance a more plausible role in transcriptional and/or epigenetic regulation [
68]. Intriguingly, their subcellular localization is modulated as embryonic development advances in a tissue-specific manner, i.e. more prominent in the embryonic ventricles as compared to the PE/EE. Thus, these observations support the notion of a dynamic lncRNA subcellular localization that might reflect changes in their maturation and/or biological function [
90], a process that seems to be distinctly regulated by transcription factors such as
Mef2c, as demonstrated herein.
Cardiac-enriched transcription factors such as
Mef2c [
91,
92,
93],
Nkx2.5 [
94,
95,
96]
Srf [
97,
98] and
Pitx2 [
99,
100,
101,
102] play pivotal roles in cardiogenesis as well as in proepicardial/embryonic epicardial development [
103,
104,
105,
106] by transcriptionally regulating multiple downstream targets, particularly in heart development. While ample evidence is reported on transcriptional regulation of protein-coding RNA [
107,
108,
109,
110], scarce evidence is available regarding non-coding RNAs, particularly lncRNAs [
62,
111,
112]. We provide herein evidence that
Mef2c exerts a robust transcriptional repression in embryonic ventricle but not in the epicardium (except for
Wt1_76127) while such repressive signals are mostly mediated by
Nkx2.5 in the epicardium. On the other hand,
Pitx2c and
Srf are essential for proper lncRNAs expression in the epicardium while mostly dispensable in the embryonic ventricle. Overall these data reinforce the notion that lncRNA transcriptional expression is modulated by distinct cardiac-enriched transcription factors, as previously reported for other lncRNAs [
62,
111] in a tissue-specific manner.
Long non-coding RNAs are poorly conserved non-coding RNA molecules [
48,
71], in contrast to microRNAs [
47]. We could only identify a mouse homologue for
Wt1_76127 in the mouse genome, given that is located in the syntenic genomic locus,
i.e. the vicinity of
Wt1 transcription factor and shares a 39,1% nucleotide homology. Analyses of its embryonic and adult tissue distribution as well as its subcellular distribution within distinct cardiovascular cellular types display a rather similar profile as its chicken homologue, further supporting its evolutionary conservation.
On the other hand, transcriptional regulation of
Gm14014 displays a mixture of similar control mechanisms as the chicken
Wt1_76127 homologue, as
Nkx2.5 and
Srf are completely dispensable for their expression in myocardial cells while
Mef2c and
Pitx2c exert repressive and activation capacities on epicardial cells. On the other hand, divergent roles for
Mef2c and
Pitx2c are observed in myocardial cells, displaying opposite transcriptional regulation in chicken vs mouse. Thus, these data support the notion that cardiac-enriched transcription factors can distinctly and selectively regulate homologue lncRNA expression in a cell-specific manner, as previously reported for other lncRNAs [
62,
111,
113]. Furthermore, transcriptional regulation of lncRNAs seems to rapidly evolved among different species, as previously reported [
48,
71], supporting the notion that they might even exert divergent functional roles in chicken
vs mice.
Cardiac injury, such as myocardial infarction, in the adult mammalian heart is unable to compensate its myocardial loss. However, during a transient neonatal period,
i.e. first week, complete healing of the heart is achieved in mice [
114]. Several injury models such as apical ventricular resection, cryoinjury and left coronary artery occlusion have been reported, displaying all of them a characteristic healing pattern [
115,
116,
117]; an early inflammatory response is firstly activated, followed by a fibrous extracellular matrix deposition phase that allows to maintain the cardiac structural and functional integrity that is finally followed by a resolution phase that promotes myocardial cell division and fibrous scar clearance (see for a recent review; [
118]). Thus, multiple cell types are required to heal the damaged heart, including the epicardial and myocardial layers, as well as fibroblasts and macrophages. Understanding of the molecular mechanisms driving cardiac regeneration has greatly advanced over the last decade [
119,
120], including the functional role of non-coding RNAs such as microRNAs [
121,
122,
123] and lncRNAs [
124,
125,
126,
127,
128], particularly with emphasis on cardiomyocyte recovery. In our study we explored the plausible contribution of
Gm14014 during cardiac injury by generating embryonic and postnatal
ex vivo cryoinjury models.
Gm14014 upregulation is observed in all embryonic stages tested as early as E10.5 when the embryonic heart is formed by epicardial, myocardial and endocardial cells in the absence of a structured fibroskeleton and resident immune cells, supporting the notion that it should be upregulated in one of these three layers. On the other hand,
Gm14014 is down-regulated at P1 stage, when all previously cited cell types are already present. Curiously, activation of distinct signaling pathways including epicardial activation (
Wt1, Tbx18, Tcf21), and epithelial-to-mesenchymal transition (EMT) (
Snail1, Snail2, Cdh5) is prominently observed at P1 cryoinjured hearts but not at early embryonic stages (E10.5 and E13.5) supporting a negative correlation between
Gm14014 up-regulation and the onset of these biological processes. Such negative correlation is also observed at P7, since those genes involved in the signaling pathways mentioned above are also down-regulated while
Gm14014 is upregulated, similarly as in the embryonic stages. Therefore, these data support a plausible role of
Gm14014 in cardiac regeneration provided its complementary pattern after cryoinjury in P1 (regenerative) and P7 (non-regenerate) stages [
82].
LncRNAs can modulate multiple biological processes, including epigenetic, transcriptional and post-transcriptional processes [
48,
71,
90]. Dissecting their subcellular localization and searching for protein partners can provide further evidence of their functional roles [
48,
71,
90]. We demonstrated herein by RT-qPCR as well as by SCRINSHOT that
Gm14014 is widely expressed in different tissues and cell types. Within the cardiovascular context, its expression is prominently observed in MEC1 epicardial and MEVEC endocardial cells as compared to HL1 myocardial cells and furthermore its subcellular localization in all these cell types is prominently nuclear. Furthermore, functional
in vitro analyses demonstrated that
Gm14014 is essential for epicardial cell migration, but not for endocardial and myocardial migration, supporting a cell-type specific role. While a large body of evidence is currently reported on the role of lncRNA in oncogenic cell migration and metastasis (see for recent reviews; [
129,
130,
131]), this is, to the best of our knowledge the first evidence of the functional role of a lncRNA in homeostatic cell migration, with implication for embryonic epicardial development.
Mechanistically, the identification of
Gm14014-associated proteins by MS revealed a similar number of nuclear, cytoplasmic and nuclear/cytoplasmic proteins. Silencing assays demonstrate that several cytoskeletal proteins such as Myh9 and Myl9 are severely impaired, supporting a plausible molecular link to impaired epicardial cell migration [
85,
132,
133,
134]. However, it remains to be reconciled why
Gm14014 is primarily localized in the nucleus but cytoskeletal reorganization is observed. A plausible explanation is that immature
Gm14014 transcripts are retained into the nucleus and only a subset is selectively translocated towards the cytoplasm exerting its functional role [
90], a process that might be developmentally and transcriptionally dynamic. This is indeed in line with the fact that chicken
Gm14014 homologue lncRNA subcellular expression,
i.e. Wt1_76127, during epicardial and myocardial maturation, progressively becomes more abundantly expressed in the cytoplasm as development proceeds. Similar observation has also been reported for other lncRNAs during development [
66,
105].
Alternatively,
Gm14014 can interact with nuclearly located proteins that transcriptionally control the expression of key master genes regulating actin cytoskeletal proteins, since also multiple nuclear proteins were detected in our MS analyses. We tested whether Rac1 might be this linking factor, but failed to demonstrate
Gm14014 physical interaction nor immunocytochemical differences. Thus, it might be plausible that
Gm14014 might be interacting with proteins that are both nuclear and cytoplasmic. We provided evidence that
Gm14014 physically interacts with Myl9. Furthermore, we demonstrated that inhibition of
Gm14014 selectively translocates Myl9 to the cytoplasm and additionally diminished
Myh9 expression, a downregulation that is selectively observed only in MEC1 epicardial cells where migration is also halted by
Gm14014 silencing but not in other cell types. Therefore, these observation supports a working model by which Gm14014 physically interacts with Myl9 in the nucleus, facilitating
Myh9 transcription and thus modulating cell migration as previously reported [
85] (
Figure 6). In absence of
Gm14014, Myl9 translocates to the cytoplasm and thus does no longer promotes
Myh9 transcription. As a consequence,
Myh9 is downregulated in the cytoplasm and furthermore, leading to cytoskeleton remodeling and thus halting cell migration (
Figure 6). Surprisingly, such molecular mechanisms are not observed in the chicken
Gm14014 homologue,
i.e. Wt1_76127, demonstrating divergent functional roles.
In summary, we provide herein evidence that chicken Wt1_76127 lncRNA and its murine homologue, i.e. Gm14014, display similar tissue distribution profile in embryonic and adult stages. However, transcriptional regulation by key cardiac enriched transcription factors display conserved and divergent profiles. Additionally, we demonstrated that murine Gm14014, but not chicken Wt1_76127 is essential for epicardial, but not myocardial, cell migration, a process that is modulated by Gm14014-Myl9 physical interaction and subsequent nuclear-to-cytoplasmic translocation and cytoskeletal re-arrangement. Thus, these data support the notion that homologues lncRNAs can exert distinct and species-specific functional capabilities.
5. Methods
5.1. Chicken Embryonic Tissues and Epicardial Explants
Fertilized eggs from white Leghorn chickens (Granja Santa Isabel, Córdoba, Spain) were incubated at 37.5°C and 50% humidity for 2–7 days. Embryos were harvested at different developmental stages (HH17, HH24 and HH32) and classified according to Hamburger and Hamilton [
56]. Embryos were removed from the egg by cutting the blastocyst margin with iridectomy scissors and placed into phosphate-buffered saline (PBS). For RT-qPCR analyses, HH17 PE were dissected out, pooled (n=10), and stored at -80°C until used. Additionally,
in vitro explants cultures, HH24 and HH32 cardiac explants were cultured for 24 hours and subsequently the ventricular and epicardial outgrowths were separated, isolated, pooled (n=10), and stored at -80°C until used.
5.2. Chicken Primary Cultures
HH36 embryonic hearts were isolated and disintegrated with iridectomy scissors and placed in PBS. The tissues were incubated with trypsin for 30 minutes at 37°C and the supernatant was collected. This trypsinization step was repeated until all the tissue was fully dissociated. The supernatant was centrifugated and the pellet was cultured in plastic Petri dishes with fibroblast culture medium (Dulbecco’s Modified Eagle’s Medium-high glucose) (Sigma) supplemented with fetal bovine serum (FBS) 10%. Two pre-plating steps were made to separate cardiac fibroblasts (CFs) and cardiomyocytes (CMs). Subsequently, CMs from supernatant was cultured in a plastic Petri dish with fresh CMs culture medium supplemented with 0,001M 5-bromo-2-deoxyuridine (Sigma) to inhibit fibroblast proliferation as previously reported [
57].
5.3. Mouse Embryonic Tissues
CD1 mice were bred and embryos were collected at different embryonic developmental stages, including embryonic day (E) E10.5, E13.5 and E19.5. Postnatal day (P1, P7 and P21) hearts were also collected. Pregnant females and neonatal mice were euthanized by cervical dislocation and by decapitation, respectively. Subsequently, embryonic and postnatal tissues were dissected, pooled and stored at -80°C until used. Approved consent of the Andalusian Ethic Committee was obtained prior to the initiation of the study.
5.4. Mouse Explants for Cryoinjury
Embryonic hearts for different developmental stages (E10.5, E13.5 and E19.5) and postnatal hearts (P1 and P7) were collected. Apical ventricular sections were isolated and dissected using dissection forceps and micro-scalpel and cultured in multi-wells plates with DMEM/F-12 + GlutaMAX 1 (1X) (Gibco) as culture medium. Apical ventricular explants were cryoinjured using a needle frozen with liquid nitrogen and were cultivated at different times (6, 12, 24 and 48 hours). The explants were collected and stored at -80°C until used.
5.5. Nucleus/Cytoplasm SUBCELLULAR isolation
Cytoplasmic and nuclear RNA fractions from HH17 PE, HH24 and HH32 epicardial/ myocardial outgrowths as well as from primary chicken cardiomyocytes, HL1 cardiomyocytes (SCC065, Sigma-Aldrich), MEVEC endocardial cells [
58], MEC1 (SCC187, Sigma-Aldrich) and EPIC epicardial cells [
59] were isolated with Cytoplasmic & Nuclear RNA Purification Kit (Norgen, Belmont, CA, USA) following the manufacturer’s instructions. After RNA isolation, RT-qPCR analysis for nuclearly enriched
Xist (isoform 2) mRNA and cytoplasmic
Gapdh mRNA were performed to validate enrichment on each subcellular fraction. RT-qPCR analysis of distinct lncRNAs was subsequently performed as detailed below.
5.6. siRNA Cell Transfections
Chicken ventricular and epicardial explants, HL1 cardiomyocytes, MEVEC endocardial cells, MEC1 and EPIC epicardial cells (6 × 10
5 cells per well) were transfected with Pitx2c-siRNA, Srf-siRNA, Nkx2.5-siRNA, Mef2c-siRNA and Myl9-siRNA (Sigma, Aldrich, Munich, Germany), respectively, as previously described [
60,
61,
62]. siRNA sequences are provided in
Supplementary Table S1.
5.7. ASO Design and Transfection
Antisense oligonucleotides (ASOs) were designed as previously reported [
63]. The structure of ASOs used in this report consist of a main backbone (10 nucleotides) with phosphorothioate groups, and on both sides 5 nucleotides with different methyl groups
. ASO transfections were carried out with Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s guidelines. A concentration of 20nM of
Wt1_76127 and
Gm14014 ASO was applied to the cells at different times, respectively while for
Myl9 ASO the concentration used was 80nM. ASOs sequences are provided in
Supplementary Table S1.
5.8. Cell Viability Assays
Cell viability was analyzed with Apoptosis/Necrosis Assay Kit (Abcam), following the manufacturer’s instructions. Cell cultures were analyzed using a Leica TCS SP5 II confocal scanning laser microscope. Viable cells were detected with Cytocalcein violet 450 reagent (blue), apoptotic cells with Apopxin Deep Red (red) and apoptotic and necrotic cells were detected with Nuclear Green DCS1 reagent (green).
5.9. Cell Migration Assays
Cell migration was analyzed by scratch assay, as described by Ascione
et al. [
64]. Primary chicken cardiomyocytes, HL1 cardiomyocytes, MEVEC endocardial cells and MEC1 epicardial cells were plated on 24-well culture dishes at a density of 6 × 10
5 cells per well and incubated until 90-100% confluence. Cell monolayers were manually scratched with a p200 pipette tip. PBS was used to wash the cells and subsequently replaced by the serum-free medium. The experimental group was treated with lipidic vesicles containing ASOs as cargo, while the control group was treated with empty lipidic vesicles, respectively. All plates were photographed after 24 h. In addition, a time lapse analysis was carried out for the two conditions. Cell monolayers were scratched, transfected and cultured for 24h taking pictures every 15 minutes. A Leica TCS SP5 II confocal microscope was used, maintaining optimal temperature and humidity conditions for cell growth.
5.10. Growth Factors and Thymosin β4 Administration
HL1 cardiomyocytes, MEVEC endocardial cells and MEC1 epicardial cells were treated for 24 h with Bmp2, Bmp4, Fgf2, Fgf8 and thymosin β4 (Prospec, East Brunswick, NJ, USA), respectively, as reported by Dueñas
et al. [
55]. Cells were collected and processed according for RT-qPCR. In all cases, 3–5 independent biological replicates were analyzed.
5.11. Immunohistochemical Analyses
Control and experimental treated cells were collected after the corresponding treatment, rinsed in PBS for 10 min and fixed with 4 % PFA for 10 minutes at room temperature. After fixation, the samples were rinsed three times (10 min each) in PBS at room temperature and then permeabilized with 1% Triton X-100 and NH4Cl 50nM in PBS for 10 min at room temperature. To block nonspecific binding sites, PBS containing gelatin solution 0,2% (Sigma) was applied twice (10 min each). After blocking, the samples were rinsed three times (10 min each) in PBS at room temperature and were immunofluorescently labeled to detect different proteins. Primary antibodies against Actn1 (ab68194, Abcam), Actn4 (ab108198, Abcam), Myh9 (ab75590, Abcam), pHH3 (CA-92590, Milipore), Myl9 (sc-28329, Santa Cruz), Rac1 (sc-514583, Santa Cruz) and MF20 (ATCC) were used, diluted (1:200) in blocking solution, and applied to each culture overnight at 4°C, respectively. Subsequently, the samples were rinsed three times (10 min each) in PBS to remove excess of primary antibody and incubated 30 min at room temperature with Alexa-Fluor 546, 488 and 633 (1:100; Invitrogen) as a secondary antibody, respectively, as corresponding. Finally, the cells were incubated with DAPI (1:2000; Sigma) for 10 min at room temperature and rinsed three times in PBS for 5 min each. Cell cultures were stored in PBS in darkness at 4°C until analyzed using a Leica TCS SP5 II confocal scanning laser microscope.
5.12. SCRINSHOT In Situ Hybridization
SCRINSHOT (
(Single-Cell Resolution IN Situ Hybridization On Tissues) assay was conducted as described by Sountoulidis
et al. [
65]. Heart tissue from mouse embryos (E14.5) and postnatal 21 days (P21) animals were collected, washed in PBS 1X (pH=7.4), embedded in OCT (FSC22 Clear, Leica), and stored at -80°C until used. 5µm thick cryostat sections were obtained (Leica CM3050S). Padlock probes (40-45 nucleotides) were designed using PrimerQuest online tool (IDT: Integrated DNA Technologies). These DNA oligos were used to design the fluorophore detection oligos, replacing 2-3 “T” nucleotides with “U” to subsequently remove these oligos after the detection cycle using the enzyme Uracil-DNA Glycosylase (Thermo, EN0362). The sequences of the padlock probes and the fluorophore detection oligos are provided in
Supplementary Table S1. Images were obtained using a Zeiss Axio Observer Z.2 fluorescent microscope. Image analysis was based on the nuclear segmentation and alignment of the different detection cycles for the genes of interest. Colocalization was obtained as the number of positive cells for the gene of interest relative to the total number of cells in each sample. Further image processing was carried out using FIJI, Cell Profiler, RStudio and TissUUmaps software.
5.13. RNA Isolation and cDNA Synthesis
Total RNA was isolated Direct-Zol RNA Miniprep Kit (Zymo Research), according to the manufacturer’s instructions. In all cases, at least three distinct pooled samples were used to perform the corresponding RT-qPCR experiments. For mRNA expression measurements, reverse transcription Maxima First Strand cDNA Synthesis Kit for RT-qPCR (Thermo Scientific) was used, according to the manufacturer’s guidelines. Negative controls to assess genomic contamination were performed for each sample, without reverse transcriptase, which resulted in all cases in no detectable amplification product.
5.14. qPCR Analyses (mRNA and lncRNA)
Real-time PCR experiments were performed with 2μL of diluted cDNA, GoTaq® qPCR Master Mix (Promega) and corresponding primer sets described in
Supplementary Table S1. All qPCRs were performed using a CFX384TM thermocycler (Bio-Rad) following the manufacturer’s recommendations. The relative level of expression of each gene was calculated as described by Livak & Schmittgen [
66] using
Gapdh as an internal control for mRNA expression analyses. All primers were designed to span exon-exon boundaries using the online Primer3 software Primer3Plus (
https://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi). Each PCR reaction was carried out in triplicate and repeated in at least three distinct biological samples to obtain representative means. No amplifications were observed in PCR control reactions containing only water as template.
5.15. LncRNA Pull Down Assays
Pull down of biotinylated RNA was carried out as described by Panda
et al. [
67]. Biotinylated RNA of exon 1 and exon 2 of
Gm14014 and
Gapdh were synthesized from PCR fragments using specific forward primers that contained the T7 RNA polymerase promoter sequence [(T7) AGTAATACGACTCACTATAGGG]. Seven fragments were obtained for the sequence of
Gm14014 and one fragment for
Gapdh. The fragments were biotinylated with Biotin-14-CTP (Thermo Fisher Scientific, Invitrogen) and transcribed with MaxiScript T7 kit (Thermo Fisher Scientific, Invitrogen). Whole-cell lysates (500 μg) from MEC1 cells were incubated with 1 μg of biotinylated RNA (biotinylated
Gm14014 and
Gapdh samples) for 2h at room temperature. An input sample incubated only with cell lysate was included as a negative control. Complexes were isolated with Streptavidin-coupled Dynabeads (Invitrogen) and analyzed by mass spectrometry (MS). DAVID database and Gene Ontology analyses were subsequently performed on the resulting proteomic data.
5.16. Western Blot
Western blot (WB) was performed with 10% of the pull-down lysate to validate the interaction between Gm14014, Myl9 and Rac1. The primary antibodies Myl9 (sc-28329, Santa Cruz) and Rac1 (sc-514583, Santa Cruz) were used at a concentration of 1:100 incubated for 5h at room temperature and the secondary antibody-HRP conjugate (170-6516, Biorad) at 1:5000 for 30 min at room temperature. Blocking was carried out with albumin and washes with PBST, according to the antibody manufacturer's recommendations.
5.17. Statistical Analyses
For statistical analyses of datasets, unpaired Student’s t-tests were used. Significance levels or P values are stated in each corresponding figure legend. P < 0.05 was considered statistically significant.
Figure 1.
Panel A. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH24 embryonic tissues, demonstrating a high expression levels in the limb buds as compared to the heart. Panel B. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH32 embryonic tissues, demonstrating a high expression levels in the limb buds, body wall and liver as compared to the heart and eye. Panel C. Schematic representation of the chicken ventricular explants. Panel D. RT-qPCR analyses of the subcellular nuclear and cytoplasmic distribution of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH17, HH24 and HH32 proepicardial (PE; HH17), embryonic epicardial (EE, HH24 and HH32) and ventricular (V, HH17, HH24 and HH32) tissues. Observed that all three lncRNAs are prominently nuclear but their cytoplasmic expression increases as development proceeds, particularly in the ventricular tissues. Panel E. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 expression in epicardial cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively. Panel F. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 expression in ventricular cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively.
Figure 1.
Panel A. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH24 embryonic tissues, demonstrating a high expression levels in the limb buds as compared to the heart. Panel B. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH32 embryonic tissues, demonstrating a high expression levels in the limb buds, body wall and liver as compared to the heart and eye. Panel C. Schematic representation of the chicken ventricular explants. Panel D. RT-qPCR analyses of the subcellular nuclear and cytoplasmic distribution of Wt1_76127, Bmp4_53170 and Fgf8_57126 in HH17, HH24 and HH32 proepicardial (PE; HH17), embryonic epicardial (EE, HH24 and HH32) and ventricular (V, HH17, HH24 and HH32) tissues. Observed that all three lncRNAs are prominently nuclear but their cytoplasmic expression increases as development proceeds, particularly in the ventricular tissues. Panel E. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 expression in epicardial cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively. Panel F. RT-qPCR analyses of Wt1_76127, Bmp4_53170 and Fgf8_57126 expression in ventricular cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively.
Figure 2.
Panel A. RT-qPCR analyses of Gm14014 expression in embryonic (E13.5) mouse tissues, demonstrating wide expression in different tissues, being higher in the lung, stomach and liver as compared to other tissues, e.g. heart. Panel B. RT-qPCR analyses of Gm14014 expression in adult (3 months old) mouse tissues, demonstrating wide expression in different tissues, being higher in the lung, stomach and liver as compared to other tissues, e.g. heart. Panel C. SCRINSHOT in situ hybridization analyses of Gm14014 expression in E14.5 mouse hearts, panel C.1 and C.2 are close ups corresponding to atrial and ventricular areas, respectively. Panel D. SCRINSHOT in situ hybridization analyses of Gm14014 expression in 21 days postnatal (P21) mouse hearts, panel D.1 and D.2 are close ups corresponding to atrial and ventricular areas, respectively. Panel E. SCRINSHOT in situ hybridization analyses of Gm14014, Acta2, Col1a2 and Sox9 expression colocalization in E14.5 and P21 mouse hearts, corresponding to ventricular closeups, respectively. Panel F. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells. Observe higher expression levels in MEVEC endocardial cells as compared to MEC1 and EPIC epicardial cells, with lower expression levels in HL1 cardiomyocytes. Panel G. RT-qPCR analyses of the subcellular nuclear and cytoplasmic distribution of Gm14014 in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells. Note that in all cases, Gm14014 is prominently nuclear. Panel H. RT-qPCR analyses of Gm14014 expression in MEC1 epicardial cells, MEVEC endocardial cells and HL1 cardiomyocytes after selective administration of Fgf2, Fgf8, Bmp2, Bmp4 and thymosin β4, respectively. Panel I. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively. Panel J. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells after selective overexpression of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors, respectively.
Figure 2.
Panel A. RT-qPCR analyses of Gm14014 expression in embryonic (E13.5) mouse tissues, demonstrating wide expression in different tissues, being higher in the lung, stomach and liver as compared to other tissues, e.g. heart. Panel B. RT-qPCR analyses of Gm14014 expression in adult (3 months old) mouse tissues, demonstrating wide expression in different tissues, being higher in the lung, stomach and liver as compared to other tissues, e.g. heart. Panel C. SCRINSHOT in situ hybridization analyses of Gm14014 expression in E14.5 mouse hearts, panel C.1 and C.2 are close ups corresponding to atrial and ventricular areas, respectively. Panel D. SCRINSHOT in situ hybridization analyses of Gm14014 expression in 21 days postnatal (P21) mouse hearts, panel D.1 and D.2 are close ups corresponding to atrial and ventricular areas, respectively. Panel E. SCRINSHOT in situ hybridization analyses of Gm14014, Acta2, Col1a2 and Sox9 expression colocalization in E14.5 and P21 mouse hearts, corresponding to ventricular closeups, respectively. Panel F. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells. Observe higher expression levels in MEVEC endocardial cells as compared to MEC1 and EPIC epicardial cells, with lower expression levels in HL1 cardiomyocytes. Panel G. RT-qPCR analyses of the subcellular nuclear and cytoplasmic distribution of Gm14014 in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells. Note that in all cases, Gm14014 is prominently nuclear. Panel H. RT-qPCR analyses of Gm14014 expression in MEC1 epicardial cells, MEVEC endocardial cells and HL1 cardiomyocytes after selective administration of Fgf2, Fgf8, Bmp2, Bmp4 and thymosin β4, respectively. Panel I. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells after selective inhibition of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors by siRNA administration, respectively. Panel J. RT-qPCR analyses of Gm14014 expression in HL1 cardiomyocytes, MEC1 and EPIC epicardial cells, and MEVEC endocardial cells after selective overexpression of Mef2c, Nkx2.5, Pitx2c and Srf transcription factors, respectively.
Figure 3.
Panel A. Schematic representation of the mouse ventricular explants for ex vivo cryoinjury. Panel B. RT-qPCR analyses of Gm14014 expression in E10.5, E13.5, E19.5 and P1 mouse ventricular explant at 6h, 12h, 24h and 48h after cryoinjury, respectively. Observe that Gm14014 progressively increases over time at embryonic stages while it decreases at postnatal stage. Panel C. Heatmap representation of the expression levels of myocardial (Tnnt2), epicardial (Wt1, Tbx18, Tcf21) and EMT (Snail1, Snail2, Cdh5) markers in E10.5, E13.5, E19.5, P1 and P7 mouse ventricular explants at 24h after cryoinjury, respectively.
Figure 3.
Panel A. Schematic representation of the mouse ventricular explants for ex vivo cryoinjury. Panel B. RT-qPCR analyses of Gm14014 expression in E10.5, E13.5, E19.5 and P1 mouse ventricular explant at 6h, 12h, 24h and 48h after cryoinjury, respectively. Observe that Gm14014 progressively increases over time at embryonic stages while it decreases at postnatal stage. Panel C. Heatmap representation of the expression levels of myocardial (Tnnt2), epicardial (Wt1, Tbx18, Tcf21) and EMT (Snail1, Snail2, Cdh5) markers in E10.5, E13.5, E19.5, P1 and P7 mouse ventricular explants at 24h after cryoinjury, respectively.
Figure 4.
Panel A. RT-qPCR analyses of Gm14014 expression in MEC1 epicardial cells after Gm14014 ASO administration at 6h, 12h, 24h and 48h after transfection. Observe that selective downregulation is achieved at 6h,12h and 48h while at 24h there is a significant upregulation. Panel B. Schematic representation of wound healing scratch assay in MEC1 epicardial cells in controls, Gm14014 ASO, siMyl9 and Myl9 ASO treated cells, respectively, HL1 cardiomyocytes and MEVEC endocardial cells at 6h, 8h, 12h and 24h in controls and Gm14014 ASO treated cells, respectively. Representative images at t=0 and t=24h in MEC1 epicardial cells. Note that Gm14014 ASO treated cells are significantly decreased in migration in MEC1 epicardial cells but not in HL1 cardiomyocytes and MEVEC endocardial cells. siMyl9 and Myl9 ASO treated cells are significantly decreased in migration in MEC1 epicardial cells. Panel C. Representative images of proliferation assays as revealed by phospho-histone 3 (pHH3) immunocytochemistry in control and scratched MEC1 epicardial cells corresponding to control and Gm14014 ASO conditions. Quantitative analyses are also shown demonstrating no significant differences in cell proliferation. Panel D. Graphical representation of lineal vs non-lineal cell migration in time-lapse confocal image analyses of control and scratched MEC1 epicardial cells corresponding to control and Gm14014 ASO conditions. Panel E. Representative images of immunocytochemical analyses of Actn1, Actn4, Myh9, Myl9, Rac1 at 6h and of Myh9, Myl9, Rac1 at 24h after administration of Gm14014 ASO as compared to controls. Observe that there is significant difference in the expression of Myh9 and Myl9 at 6h but not at 24h. Note also that Myl9 displays both nuclear and cytoplasmic distribution in controls while in Gm14014 treated cell is exclusively cytoplasmic. Panel F. Quantitative analysis and representative blot of Myl9 after Gm14014 PD. Observe that Myl9 protein interacts with Gm14014.
Figure 4.
Panel A. RT-qPCR analyses of Gm14014 expression in MEC1 epicardial cells after Gm14014 ASO administration at 6h, 12h, 24h and 48h after transfection. Observe that selective downregulation is achieved at 6h,12h and 48h while at 24h there is a significant upregulation. Panel B. Schematic representation of wound healing scratch assay in MEC1 epicardial cells in controls, Gm14014 ASO, siMyl9 and Myl9 ASO treated cells, respectively, HL1 cardiomyocytes and MEVEC endocardial cells at 6h, 8h, 12h and 24h in controls and Gm14014 ASO treated cells, respectively. Representative images at t=0 and t=24h in MEC1 epicardial cells. Note that Gm14014 ASO treated cells are significantly decreased in migration in MEC1 epicardial cells but not in HL1 cardiomyocytes and MEVEC endocardial cells. siMyl9 and Myl9 ASO treated cells are significantly decreased in migration in MEC1 epicardial cells. Panel C. Representative images of proliferation assays as revealed by phospho-histone 3 (pHH3) immunocytochemistry in control and scratched MEC1 epicardial cells corresponding to control and Gm14014 ASO conditions. Quantitative analyses are also shown demonstrating no significant differences in cell proliferation. Panel D. Graphical representation of lineal vs non-lineal cell migration in time-lapse confocal image analyses of control and scratched MEC1 epicardial cells corresponding to control and Gm14014 ASO conditions. Panel E. Representative images of immunocytochemical analyses of Actn1, Actn4, Myh9, Myl9, Rac1 at 6h and of Myh9, Myl9, Rac1 at 24h after administration of Gm14014 ASO as compared to controls. Observe that there is significant difference in the expression of Myh9 and Myl9 at 6h but not at 24h. Note also that Myl9 displays both nuclear and cytoplasmic distribution in controls while in Gm14014 treated cell is exclusively cytoplasmic. Panel F. Quantitative analysis and representative blot of Myl9 after Gm14014 PD. Observe that Myl9 protein interacts with Gm14014.
Figure 5.
Panel A. RT-qPCR expression analyses of Wt1_76127 in primary cultures of chicken embryonic cardiomyocytes (CMs) and cardiac fibroblasts (CFs). Panel B. Subcellular distribution of Wt1_76127 in the nucleus and cytoplasm primary cultures of chicken embryonic cardiomyocytes. Panel C. RT-qPCR expression analyses of Wt1_76127 in embryonic cardiomyocytes and embryonic epicardium (EE) in control and Wt1_76127 ASO treated conditions. Observe a significant downregulation in Wt1_76127 after ASO administration. Panel D. Quantitative analyses of cell viability, apoptosis and necrosis in primary cultures of chicken embryonic cardiomyocytes (CMs) of control and Wt1_76127 ASO treated conditions. Representative images are also depicted. Note not significant differences for cell viability and necrosis, whereas there are significant differences in apoptosis. Panel E Immunocytochemical characterization of primary culture of embryonic cardiomyocytes with MF20 antibody, demonstrating that most of these cell display sarcomeric myosins. Panel F. Immunocytochemical characterization of proliferation in scratched primary culture of embryonic cardiomyocytes. Observe that there are no significant differences. Panel G. Schematic representation of wound healing scratch assay in primary culture of embryonic cardiomyocytes at 6h, 8h, 12h and 24h in controls and Wt1_76127 ASO treated cells, respectively. Panel H. Quantitative analyses of epicardial cell migration in ventricular explant in control and Wt1_76127 ASO treated conditions. Observe that there are no significant differences. Panel I. Representative images of ventricular explants immunostained against Actn1, Actn4 and Myh9 in control and Wt1_76127 ASO treated conditions. Observe that there are no significant differences. Panel J. Representative images of primary culture of embryonic cardiomyocytes immunostained against Actn1, Actn4 and Myh9 in control and Wt1_76127 ASO treated conditions. Observe that and Myh9, but not Actn4 display significant differences.
Figure 5.
Panel A. RT-qPCR expression analyses of Wt1_76127 in primary cultures of chicken embryonic cardiomyocytes (CMs) and cardiac fibroblasts (CFs). Panel B. Subcellular distribution of Wt1_76127 in the nucleus and cytoplasm primary cultures of chicken embryonic cardiomyocytes. Panel C. RT-qPCR expression analyses of Wt1_76127 in embryonic cardiomyocytes and embryonic epicardium (EE) in control and Wt1_76127 ASO treated conditions. Observe a significant downregulation in Wt1_76127 after ASO administration. Panel D. Quantitative analyses of cell viability, apoptosis and necrosis in primary cultures of chicken embryonic cardiomyocytes (CMs) of control and Wt1_76127 ASO treated conditions. Representative images are also depicted. Note not significant differences for cell viability and necrosis, whereas there are significant differences in apoptosis. Panel E Immunocytochemical characterization of primary culture of embryonic cardiomyocytes with MF20 antibody, demonstrating that most of these cell display sarcomeric myosins. Panel F. Immunocytochemical characterization of proliferation in scratched primary culture of embryonic cardiomyocytes. Observe that there are no significant differences. Panel G. Schematic representation of wound healing scratch assay in primary culture of embryonic cardiomyocytes at 6h, 8h, 12h and 24h in controls and Wt1_76127 ASO treated cells, respectively. Panel H. Quantitative analyses of epicardial cell migration in ventricular explant in control and Wt1_76127 ASO treated conditions. Observe that there are no significant differences. Panel I. Representative images of ventricular explants immunostained against Actn1, Actn4 and Myh9 in control and Wt1_76127 ASO treated conditions. Observe that there are no significant differences. Panel J. Representative images of primary culture of embryonic cardiomyocytes immunostained against Actn1, Actn4 and Myh9 in control and Wt1_76127 ASO treated conditions. Observe that and Myh9, but not Actn4 display significant differences.
Figure 6.
Schematic representation of the working model of Gm14014 during epicardial cell migration. Panel A illustrate the role of Gm14014 in homeostasis and Panel B in Gm14014 silencing assays. Note that Gm14014 is bound to Myl9 in the nucleus, providing clues to cell migration by interacting with Myh9 in the cytoplasm, while if Gm14014 is knockdown, no Myl9 nuclear binding occurs and thus is completely translocated to the cytoplasm.
Figure 6.
Schematic representation of the working model of Gm14014 during epicardial cell migration. Panel A illustrate the role of Gm14014 in homeostasis and Panel B in Gm14014 silencing assays. Note that Gm14014 is bound to Myl9 in the nucleus, providing clues to cell migration by interacting with Myh9 in the cytoplasm, while if Gm14014 is knockdown, no Myl9 nuclear binding occurs and thus is completely translocated to the cytoplasm.