Analysis of characteristic genes and ceRNA regulation mechanism of endometriosis based on full transcriptional sequencing.

Xie C ，Yin Z ，Liu Y 《Frontiers in Genetics》

Gene expression analysis in endometriosis: Immunopathology insights, transcription factors and therapeutic targets.

Endometriosis is recognized as an estrogen-dependent inflammation disorder, estimated to affect 8%-15% of women of childbearing age. Currently, the etiology and pathogenesis of endometriosis are not completely clear. Underlying mechanism for endometriosis is still under debate and needs further exploration. The involvement of transcription factors and immune mediations may be involved in the pathophysiological process of endometriosis, but the specific mechanism remains to be explored. This study aims to investigate the underlying molecular mechanisms in endometriosis. The gene expression profile of endometriosis was obtained from the gene expression omnibus (GEO) database. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied to the endometriosis GSE7305 datasets. Cibersort and MCP-counter were used to explore the immune response gene sets, immune response pathway, and immune environment. Differentially expressed genes (DEGs) were identified and screened. Common biological pathways were being investigated using the kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis. Transcription factors were from The Human Transcription Factors. The least absolute shrinkage and selection operator (Lasso) model identified four differential expressions of transcription factors (AEBP1, HOXB6, KLF2, and RORB). Their diagnostic value was calculated by receiver operating characteristic (ROC) curve analysis and validated in the validation cohort (GSE11691, GSE23339). By constructing the interaction network of crucial transcription factors, weighted gene coexpression network analysis (WGCNA) was used to search for key module genes. Metascape was used for enrichment analysis of essential module genes and obtained HOXB6, KLF2. The HOXB6 and KLF2 were further verified as the only two intersection genes according to Support Vector Machine Recursive Feature Elimination (SVM-RFE) and random forest models. We constructed ceRNA (lncRNA-miRNA-mRNA) networks with four potential transcription factors. Finally, we performed molecular docking for goserelin and dienogest with four transcription factors (AEBP1, HOXB6, KLF2, and RORB) to screen potential drug targets. Immune and metabolic pathways were enriched in GSVA and GSEA. In single sample gene set enrichment analysis (ssGSEA), most immune infiltrating cells, immune response gene sets, and immune response pathways are differentially expressed between endometriosis and non-endometriosis. Twenty-seven transcription factors were screened from differentially expressed genes. Most of the twenty-seven transcription factors were correlated with immune infiltrating cells, immune response gene sets and immune response pathways. Furthermore, Adipocyte enhancer binding protein 1 (AEBP1), Homeobox B6 (HOXB6), Kruppel Like Factor 2 (KLF2) and RAR Related Orphan Receptor B (RORB) were selected out from twenty-seven transcription factors. ROC analysis showed that the four genes had a high diagnostic value for endometriosis. In addition, KLF2 and HOXB6 were found to play particularly important roles in multiple modules (String, WGCNA, SVM-RFE, random forest) on the gene interaction network. Using the ceRNA network, we found that NEAT1 may regulate the expressions of AEBP1, HOXB6 and RORB, while X Inactive Specific Transcript (XIST) may control the expressions of HOXB6, RORB and KLF2. Finally, we found that goserelin and dienogest may be potential drugs to regulate AEBP1, HOXB6, KLF2 and RORB through molecular docking. AEBP1, HOXB6, KLF2, and RORB may be potential biomarkers for endometriosis. Two of them, KLF2 and HOXB6, are critical molecules in the gene interaction network of endometriosis. Discovered by molecular docking, AEBP1, HOXB6, KLF2, and RORB are targets for goserelin and dienogest.

Geng R ，Huang X ，Li L ，Guo X ，Wang Q ，Zheng Y ，Guo X ... -  《Frontiers in Immunology》

Bioinformatics identification and validation of biomarkers and infiltrating immune cells in endometriosis.

Endometriosis (EM) is a common gynecological disorder that often leads to irregular menstruation and infertility. The pathogenesis of EM remains unclear and delays in diagnosis are common. Thus, it is urgent to explore potential biomarkers and underlying molecular mechanisms for EM diagnosis and therapies. Three EM-related datasets (GSE11691, GSE25628, and GSE86534) were downloaded from the Gene Expression Omnibus (GEO) which were integrated into a combined dataset after removing batch effect. Differentially expressed immune cell-related genes were obtained by CIBERSORT, WGCNA, and the identification of differentially expressed genes. Random forest model (RF), support vector machine model (SVM), and generalized linear model (GLM) were then constructed and the biomarkers for EM were determined. A nomogram evaluating the risk of disease was constructed and the validity was assessed by the calibration curve, DCA curve, and clinical impact curve. Single-gene Gene Set Enrichment Analysis (GSEA)was performed to explore the molecular mechanisms of biomarkers. The ceRNA regulatory network of biomarkers was created by Cytoscape and potential target drugs were obtained in the DGIdb database (Drug-Gene Interaction database).The expression levels of biomarkers from clinical samples was quantified by RT-qPCR. The ratio of eight immune cells was significantly different between the eutopic and ectopic endometrium samples. A total of eight differentially expressed immune cell-related genes were investigated. The SVM model was a relatively suitable model for the prediction of EM and five genes (CXCL12, PDGFRL, AGTR1, PTGER3, and S1PR1) were selected from the model as biomarkers. The calibration curve, DCA curve, and clinical impact curve indicated that the nomogram based on the five biomarkers had a robust ability to predict disease. Single gene GSEA result suggested that all five biomarkers were involved in labyrinthine layer morphogenesis and transmembrane transport-related biological processes in EM. A ceRNA regulatory network containing 184 nodes and 251 edges was constructed. Seven drugs targeting CXCL12, 49 drugs targeting AGTR1, 16 drugs targeting PTGER3, and 21 drugs targeting S1PR1 were extracted as potential drugs for EM therapy. Finally, the expression of PDGFRL and S1PR1 in clinical samples was validated by RT-qPCR, which was consistent with the result of public database. In summary, we identified five biomarkers (CXCL12, PDGFRL, AGTR1, PTGER3, and S1PR1) and constructed diagnostic model, furthermore predicted the potential therapeutic drugs for EM. Collectively, these findings provide new insights into EM diagnosis and treatment.

Jiang H ，Zhang X ，Wu Y ，Zhang B ，Wei J ，Li J ，Huang Y ，Chen L ，He X ... -  《Frontiers in Immunology》

Identification of characteristic genes and construction of regulatory network in gallbladder carcinoma.

Gallbladder carcinoma (GBC) is a highly malignant tumor with a poor overall prognosis. This study aimed to identify the characteristic microRNAs (miRNAs) of GBC and the competing endogenous RNA (ceRNA) regulatory mechanisms. The microarray data of GBC tissue samples and normal gallbladder (NGB) tissue samples from the Gene Expression Omnibus (GEO) database was downloaded. GBC-related differentially expressed miRNAs (DE-miRNAs) were identified by inter-group differential expression analysis and weighted gene co-expression network analysis (WGCNA). Machine learning algorithms were used to screen the characteristic miRNA based on the intersect between least absolute shrinkage and selection operator (LASSO) and Support vector machine-recursive feature elimination (SVM-RFE). Based on the differential expression analysis of GEO database, the ceRNA network of characteristic miRNA was predicted and constructed. The biological functions of the ceRNA network were revealed by carrying out the gene enrichment analysis was implemented. We further screened the key genes of ceRNA network and constructed a protein-protein interaction (PPI) network, and predicted and generated the transcription factors (TFs) network of signature miRNAs. The expression of characteristic miRNA in clinical samples was verified by quantitative real-time polymerase chain reaction (qRT-PCR). A total of 131 GBC-related DE-miRNAs were obtained. The hsa-miR-4770 was defined as characteristic miRNA for GBC. The ceRNA network containing 211 mRNAs, one miRNA, two lncRNAs, and 48 circRNAs was created. Gene enrichment analysis suggested that the downstream genes were mainly involved in actin filament organization, cell-substrate adhesion, cell-matrix adhesion, reactive oxygen species metabolic process, glutamine metabolic process and extracellular matrix (ECM)-receptor interaction pathway. 10 key genes in the network were found to be most correlated with disease, and involved in cell cycle-related processes, p53, and extrinsic apoptotic signaling pathways. The qRT-PCR result demonstrated that hsa-miR-4770 is down-regulated in GBC, and the expression trend is consistent with the public database. We identified hsa-miR-4770 as the characteristic miRNA for GBC. The ceRNA network of hsa-miR-4770 may play key roles in GBC. This study provided some basis for potential pathogenesis of GBC.

Shao H ，Zhu J ，Zhu Y ，Liu L ，Zhao S ，Kang Q ，Liu Y ，Zou H ... -  《BMC Medical Genomics》

Identification of immune- and autophagy-related genes and effective diagnostic biomarkers in endometriosis: a bioinformatics analysis.

To identify autophagy- and immune-related hub genes affecting the diagnosis and treatment of endometriosis. Gene expression data were downloaded from the Gene Expression Omnibus (GEO) (GSE11691 and GSE120103 for training, and GSE7305 for validation). By overlapping the differentially expressed genes (DEGs), Weighted gene co-expression network analysis (WGCNA) module genes, and autophagy-related genes (ARGs), and immune-related genes (IRGs) separately, hub genes were identified using the least absolute shrinkage and selection operator (LASSO)and support vector machine recursive feature elimination (SVM-RFE). The hub genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. A hub gene-prediction model was constructed and assessed using five-fold cross-validation via five supervised machine-learning algorithms: random forest, the sequential minimal optimization (SMO), K-nearest neighbours (IBK), C4.5 decision tree (J48), and logistics regression. The area under the receiver operating characteristic curve (AUC) was adopted to assess the identification ability of characteristic genes. 1,116 DEGs were obtained from the training cohort, and 22 endometriosis-related IRGs were identified by overlapping the 1,116 DEGs, 3,222 module genes, and 1,793 IRGs. Meanwhile, 45 endometriosis-related ARGs were obtained (1,928 ARGs). Subsequently, nine IRG hub genes (BST2, CCL13, CD86, CSF1, FAM3C, GREM1, ISG20, PSMB8, and S100A11) and nine ARG hub genes (GSK3A, HTR2B, RAB3GAP1, ARFIP2, BNIP3, CSF1, MAOA, PPP1R13L, and SH3GLB2) were obtained by LASSO and SVM-RFE. GO analysis indicated that the ARG hub genes responded to the regulation of autophagy and mitochondrial outer membrane permeabilization, and KEGG enrichment analysis involved serotonergic and dopaminergic synapses. GO analysis also indicated that the IRG hub genes responded to the regulation of leukocyte proliferation and mononuclear cell migration, and KEGG analysis showed enrichment involved in viral protein interaction with cytokines and cytokine receptors. The AUC of the random-forest algorithm of ARGs was 0.975 in the training cohort and 0.940 in the validation cohort, and the AUC of the SMO algorithm of IRGs was 0.907 in the training cohort and 0.8 in the validation cohort. Seventeen hub genes are closely associated with endometriosis. These genes are potential autophagy- and immune-related biomarkers for diagnosis and treatment of endometriosis.

Ji X ，Huang C ，Mao H ，Zhang Z ，Zhang X ，Yue B ，Li X ，Wu Q ... -  《-》