Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 10 November 2021
Sec. Cardiovascular Genetics and Systems Medicine
This article is part of the Research Topic Insights in Cardiovascular Genetics and Systems Medicine: 2021 View all 8 articles

Expression Signatures of Long Noncoding RNAs in Left Ventricular Noncompaction

\nQingshan TianQingshan Tian1Hanxiao Niu,Hanxiao Niu1,2Dingyang Liu,Dingyang Liu1,3Na Ta,Na Ta1,4Qing YangQing Yang1Vikram NortonVikram Norton5Yujing Wu,Yujing Wu1,6Amit K. MaitiAmit K. Maiti7Hao Wu
Hao Wu5*Zhenzhong Zheng
Zhenzhong Zheng1*
  • 1Department of Cardiology, The First Affiliated Hospital of Nanchang University, Nanchang, China
  • 2Department of Functional Examination, Shaanxi Provincial People's Hospital, Xi'an, China
  • 3Department of Cardiology, Jiangxi Provincial People's Hospital, Nanchang, China
  • 4Postoperative Cardiac Intensive Care Unit, Fuwai Hospital Chinese Academy of Medical Sciences, Shenzhen, China
  • 5Vascular Biology Program, Department of Surgery, Harvard Medical School, Boston Children's Hospital, Boston, MA, United States
  • 6Department of Emergency, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China
  • 7Mydnavar, Department of Genetics and Genomics, Troy, MI, United States

Long noncoding RNAs have gained widespread attention in recent years for their crucial role in biological regulation. They have been implicated in a range of developmental processes and diseases including cancer, cardiovascular, and neuronal diseases. However, the role of long noncoding RNAs (lncRNAs) in left ventricular noncompaction (LVNC) has not been explored. In this study, we investigated the expression levels of lncRNAs in the blood of LVNC patients and healthy subjects to identify differentially expressed lncRNA that develop LVNC specific biomarkers and targets for developing therapies using biological pathways. We used Agilent Human lncRNA array that contains both updated lncRNAs and mRNAs probes. We identified 1,568 upregulated and 1,141 downregulated (log fold-change > 2.0) lncRNAs that are differentially expressed between LVNC and the control group. Among them, RP11-1100L3.7 and XLOC_002730 are the most upregulated and downregulated lncRNAs. Using quantitative real-time reverse transcription polymerase chain reaction (RT-QPCR), we confirmed the differential expression of three top upregulated and downregulated lncRNAs along with two other randomly picked lncRNAs. Gene Ontology (GO) and KEGG pathways analysis with these differentially expressed lncRNAs provide insight into the cellular pathway leading to LVNC pathogenesis. We also identified 1,066 upregulated and 1,017 downregulated mRNAs. Gene set enrichment analysis (GSEA) showed that G2M, Estrogen, and inflammatory pathways are enriched in differentially expressed genes (DEG). We also identified miRNA targets for these differentially expressed genes. In this study, we first report the use of LncRNA microarray to understand the pathogenesis of LVNC and to identify several lncRNA and genes and their targets as potential biomarkers.

Introduction

Left ventricular noncompaction is a rare disease that is characterized by the failure of densification in the normal embryonic cardiac tissue. It can occur in isolation or association with congenital heart defects (CHDs), neuromuscular disorders, and systemic heart anomalies (1, 2). Morphological features of left ventricular noncompaction (LVNC) consist of deep trabecular recesses in the myocardial wall with noncompaction of the loosely interwoven meshwork in the left ventricular cavity (3). The typical clinical result is a triad of heart failure, arrhythmias, and systemic embolism (4).

Left ventricular noncompaction is a genetically heterogeneous disease (5) and is associated with genes and/or proteins that are involved in sarcomere (MYH7, ACTC, MYBPC3, TNNT2), cytoskeletal (ZASP, LMNA), and mitochondrial structure or function (TAZ) (69). However, using exome and mitochondrial DNA sequencing from myocardial tissue samples, Liu Z et al. identified mutations in 16S rRNA 2336T>C mitochondrial mutation but did not detect any pathogenic mutations in TNNT2 and MYBPC3 genes (10).

Long noncoding RNA molecules are >200 nt in length, similar to mRNA, but do not code for structural proteins. Rather, they modulate the expression of protein-coding genes (11). Particularly, relating lncRNA with disease pathogenesis has widespread implications. Long noncoding RNAs (LncRNAs) are easily targetable, and their expression can be controlled from external oligo DNA or RNA that can be administered through blood or tissue specifically (12). A disease development could be delayed or cured by externally manipulating their expressions. Such advantages recently gained specific interest in identifying disease-specific lncRNA (13, 14).

Numerous studies show that lncRNAs play a positive and negative role in the differentiation, development, and progression of many diseases such as cancer, neurodegenerative diseases, and heart diseases (13, 1517). The role of lnCRNA in cardiac dysfunction is especially beginning to be elucidated (18, 19). However, the role of lncRNA in LVNC has not been reported. This study compared the lncRNA and mRNA expression in the blood of six LVNC patients with six healthy subjects as controls and correlated their association with LVNC. Identification of these differentially expressed lncRNAs and mRNA may provide essential targets to develop blood-based biomarkers and potential therapeutic interventions for treating patients with LVNC.

Methods

Patients and Healthy Subjects

The blood samples of LVNC patients and healthy subjects were obtained from May 2015 to October 2015 in the Department of Cardiology in the First Affiliated Hospital of Nanchang University. Six patients with LVNC were selected for the case-study group and six healthy subjects were selected as the control group. Clinical characteristics of LVNC patients and control subjects are presented in Table 1. Written informed consents were obtained from all patients and the study was approved by the Institutional Review Board (IRB) of the First Affiliated Hospital of Nanchang University. We followed all ethical compliance for the study.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of study populations.

RNA Isolation and Array Hybridization

Blood samples from each patient and healthy subject were centrifuged at 1,000 rpm for 10 min to separate plasma from cells. They were then stored at −80°C. Total RNAs from each sample were prepared from blood cells with the RNA isolation kit (Qiagen, Canada) after removal of rRNA from a Eukaryotic mRNA Isolation Kit, Epicenter. Then, each sample was amplified and transcribed into fluorescent cRNA along the entire length of the transcripts, without 3' bias, utilizing a random priming method (Arraystar Flash RNA Labeling Kit, Arraystar). The labeled cRNAs were purified by RNeasy Mini Kit (Qiagen). The concentration and specific activity of the labeled cRNAs (pmol Cy3/μg cRNA) were measured by NanoDrop ND-1000 (Table 2). One microgram and each labeled cRNA was fragmented by adding 5 μl 10 × Blocking Agent and 1 μl of 25 × Fragmentation Buffer. The mixture was then heated at 60°C for 30 min. Finally, 25 μl 2 × GE Hybridization buffer was added to dilute the labeled cRNA. Fifty microliter of hybridization solution was dispensed into the gasket slide and assembled to the lncRNA expression microarray slide (Agilent lncRNA microarray chip (catalog no.: G4851C, Agilent technology) that contained 26,803 Entrez unique genes and 30,606 unique lncRNAs. Microarray hybridization was performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology) with minor modifications. The slides were incubated for 17 h at 65°C in an Agilent Hybridization Oven. The hybridized arrays were washed, fixed, and scanned using the Agilent DNA Microarray Scanner (part number G2505C).

TABLE 2
www.frontiersin.org

Table 2. The primers for RT-qPCR.

Validation of the Differentially Expressed LncRNAs by Quantitative Real-Time PCR (qRT-PCR)

RT-QPCR was used to validate the microarray lncRNA expression among the six LVNC and six healthy patients. The primer sequences of β-actin (H, ENST00000423402, uc021ybu.1, uc002nqf.1, ENST00000550301, ENST00000417844, and ENST00000554022 were shown in Table 2. The primers used for RT-QPCR were designed and synthesized by Invitrogen (Shanghai, China). β-actin was used as an internal control (housekeeping gene). RT-QPCR was performed using GeneAmp PCR System 9700 (Applied Biosystems). The PCR reaction conditions were an initial denaturation at 95°C for 10 min, followed by 40 PCR cycles at 95°C for 10 s, and 60°C for 60 s. The products were amplified for 10 min by slowly heating them from 60 to 99°C after the extension. The target (lncRNA or gene) and housekeeping genes for each sample were used for Real-time PCR reactions. The target gene/lncRNA concentration of each sample was divided by the concentration of the housekeeping gene to obtain the relative concentration of the target gene for each sample. The amount of target sequence for the sample was derived from the standard curve throughout the experiment.

LncRNA Microarray Statistical Analysis

Agilent Feature Extraction software (version 11.0.1.1) was used to analyze acquired array images. Quantile normalization and subsequent data processing were performed by using the GeneSpring GX v12.1 software package (Agilent Technologies). After quantile normalization of the raw data, lncRNAs and mRNAs which had flags in at least six out of 12 samples in Present or Marginal (“All Targets Value”) were chosen for further data analysis. Differentially expressed lncRNAs and mRNAs with statistical significance between the two groups of healthy and LVNC patients were identified through p-value or FDR filtering. Differentially expressed lncRNAs and mRNAs between the two samples were identified through Fold Change filtering. Hierarchical Clustering and combined analysis were performed and the differences between mRNAs were analyzed to infer their role in cellular pathway using Gene Ontology (GO) and KEGG.

Gene Set Enrichment Analysis (GSEA)

Gene set enrichment analysis (GSEA) is essentially done as described by Subramanian et al. (20). The top 1,500 differentially expressed genes (DEG) are used with the Hallmark Gene Set database. Analysis was done with a “weighted” scoring scheme with “no collapse” as the main parameters. This analysis picked 13 enriched pathways. Among these pathways, the G2M transition, estrogen response, and inflammatory pathways had the highest enrichment score and were studied further. Individual gene expression of GSEA significant genes in these pathways across the LVNC and healthy subjects are performed in R package with “plot” and “bar plot” function.

Determination of Biological Functions and Identification of Targeted MiRNA

Various biological functions and interactome analysis of DEG are determined by the use of FUNRICH software (21). The top 500 DEGs are used as input to carry out the analysis.

Information about biological pathways and target miRNA for these DEGs are deduced by using EnrichR software (2224). The top 500 DEGs are also used as input to identify pathways and miRNA.

Result

Microarray Expression Analysis

Microarray expression analysis identified numerous lncRNAs and mRNAs whose expression varies between LVNC and the control group. In this study, we used a t-test to identify differentially expressed RNAs. According to the Benjamini Hochberg FDR method and the corrected p-value (FDR), differentially expressed lncRNAs and mRNAs were screened with Fold change ≥ 2.0, FDR <0.05. Compared with the p-value filter, this method can better reduce the false positive rate. Two thousand six hundred eighty-two lncRNA (1,568 up and, 1,114 downregulated) and 2,083 mRNA (1,066 up and, 1,017 down) transcripts were identified in serum compared with the control group (log fold-change > 2.0). The top 10 upregulated lncRNA that were significantly downregulated were ENST00000550301, ENST00000433152, uc002nqf.1, ENST00000415393, ENST00000584722, TCONS_00029778, ENST00000443008, TCONS_00024700, NR_003256, ENST00000485225, TCONS_00006918, ENST00000423402, TCONS_00006919, uc022bqs.1, uc021ybu.1, ENST00000537478, and ENST00000447277. The differentially expressed top twenty lncRNA and mRNA were listed in Tables 3, 4. A Scatterplot is used for assessing the lncRNA expression variation between the two groups of LVNC and the control (Figures 1A,B). Hierarchical Cluster shows lncRNA expression patterns (Figure 1C).

TABLE 3
www.frontiersin.org

Table 3. Differentially expressed long noncoding RNA (LncRNA).

TABLE 4
www.frontiersin.org

Table 4. Differentially expressed mRNA.

FIGURE 1
www.frontiersin.org

Figure 1. Using arraystar human long noncoding RNAs (lncRNA) microarray V3.0. Two thousand six hundred eighty-two lncRNAs (log fold- changes > 2.3) were detected. The scatter plot is a visualization method used for assessing the lncRNAs (A) and mRNAs (B) expression variations between left ventricular noncompaction (LVNC) and the control group. The values of the X and Y axes in the scatter plot are the averaged normalized signal values of the group (log2 scaled). The green lines are fold-change lines. (C) Differentially expressed lncRNAs and mRNAs in plasma of LVNC and the control were analyzed using hierarchical clustering. Hierarchical clustering analysis arranges samples into groups by expression level, “red” indicates high relative expression, and “green” indicates low relative expression.

Real-Time Quantitative PCR Validation

By using quantitative real-time reverse-transcription PCR (RT-QPCR), the top three up-regulated lncRNAs and top three downregulated lncRNAs with log fold-changes > 2 were selected to test and verify the microarray data in different blood samples of LVNC patients and paired healthy people. Due to the failure of TCONS_00006918 and TCONS_00006919 primer design, we selected the downregulated expression ENST00000423402, uc022bqs.1, and uc021ybu.1 for PCR verification. The downregulated expression of ENST00000423402, uc021ybu.1, and the upregulated expression of uc002nqf.1 and ENST00000550301 were statistically significant (p <0.05). Two lncRNAs (ENST00000417844 and, ENST00000554022) were randomly selected to prove the consistency of microarray and qPCR. The expression of lncRNA uc002nqf.1, ENST00000550301, ENST00000417844 and, ENST00000554022 was up-regulated, and ENST00000423402 and uc021ybu.1 was downregulated in the LVNC samples vs. the control samples (Figure 2), consistent with the microarray results. Thus, microarray data profiling indicated a series of lncRNAs that were constantly differentially expressed between LVNC blood and healthy blood.

FIGURE 2
www.frontiersin.org

Figure 2. RT-QPCR validation. RT-QPCR was used to confirm the expression of six long non-coding RNAs, ENST00000423402, uc021ybu.1, uc002nqf.1, ENST00000550301, ENST00000417844, and ENST00000554022, in the LVNC and control groups. β-actin was used as an internal control. Values are presented as the mean ± standard error of the mean. *p <0.05, **p <0.01, ***p <0.001.

GO Analysis and Pathway Analysis

To establish the involvement of these RNAs in cellular pathways, we used GO databases to analyze their potential biological functions. As previously reported, GO analysis was performed to determine the gene and gene product enrichment, which covered three domains: biological processes (BP), cellular components (CC), and molecular functions (MF) (25).

We found that the highest GO classification targeted by underregulated transcripts were intracellular signal transduction (BP), plasma membrane (CC), and enzyme binding (MF) (Figures 3A–C). Meanwhile, the highest GO classification targeted by the overregulated transcripts were protein modification process (BP), an integral component of the plasma membrane (CC), and protein binding (MF), (Figures 3D–F). In addition, the KEGG pathway analysis demonstrated that the “bile secretion” and “thyroid hormone synthesis” were the most significant pathways containing differentially upregulated mRNAs. Increased bile secretion is shown to have an adverse effect and promotes heart failure in the long run and thus, may be relevant for developing LVNC (26, 27). Thyroid hormone disbalance has long been known to be a risk factor for heart diseases although its role in LVNC is not demonstrated (28, 29). The “neurotrophin signaling pathways” and “GnRH signaling pathways” were also important pathways containing differentially downregulated mRNAs (Figures 4A,B).

FIGURE 3
www.frontiersin.org

Figure 3. Gene ontology (GO) analysis. The top 10 GO terms associated with coding gene functions of upregulated lncRNAs (A–C) and downregulated lncRNAs (D–F) are listed.

FIGURE 4
www.frontiersin.org

Figure 4. Pathway analysis. The top 10 pathways that associated coding gene of up-regulated lncRNAs (A) and downregulated lncRNAs (B) are listed.

GSEA Enriched Pathway

The Hallmark gene set with which GSEA has analyzed covers a wide range of pathways that would capture a set of genes represented in DEG. Our analysis with the top 1,500 genes identified 13 pathways that are enriched significantly. The rank order of genes with top enrichment score and their heatmap is shown in Supplementary Table 1 and Supplementary Figure 1. Among the pathways with the highest enrichment scores are the G2M transition pathway, Estrogen response pathway, and inflammatory pathway (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Gene set enrichment analysis (GSEA) pathways for differentially expressed mRNA in LVNC. (A) Heatmap of the three most enriched pathways showing genes. (B) Enrichment score of these pathways.

The significance of the G2M pathway is relevant to the LVNC since many stress-responsive genes often play a role in heart development and pathogenesis (30). The glutathione-S-Transferase pi 1 (GSTP1) gene is highly polymorphic and involved in many diseases including cancer and heart failure (31). Its role in LVNC could be further investigated. HSPA8 (Hsc-70) is reported to be a titin interacting protein in an LVNC gene-edited mice model. Therefore, this gene could have a critical role in pathogenesis (32). Trophinin associated protein (TROAP) is shown to be involved in myocardial infarction (33), thus, its role in LVNC could also be important. Dystrophin (DMD) is known to play a role in LVNC development and pathogenesis (34), justifying its higher expression in patients (Figure 6). MYB is a nuclear protein and acts as a transcription factor for many cellular pathways in cardiac development (35). Adenylate cyclase (ADCY1) is involved in cardiac development and cardiomyopathy (36). HSPB8 is a heat shock protein and is also involved in cardiac development and cardiomyopathy (37). IFITM1 is an interferon-induced transmembrane protein that helps with the proliferation of neonatal cardiomyocytes in rats (38). TNF-inducing gene (TNFSF15) plays a critical role in neovascularization and inflammation in the heart. Thus, it could also play a critical role in developing LVNC (39).

FIGURE 6
www.frontiersin.org

Figure 6. Differential expression of genes in LVNC and control in GSEA enriched pathways. Expression differences of top three genes in LVNC and healthy subjects in three enriched pathways.

Biological Functions of DEGs

The biological functions of these DEGs are deduced as components, such as biological processes, transcription factors, cellular components, and molecular functions. Major biological processes include ion transport, transcription, purine cleavage, neurogenesis, and cell migration (Figures 7A–D). Interactome analysis identified a major pathway (Figure 8A) that includes CTLA4, CDKN2A, and HSPA8 which denotes cell cycle, apoptosis, autophagy, and cellular remodeling (40).

FIGURE 7
www.frontiersin.org

Figure 7. Identification of biological functions of differentially expressed genes (DEGs). (A–D) FUNRICH identified several critical biological functions ranging from cellular components to transcription factors.

FIGURE 8
www.frontiersin.org

Figure 8. Interactome and miRNA target analysis using FUNRICH and EnrichR. (A) Major pathway includes CTLA4, CDKN2A, and HSPA8 involves cell cycle, apoptosis, and autophagy pathway. (B) Target miRNA analysis by miRTarbase that identifies CISD2 and CDKN2A. (C) Target miRNA analysis by targetmiRNAScan identifies RUNXIT1 and KSR2.

DEGs Modulate a Group of Targeted MiRNAs

Analysis with EnrichR identified several important pathways (Supplementary Figure 2) such as transcription factor (TF) perturbation, epigenetic genes, and more importantly, miRNA. It used the miRtargetBase 2017 and targetScanmicroRNA2017 database and identified miRNAs that are linked to these DEGs (Tables 5, 6). miRTargetBase miRNAs identified CISD2, CDKN2A, etc. (Figure 8B) as target genes, and miRTargetScan identified RUNXT1 and KSR2 (Figure 8C) as target genes. Their role in developing LVNC could be further established.

TABLE 5
www.frontiersin.org

Table 5. Identification of target miRNA of differentially expressed genes (DEGs) (miRtargetBase 2017).

TABLE 6
www.frontiersin.org

Table 6. Identification of target miRNA of DEGs (targetScanmicroRNA2017).

The top significant miRNAs, hsa-miR-3197 and hsa-miR-6780-3p are identified by miRTargetbase and miRTargetScan, respectively. The function of Has-miR-3197 has not been researched extensively though it has been reported to be a potential biomarker for diabetic retinopathy (41). However, its role in cardiovascular disease and LVNC pathogenesis could be investigated. The has-miR-6780-3p is shown to be downregulated in the plasma of patients with increased low-density lipoprotein (LDL) which could have an impact on developing cardiac dysfunction (42).

Discussion

Although noncompaction of the myocardium is seen in association with other congenital cardiac abnormalities, it can also occur as a primary disorder in the absence of other structural heart diseases (3). It remains unclear whether it represents a distinct disease or a morphological trait shared by many phenotypically different cardiomyopathies (43). Samples selected for this study were isolated cases of left ventricular myocardium which appears in a variety of clinical manifestations, including heart failure, thromboembolic events, and tachyarrhythmias. The aim of this study was to determine the differently expressed lncRNA in the blood of LVNC and to provide the detectable biomarkers for clinical diagnosis. These detectable biomarkers may provide information for further research on regulation during the myocardium noncompaction process through LncRNA and coding gene analysis.

Long noncoding RNA (LncRNA) and coding gene analysis showed that one upregulated lncRNA AC068542.1 upregulates the expression of the mRNA-DPP10. DPP10 encodes a single-pass type II membrane protein that is a member of the S9B family in clan SC of the serine proteases. It binds specific voltage-gated potassium channels and alters their expression with biophysical properties. Turnow K. et al. reported that the interaction of DPP10a, expressed in the human atrium with Kv4.3 channels, generates a sustained current component of Ito, which may affect the late repolarization phase of atrial action potentials (44). Thus, the modulated expression of this gene could have an impact on LVNC pathogenesis.

AC068542.1 lncRNA has been implicated in modulating transcription factors and has a distinct epigenetic mark at its transcription start site (TSS). The author suggested an alternative role of its regulation (45). This lncRNA is also implicated in the modulation of miRNA in breast cancer (46), however, its role in LVNC is not clear.

Another upregulated lncRNA RP11-973D8.4 downregulates the adjacent gene CDK2, which encodes a member of a family of serine/threonine protein kinases that participates in the cell cycle regulation. The activity of this protein is especially critical during the G1 to S phase transition. Spencer SL et al. found that cells at the end of mitosis decide whether or not to start the next cell cycle by immediately activating CDK2 or by entering a transient G0-like state by suppressing CDK2 activity (47).

A differentially expressed top mRNA is the METTL1 (methyltransferase 1) gene which has a binding motif of S-S-adenosylmethionine binding motif. Methyltransferase mediates the formation of N(7)-methylguanine in a subset of RNA species, such as tRNAs, mRNAs, and microRNAs (miRNAs) (PubMed:12403464, PubMed:31031084, PubMed:31031083) and catalyzes the formation of N(7)-methylguanine at position 46 (m7G46) in tRNA (48). METTL1 gene promotes VEGFA mRNA translation that has implications in heart development (49).

Gene Ontology and KEGG pathway analysis identified the “bile acid secretion pathway” and “Thyroid hormone synthesis pathway” that are known to have a profound impact on heart diseases. However, the roles of these pathways have not been extensively studied in relation to the development of LVNC (27, 28). Furthermore, GSEA analysis identified HSPA1, DMD, and ADCY1 as having significant relevance in developing heart ailments and which presumably play a critical role in developing LVNC (32, 34, 36).

Finally, we identified DEG-specific target miRNAs that may play a critical role in the pathogenesis of LVNC but they have not been investigated extensively. It is worthwhile to further investigate these miRNAs and their target genes to explore their role in LVNC pathogenesis. To further confirm the roles of these genes in patients with LVNC, animal models and gene editing techniques, such as CRISPR-Cas9, could effectively discover mechanisms underlying this disease and, thus, shed light on potential novel targeting strategies for treating patients with LVNC.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board (IRB) of First Affiliated Hospital of Nanchang University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ZZ and HW: conceptualization and funding acquisition. QT, HN, DL, AM, ZZ, and HW: methodology. QT, HN, DL, NT, QY, VN, YW, AM, ZZ, and HW: validation. QT, HN, DL, NT, VN, AM, ZZ, and HW: writing and editing. All authors contributed to the article and approved the final manuscript.

Funding

This study was supported by a grant from the National Nature Science Foundation of China (No. 30960119, 81260044), Science and Technology Planning Project of Jiangxi Provincial Health Commission (No. 20195129), and Natural Science Foundation of Jiangxi Province (No. 20202ACBL206002).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We gratefully acknowledge the nurses and patients of the First Affiliated Hospital of Nanchang University for their support during the study period.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.763858/full#supplementary-material

References

1. Oechslin E, Jenni R. Left ventricular non-compaction revisited: a distinct phenotype with genetic heterogeneity? Eur Heart J. (2011) 32:1446–56. doi: 10.1093/eurheartj/ehq508

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Stöllberger C, Finsterer J, Blazek G. Left ventricular hypertrabeculation/noncompaction and association with additional cardiac abnormalities and neuromuscular disorders. Am J Cardiol. (2002) 90:899–902. doi: 10.1016/s0002-9149(02)02723-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Jenni R, Oechslin E, Schneider J, Attenhofer Jost C, Kaufmann PA. Echocardiographic and pathoanatomical characteristics of isolated left ventricular non-compaction: a step towards classification as a distinct cardiomyopathy. Heart. (2001) 86:666–71. doi: 10.1136/heart.86.6.666

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Paterick TE, Tajik AJ. Left ventricular noncompaction: a diagnostically challenging cardiomyopathy. Circ J. (2012) 76:1556–62. doi: 10.1253/circj.cj-12-0666

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Captur G, Nihoyannopoulos P. Left ventricular non-compaction: genetic heterogeneity, diagnosis and clinical course. Int J Cardiol. (2010) 140:145–53. doi: 10.1016/j.ijcard.2009.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vatta M, Mohapatra B, Jimenez S, Sanchez X, Faulkner G, Perles Z, et al. Mutations in Cypher/ZASP in patients with dilated cardiomyopathy and left ventricular non-compaction. J Am Coll Cardiol. (2003) 42:2014–27. doi: 10.1016/j.jacc.2003.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Klaassen S, Probst S, Oechslin E, Gerull B, Krings G, Schuler P, et al. Mutations in sarcomere protein genes in left ventricular noncompaction. Circulation. (2008) 117:2893–901. doi: 10.1161/CIRCULATIONAHA.107.746164

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Vermeer AM, van Engelen K, Postma AV, Baars MJ, Christiaans I, De Haij S, et al. Ebstein anomaly associated with left ventricular noncompaction: an autosomal dominant condition that can be caused by mutations in MYH7. Am J Med Genet C Semin Med Genet. (2013) 163C:178–84. doi: 10.1002/ajmg.c.31365

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Monserrat L, Hermida-Prieto M, Fernandez X, Rodríguez I, Dumont C, Cazón L, et al. Mutation in the alpha-cardiac actin gene associated with apical hypertrophic cardiomyopathy, left ventricular non-compaction, and septal defects. Eur Heart J. (2007) 28:1953–61. doi: 10.1093/eurheartj/ehm239

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Liu Z, Song Y, Li D, He X, Li S, Wu B, et al. The novel mitochondrial 16S rRNA 2336T>C mutation is associated with hypertrophic cardiomyopathy. J Med Genet. (2014) 51:176–84. doi: 10.1136/jmedgenet-2013-101818

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kung JT, Colognori D, Lee JT. Long noncoding RNAs: past, present, and future. Genetics. (2013) 193:651–69. doi: 10.1534/genetics.112.146704

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Winkle M, El-Daly SM, Fabbri M, Calin GA. Noncoding RNA therapeutics - challenges and potential solutions. Nat Rev Drug Discov. (2021) 20:629–51. doi: 10.1038/s41573-021-00219-z

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang Y, Sun X. The functions of LncRNA in the heart. Diabetes Res Clin Pract. (2020) 168:108249. doi: 10.1016/j.diabres.2020.108249

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Batista PJ, Chang HY. Long noncoding RNAs: cellular address codes in development and disease. Cell. (2013) 152:1298–307. doi: 10.1016/j.cell.2013.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Peng WX, Koirala P, Mo YY. LncRNA-mediated regulation of cell signaling in cancer. Oncogene. (2017) 36:5661–67. doi: 10.1038/onc.2017.184

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Bhan A, Soleimani M, Mandal SS. Long Noncoding RNA and Cancer: A New Paradigm. Cancer Res. (2017) 77:3965–81. doi: 10.1158/0008-5472.CAN-16-2634

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Riva P, Ratti A, Venturin M. The long non-coding RNAs in neurodegenerative diseases: novel mechanisms of pathogenesis. Curr Alzheimer Res. (2016) 13:1219–31. doi: 10.2174/1567205013666160622112234

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Hobuß L, Bär C, Thum T. Long non-coding RNAs: at the heart of cardiac dysfunction? Front Physiol. (2019) 10:30. doi: 10.3389/fphys.2019.00030

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Das S, Shah R, Dimmeler S, Freedman JE, Holley C, Lee JM, et al. Noncoding RNAs in cardiovascular disease: current knowledge, tools and technologies for investigation, and future directions: a scientific statement from the american heart association. Circ Genom Precis Med. (2020) 13:e000062. doi: 10.1161/HCG.0000000000000062

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, et al. FunRich: An open access standalone functional enrichment and interaction network analysis tool. Proteomics. (2015) 15:2597–601. doi: 10.1002/pmic.201400515

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. (2013) 14:128. doi: 10.1186/1471-2105-14-128

PubMed Abstract | CrossRef Full Text

23. Xie Z, Bailey A, Kuleshov MV, Clarke DJB, Evangelista JE, Jenkins SL, et al. Gene Set Knowledge Discovery with Enrichr. Curr Protoc. (2021) 1:e90. doi: 10.1002/cpz1.90

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. (2016) 44:W90–7. doi: 10.1093/nar/gkw377

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium Nat Genet. (2000) 25:25–9. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hanafi NI, Mohamed AS, Sheikh Abdul Kadir SH, Othman MHD. Overview of bile acids signaling and perspective on the signal of ursodeoxycholic acid, the most hydrophilic bile acid, in the heart. Biomolecules. (2018) 8:159. doi: 10.3390/biom8040159

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mayerhofer CCK, Ueland T, Broch K, Vincent RP, Cross GF, Dahl CP, et al. Increased secondary/primary bile acid ratio in chronic heart failure. J Card Fail. (2017) 23:666–71. doi: 10.1016/j.cardfail.2017.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Klein I, Danzi S. Thyroid disease and the heart. Curr Probl Cardiol. (2016) 41:65–92. doi: 10.1016/j.cpcardiol.2015.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Osuna PM, Udovcic M, Sharma MD. Hyperthyroidism and the heart. Methodist Debakey Cardiovasc J. (2017) 13:60–3. doi: 10.14797/mdcj-13-2-60

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Bainbridge MN, Davis EE, Choi WY, Dickson A, Martinez HR, Wang M, et al. Loss of function mutations in NNT are associated with left ventricular noncompaction. Circ Cardiovasc Genet. (2015) 8:544–52. doi: 10.1161/CIRCGENETICS.115.001026

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Simeunovic D, Odanovic N, Pljesa-Ercegovac M, Radic T, Radovanovic S, Coric V, et al. Glutathione Transferase P1 Polymorphism Might Be a Risk Determinant in Heart Failure. Dis Markers. (2019) 2019:6984845. doi: 10.1155/2019/6984845

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jiang H, Hooper C, Kelly M, Steeples V, Simon JN, Beglov J, et al. Functional analysis of a gene-edited mouse model to gain insights into the disease mechanisms of a titin missense variant. Basic Res Cardiol. (2021) 116:14. doi: 10.1007/s00395-021-00853-z

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Deng B, Zhang X, Liang Y, Jiang H, Huang W, Wu Y, et al. Nonadherent culture method promotes MSC-mediated vascularization in myocardial infarction via miR-519d/VEGFA pathway. Stem Cell Res Ther. (2020) 11:266. doi: 10.1186/s13287-020-01780-x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Statile CJ, Taylor MD, Mazur W, Cripe LH, King E, Pratt J, et al. Left ventricular noncompaction in Duchenne muscular dystrophy. J Cardiovasc Magn Reson. (2013) 15:67. doi: 10.1186/1532-429X-15-67

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Sun C, Kontaridis MI. Physiology of cardiac development: from genetics to signaling to therapeutic strategies. Curr Opin Physiol. (2018) 1:123–39. doi: 10.1016/j.cophys.2017.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zou H, Wang K, Gao Y, Song H, Xie Q, Jin M, et al. Chronic alcohol consumption from adolescence-to-adulthood in mice–hypothalamic gene expression changes in the dilated cardiomyopathy signaling pathway. BMC Neurosci. (2014) 15:61. doi: 10.1186/1471-2202-15-61

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wadmore K, Azad AJ, Gehmlich K. The Role of Z-disc Proteins in Myopathy and Cardiomyopathy. Int J Mol Sci. (2021) 22:3058. doi: 10.3390/ijms22063058

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Lau SL, Yuen ML, Kou CY, Au KW, Zhou J, Tsui SK. Interferons induce the expression of IFITM1 and IFITM3 and suppress the proliferation of rat neonatal cardiomyocytes. J Cell Biochem. (2012) 113:841–7. doi: 10.1002/jcb.23412

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zhang Z, Li LY. TNFSF15 modulates neovascularization and inflammation. Cancer Microenviron. (2012) 5:237–47. doi: 10.1007/s12307-012-0117-8

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Qiu JX, Zhou ZW, He ZX, Zhao RJ, Zhang X, Yang L, et al. Plumbagin elicits differential proteomic responses mainly involving cell cycle, apoptosis, autophagy, and epithelial-to-mesenchymal transition pathways in human prostate cancer PC-3 and DU145 cells. Drug Des Devel Ther. (2015) 9:349–417. doi: 10.2147/DDDT.S71677

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Ji H, Yi Q, Chen L, Wong L, Liu Y, Xu G, et al. Circulating miR-3197 and miR-2116-5p as novel biomarkers for diabetic retinopathy. Clin Chim Acta. (2020) 501:147–53. doi: 10.1016/j.cca.2019.10.036

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Liu C, Lai Y, Ying S, Zhan J, Shen Y. Plasma exosome-derived microRNAs expression profiling and bioinformatics analysis under cross-talk between increased low-density lipoprotein cholesterol level and ATP-sensitive potassium channels variant rs1799858. J Transl Med. (2020) 18:459. doi: 10.1186/s12967-020-02639-8

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Hussein A, Karimianpour A, Collier P, Krasuski RA. Isolated noncompaction of the left ventricle in adults. J Am Coll Cardiol. (2015) 66:578–85. doi: 10.1016/j.jacc.2015.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Turnow K, Metzner K, Cotella D, Morales MJ, Schaefer M, Christ T, et al. Interaction of DPP10a with Kv43 channel complex results in a sustained current component of human transient outward current. Ito Basic Res Cardiol. (2015) 110:5. doi: 10.1007/s00395-014-0457-2

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Ghosh S, Sati S, Sengupta S, Scaria V. Distinct patterns of epigenetic marks and transcription factor binding sites across promoters of sense-intronic long noncoding RNAs. J Genet. (2015) 94:17–25. doi: 10.1007/s12041-015-0484-2

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zhang J, Xu T, Liu L, Zhang W, Zhao C, Li S, et al. LMSM: A modular approach for identifying lncRNA related miRNA sponge modules in breast cancer. PLoS Comput Biol. (2020) 16:e1007851. doi: 10.1371/journal.pcbi.1007851

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Spencer SL, Cappell SD, Tsai FC, Overton KW, Wang CL, Meyer T. The proliferation-quiescence decision is controlled by a bifurcation in CDK2 activity at mitotic exit. Cell. (2013) 155:369–83. doi: 10.1016/j.cell.2013.08.062

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Alexandrov A, Martzen MR, Phizicky EM. Two proteins that form a complex are required for 7-methylguanosine modification of yeast tRNA. RNA. (2002) 8:1253–66. doi: 10.1017/s1355838202024019

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Zhao Y, Kong L, Pei Z, Li F, Li C, Sun X, et al. m7G methyltransferase METTL1 promotes post-ischemic angiogenesis via promoting VEGFA mRNA translation. Front Cell Dev Biol. (2021) 9:642080. doi: 10.3389/fcell.2021.642080

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: lncRNAs, LVNC, mRNAs, miRNA targets, biomarker

Citation: Tian Q, Niu H, Liu D, Ta N, Yang Q, Norton V, Wu Y, Maiti AK, Wu H and Zheng Z (2021) Expression Signatures of Long Noncoding RNAs in Left Ventricular Noncompaction. Front. Cardiovasc. Med. 8:763858. doi: 10.3389/fcvm.2021.763858

Received: 25 August 2021; Accepted: 05 October 2021;
Published: 10 November 2021.

Edited by:

Mohammad Taheri, Jena University Hospital, Germany

Reviewed by:

Tahereh Azimi, Shahid Beheshti University of Medical Sciences, Iran
Bashdar Mahmud Hussen, Hawler Medical University, Iraq

Copyright © 2021 Tian, Niu, Liu, Ta, Yang, Norton, Wu, Maiti, Wu and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhenzhong Zheng, Z3JlYXRlZnVsJiN4MDAwNDA7MTYzLmNvbQ==; Hao Wu, aGFvLnd1MyYjeDAwMDQwO2NoaWxkcmVucy5oYXJ2YXJkLmVkdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.