Abstract:Objective: Hypertrophic cardiomyopathy (HCM) is a complex heterogeneous heart disease. Recent reports found that long non-coding RNAs (lncRNAs) play an important role in the progression of cardiovascular diseases. The present study aimed to identify the novel lncRNAs and messenger RNAs (mRNAs) and determine the key pathways involved in HCM. Methods: The lncRNA and mRNA sequencing datasets of GSE68316 and GSE130036 were downloaded from the Gene Expression Omnibus (GEO) database. An integrated co-expression network analysis was conducted to identify differentially expressed lncRNAs and differentially expressed mRNAs in patients with HCM. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were explored to identify the biological functions and signaling pathways of the co-expression network. Protein–protein interaction (PPI) and hub gene networks were constructed by using Cytoscape software. Plasma samples of patients with HCM and the GSE89714 dataset were used to validate the bioinformatics results. Results: A total of 1,426 differentially expressed long non-coding RNAs (lncRNAs) and 1,715 differentially expressed mRNAs were obtained from GSE68316, of which 965 lncRNAs and 896 mRNAs were upregulated and 461 lncRNAs and 819 mRNAs were downregulated. A total of 469 differentially expressed lncRNAs and 2,407 differentially expressed mRNAs were screened from GSE130036, of which 183 lncRNAs and 1,283 mRNAs were upregulated and 286 lncRNAs and 1,124 mRNAs were downregulated. A co-expression network was constructed and contained 30 differentially expressed lncRNAs and 63 differentially expressed mRNAs, which were primarily involved in 'G-protein beta/gamma-subunit complex binding,' 'polyubiquitin modification-dependent protein binding,' 'Apelin signaling pathway,' and 'Wnt signaling pathway.' The 10 hub genes in the upregulated network [G Protein Subunit Alpha I2 (GNAI2), G Protein Subunit Alpha I1 (GNAI1), G Protein Subunit Alpha I3 (GNAI3), G Protein Subunit Gamma 2 (GNG2), G Protein Subunit Beta 1 (GNB1), G Protein Subunit Gamma 13 (GNG13), G Protein Subunit Gamma Transducin 1 (GNGT1), G Protein Subunit Gamma 12 (GNG12), AKT Serine/Threonine Kinase 1 (AKT1) and GNAS Complex Locus (GNAS)] and the 10 hub genes in the downregulated network [Nucleotide-Binding Oligomerization Domain Containing Protein 2 (NOD2), Receptor-Interacting Serine/Threonine Kinase 2 (RIPK2), Nucleotide-Binding Oligomerization Domain Containing Protein 1 (NOD1), Mitochondrial Antiviral Signaling Protein (MAVS), Autophagy Related 16-Like 1 (ATG16L1), Interferon Induced With Helicase C Domain 1 (IFIH1), Autophagy Related 5 (ATG5), TANK-Binding Kinase 1 (TBK1), Caspase Recruitment Domain Family Member 9 (CARD9), and von Willebrand factor (VWF)] were screened using cytoHubba. The expression of LA16c-312E8.2 and RP5-1160K1.3 in the plasma of patients with HCM was elevated, and the expression of the MIR22 host gene (MIR22HG) was decreased, which was consistent with our analysis, while the expression of LINC00324 and Small Nucleolar RNA Host Gene 12 (SNHG12) was not significantly different between the two groups. Verification analyses performed on GSE89714 showed the upregulated mRNAs of Chloride Voltage-Gated Channel 7 (CLCN7), N-Acetylglucosamine-1-Phosphate Transferase Subunit Gamma (GNPTG), Unk Like Zinc Finger (UNKL), Adenosine Monophosphate Deaminase 2 (AMPD2), GNAI3, WD Repeat Domain 81 (WDR81), and Serpin Family F Member 1 (SERPINF1) and downregulated mRNAs of TATA-Box Binding Protein Associated Factor 12 (TAF12) co-expressed with five crucial lncRNAs. Moreover, GNAI2, GNAI3, GNG12, and vWF were upregulated and GNAS was downregulated in the top 10 hub genes of upregulated and downregulated PPI networks. Conclusion: These findings from integrative biological analysis of lncRNA-mRNA co-expression networks explored the key genes and pathways and provide new insights into the understanding of the mechanism and discovering new therapeutic targets for HCM. Three differentially expressed pivotal lncRNAs (LA16c-312E8.2, RP5-1160K1.3, and MIR22HG) in the co-expression network may serve as biomarkers and intervention targets for the diagnosis and treatment of HCM.

Identification of a Lncrna-Mirna-mrna Network Based on Competitive Endogenous RNA Theory Reveals Functional Lncrnas in Hypertrophic Cardiomyopathy.

Identification and Functional Prediction of Long Non-Coding RNAs in Dilated Cardiomyopathy by Bioinformatics Analysis

Analysis of Changes in Circular RNA Expression and Construction of Cerna Networks in Human Dilated Cardiomyopathy.

Pgpm_a_314880 823..837

Identification of Underlying Hub Genes Associated with Hypertrophic Cardiomyopathy by Integrated Bioinformatics Analysis

SNHG 12 and Hsa-Mir-140-5p May Play an Important Role in the Cerna Network Related to Hypertrophic Cardiomyopathy

Identification of Circrna-Mirna-mrna Regulatory Network and Its Role in Cardiac Hypertrophy

Identification of key genes for hypertrophic cardiomyopathy using integrated network analysis of differential lncRNA and gene expression

Identification of Potential miRNA-mRNA Regulatory Network Contributing to Hypertrophic Cardiomyopathy (HCM)

Identification of Circulating Hub Long Noncoding RNAs Associated with Hypertrophic Cardiomyopathy Using Weighted Correlation Network Analysis

Construction and Analysis of Cardiac Hypertrophy-Associated Lncrna-Mrna Network Based on Competitive Endogenous RNA Reveal Functional Lncrnas in Cardiac Hypertrophy

Comprehensive Construction of a Circular RNA-Associated Competing Endogenous RNA Network Identified Novel Circular RNAs in Hypertrophic Cardiomyopathy by Integrated Analysis

Transcriptomic Analyses and Experimental Validation Identified Immune-Related lncRNA–mRNA Pair MIR210HG–BPIFC Regulating the Progression of Hypertrophic Cardiomyopathy

Identification of the Potential Function of Circrna in Hypertrophic Cardiomyopathy Based on Mutual RNA-RNA and RNA-RBP Relationships Shown by Microarray Data

Distinct Circulating Expression Profiles of Long Noncoding RNAs in Heart Failure Patients with Ischemic and Nonischemic Dilated Cardiomyopathy.

Construction and analysis of the lncRNA‑miRNA‑mRNA network based on competitive endogenous RNA reveals functional genes in heart failure.

The Role of Competing Endogenous RNA (cerna) Network in Ischemic Cardiomyopathy by Bioinformatic Analysis

Identification and analysis of a key long non-coding RNAs (lncRNAs)-associated module reveal functional lncRNAs in cardiac hypertrophy.

Genetic Dissection of Hypertrophic Cardiomyopathy with Myocardial RNA-Seq

Identification of Biomarkers Correlated with Hypertrophic Cardiomyopathy with Co-Expression Analysis.

Long Non-Coding and Coding RNA Profiling Using Strand-Specific RNA-seq in Human Hypertrophic Cardiomyopathy