Transcriptome Sequencing Data Reveal LncRNA-miRNA-mRNA Regulatory Network in Calcified Aortic Valve Disease

Calcified aortic valve disease (CAVD) is one of the most common valvular heart diseases in the elderly population. However, no effective medical treatments have been found to interfere with the progression of CAVD, and specific molecular mechanisms of CAVD remain unclear.

Transcriptome sequencing data of GSE55492 and GSE148219 were downloaded from the European Nucleotide Archive, and the microarray dataset, GSE12644 was acquired from the Gene Expression Omnibus database. Software, including FastQC, HISAT2, samtools, and featureCounts was applied to generate the read count matrix. The "Limma" package in R was utilized to analyze differentially expressed genes (DEGs). Thereafter, weighted gene co-expression network analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and the protein-protein interaction (PPI) network were used to identify hub genes associated with CAVD, which were further validated by receiver operating characteristic curve (ROC) analysis using GSE12644. The long non-coding RNA (LncRNA)-mediated regulatory network was established based on the differentially expressed LncRNAs and hub genes, which were detected using quantitative real-time PCR (qRT-PCR) in clinical samples and valve interstitial cells. Moreover, CIBERSORT was used to calculate the expression distribution of immune cell infiltration in CAVD.

A total of 126 DEGs were included in the PPI network. PI3K-Akt signaling pathway, ECM-receptor interaction, hematopoietic cell lineage, cell adhesion molecules, and focal adhesion were the most enriched pathways revealed by KEGG. Four LncRNAs, including TRHDE-AS1, LINC00092, LINC01094, and LINC00702 were considered the differentially expressed LncRNA. SPP1, TREM1, GPM6A, CCL19, CR1, NCAM1, CNTN1, TLR8, SDC1, and COL6A6 were the 10 hub genes identified to be associated with CAVD. Moreover, the calcified aortic valve samples had a greater level of Tregs, naïve B cells, and M0 macrophages than the noncalcified ones, whereas CAVD samples had a lower M2 macrophage expression compared to the noncalcified valve tissues.

The current study identified SPP1, TREM1, TLR8, SDC1, GPM6A, and CNTN1 as hub genes that could potentially be associated with CAVD. The LINC00702-miR-181b-5p-SPP1 axis might participate in the development of CAVD. Additionally, M2 macrophages, Tregs, naïve B cells, and M0 macrophages might possibly play a role in the initiation of CAVD.