Bacillus cereus NJ01 induces SA- and ABA-mediated immunity against bacterial pathogens through the EDS1-WRKY18 module.

Emerging evidence suggests a beneficial role of rhizobacteria in ameliorating plant disease resistance in an environment-friendly way. In this study, we characterize a rhizobacterium, Bacillus cereus NJ01, that enhances bacterial pathogen resistance in rice and Arabidopsis. Transcriptome analyses show that root inoculation of NJ01 induces the expression of salicylic acid (SA)- and abscisic acid (ABA)-related genes in Arabidopsis leaves. Genetic evidence showed that EDS1, PAD4, and WRKY18 are required for B. cereus NJ01-induced bacterial resistance. An EDS1-PAD4 complex interacts with WRKY18 and enhances its DNA binding activity. WRKY18 directly binds to the W box in the promoter region of the SA biosynthesis gene ICS1 and ABA biosynthesis genes NCED3 and NCED5 and contributes to the NJ01-induced bacterial resistance. Taken together, our findings indicate a role of the EDS1/PAD4-WRKY18 complex in rhizobacteria-induced disease resistance.

Fusarium graminearum rapid alkalinization factor peptide negatively regulates plant immunity and cell growth via the FERONIA receptor kinase.

The plant rapid alkalinization factor (RALF) peptides function as key regulators in cell growth and immune responses through the receptor kinase FERONIA (FER). In this study, we report that the transcription factor FgPacC binds directly to the promoter of FgRALF gene, which encodes a functional homologue of the plant RALF peptides from the wheat head blight fungus Fusarium graminearum (FgRALF). More importantly, FgPacC promotes fungal infection via host immune suppression by activating the expression of FgRALF. The FgRALF peptide also exhibited typical activities of plant RALF functions, such as inducing plant alkalinization and inhibiting cell growth, including wheat (Triticum aestivum), tomato (Solanum lycopersicum) and Arabidopsis thaliana. We further identified the wheat receptor kinase FERONIA (TaFER), which is capable of restoring the defects of the A. thaliana FER mutant. In addition, we found that FgRALF peptide binds to the extracellular malectin-like domain (ECD) of TaFER (TaFERECD) to suppress the PAMP-triggered immunity (PTI) and cell growth. Overexpression of TaFERECD in A. thaliana confers plant resistance to F. graminearum and protects from FgRALF-induced cell growth inhibition. Collectively, our results demonstrate that the fungal pathogen-secreted RALF mimic suppresses host immunity and inhibits cell growth via plant FER receptor. This establishes a novel pathway for the development of disease-resistant crops in the future without compromising their yield potential.

Nanopore long-read RNA sequencing reveals functional alternative splicing variants in human vascular smooth muscle cells.

Vascular smooth muscle cells (VSMCs) are the major contributor to vascular repair and remodeling, which showed high level of phenotypic plasticity. Abnormalities in VSMC plasticity can lead to multiple cardiovascular diseases, wherein alternative splicing plays important roles. However, alternative splicing variants in VSMC plasticity are not fully understood. Here we systematically characterized the long-read transcriptome and their dysregulation in  human aortic smooth muscle cells (HASMCs) by employing the Oxford Nanopore Technologies long-read RNA sequencing in HASMCs that are separately treated with platelet-derived growth factor, transforming growth factor, and hsa-miR-221-3P transfection. Our analysis reveals frequent alternative splicing events and thousands of unannotated transcripts generated from alternative splicing. HASMCs treated with different factors exhibit distinct transcriptional reprogramming modulated by alternative splicing. We also found that unannotated transcripts produce different open reading frames compared to the annotated transcripts. Finally, we experimentally validated the unannotated transcript derived from gene CISD1, namely CISD1-u, which plays a role in the phenotypic switch of HASMCs. Our study characterizes the phenotypic modulation of HASMCs from an insight of long-read transcriptome, which would promote the understanding and the manipulation of HASMC plasticity in cardiovascular diseases.

AUF1-induced circular RNA hsa_circ_0010467 promotes platinum resistance of ovarian cancer through miR-637/LIF/STAT3 axis.

BACKGROUND: Increasing evidences has indicated that primary and acquired resistance of ovarian cancer (OC) to platinum is mediated by multiple molecular and cellular factors. Understanding these mechanisms could promote the therapeutic efficiency for patients with OC. METHODS: Here, we screened the expression pattern of circRNAs in samples derived from platinum-resistant and platinum-sensitive OC patients using RNA-sequencing (RNA-seq). The expression of hsa_circ_0010467 was validated by Sanger sequencing, RT-qPCR, and fluorescence in situ hybridization (FISH) assays. Overexpression and knockdown experiments were performed to explore the function of hsa_circ_0010467. The effects of hsa_circ_0010467 on enhancing platinum treatment were validated in OC cells, mouse model and patient-derived organoid (PDO). RNA pull-down, RNA immunoprecipitation (RIP), and dual-luciferase reporter assays were performed to investigate the interaction between hsa_circ_0010467 and proteins. RESULTS: Increased expression of hsa_circ_0010467 is observed in platinum-resistant OC cells, tissues and serum exosomes, which is positively correlated with advanced tumor stage and poor prognosis of OC patients. Hsa_circ_0010467 is found to maintain the platinum resistance via inducing tumor cell stemness, and silencing hsa_circ_0010467 substantially increases the efficacy of platinum treatment on inhibiting OC cell proliferation. Further investigation reveals that hsa_circ_0010467 acts as a miR-637 sponge to mediate the repressive effect of miR-637 on leukemia inhibitory factor (LIF) and activates the LIF/STAT3 signaling pathway. We further discover that AUF1 could promote the biogenesis of hsa_circ_0010467 in OC. CONCLUSION: Our study uncovers the mechanism that hsa_circ_0010467 mediates the platinum resistance of OC through AUF1/hsa_circ_0010467/miR-637/LIF/STAT3 axis, and provides potential targets for the treatment of platinum-resistant OC patients.

CURLY LEAF modulates apoplast liquid water status in Arabidopsis leaves.

The apoplast of plant leaves, the intercellular space between mesophyll cells, is normally largely filled with air with a minimal amount of liquid water in it, which is essential for key physiological processes such as gas exchange to occur. Phytopathogens exploit virulence factors to induce a water-rich environment, or "water-soaked" area, in the apoplast of the infected leaf tissue to promote disease. We propose that plants evolved a "water soaking" pathway, which normally keeps a nonflooded leaf apoplast for plant growth but is disturbed by microbial pathogens to facilitate infection. Investigation of the "water soaking" pathway and leaf water control mechanisms is a fundamental, yet previously overlooked, aspect of plant physiology. To identify key components in the "water soaking" pathway, we performed a genetic screen to isolate Arabidopsis (Arabidopsis thaliana) severe water soaking (sws) mutants that show liquid water overaccumulation in the leaf under high air humidity, a condition required for visible water soaking. Here, we report the sws1 mutant, which displays rapid water soaking upon high humidity treatment due to a loss-of-function mutation in CURLY LEAF (CLF), encoding a histone methyltransferase in the POLYCOMB REPRESSIVE COMPLEX 2 (PRC2). We found that the sws1 (clf) mutant exhibits enhanced abscisic acid (ABA) levels and stomatal closure, which are indispensable for its water soaking phenotype and mediated by CLF's epigenetic regulation of a group of ABA-associated NAM, ATAF, and CUC (NAC) transcription factor genes, NAC019/055/072. The clf mutant showed weakened immunity, which likely also contributes to the water soaking phenotype. In addition, the clf plant supports a substantially higher level of Pseudomonas syringae pathogen-induced water soaking and bacterial multiplication, in an ABA pathway and NAC019/055/072-dependent manner. Collectively, our study sheds light on an important question in plant biology and demonstrates CLF as a key modulator of leaf liquid water status via epigenetic regulation of the ABA pathway and stomatal movement.

Long-read transcriptome sequencing reveals allele-specific variants at high resolution.

Disturbance in the regulation of transcript structure plays a crucial role in human disease. In a recent study, Glinos et al. characterized allele-specific transcript alterations in long-read RNA sequencing (RNA-seq) data derived from multiple human tissues and provide a high-resolution view of how disease-associated genetic variants affect transcript structure.

Systematic characterization of cancer transcriptome at transcript resolution.

Transcribed RNAs undergo various regulation and modification to become functional transcripts. Notably, cancer transcriptome has not been fully characterized at transcript resolution. Herein, we carry out a reference-based transcript assembly across >1000 cancer cell lines. We identify 498,255 transcripts, approximately half of which are unannotated. Unannotated transcripts are closely associated with cancer-related hallmarks and show clinical significance. We build a high-confidence RNA binding protein (RBP)-transcript regulatory network, wherein most RBPs tend to regulate transcripts involved in cell proliferation. We identify numerous transcripts that are highly associated with anti-cancer drug sensitivity. Furthermore, we establish RBP-transcript-drug axes, wherein PTBP1 is experimentally validated to affect the sensitivity to decitabine by regulating KIAA1522-a6 transcript. Finally, we establish a user-friendly data portal to serve as a valuable resource for understanding cancer transcriptome diversity and its potential clinical utility at transcript level. Our study substantially extends cancer RNA repository and will facilitate anti-cancer drug discovery.

FOXC2 Promotes Vasculogenic Mimicry in Ovarian Cancer.

FOXC2 is a forkhead family transcription factor that plays a critical role in specifying mesenchymal cell fate during embryogenesis. FOXC2 expression is associated with increased metastasis and poor survival in various solid malignancies. Using in vitro and in vivo assays in mouse ovarian cancer cell lines, we confirmed the previously reported mechanisms by which FOXC2 could promote cancer growth, metastasis, and drug resistance, including epithelial-mesenchymal transition, stem cell-like differentiation, and resistance to anoikis. In addition, we showed that FOXC2 expression is associated with vasculogenic mimicry in mouse and human ovarian cancers. FOXC2 overexpression increased the ability of human ovarian cancer cells to form vascular-like structures in vitro, while inhibition of FOXC2 had the opposite effect. Thus, we present a novel mechanism by which FOXC2 might contribute to cancer aggressiveness and poor patient survival.

Focal Serous Tubal Intra-Epithelial Carcinoma Lesions Are Associated With Global Changes in the Fallopian Tube Epithelia and Stroma.

Objective: Serous tubal intra-epithelial carcinoma (STIC) lesions are thought to be precursors to high-grade serous ovarian cancer (HGSOC), but HGSOC is not always accompanied by STIC. Our study was designed to determine if there are global visual and subvisual microenvironmental differences between fallopian tubes with and without STIC lesions. Methods: Computational image analyses were used to identify potential morphometric and topologic differences in stromal and epithelial cells in samples from three age-matched groups of fallopian tubes. The Benign group comprised normal fallopian tubes from women with benign conditions while the STIC and NoSTIC groups consisted of fallopian tubes from women with HGSOC, with and without STIC lesions, respectively. For the morphometric feature extraction and analysis of the stromal architecture, the image tiles in the STIC group were further divided into the stroma away from the STIC (AwaySTIC) and the stroma near the STIC (NearSTIC). QuPath software was used to identify and quantitate secretory and ciliated epithelial cells. A secretory cell expansion (SCE) or a ciliated cell expansion (CCE) was defined as a monolayered contiguous run of >10 secretory or ciliated cells uninterrupted by the other cell type. Results: Image analyses of the tubal stroma revealed gradual architectural differences from the Benign to NoSTIC to AwaySTIC to NearSTIC groups. In the epithelial topology analysis, the relative number of SCE and the average number of cells within SCE were higher in the STIC group than in the Benign and NoSTIC groups. In addition, aging was associated with an increased relative number of SCE and a decreased relative number of CCE. ROC analysis determined that an average of 15 cells within SCE was the optimal cutoff value indicating the presence of a STIC lesion in the tubal epithelium. Conclusions: Our findings suggest that global stromal alterations and age-associated reorganization of tubal secretory and ciliated cells are associated with STIC lesions. Further studies will need to determine if these alterations precede STIC lesions and provide permissible conditions for the formation of STIC.

Bacterial effectors manipulate plant abscisic acid signaling for creation of an aqueous apoplast.

Phytopathogens like Pseudomonas syringae induce "water soaking" in the apoplastic space of plant leaf tissue as a key virulence mechanism. Water soaking is commonly observed in diverse pathosystems, yet the underlying physiological basis remains largely elusive. Here, we show that one of the strong P. syringae water-soaking inducers, AvrE, alters the regulation of abscisic acid (ABA) to induce ABA signaling, stomatal closure, and, thus, water soaking. AvrE binds and inhibits the function of Arabidopsis type one protein phosphatases (TOPPs), which negatively regulate ABA by suppressing SnRK2s, a key node of the ABA signaling pathway. The topp12537 quintuple mutants display significantly enhanced water soaking after P. syringae inoculation, whereas the loss of the ABA pathway dampens P. syringae-induced water soaking and disease. Our study uncovers the hijacking of ABA signaling and stomatal closure by P. syringae effectors as key mechanisms of disease susceptibility.

LncPep: A Resource of Translational Evidences for lncRNAs.

Long noncoding RNAs (lncRNAs) are a type of transcript that is >200 nucleotides long with no protein-coding capacity. Accumulating studies have suggested that lncRNAs contain open reading frames (ORFs) that encode peptides. Although several noncoding RNA-encoded peptide-related databases have been developed, most of them display only a small number of experimentally validated peptides, and resources focused on lncRNA-encoded peptides are still lacking. We used six types of evidence, coding potential assessment tool (CPAT), coding potential calculator v2.0 (CPC2), N6-methyladenosine modification of RNA sites (m6A), Pfam, ribosome profiling (Ribo-seq), and translation initiation sites (TISs), to evaluate the coding potential of 883,804 lncRNAs across 39 species. We constructed a comprehensive database of lncRNA-encoded peptides, LncPep (http://www.shenglilabs.com/LncPep/). LncPep provides three major functional modules: 1) user-friendly searching/browsing interface, 2) prediction and BLAST modules for exploring novel lncRNAs and peptides, and 3) annotations for lncRNAs, peptides and supporting evidence. Taken together, LncPep is a user-friendly and convenient platform for discovering and investigating peptides encoded by lncRNAs.

Dissecting the single-cell transcriptomeunderlying chronic liver injury.

Chronic liver disease (CLD) is currently a major health problem worldwide, which is accompanied by chronic liver injury and lack of clinically effective treatment; however, systematic characterization of chronic liver injury procedures at single-cell resolution is lacking. In the present study, we established chronic liver injury mouse models and conducted single-cell RNA sequencing (scRNA-seq), including choline-deficient, ethionine-supplemented (CDE) and 3,5-diethoxycarbonyl 1,4-dihydrocollidinen (DDC) mouse models. We captured in total 16,389 high-quality cells and identified 12 main cell types in scRNA-seq data. Macrophages and endothelial cells are the largest cell populations in our dataset. Transcriptional trajectory analysis revealed different expression patterns of cells between CDE and DDC models and identified potential liver injury markers, such as Ets1, Gda, Itgam, and Sparc. Differential analysis identified 25 and 152 differentially expressed genes in CDE and DDC macrophages, respectively. In addition, 413 genes were detected to exclusively express in specific pseudotime states of macrophages. These genes were found to participate in immune-related biological processes. Further cell-cell communication analysis found extensive receding of cell-cell interactions between different cell types in the liver injury process, especially in the DDC model. Our study characterized the single-cell transcriptional landscape in the process of chronic liver injury, promoting the understanding of the underlying molecular mechanisms and providing candidate clinical strategy for effective intervention of chronic liver diseases.

Maximizing the Utility of Transcriptomics Data in Inflammatory Skin Diseases.

Inflammatory skin diseases are induced by disorders of the host defense system of the skin, which is composed of a barrier, innate and acquired immunity, as well as the cutaneous microbiome. These disorders are characterized by recurrent cutaneous lesions and intense itch, which seriously affecting life quality of people across all ages and ethnicities. To elucidate molecular factors for typical inflammatory skin diseases (such as psoriasis and atopic dermatitis), transcriptomic profiling assays have been largely performed. Additionally, single-cell RNA sequencing (scRNA-seq) as well as spatial transcriptomic profiling have revealed multiple potential translational targets and offered guides to improve diagnosis and treatment strategies for inflammatory skin diseases. High-throughput transcriptomics data has shown unprecedented power to disclose the complex pathophysiology of inflammatory skin diseases. Here, we will summarize discoveries from transcriptomics data and discuss how to maximize the transcriptomics data to propel the development of diagnostic biomarkers and therapeutic targets in inflammatory skin diseases.

Integrative analysis reveals clinically relevant molecular fingerprints in pancreatic cancer.

Pancreatic cancer is a highly aggressive cancer with an exceedingly low rate of response to treatments, which calls for comprehensive molecular characterization of pancreatic cancer cell lines (PCCLs). We screened multi-layer molecular data of 36 PCCLs, including gene mutation, gene expression, microRNA (miRNA) expression, and protein profiles. Our comparative analysis of genomic mutations found that PCCLs recapitulated genomic alterations of the primary tumor and suggested potential therapeutic strategies for clinical interventions. The panel of 36 PCCLs was classified into 2 subgroups based on transcriptomic mRNA expression, wherein the C1 subgroup was characterized with differentiation, whereas C2 cell lines were featured with immunity, angiogenesis, epidermis, and proliferation. Transcriptomic classification was further recapitulated by miRNA and protein expression. Additionally, the differential proteins between C1 and C2 subgroups were prominently involved in epidermal growth factor receptor (EGFR) signaling, phosphatidylinositol 3-kinase (PI3K) signaling, and mitogen-activated protein kinase (MAPK) signaling pathways. Tumor samples from different subgroups exhibited distinct infiltration of CD4 naive cells and monocytes. Remarkably, patients in subgroups C1 showed longer survival, whereas those in C2 had worse clinical outcome. Further integrative analysis revealed that temozolomide and NVP-TAE684 showed higher sensitivity in the C1 subgroup, whereas the C2 cell lines were more sensitive to SR1001 and SRT-1720. Our results also showed that PCCLs with mutations in CDKN2A, TP53, and SMAD4 were more sensitive to certain anti-cancer drugs. Our integrative analysis identified molecular features of pancreatic cancer that were associated with clinical significance and drug sensitivity, providing potentially effective strategies for precision treatments of patients with pancreatic cancer.

Identification of Key Genes and Pathways associated with Endometriosis by Weighted Gene Co-expression Network Analysis.

Background: Endometriosis is a common gynecological disorder with high rates of infertility and pelvic pain. However, its pathogenesis and diagnostic biomarkers remain unclear. This study aimed to elucidate potential hub genes and key pathways associated with endometriosis in ectopic endometrium (EC) and eutopic endometrium (EU). Material and Method: EC and EU-associated microarray datasets were obtained from the gene expression omnibus (GEO) database. Gene set enrichment analysis was performed to obtain further biological insight into the EU and EC-associated genes. Weighted gene co-expression network analysis (WGCNA) was performed to find clinically significant modules of highly-correlated genes. The hub genes that belong to both the weighted gene co-expression network and protein-protein interaction (PPI) network were identified using a Venn diagram. Results: We obtained EC and EU-associated microarray datasets GSE7305 and GSE120103. Genes in the EC were mainly enriched in the immune response and immune cell trafficking, and genes in the EU were mainly enriched in stress response and steroid hormone biosynthesis. PPI networks and weighted gene co-expression networks were constructed. An EC-associated blue module and an EU-associated magenta module were identified, and their function annotations revealed that hormone receptor signaling or inflammatory microenvironments may promote EU passing through the oviducts and migrating to the ovarian surfaces, and adhesion and immune correlated genes may induce the successful ectopic implantation of the endometrium (EC). Twelve hub genes in the EC and sixteen hub genes in the EU were recognized and further validated in independent datasets. Conclusion: Our study identified, for the first time, the hub genes and enrichment pathways in the EC and EU using WGCNA, which may provide a comprehensive understanding of the pathogenesis of endometriosis and have important clinical implications for the treatment and diagnosis of endometriosis.

The predictive and prognostic effects of PD-L1 expression on TKI treatment and survival of EGFR-mutant NSCLC: A meta-analysis.

ABSTRACT: Whether programmed death-ligand 1 (PD-L1) expression could predict the outcome of tyrosine kinase inhibitor (TKI) treatment and prognosis of epidermal growth factor receptor (EGFR)-mutant nonsmall cell lung cancer (NSCLC) is remaining controversial.Potential studies were search from PubMed, Embase, and Web of Science databases. Pooled odds ratio of objective response rate was used to describe the relationship between PD-L1 expression and primary resistance to EGFR-TKIs. Pooled hazard ratios (HRs) of progression-free survival (PFS) and overall survival (OS) were included to assess the effects of PD-L1 status on the outcome of EGFR-TKI treatment and survival of EGFR-mutant NSCLCs.Eighteen eligible studies (1986 EGFR-mutant NSCLCs) were included in this meta-analysis. Positive PD-L1 expression correlated with lower objective response rate of EGFR-TKI treatment (odds ratio [95% confidence interval {CI}] = 0.52 [0.28-0.98], P = .043), while PFS (adjusted HR [95% CI] = 1.49 [0.96-1.89], P = .332) and OS (HR [95% CI] = 1.24 [0.70-2.20], P = .456) of EGFR-TKI treatment did not correlated with PD-L1 status. Furthermore, PD-L1 expression was not a predictive biomarker for the OS (HR [95% CI] = 1.43 [0.98-2.08], P = .062) in overall EGFR-mutant cohort.Positive PD-L1 expression indicated a higher incidence of primary resistance, but did not correlate with the PFS or OS of EGFR-TKI therapy. In addition, PD-L1 expression was unlikely a predictive biomarker for prognosis of EGFR-mutant NSCLCs.

ANKRD22 Drives Rapid Proliferation of Lgr5+ Cells and Acts as a Promising Therapeutic Target in Gastric Mucosal Injury.

BACKGROUND & AIMS: Rapid gastric epithelial progenitor cell (EPC) proliferation and inflammatory response inhibition play key roles in promoting the repair of gastric mucosal damage. However, specific targets inducing these effects are unknown. In this study, we explored the effects of a potential target, Ankyrin repeat domain 22 (ANKRD22). METHODS: An acute gastric mucosal injury model was established with Ankrd22-/- and Ankrd22+/+ mice by intragastric administration of acidified ethanol. Organoid culture and flow cytometry were performed to evaluate the effects of ANKRD22 on leucine-rich repeat-containing G-protein-coupled receptor 5-positive (Lgr5+) gastric EPC proliferation. The mechanisms by which ANKRD22 affects gastric EPC proliferation and inflammatory responses were explored by mitochondrial Ca2+ influx and immunoblotting. Candidate ANKRD22 inhibitors then were screened virtually and validated in vitro and in vivo. RESULTS: After acute gastric mucosal injury, the number of Lgr5+ gastric EPCs was increased significantly in Ankrd22-/- mice compared with that in Ankrd22+/+ mice. Moreover, Ankrd22 knockout attenuated inflammatory cell infiltration into damaged gastric tissues. ANKRD22 deletion also reduced mitochondrial Ca2+ influx and cytoplasmic nuclear factor of activated T cells in gastric epithelial cells and macrophages, which further induced Lgr5+ gastric EPC proliferation and decreased macrophage release of tumor necrosis factor-α and interleukin 1α. In addition, a small molecule, AV023, was found to show similar effects to those produced by ANKRD22 deletion in vitro. Intraperitoneal injection of AV023 into the mouse model promoted the repair of gastric mucosal damage, with increased proliferation of Lgr5+ gastric EPCs and visible relief of inflammation. CONCLUSIONS: ANKRD22 inhibition is a potential target-based therapeutic approach for promoting the repair of gastric mucosal damage.

Construction and topological analysis of an endometriosis-related exosomal circRNA-miRNA-mRNA regulatory network.

Novel biomarkers are needed to accelerate the diagnosis and treatment of endometriosis. We performed RNA sequencing to explore the expression profiles of exosomal circular RNAs (circRNAs), microRNAs (miRNAs) and mRNAs in patients with ovarian endometriomas, eutopic endometria and normal endometria. Differentially expressed genes between the different pairs of groups were analyzed and functionally annotated. Then, miRNA-target RNA pairs were identified, competing endogenous RNA (ceRNA) scores were calculated, gene expression characteristics were determined, and these parameters were used to construct an exosomal ceRNA network. We identified 36 candidate hub genes with high degrees of gene connectivity. We also topologically analyzed the ceRNA network to obtain a hub ceRNA network of circRNAs with the highest closeness and ceRNA efficiency. Twelve genes overlapped between the 36 candidate hub genes and the genes in the hub ceRNA network. These 12 genes were considered to be exosomal RNA-based biomarkers, and circ_0026129/miRNA-15a-5p/ATPase H+ transporting V1 subunit A (ATP6V1A) were at the center of the ceRNA network. By determining the exosomal RNA expression profiles of endometriosis patients and constructing a circRNA-associated ceRNA network, these findings provide insight into the molecular pathways of endometriosis and new resources for its diagnosis and treatment.

Identification of Infertility-Associated Topologically Important Genes Using Weighted Co-expression Network Analysis.

Endometriosis has been associated with a high risk of infertility. However, the underlying molecular mechanism of infertility in endometriosis remains poorly understood. In our study, we aimed to discover topologically important genes related to infertility in endometriosis, based on the structure network mining. We used microarray data from the Gene Expression Omnibus (GEO) database to construct a weighted gene co-expression network for fertile and infertile women with endometriosis and to identify gene modules highly correlated with clinical features of infertility in endometriosis. Additionally, the protein-protein interaction network analysis was used to identify the potential 20 hub messenger RNAs (mRNAs) while the network topological analysis was used to identify nine candidate long non-coding RNAs (lncRNAs). Functional annotations of clinically significant modules and lncRNAs revealed that hub genes might be involved in infertility in endometriosis by regulating G protein-coupled receptor signaling (GPCR) activity. Gene Set Enrichment Analysis showed that the phospholipase C-activating GPCR signaling pathway is correlated with infertility in patients with endometriosis. Taken together, our analysis has identified 29 hub genes which might lead to infertility in endometriosis through the regulation of the GPCR network.

LINC01018 and SMIM25 sponged miR-182-5p in endometriosis revealed by the ceRNA network construction.

The current study intended to explore the interaction of the long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) under the background of competitive endogenous RNA (ceRNA) network in endometriosis (EMs). The differentially expressed miRNAs (DEmiRs), differentially expressed lncRNA (DELs), and differentially expressed genes (DEGs) between EMs ectopic (EC) and eutopic (EU) endometrium based on three RNA-sequencing datasets (GSE105765, GSE121406, and GSE105764) were identified, which were used for the construction of ceRNA network. Then, DEGs in the ceRNA network were performed with Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analysis. Besides, the DEmiRs in the ceRNA network were validated in GSE124010. And the target DELs and DEGs of verified DEmiRs were validated in GSE86534. The correlation of verified DEmiRs, DEGs, and DELs was explored. Moreover, gene set enrichment analysis (GSEA) was applied to investigate the function of verified DEmiRs, DEGs, and DELs. Overall, 1352 DEGs and 595 DELs from GSE105764, along with 27 overlapped DEmiRs between GSE105765 and GSE121406, were obtained. Subsequently, a ceRNA network, including 11 upregulated and 16 downregulated DEmiRs, 7 upregulated and 13 downregulated DELs, 48 upregulated and 46 downregulated DEGs, was constructed. The GO and KEGG pathway analysis showed that this ceRNA network probably was associated with inflammation-related pathways. Furthermore, hsa-miR-182-5p and its target DELs (LINC01018 and SMIM25) and DEGs (BNC2, CHL1, HMCN1, PRDM16) were successfully verified in the validation analysis. Besides, hsa-miR-182-5p was significantly negatively correlated with these target DELs and DEGs. The GSEA analysis implied that high expression of LINC01018, SMIM25, and CHL1, and low expression of hsa-miR-182-5p would activate inflammation-related pathways in endometriosis EU samples.LINC01018 and SMIM25 might sponge hsa-miR-182-5p to upregulate downstream genes such as CHL1 to promote the development of endometriosis.