Identification of a potentially novel LncRNA-miRNA-mRNA competing endogenous RNA network in pulmonary arterial hypertension via integrated bioinformatic analysis

Pathology and pathobiology of pulmonary hypertension: current insights and future directions

In recent years, major advances have been made in the understanding of the cellular and molecular mechanisms driving pulmonary vascular remodelling in various forms of pulmonary hypertension, including pulmonary arterial hypertension, pulmonary hypertension associated with left heart disease, pulmonary hypertension associated with chronic lung disease and hypoxia, and chronic thromboembolic pulmonary hypertension. However, the survival rates for these different forms of pulmonary hypertension remain unsatisfactory, underscoring the crucial need to more effectively translate innovative scientific knowledge into healthcare interventions. In these proceedings of the 7th World Symposium on Pulmonary Hypertension, we delve into recent developments in the field of pathology and pathophysiology, prioritising them while questioning their relevance to different subsets of pulmonary hypertension. In addition, we explore how the latest omics and other technological advances can help us better and more rapidly understand the myriad basic mechanisms contributing to the initiation and progression of pulmonary vascular remodelling. Finally, we discuss strategies aimed at improving patient care, optimising drug development, and providing essential support to advance research in this field.

Hsa_circ_0014606 Derived from Exosomes Promotes Gastric Carcinoma Tumorigenesis and Proliferation by Sponging miR-514b-3p to Upregulate HNRNPC

Gastric cancer is a common malignant tumor, and due to its insidious onset and limited screening methods, most patients are diagnosed with advanced disease and have a poor prognosis. The circRNA in exosomes has an essential role in cancer diagnosis and treatment. However, the part of hsa_circ_0014606 within exosomes in gastric cancer progression is unclear. Firstly, we extracted exosomes from the serum of gastric cancer patients and healthy individuals by ultracentrifugation and analyzed the expression of hsa_circ_0014606 in both exosomes; then knocked down hsa_circ_0014606 in vivo and in vitro, respectively, to observe its effect on the physiological function of gastric cancer cells; finally, we used bioinformatics to screen hsa_circ_0014606 targeting miRNAs and mRNAs, and experiments were performed to verify the interrelationship between the three. The results showed that the level of hsa_circ_0014606 in the serum exosomes of gastric cancer patients was significantly higher than that of the healthy population. The knockdown of hsa_circ_0014606 slowed the proliferation of gastric cancer cells, significantly reduced migration and invasion ability, accelerated apoptosis, and reduced tumor size in mice. In addition, the expression of hsa_circ_0014606 was negatively correlated with the expression of miR-514b-3p and positively correlated with the expression of heterogeneous nuclear ribonucleoprotein C (HNRNPC). In conclusion, hsa_circ_0014606 exerted a pro-cancer effect indirectly through miR-514b-3p targeting gene HNRNPC, and this study provides a new potential target for treating gastric cancer.

Endothelial downregulation of nuclear m6A reader YTHDC1 promotes pulmonary vascular remodeling in sugen hypoxia model of pulmonary hypertension

Background

Pulmonary hypertension (PH) is characterized with vascular remodeling, which is intiated by vascular endothelial dysfunction. N6-methyladenosine (m6A) modification mediates gene expression in many ways including mediating RNA degradation, splicing, nuclear export et al. m6A modification have been found to be associated with the development of PH. However, the role of m6A regulators in pulmonary artery endothelial cells (PAECs) dysfunction of PH is still under research.

Methods

The expression levels of m6A regulators in PAECs were analyzed with the single-cell sequencing Data(scRNA). Next, the target differentially expressed genes (DEGs) of m6A regulators in PAECs were functionally annotated. The analysis of cellular interactions included the examination of receptor-ligand pairs regulated by m6A regulators. Pseudo-time trajectory analyses and a ceRNA network involving lncRNAs, miRNAs, and mRNAs were conducted in PAECs. Furthermore, microarray data (GSE180169) for Sugen Hypoxia PH (SuHx PH) mouse models was screened for DEGs and m6A regulators in PAECs. Moreover, the expression of YTHDC1 in the lung samples of SuHx PH models was determined using immunofluorescence. In vitro, the mRNA expression of YTHDC1 in HPAECs under hypoxia conditions was detected. The effect of YTHDC1 recombinant protein on HPAEC proliferation was detected by Cell Counting Kit-8 (CCK8).

Results

Dysregulation of m6A regulators was observed in mouse PAECs. The m6A reader of YTHDC1 was decreased in PAECs in scRNA data and RNAseq data of isolated PAECs of SuHx PH models. Downregulation of YTHDC1 was caused by hypoxia in PAECs in vitro and similar results was observed in PAECs of SuHx PH mouse models. Next, YTHDC1 recombinant protein was found to inhibit HPAECs proliferation. The DEGs targeted by YTHDC1 were enriched in angiogenesis, endothelial cell migration, fluid shear stress, and stem cell maintenance. Analysis indicates that interactions among endothelial cells, smooth muscle cells, fibroblasts, and immune cells, mediated by specific YTHDC1 target genes (e.g., PTPRC-MRC1, ITBG2-ICAM1, COL4A1-CD44), contribute to PH development. Also, the YTHDC1 expression were consistent with Thioredoxin interacting protein (TXNIP). What's more, the predicted transcription factors showed that NFKB1, Foxd3 may be involved in the regulation of YTHDC1. Lastly, our data suggest that YTHDC1 may be involved in regulating PAECs dysfunction through lncRNA/miRNA/mRNA network.

Conclusion

For the first time, we analyzed changes in the expression and biological functions of m6A regulators in SuHx PH mouse models. We causatively linked YTHDC1 to PAECs dysfunction, providing novel insight into and opportunities to diagnose and treat PH.

O-GlcNAcylation of RAB10 promotes hepatocellular carcinoma progression

Ras-related protein Rab-10 (RAB10) is involved in tumorigenesis and progression of hepatocellular carcinoma (HCC). Here, we found RAB10, O-GlcNAc transferase (OGT), and O-GlcNAcylation were upregulated in HCC. In addition, RAB10 protein level was prominently positively correlated with the expression of OGT. O-GlcNAcylation modification of RAB10 was then investigated. Here we showed that RAB10 interacts directly with OGT in HCC cell lines, Meanwhile, O-GlcNAcylation enhanced RAB10 protein stability. Furthermore, knockdown of OGT suppressed aggressive behaviors of HCC in vitro and in vivo, while elevated RAB10 reversed these. Taken together, these results indicated that OGT mediated O-GlcNAcylation stabilized RAB10, thus accelerating HCC progression.

Integrated analysis of the lncRNA-associated competing endogenous RNA network in salt sensitivity of blood pressure

Accumulating evidence showed that competing endogenous RNA (ceRNA) mechanism plays a pivotal role in salt sensitivity of blood pressure (SSBP). We constructed a ceRNA network based on SSBP-related differently expressed lncRNAs (2), mRNAs (73) and miRNAs (18). Bioinformatic analyses were utilized to analyze network and found network genes participate in biological pathways related to SSBP pathogenesis such as regulation of nitric oxide biosynthetic process (GO:0045,428) and cellular response to cytokine stimulus (GO:0071,345). Fourteen candidate ceRNA pathways were selected from network to perform qRT-PCR validation and found nine RNAs (KCNQ1OT1, SLC8A1-AS1, IL1B, BCL2L11, KCNJ15, CX3CR1, KLF2, hsa-miR-362-5p and hsa-miR-423-5p) differently expressed between salt-sensitive (SS) and salt-resistant (SR) groups (P < 0.05). Four ceRNA pathways were further validated by luciferase reporter assay and found KCNQ1OT1→hsa-miR-362-5p/hsa-miR-423-5p→IL1B pathways may influence the pathogenic mechanism of SS. Our findings suggested the ceRNA pathway and network may affect SS occurrence mainly through endothelial dysfunction and inflammatory activation.

Whole-genome microRNA sequencing analysis in patients with pulmonary hypertension

Pulmonary hypertension (PH) is a pathological disorder with multiple clinical manifestations that lead to cardiovascular and respiratory diseases in most patients. Recent studies have revealed that microRNAs (miRNAs) play important roles as upstream signaling molecules in several diseases, including PH. However, miRNAs that can be used as diagnostic or prognostic biomarkers for PH have not been identified. Thus, in this study, peripheral blood samples obtained from patients with PH and healthy individuals were subjected to genome-wide miRNA sequencing and transcriptome analysis. We screened 136 differentially expressed miRNAs in patients with PH and verified that four differentially expressed miRNAs, namely, hsa-miR-1304-3p, hsa-miR-490-3p, hsa-miR-11400, and hsa-miR-31-5p, could be used as clinical diagnostic biomarkers for pulmonary arterial hypertension. Our findings provide a basis for further in-depth investigations of the specific mechanisms of miRNAs in PH.

Competing endogenous RNA network analysis of the molecular mechanisms of ischemic stroke

BACKGROUND

Ischemic stroke (IS) is a serious neurological disease that largely results in long-term disability and death. Extensive evidence has indicated that the activation of inflammation and ferroptosis significantly contribute to the development of IS pathology. However, the underlying molecular mechanism remains unclear. In this study, we aimed to identify potential biomarkers associated with IS through the construction of a competing endogenous RNA (ceRNA) network and to investigate the possible inflammatory and ferroptosis-related molecular mechanisms.

RESULTS

We identified 178 differentially expressed target messenger RNAs (DETmRNAs) associated with IS. As revealed through enrichment analysis, the DEmRNAs were mainly enriched in the inflammatory signaling pathways and also related to ferroptosis mechanism. The CIBERSORT algorithm showed immune infiltration landscapes in which the naïve B cells, naïve T cells, and monocytes had statistically different numbers in the cerebral infarction group compared with the control group. A ceRNA network was constructed in this study involving 44 long non-coding RNAs (lncRNAs), 15 microRNAs (miRNAs), and 160 messenger RNAs (mRNAs). We used the receiver operating characteristic (ROC) analysis to identify three miRNAs (miR-103a-3p, miR-140-3p, and miR-17-5p), one mRNA (TLR4), and one lncRNA (NEAT1) as the potential key biomarkers of the ceRNA network. The key mRNA and lncRNA were shown to be highly related to the ferroptosis mechanism of IS. The expression of these key biomarkers was also further validated by a method of quantitative real-time polymerase chain reaction in SH-SY5Y cells, and the validated results were consistent with the findings predicted by bioinformatics.

CONCLUSION

Our results suggest that the ceRNA network may exert an important role in the inflammatory and ferroptosis molecular mechanisms of IS, providing new insight into therapeutic IS targets.

Hub Genes and Immune Cell Infiltration in Hypoxia-Induced Pulmonary Hypertension: Bioinformatics Analysis and In Vivo Validation

BACKGROUND

Hypoxia-induced pulmonary hypertension (HPH) represents a severe pulmonary disorder with high morbidity and mortality, which necessitates identifying the critical molecular mechanisms underlying HPH pathogenesis.

METHODS

The mRNA expression microarray GSE15197 (containing 8 pulmonary tissues from HPH and 13 normal controls) was downloaded from Gene Expression Omnibus (GEO). Gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) were executed by RStudio software. The Protein-Protein Interaction (PPI) network was visualized and established using Cytoscape, and the cytoHubba app from Cytoscape was used to pick out the hub modules. The infiltration of immune cells in HPH was analyzed using the CIBERSORTx. To confirm the potential hub genes, real-time quantitative reverse transcription PCR (qRT-PCR) was conducted using lung tissues of rat HPH models and controls.

RESULTS

A total of 852 upregulated and 547 downregulated genes were identified. The top terms in biological processes were apoptosis, proliferation, and regulation of the MAPK cascade, including ERK1/2. Cytoplasm, cytosol, and membrane were enriched in cellular component groups. Molecular functions mainly focus on protein binding, protein serine/threonine kinase activity and identical protein binding. KEGG analysis identified pathways in cancer, regulation of actin cytoskeleton and rap1 signaling pathway. There was significantly different immune cell infiltration between HPH and normal control samples. High proportions of the memory subsets of B cells and CD4 cells, Macrophages M2 subtype, and resting Dendritic cells were found in HPH samples, while high proportions of naive CD4 cells and resting mast cells were found in normal control samples. The qRT-PCR results showed that among the ten identified hub modules, FBXL3, FBXL13 and XCL1 mRNA levels were upregulated, while NEDD4L, NPFFR2 and EDN3 were downregulated in HPH rats compared with control rats.

CONCLUSION

Our study revealed the key genes and the involvement of immune cell infiltration in HPH, thus providing new insight into the pathogenesis of HPH and potential treatment targets for patients with HPH.

The identification and verification of hub genes associated with pulmonary arterial hypertension using weighted gene co-expression network analysis

BACKGROUND

Pulmonary arterial hypertension (PAH) is characterized by a progressive increase in pulmonary vascular resistance and pulmonary arterial pressure, with complex etiology, difficult treatment and poor prognosis. The objective of this study was to investigate the potential biomarkers for PAH based on bioinformatics analysis.

METHODS

The GSE117261 datasets were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified by screening PAH patients and controls. Then the DEGs were analyzed using a Weighted Gene Co-expression Network Analysis (WGCNA) and the key modules were determined, and to further explore their potential biological functions via Gene Ontology analysis (GO), Kyoto Encyclopedia of Genes and Genomes Pathway analysis (KEGG), and Gene Set Enrichment Analysis (GSEA). Moreover, Protein-protein interaction (PPI) networks were constructed to identify hub gene candidates in the key modules. Finally, real-time quantitative polymerase chain reaction was supplied to detect the expressions of hub genes in human pulmonary arterial smooth cells treated with cobalt chloride (COCl₂) which was used to mimic hypoxia.

RESULTS

There were 2299 DEGs identified. WGCNA indicated that yellow module was the key one correlated with PAH. GO and KEGG analysis demonstrated that genes in the yellow module were mainly enriched in 'Pathways in cancer'. GSEA revealed that 'HALLMARK_MYC_TARGETS_V1' was remarkably enriched in PAH. Based on the PPI network, vascular endothelial growth factor A, proto-oncogene receptor tyrosine kinase (KIT), PNN interacting serine and arginine rich protein (PNISR) and heterogeneous nuclear ribonucleoprotein H1 (HNRNPH1) were identified as the hub genes. Additionally, the PCR indicated that the elevated expressions of PNISR and HNRNPH1 were in line with the bioinformatics analysis. ROC analysis determined that PNISR and HNRNPH1 may be potential biomarkers to provide better diagnosis of PAH.

CONCLUSION

PNISR and HNRNPH1 were potential biomarkers to diagnosis PAH. In summary, the identified DEGs, modules, pathways, and hub genes provide clues and shed light on the potential molecular mechanisms of PAH.

Comprehensive analysis of the expression of N6-methyladenosine RNA methylation regulators in pulmonary artery hypertension

Background: Pulmonary arterial hypertension (PAH) is a progressive disease characterized by pulmonary vascular remodeling. The development of PAH involves N6-methyladenosine (m6A) modification. However, the functional role of m6A regulators in PAH and the underlying regulatory mechanisms remain unknown so far.

Methods: Microarray data (GSE149713) for monocrotaline induced PAH (MCT-PAH) rat models were downloaded and screened for differentially expressed genes (DEGs) and m6A regulators. Next, we screened for differentially expressed m6A regulators in endothelial cells (ECs), smooth muscle cells (SMCs), fibroblasts, interstitial macrophages, NK cells, B cells, T cells, regulatory T cells (Tregs) using scRNA sequencing data. The target DEGs of m6A regulators in ECs, SMCs, fibroblasts, and Tregs were functionally annotated using the Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. In addition, the cellular interaction analysis was performed to reveal the receptor—ligand pairs regulated by m6A regulators. Pseudo-time trajectory analyses were performed and a ceRNA network of lncRNAs-miRNAs-mRNAs was constructed in SMCs. Furthermore, the RNA transcriptome sequencing data for the SMCs isolated from idiopathic PAH (IPAH) patients (GSE144274) were validated for differentially expressed m6A regulators. Moreover, the HNRNPA2B1 levels in the lung samples from PAH patients and MCT-PAH were determined using immunohistochemistry.

Results: The m6A regulators were observed to be dysregulated in PAH. HNRNPA2B1expression level was increased in the PASMCs of scRNAs and IPAH patients. The target DEGs of HNRNPA2B1 were enriched in the regulation of muscle cell differentiation and vasculature development in PASMCs. The HNRNPA2B1 expression levels determined were consistent with the proliferation-related and collagen synthesis-related gene COL4A1. Moreover, the predicted transcription factors (TFs) foxd2/3 and NFκB could be involved in the regulation of HNRNPA2B1. HNRNPA2B1 might be regulating SMCs proliferation and phenotypic transition via rno-miR-330–3p/TGFβR3 and rno-miR-125a-3p/slc39a1. In addition, HNRNPA2B1 was observed to be highly expressed in the lung samples from MCT-PAH rat models and patients with PAH.

Conclusion: In summary, the present study identified certain key functional m6A regulators that are involved in pulmonary vascular remodeling. The investigation of m6A patterns might be promising and provide biomarkers for diagnosis and treatment of PAH in the future.

THUMPD3-AS1 facilitates cell growth and aggressiveness by the miR-218-5p/SKAP1 axis in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a malignant cancer with a high mortality. Accumulating studies have revealed that mRNAs involved in ceRNA (competing endogenous RNA) network are implicated in the tumorigenesis and development of CRC. Here, we aimed to elucidate the ceRNA network involving Src kinase associated phosphoprotein 1 (SKAP1) in the biological characteristics of CRC.

METHODS

Expression levels of genes in colon adenocarcinoma (COAD) samples and prognosis of COAD patients were predicted using publicly available online tool. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), clony formation and Transwell assays were conducted to test the biological functions of SKAP1 and THUMPD3 antisense RNA 1 (THUMPD3-AS1) in CRC cells. Western blot was used to measure the protein levels of SKAP1. Gene expression in CRC cells was detected by reverse transcription quantitative polymerase chain reaction (RT-qPCR). The interaction between miR-218-5p and THUMPD3-AS1 (or SKAP1) was verified by RNA pulldown and luciferase reporter assays.

RESULTS

SKAP1 was upregulated in COAD tissues and CRC cells and it reflected a poor prognosis in patients with COAD. SKAP1 knockdown inhibited CRC (HT-29 and HCT-116) cell proliferation, migration and invasion. Mechanistically, THUMPD3-AS1 acted as a ceRNA to sponge miR-218-5p and subsequently upregulated SKAP1 expression in CRC cells. SKAP1 overexpression reversed the suppressive effect of THUMPD3-AS1 knockdown on proliferation, migration and invision of CRC cells.

CONCLUSIONS

THUMPD3-AS1 promotes CRC cell growth and aggressiveness by regulating the miR-218-5p/SKAP1 axis.