Comprehensive Analysis of lncRNA Expression Pattern and lncRNA-miRNA-mRNA Network in a Rat Model With Cavernous Nerve Injury Erectile Dysfunction

Comprehensive analysis of biological landscape of oxidative stress-related genes in diabetic erectile dysfunction

Oxidative stress plays a pivotal role in the pathogenesis of diabetic erectile dysfunction, while specific mechanisms have not been illuminated. The study aims to reveal the genetic expression patterns of oxidative stress in diabetic erectile dysfunction. Transcriptome data of diabetic erectile dysfunction and oxidative stress-related genes (OSRGs) in the Gene Expression Omnibus database were downloaded and analyzed based on differential expression. Functional enrichment analyses were conducted to clarify the biological functions. A protein interaction framework was established, and significant gene profiles were validated in the cavernous endothelial cells, clinical patients, and rat models. A miRNA-OSRGs network was predicted and validated. The results were analyzed using Student's t-test. The analysis screened 203 differentially expressed OSRGs (p < 0.05), which had a close association with oxidoreductase activities, glutathione metabolism, and autophagy. A four-gene signature comprised of EPS8L2 (p = 0.044), GSTA3 (p = 0.015), LOX (p < 0.001) and MGST1 (p = 0.002) was well validated and regarded as the hub OSRGs. Compared with the control group, notable increases and decreases were observed in the expressions of GSTA3 (3.683 ± 0.636 vs. 0.416 ± 0.507) and LOX (2.104 ± 1.895 vs. 18.804 ± 2.751) in the validated diabetic erectile dysfunction group. The hub OSRGs-related miRNAs participated in smooth muscle cell proliferation. Besides, miR-125a-3p (p = 0.034) and miR-138-2-3p (p = 0.012) were validated as promising oxidative stress-related miRNA biomarkers. Our findings revealed the genetic alternations of oxidative stress in diabetic erectile dysfunction. These results will be instructive to explore the molecular landscape and the potential treatment for diabetic erectile dysfunction.

Extracellular Matrix and Protein Phosphorylation Dysregulation Related to Diabetes-Induced Erectile Dysfunction

Diabetes can cause erectile dysfunction (ED) in more than half of male patients. However, the mechanisms underlying diabetes-induced erectile dysfunction (DED) remain unknown. This study is aimed at systematically analyzing the cellular and molecular mechanisms leading to DED using bioinformatic analysis and providing molecular targets for predicting and treating DED. In total, we identified 800 DEGs in the DED samples compared with those in the control group. The 407 upregulated DEGs were mainly enriched in glucose and lipid metabolism-related pathways, and the 393 downregulated DEGs were primarily enriched in tissue development and structure. Dysregulated extracellular matrix genes (especially collagen and elastin) may be closely related to damage to the erectile function of the corpus cavernosum. Sixteen hub genes and 24 modules were detected with hub genes and MCODE analysis. The consensus sequence AAA (G/C) AAA was observed at the promoter sites of most genes that were enriched in the “posttranslational protein phosphorylation” pathway. These genes had abundant phosphorylation sites. Furthermore, 20 TFs targeting DEGs were identified using ChEA3 tool. In conclusion, our research comprehensively and systematically describes the molecular characteristics of DED and suggests that dysregulated extracellular matrix genes and protein phosphorylation may play critical roles in DED. Therefore, they may be potential markers for diagnosing and treating DED.

Pathway analysis of microarray data from corpora cavernosal tissue of patients with a prostatectomy or Peyronie disease in comparison with a cavernous nerve-injured rat model of erectile dysfunction

INTRODUCTION

Patients with a prostatectomy are at high risk of developing erectile dysfunction (ED) that is refractory to phosphodiesterase type 5 inhibitors. The cavernous nerve (CN) is frequently damaged during prostatectomy, causing loss of innervation to the penis. This initiates corpora cavernosal remodeling (apoptosis and fibrosis) and results in ED.

AIM

To aid in the development of novel ED therapies, the current aim was to obtain a global understanding of how signaling mechanisms alter in the corpora cavernosa with loss of CN innervation that results in ED.

METHODS

Microarray and pathway analysis were performed on the corpora cavernosal tissue of patients with a prostatectomy (n = 3) or Peyronie disease (control, n = 3). Results were compared with an analysis of a Sprague-Dawley rat CN injury model (n = 10). RNA was extracted by TRIzol, DNase treated, and purified by a Qiagen Mini Kit. Microarray was performed with the Human Gene 2.0 ST Array and the RU34 rat array. Differentially expressed genes were identified through several analytic tools (ShinyGO, Ingenuity, WebGestalt) and databases (GO, Reactome). A 2-fold change was used as the threshold for differential expression.

OUTCOMES

Pathways that were altered (up- or downregulated) in response to CN injury in the prostatectomy patients and a rat CN injury model were determined.

RESULTS

Microarray identified 197 differentially expressed protein-coding genes in the corpora cavernosa from the prostatectomy cohort, with 100 genes upregulated and 97 genes downregulated. Altered signaling pathways that were identified that affect tissue morphology included the following: neurologic disease, cell death and survival, tissue and cellular development, skeletal and muscle development and disorders, connective tissue development and function, tissue morphology, embryonic development, growth and proliferation, cell-to-cell signaling, and cell function and maintenance. These human pathways have high similarity to those observed in the CN-injured rat ED model.

CLINICAL IMPLICATIONS

Significant penile remodeling continues in patients long after the acute surgical injury to the CN takes place, offering the opportunity for clinical intervention to reverse penile remodeling and improve erectile function.

STRENGTHS AND LIMITATIONS

Understanding how signaling pathways change in response to CN injury and how these changes translate to altered morphology of the corpora cavernosa and ensuing ED is critical to identify strategic targets for therapy development.

CONCLUSION

Altered signaling in pathways that regulate tissue homeostasis, morphogenesis, and development was identified in penes of patients with a prostatectomy, and competitive forces of apoptosis and proliferation/regeneration were found to compete to establish dominance after CN injury. How these pathways interact to regulate penis tissue homeostasis is a complex process that requires further investigation.

A comprehensive analysis on the relationship between BDE-209 exposure and erectile dysfunction

Decabromodiphenyl ether (mainly BDE-209) is a commonly used brominated flame retardant in various industrial products. Although its damage to the reproduction system has been established, its effect on erectile function remains unclear. The present study investigated whether BDE-209 induced erectile dysfunction in male SD rats and the underlying mechanisms. Pubertal male rats were exposed to BDE-209 orally (0, 5, 50, and 500 mg/kg/day) for 28 days and the ICP (intracavernous pressure) and MAP (mean arterial pressure) were measured. After the rats were euthanized, the fibrosis and apoptosis levels were evaluated. Additionally, the endothelial function of the rat vascular endothelium cells and the human umbilical vein endothelial cells were impaired after treatment with 50 μM and 100 μM BDE-209. Moreover, the bioinformatics based on CTD database and ChIP-X Enrichment Analysis, version 3 (ChEA3) and molecular docking analysis demonstrated that 5 transcription factors (NFKB1, NR3C1, E2F5, REL, IRF4) might regulate endothelial function by affecting the expression of interactive genes (BCL-2, CAP3, CAT, TNF, MAPK1, and MAPK3). In summary, the present study demonstrated that BDE-209 might affect downstream interactive genes by binding to transcription factors, leading to corpus cavernosum endothelial dysfunction, thus contributing to erectile dysfunction in rats.

Comparative Transcriptome Analyses of Geriatric Rats Associate Age-Related Erectile Dysfunction With a lncRNA-miRNA-mRNA Regulatory Network

BACKGROUND

The key regulatory roles of long non-coding RNAs (lncRNAs) in age-related erectile dysfunction (A-ED) are unknown.

AIM

This study aimed to identify putative lncRNAs that regulate age-related erectile dysfunction via transcriptome analyses, and to predict their specific regulatory routes via bioinformatics methods.

METHODS

22 geriatric male SD rats were divided into age-related erectile dysfunction (A-ED) and negative control (NC) groups after evaluations of intracavernous pressure (ICP). By comparative analysis of transcriptomes of cavernosal tissues from both groups, we identified differentially expressed lncRNAs, miRNAs, and mRNAs. Seven differentially expressed lncRNAs were selected and further verified by quantitative real-time polymerase chain reactions (RT-qPCR). The construction of the lncRNA-miRNA-mRNA network, the Gene Ontology (GO) term enrichment, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed in Cytoscape.

RESULTS

From comparative transcriptome analyses of A-ED and NC groups, 69, 29, and 364 differentially expressed lncRNAs, miRNAs, and mRNAs were identified respectively. Differentially expressed lncRNAs were culled to seven, which were all verified by qPCR. Three of these lncRNAs (ENSRNOT00000090050, ENSRNOT00000076482, and ENSRNOT00000029245) were used to build regulatory networks, of which only ENSRNOT00000029245 was successful. Moreover, GO and KEGG analyses demonstrated that these lncRNAs possibly regulated muscle myosin complex, muscle cell cellular homeostasis, and ultimately erectile function in rats through PI3K-Akt, fluid shear stress, and atherosclerosis pathways.

CONCLUSION

Our study identified differentially expressed lncRNAs, miRNAs, and mRNAs through comparisons of transcriptomes of geriatric rats. An identified lncRNA verified by qPCR, was used to construct a lncRNA-miRNA-mRNA regulatory network. LncRNA ENSRNOT00000029245 possibly regulated downstream mRNAs through this regulatory network, leading to apoptosis in the cavernous tissue, fibrosis, and endothelial dysfunction, which ultimately caused ED. These findings provide seminal insights into the molecular biology of aging-related ED, which could spur the development of effective therapeutics.

Long non-coding RNA SNHG7 upregulates FGF9 to alleviate oxygen and glucose deprivation-induced neuron cell injury in a miR-134-5p-dependent manner

Long non-coding RNA small nucleolar RNA host gene 7 (SNHG7) was reported to regulate the pathogenesis of ischemic stroke. The study aimed to disclose SNHG7 role in oxygen and glucose deprivation (OGD)-induced Neuro-2a (N2a) cell disorders. An OGD injury cell model was established using N2a cells. The expression of SNHG7, microRNA-134-5p (miR-134-5p) and fibroblast growth factor 9 (FGF9) was determined by quantitative real-time polymerase chain reaction. Protein expression was detected by western blot. Cell viability and Lactate Dehydrogenase (LDH) leakage were determined by cell counting kit-8 and LDH activity detection assays. Oxidative stress was investigated by Superoxide Dismutase and Catalase activity assays as well as Malondialdehyde and Reactive Oxygen Species detection kits. Cell apoptosis and caspase-3 activity were severally demonstrated by flow cytometry and caspase-3 activity assays. The interaction between miR-134-5p and SNHG7 or FGF9 was predicted by online databases, and identified by mechanism assays. OGD treatment decreased SNHG7 and FGF9 expression, but increased miR-134-5p expression. OGD treatment repressed cell viability, promoted LDH leakage and induced oxidative stress and apoptosis in N2a cells, which was rescued by SNHG7 overexpression. SNHG7 acted as a sponge for miR-134-5p, and regulated OGD-triggered cell damage by associating with miR-134-5p. Additionally, miR-134-5p depletion protected N2a cells from OGD-induced injury by targeting FGF9. Ectopic SNHG7 expression protected against OGD-induced neuronal cell injury by inducing FGF9 through sponging miR-134-5p, providing a novel therapeutic target for ischemic stroke.

Transcriptome profiling of lncRNA and co-expression network in the vaginal epithelial tissue of women with lubrication disorders

Background

Vaginal lubrication is a crucial physiological response that occurs at the beginning of sexual arousal. However, research on lubrication disorders (LD) is still in its infancy, and the role of long non-coding RNAs (lncRNAs) in LD remains unclear. This study aimed to explore the function of lncRNAs in the pathogenesis of vaginal LD.

Methods

The expression profiles of LD and normal control (NC) lncRNAs were examined using next-generation sequencing (NGS), and eight selected differentially expressed lncRNAs were verified by quantitative real-time PCR. We conducted GO annotation and KEGG pathway enrichment analyses to determine the principal functions of significantly deregulated genes. LncRNA-mRNA co-expression and protein-protein interaction (PPI) networks were constructed and the lncRNA transcription factors (TFs) were predicted.

Results

From the results, we identified 181,631 lncRNAs and 145,224 mRNAs in vaginal epithelial tissue. Subsequently, our preliminary judgment revealed a total of 499 up-regulated and 337 down-regulated lncRNAs in LD. The top three enriched GO items of the dysregulated lncRNAs included the following significant terms: “contractile fiber part,” “actin filament-based process,” and “contractile fiber”. The most enriched pathways were “cell-extracellular matrix interactions,” “muscle contraction,” “cell-cell communication,” and “cGMP-PKG signaling pathway”. Our results also showed that the lncRNA-mRNA co-expression network was a powerful platform for predicting lncRNA functions. We determined the three hub genes, ADCY5, CXCL12, and NMU, using PPI network construction and analysis. A total of 231 TFs were predicted with RHOXF1, SNAI2, ZNF354C and TBX15 were suspected to be involved in the mechanism of LD.

Conclusion

In this study, we constructed the lncRNA-mRNA co-expression network, predicted the lncRNA TFs, and comprehensively analyzed lncRNA expression profiles in LD, providing a basis for future studies on LD clinical biomarkers and therapeutic targets. Further research is also needed to fully determine lncRNA’s role in LD development.

Identification of miRNA Signature Associated With Erectile Dysfunction in Type 2 Diabetes Mellitus by Support Vector Machine-Recursive Feature Elimination

Diabetic mellitus erectile dysfunction (DMED) is one of the most common complications of diabetes mellitus (DM), which seriously affects the self-esteem and quality of life of diabetics. MicroRNAs (miRNAs) are endogenous non-coding RNAs whose expression levels can affect multiple cellular processes. Many pieces of studies have demonstrated that miRNA plays a role in the occurrence and development of DMED. However, the exact mechanism of this process is unclear. Hence, we apply miRNA sequencing from blood samples of 10 DMED patients and 10 DM controls to study the mechanisms of miRNA interactions in DMED patients. Firstly, we found four characteristic miRNAs as signature by the SVM-RFE method (hsa-let-7E-5p, hsa-miR-30 days-5p, hsa-miR-199b-5p, and hsa-miR-342–3p), called DMEDSig-4. Subsequently, we correlated DMEDSig-4 with clinical factors and further verified the ability of these miRNAs to classify samples. Finally, we functionally verified the relationship between DMEDSig-4 and DMED by pathway enrichment analysis of miRNA and its target genes. In brief, our study found four key miRNAs, which may be the key influencing factors of DMED. Meanwhile, the DMEDSig-4 could help in the development of new therapies for DMED.

Identification of Regeneration and Hub Genes and Pathways at Different Time Points after Spinal Cord Injury

Spinal cord injury (SCI) is a neurological injury that can cause neuronal loss around the lesion site and leads to locomotive and sensory deficits. However, the underlying molecular mechanisms remain unclear. This study aimed to verify differential gene time-course expression in SCI and provide new insights for gene-level studies. We downloaded two rat expression profiles (GSE464 and GSE45006) from the Gene Expression Omnibus database, including 1 day, 3 days, 7 days, and 14 days post-SCI, along with thoracic spinal cord data for analysis. At each time point, gene integration was performed using "batch normalization." The raw data were standardized, and differentially expressed genes at the different time points versus the control were analyzed by Gene Ontology enrichment analysis, the Kyoto Encyclopedia of Genes and Genomes pathway analysis, and gene set enrichment analysis. A protein-protein interaction network was then built and visualized. In addition, ten hub genes were identified at each time point. Among them, Gnb5, Gng8, Agt, Gnai1, and Psap lack correlation studies in SCI and deserve further investigation. Finally, we screened and analyzed genes for tissue repair, reconstruction, and regeneration and found that Anxa1, Snap25, and Spp1 were closely related to repair and regeneration after SCI. In conclusion, hub genes, signaling pathways, and regeneration genes involved in secondary SCI were identified in our study. These results may be useful for understanding SCI-related biological processes and the development of targeted intervention strategies.