Integrative Bioinformatics Approaches to Map Potential Novel Genes and Pathways Involved in Ovarian Cancer

Identification and characterization of intermediate states in mammalian neural crest cell epithelial to mesenchymal transition and delamination

Epithelial to mesenchymal transition (EMT) is a cellular process that converts epithelial cells to mesenchymal cells with migratory potential in both developmental and pathological processes. Although originally considered a binary event, EMT in cancer progression involves intermediate states between a fully epithelial and a fully mesenchymal phenotype, which are characterized by distinct combinations of epithelial and mesenchymal markers. This phenomenon has been termed epithelial to mesenchymal plasticity (EMP), however, the intermediate states remain poorly described and it's unclear whether they exist during developmental EMT. Neural crest cells (NCC) are an embryonic progenitor cell population that gives rise to numerous cell types and tissues in vertebrates, and their formation is a classic example of developmental EMT. An important feature of NCC development is their delamination from the neuroepithelium via EMT, following which NCC migrate throughout the embryo and undergo differentiation. NCC delamination shares similar changes in cellular state and structure with cancer cell invasion. However, whether intermediate states also exist during NCC EMT and delamination remains unknown. Through single cell RNA sequencing, we identified intermediate NCC states based on their transcriptional signature and then spatially defined their locations in situ in the dorsolateral neuroepithelium. Our results illustrate the progressive transcriptional and spatial transitions from premigratory to migratory cranial NCC during EMT and delamination. Of note gene expression and trajectory analysis indicate that distinct intermediate populations of NCC delaminate in either S phase or G2/M phase of the cell cycle, and the importance of cell cycle regulation in facilitating mammalian cranial NCC delamination was confirmed through cell cycle inhibition studies. Additionally, transcriptional knockdown revealed a functional role for the intermediate stage marker Dlc1 in regulating NCC delamination and migration. Overall, our work identifying and characterizing the intermediate cellular states, processes, and molecular signals that regulate mammalian NCC EMT and delamination furthers our understanding of developmental EMP and may provide new insights into mechanisms regulating pathological EMP.

Fibronectin Type III Domain Containing 3B as a Potential Prognostic and Therapeutic Biomarker for Glioblastoma

Glioblastoma (GBM) is a representative malignant brain tumor characterized by a dismal prognosis, with survival rates of less than 2 years and high recurrence rates. Despite surgical resection and several alternative treatments, GBM remains a refractory disease due to its aggressive invasiveness and resistance to anticancer therapy. In this report, we explore the role of fibronectin type III domain containing 3B (FNDC3B) and its potential as a prognostic and therapeutic biomarker in GBM. GBM exhibited a significantly higher cancer-to-normal ratio compared to other organs, and patients with high FNDC3B expression had a poor prognosis (p < 0.01). In vitro studies revealed that silencing FNDC3B significantly reduced the expression of Survivin, an apoptosis inhibitor, and also reduced cell migration, invasion, extracellular matrix adhesion ability, and stem cell properties in GBM cells. Furthermore, we identified that FNDC3B regulates PTEN/PI3K/Akt signaling in GBM cells using MetaCore integrated pathway bioinformatics analysis and a proteome profiler phospho-kinase array with sequential western blot analysis. Collectively, our findings suggest FNDC3B as a potential biomarker for predicting GBM patient survival and for the development of treatment strategies for GBM.

Rare variant burden analysis from exomes of three consanguineous families reveals LILRB1 and PRSS3 as potential key proteins in inflammatory bowel disease pathogenesis

Background

Inflammatory bowel disease (IBD) is a chronic autoimmune disorder characterized by severe inflammation and mucosal destruction of the intestine. The specific, complex molecular processes underlying IBD pathogenesis are not well understood. Therefore, this study is aimed at identifying and uncovering the role of key genetic factors in IBD.

Method

The whole exome sequences (WESs) of three consanguineous Saudi families having many siblings with IBD were analyzed to discover the causal genetic defect. Then, we used a combination of artificial intelligence approaches, such as functional enrichment analysis using immune pathways and a set of computational functional validation tools for gene expression, immune cell expression analyses, phenotype aggregation, and the system biology of innate immunity, to highlight potential IBD genes that play an important role in its pathobiology.

Results

Our findings have shown a causal group of extremely rare variants in the LILRB1 (Q53L, Y99N, W351G, D365A, and Q376H) and PRSS3 (F4L and V25I) genes in IBD-affected siblings. Findings from amino acids in conserved domains, tertiary-level structural deviations, and stability analysis have confirmed that these variants have a negative impact on structural features in the corresponding proteins. Intensive computational structural analysis shows that both genes have very high expression in the gastrointestinal tract and immune organs and are involved in a variety of innate immune system pathways. Since the innate immune system detects microbial infections, any defect in this system could lead to immune functional impairment contributing to IBD.

Conclusion

The present study proposes a novel strategy for unraveling the complex genetic architecture of IBD by integrating WES data of familial cases, with computational analysis.

Bioinformatics and Next-Generation Data Analysis for Identification of Genes and Molecular Pathways Involved in Subjects with Diabetes and Obesity

Background and Objectives: A subject with diabetes and obesity is a class of the metabolic disorder. The current investigation aimed to elucidate the potential biomarker and prognostic targets in subjects with diabetes and obesity. Materials and Methods: The next-generation sequencing (NGS) data of GSE132831 was downloaded from Gene Expression Omnibus (GEO) database. Functional enrichment analysis of DEGs was conducted with ToppGene. The protein-protein interactions network, module analysis, target gene-miRNA regulatory network and target gene-TF regulatory network were constructed and analyzed. Furthermore, hub genes were validated by receiver operating characteristic (ROC) analysis. A total of 872 DEGs, including 439 up-regulated genes and 433 down-regulated genes were observed. Results: Second, functional enrichment analysis showed that these DEGs are mainly involved in the axon guidance, neutrophil degranulation, plasma membrane bounded cell projection organization and cell activation. The top ten hub genes (MYH9, FLNA, DCTN1, CLTC, ERBB2, TCF4, VIM, LRRK2, IFI16 and CAV1) could be utilized as potential diagnostic indicators for subjects with diabetes and obesity. The hub genes were validated in subjects with diabetes and obesity. Conclusion: This investigation found effective and reliable molecular biomarkers for diagnosis and prognosis by integrated bioinformatics analysis, suggesting new and key therapeutic targets for subjects with diabetes and obesity.

Whole-genome sequencing identifies potential candidate genes for egg production traits in laying ducks (Anas platyrhynchos)

Egg production traits are economically important in laying ducks. Genetic molecular mechanisms and candidate genes underlying these traits remain unclear. In this study, whole genome variants were identified through whole-genome resequencing using three high-egg producing (HEN) and three low-egg producing (LEN) laying ducks. The gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genome (KEGG) pathways for the genes of common differential variants between HEN and LEN ducks were determined. Frizzled class receptor 6 (FZD6) was further genotyped using the Sequenom MassARRAY iPLEX platform. The association of FZD6 gene polymorphisms with 73 egg production and weight traits in 329 female ducks were estimated. A total of 65,535 single nucleotide polymorphisms (SNPs) and 4,702 indels were identified across the genome. Fourteen GO terms and 14 KEGG pathways were determined for the genes of common differential variants, including MAPK signaling, Wnt signaling, melanogenesis and calcium signaling pathways, which are key functional pathways for poultry egg production reported in previous reports. Further analysis showed that 27 SNPs of FZD6 were associated with three early egg production of duck and egg weight traits, including egg production at 17 weeks (EP17), 18 weeks (EP18) and 19 weeks (EP19) and egg weight at 59 weeks (EW59). The FZD6 should be considered a novel candidate gene for egg production traits in laying ducks.

Identification of dysregulated canonical pathways associated with pathogenesis and progression of Amyotrophic Lateral Sclerosis-An integrated bioinformatics approach

The mechanisms responsible for the pathogenesis and progression of Amyotrophic Lateral Sclerosis (ALS) remain poorly understood, making the diagnosis of ALS challenging. We aimed to find the novel gene biomarkers via computationally analyzing microarray expression studies, in three different cell lineages, namely myotube cells, astrocyte cells and oligodendrocyte cells. Microarray gene expression profiles were obtained and analyzed for three cell types: myotube cell lineage (GSE122261), astrocyte, and oligodendrocyte cell lineage (GSE87385). A comprehensive computational pipeline, tailored explicitly for microarray gene expression profiling studies, was devised to analyze the sample groups, wherein the myotube sample group comprised of six control (GSM3462697, GSM3462698, GSM3462699, GSM3462700, GSM3462701, GSM3462702) & six diseased (GSM3462691, GSM3462692, GSM3462693, GSM3462694, GSM3462695, GSM3462696) samples were considered. Similarly, for the astrocyte sample group two samples each for the control (GSM2330040, GSM2330042) and the diseased (GSM2330039, GSM2330041), and for the oligodendrocyte sample group, 2 control (GSM2330043, GSM2330045) samples and two diseased (GSM2330044, GSM2330046) samples were considered for the current study. The in-depth interaction of these DEGs was studied using MCODE and subjected to preliminary functional analysis using ClueGO/CluePedia plug-in. Qiagen's IPA software was employed for enrichment analysis, which generated the key canonical pathways and a list of potential biomarker molecules specific to each sample group. The preliminary analysis yielded 512 DEGs across all 3-sample groups, wherein 139 DEGs belonged to the myotube sample group, 216 DEGs for the astrocyte sample group, and 157 DEGs for the oligodendrocytes sample group. The data suggests growth hormone signaling and its activity, ErbB signaling pathway, and JAK/STAT signaling pathway are some of the pathways that are significantly dysregulated and play a crucial role in the development and progression of ALS. KISS1R and CSHL1 are potential genes that could act as diagnostic biomarkers in myotube cell types. Also, KRAS, TGFB2, JUN, and SMAD6 genes may be used as prognostic biomarkers to differentiate between early and late-stage ALS-diseased patients.

Analysis of signaling cascades from myeloma cells treated with pristimerin

Multiple myeloma (MM) is the 2nd most frequently diagnosed blood cancer after non-Hodgkin's lymphoma. The present study aimed to identify the differentially expressed genes (DEGs) between the control and pristimerin-treated MM cell lines. We examined the GSE14011 microarray dataset and screened DEGs with GEO2R statistical tool using the inbuilt limma package. We used a bioinformatics pipeline to identify the differential networks, signaling cascades, and the survival of the hub genes. We implemented two different enrichment analysis including ClueGO and Metacore™, to get accurate annotation for most significant DEGs. We screened the most significant 408 DEGs from the dataset based on p-values and logFC values. Using protein network analysis, we found the genes UBC, HSP90AB1, HSPH1, HSPA1B, HSPA1L, HSPA6, HSPD1, DNAJB1, HSPE1, DNAJC10, BAG3, and DNAJC7 had higher node degree distribution. In contrast, the functional annotation provided that the DEGs were predominantly enriched in B-cell receptor signaling, unfolded protein response, positive regulation of phagocytosis, HSP70, and HSP40-dependent folding, and ubiquitin-proteasomal proteolysis. Using network algorithms, and comparing enrichment analysis, we found the hub genes enriched were INHBE, UBC, HSPA1A, HSP90AB1, IKBKB, and BAG3. These DEGs were further validated with overall survival and gene expression analysis between the tumor and control groups. Finally, pristimerin effects were validated independently in a cell line model consisting of IM9 and U266 MM cells. Pristimerin induced in vitro cytotoxicity in MM cells in a dose-dependent manner. Pristimerin inhibited NF-κB, induced accumulation of ubiquitinated proteins and inhibited HSP60 in the validation of bioinformatics findings, while pristimerin-induced caspase-3 and PARP cleavage confirmed cell death. Taken together, we found that the identified DEGs were strongly associated with the apoptosis induced in MM cell lines due to pristimerin treatment, and combinatorial therapy derived from pristimerin could act as novel anti-myeloma multifunctional agents.

Genomic analysis of matrix metalloproteinases affecting the prognosis and immunogenic profile of gastric cancer

This study systematically and comprehensively analyzed the characteristics of matrix metalloproteinases (MMPs) in gastric cancer (GC) and revealed the relationship between MMPs and prognoses, clinicopathological features, tumor microenvironment, gene mutations, and drug therapy response in patients with GC. Based on the mRNA expression profiles of 45 MMP-related genes in GC, we established a model that classified GC patients into three groups based on cluster analysis of the mRNA expression profiles. The 3 groups of GC patients showed significantly different prognoses as well as tumor microenvironmental characteristics. Next, we used Boruta's algorithm and PCA method to establish an MMP scoring system and found that lower MMP scores were associated with better prognoses, lower clinical stages, better immune cell infiltration, lower degrees of immune dysfunction and rejection, and more genetic mutations. Whereas a high MMP score was the opposite. These observations were further validated with data from other datasets, showing the robustness of our MMP scoring system. Overall, MMP could be involved in the tumor microenvironment (TME), clinical features, and prognosis of GC. An in-depth study of MMP patterns can better understand the indispensable role of MMP in the development of GC and reasonably assess the survival prognosis, clinicopathological features, and drug efficacy of different patients, thus providing clinicians with a broader vision of GC progression and treatment.

Liver transcriptomics reveals microRNA features of the host response in a mouse model of dengue virus infection

BACKGROUND

Organ dysfunction, especially liver injury, caused by dengue virus (DENV) infection has been associated with fatal cases in dengue patients around the world. However, the pathophysiological mechanisms of liver involvement in dengue remain unclear. There is accumulating evidence that miRNAs are playing an important role in regulating viral pathogenesis, and it can help in diagnostic and anti-viral therapies development.

METHODS

We collected liver tissues of DENV-infected for small RNA sequencing to identify significantly different express miRNAs during dengue virus infection, and the identified target genes of these miRNAs were annotated by biological function and pathway enrichment.

RESULTS

31 significantly altered miRNAs were identified, including 16 up-regulated and 15 down-regulated miRNAs. By performing a series of miRNA prediction and signaling pathway enrichment analyses, the down-regulated miRNAs of mmu-miR-484, mmu-miR-1247-5p and mmu-miR-6538 were identified to be the crucial miRNAs. Further analysis revealed that the inflammation and immune responses involving Hippo, PI3K-Akt, MAPK, Wnt, mTOR, TGF-beta, Tight junction, and Platelet activation were modulated collectively by these three key miRNAs during DENV infection. These pathways are considered to be closely associated with the pathogenic mechanism and treatment strategy of dengue patients.

CONCLUSION

The miRNAs identified by sequencing, especially miR-484 may be the potential therapeutic targets for liver involvement in dengue patients which involves the regulation of vascular permeability and expression of inflammatory cytokines.

Transcriptomic discovery of a theranostic signature (SERPINE1/MMP3/COL1A1/SPP1) for head and neck squamous cell carcinomas and identification of antrocinol as a candidate drug

Head and neck squamous cell carcinomas (HNSCC) are prevalent malignancies with a disappointing prognosis, necessitating the search for theranostic biomarkers for better management. Based on a meta-analysis of transcriptomic data containing ten clinical datasets of HNSCC and matched nonmalignant samples, we identified SERPINE1/MMP3/COL1A1/SPP1 as essential hub genes as the potential theranostic biomarkers. Our analysis suggests these hub genes are associated with the extracellular matrix, peptidoglycans, cell migration, wound-healing processes, complement and coagulation cascades, and the AGE-RAGE signaling pathway within the tumor microenvironment. Also, these hub genes were associated with tumor-immune infiltrating cells and immunosuppressive phenotypes of HNSCC. Further investigation of The Cancer Genome Atlas (TCGA) cohorts revealed that these hub genes were associated with staging, metastasis, and poor survival in HNSCC patients. Molecular docking simulations were performed to evaluate binding activities between the hub genes and antrocinol, a novel small-molecule derivative of an anticancer phytochemical antrocin previously discovered by our group. Antrocinol showed high affinities to MMP3 and COL1A1. Notably, antrocinol presented satisfactory drug-like and ADMET properties for therapeutic applications. These results hinted at the potential of antrocinol as an anti-HNSCC candidate via targeting MMP3 and COL1A1. In conclusion, we identified hub genes: SERPINE1/MMP3/COL1A1/SPP1 as potential diagnostic biomarkers and antrocinol as a potential new drug for HNSCC.

CKAP2L Promotes Esophageal Squamous Cell Carcinoma Progression and Drug-Resistance by Modulating Cell Cycle

Esophageal squamous cell carcinoma (ESCC) is one of the most common types of cancer and the leading cause of cancer-related mortality worldwide, especially in Asia. In this study, the gene CKAP2L was selected by GEO, TCGA, and GTEx database analysis. The high expression of CKAP2L is related to the occurrence and development of ESCC. In addition, CKAP2L knockdown can inhibit the growth and migration of ESCC cells, while CKAP2L overexpression has the opposite effect. Furthermore, in vivo experiments indicated that down-regulation of CKAP2L can inhibit the tumorigenesis of ESCC cells. KEGG pathway analysis and the STRING database explored the relationship between cell cycle and CKAP2L and verified that depletion of CKAP2L markedly arrested cell cycle in the G2/M phase. Meanwhile, CKAP2L knockdown increased the sensitivity of ESCC cells to flavopiridol, the first CDK inhibitor to be tested in clinical trials, leading to an observable reduction in cell proliferation and an increase in cellular apoptosis. In brief, we identified CKAP2L as a tumor promoter, potential prognostic indicator, and therapeutic target of ESCC, which may play a role in regulating cell cycle progression.

CDC20 is a novel biomarker for improved clinical predictions in epithelial ovarian cancer

Epithelial ovarian cancer (EOC), a common tumor of the female reproductive system, ranks first in fatalities among gynecological malignancies. Most patients find tumors at late stage and have extremely poor prognoses, which necessitates improvements in early detection. This study applied bioinformatic methods to identify potential biomarkers of EOC. First, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on differentially expressed genes (DEGs) and hub genes, and a protein-protein interaction (PPI) network was constructed. The network of hub genes was analyzed using GeneMANIA, and an analysis of biological processes was constructed with BINGO. Lastly, hub genes were analyzed for EOC-related oncology using the Oncomine and TCGA databases, and the cBioPortal online platform. Overall, cell division cycle 20 (CDC20) was identified as a key gene in EOC. Short hairpin RNA (shRNA) was used to silence CDC20 to explore its effects on EOC cell proliferation, apoptosis and SRY-related HMG-box 2 (SOX2) expression. DEGs were enriched in pathways related to cell cycle signaling, cancer, progesterone-mediated oocyte maturation, Wnt signaling and P53 signaling. Analysis revealed high expression of CDC20 in EOC tissues and a correlation with histology and tumor grade. CDC20 levels are highest in serous adenocarcinoma, when compared to ovarian clear cell carcinoma, ovarian endometrioid carcinoma and mucinous adenocarcinoma. High CDC20 expression within the tumor is associated with poor EOC prognosis. After silencing CDC20, EOC cell proliferation and migration decreased, apoptosis increased, and SOX2 expression decreased. In conclusion, CDC20 is likely a key biomarker of EOC and may act as an upstream regulator of SOX2 to mediate the SOX2 signaling in the progression of EOC. Future application of CDC20 analysis to early detection may improve prognosis, and it has the potential to be a therapeutic target.

Investigation of differentially expressed genes and dysregulated pathways involved in multiple sclerosis

Multiple Sclerosis (MS) is a neurodegenerative autoimmune and organ-specific demyelinating disorder, known to affect the central nervous system (CNS). While genetic studies have revealed several critical genes and diagnostic biomarkers associated with MS, the etiology of the disease remains poorly understood. This study is aimed at screening and identifying the key genes and canonical pathways associated with MS. Gene expression profiling of the microarray dataset GSE38010 was used to analyze two control brain samples (control 1; GSM931812, control 2; GSM931813), active inflammation stage samples (CAP1; GSM931815, CAP2; GSM931816) and late subsided stage samples (CP1; GSM931817, CP2; GSM931818) collected from patients ranging between 23 and 54years and both genders. This analysis yielded a list of 58,866 DEGs (29,433 for active-inflammation stage and 29,433 for late-subsided Stage). The interactions between the DEGs were then studied using STRING, Cytoscape software, and MCODE was employed to find the genes that form clusters. Functional enrichment and integrative analysis were performed using ClueGO/CluePedia and MetaCore™. Our data revealed dysregulated key canonical pathways in MS patients. In addition, we identified three hub genes (SCN2A, HTR2A, and HCN1) that may serve as potential biomarkers for the prognosis of MS. Furthermore, the expression patterns of HPCA and PLCB1 provide insights into the progressive stages of MS, indicating that these genes could be used in predicting MS progression. We were able to map potential biomarkers that could be used for the prognosis and diagnosis of MS.

Integrative ontology and pathway-based approach identifies distinct molecular signatures in transcriptomes of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) remains a serious concern globally due to many factors that including late diagnosis, lack of an ideal biomarker for diagnosis and prognosis, and high rate of mortality. In this study, we aimed to identify the essential dysregulated genes and molecular signatures associated with the progression and development of ESCC. The dataset with 15 ESCCs and the 15 adjacent normal tissue samples from the surrounding histopathologically tumor-free mucosa was selected. We applied bioinformatics pipelines including various topological parameters from MCODE, CytoNCA, and cytoHubba to prioritize the most significantly associated DEGs with ESCC. We performed functional enrichment annotation for the identified DEGs using DAVID and MetaCore™ GeneGo platforms. Furthermore, we validated the essential core genes in TCGA and GTEx datasets between the normal mucosa and ESCC for their expression levels. These DEGs were primarily enriched in positive regulation of transferase activity, negative regulation of organelle organization, cell cycle mitosis/S-phase transition, spindle organization/assembly, development, and regulation of angiogenesis. Subsequently, the DEGs were associated with the pathways such as oocyte meiosis, cell cycle, and DNA replication. Our study identified the eight-core genes (AURKA, AURKB, MCM2, CDC20, TPX2, PLK1, FOXM1, and MCM7) that are highly expressed among the ESCC, and TCGA dataset. The multigene comparison and principal component analysis resulted in elevated signals for the AURKA, MCM2, CDC20, TPX2, PLK1, and FOXM1. Overall, our study reported GO profiles and molecular signatures that might help researchers to grasp the pathological mechanisms underlying ESCC development and eventually provide novel therapeutic and diagnostic strategies.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Whole-exome sequencing analysis of NSCLC reveals the pathogenic missense variants from cancer-associated genes

BACKGROUND

Non-small-cell lung cancer (NSCLC) is the most common type of lung cancer. NSCLC accounts for 84% of all lung cancer cases. In recent years, advances in pathway understanding, methods for discovering novel genetic biomarkers, and new drugs designed to inhibit the signaling cascades have enabled clinicians to personalize therapy for NSCLC.

OBJECTIVES

The primary aim of this study is to identify the genes associated with NSCLC that harbor pathogenic variants that could be causative for NSCLC. The second aim is to investigate their roles in different pathways that lead to NSCLC.

METHODS

We examined exome-sequencing datasets from 54 NSCLC patients to characterize the variants associated with NSCLC.

RESULTS

Our findings revealed that 17 variants in 14 genes were considered highly pathogenic, including CDKN2A, ERBB2, FOXP1, IDH1, JAK3, KMT2D, K-Ras, MSH3, MSH6, POLE, RNF43, TCF7L2, TP53, and TSC1. Gene set enrichment analysis revealed the involvement of transmembrane receptor protein tyrosine kinase activity, protein binding, ATP binding, phosphatidylinositol-4,5-bisphosphate 3-kinase, and Ras guanyl-nucleotide exchange factor activity. Pathway analysis of these genes yielded different cancer-related pathways, including colorectal, prostate, endometrial, pancreatic, PI3K-Akt signaling pathways, and signaling pathways regulating pluripotency of stem cells. Module 1 from protein-protein interactions (PPIs) identified genes that harbor pathogenic SNPs. Three of the most deleterious SNPs are ERBB2 (rs1196929947), K-Ras (rs121913529), and POLE (rs751425952). Interestingly, one patient has a pathogenic K-Ras variant (rs121913529) co-occurred with the missense variant (rs752054698) inTSC1 gene.

CONCLUSION

This study maps highly pathogenic variants associated with NSCLC and investigates their contributions to the pathogenesis of NSCLC. This study sheds light on the potential applications of precision medicine in patients with NSCLC.

GBP5 and ACSS3: two potential biomarkers of high-grade ovarian cancer identified through downstream analysis of microarray data

Among all malignancies of the reproductive organs, ovarian cancer is the sixth leading cause of death for women. Several factors contribute to the uncontrolled expression of certain genes in cancer thus making them act as oncogenes or tumour suppressors. In this study, we have examined four microarray datasets of high-grade ovarian cancer cells to identify differentially expressed genes (DEGs). 362 and 94 common DEGs were identified as up-regulated and down-regulated, respectively from 119 disease and 31 control samples. The DEGs were further analysed for their gene ontologies (GO), pathway, protein-protein interactions and co-expression. Most of the biological processes were associated with cellular processes, biological regulation, metabolic processes, and developmental processes. Further, regulatory networks were constructed by the DEGs which are also co-expressed and the hub genes were identified. The hub genes targeted by a large number of microRNAs (miRNAs) were further analyzed to reveal their role in the overall survival of cancer patients. Finally, GBP5 and ACSS3 were highlighted as potential biomarkers for ovarian cancer research.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation and Essential Dynamics of Deleterious Proline 12 Alanine Single-Nucleotide Polymorphism in PPARγ2 Associated with Type 2 Diabetes, Cardiovascular Disease, and Nonalcoholic Fatty Liver Disease

Background. Peroxisome proliferator-activated receptor-γ (PPARγ) gene is located at 3p25 position. PPARγ functions as the master regulator of glucose homeostasis and lipoprotein metabolism, and recent studies have reported that it is involved in various metabolic diseases such as diabetes mellitus, hyperlipidemia, coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). PPARγ1 and PPARγ2 are necessary for the development of adipose tissue and insulin sensitivity regulation. But PPARγ2 is the isoform that was controlled in response to nutrient intake and obesity. Methodology. In this study, we used computational techniques to show the impact of Pro12Ala polymorphism on PPARγ2. The 3-D structure of PPARγ2 was modeled using I-TASSER server. The modeled structure was validated with the ZLab server, and the mutation was created with SPDB viewer. Stability prediction tools were used. Molecular dynamics simulation (MDS) was used to understand the structural and functional behavior of the wild type and mutant. Essential dynamics was also applied. Results and Conclusions. Stability prediction tools were showed that this mutation has a destabilizing effect on the PPARγ2 structure. The RMSD, RMSF, Rg, SASA, and DSSP were in line with H-bond results that showed less flexibility in the mutant structure. Essential dynamics was used to verify MDS results. Pro12Ala polymorphism leads to rigidity of the PPARγ2 protein and might disturb the conformational changes and interactions of PPARγ2 and results in type 2 diabetes mellitus (T2DM), CAD, and NAFLD. This study can help scientists to develop a drug therapy against these diseases.

Functional network analysis of p85 and PI3K as potential gene targets and mechanism of oleanolic acid in overcoming breast cancer resistance to tamoxifen

Background

Tamoxifen resistance in estrogen receptor positive (ER+) breast cancer therapy increases, which is the leading cause of cancer treatment failure, as it can impair patients’ prognoses, cause cancer recurrence, metastasis, and death. Combination therapy with compounds is needed to overcome tamoxifen resistance. Oleanolic acid (OA) was known to increase tamoxifen sensitivity in tamoxifen-resistant breast cancer; however, the molecular mechanism of OA and its involvement in overcoming tamoxifen resistance remain unknown and need further investigation. This study was conducted to identify the potential gene targets and molecular mechanisms of OA in overcoming tamoxifen resistance.

Results

A bioinformatic approach for functional network analysis was used in silico by utilizing secondary data in the Gene Expression Omnibus (GEO) database and analyzing them with GEO2R to obtain data on differentially expressed genes (DEGs). The DEG data were further examined with Database for Annotation, Visualization, and Integrated Discovery (DAVID), STRING, cBioPortal website, and Cytoscape with its plugin CytoHubba. Molecular docking was performed to predict the binding properties of OA on the protein encoded by the potential gene. CD44, FGFR2, PIK3R1, and MDM2 were designated as potential target genes (PTGs), and PIK3R1 was suspected as the potential gene for OA to overcome tamoxifen resistance. Molecular docking confirms that OA can inhibit p85 activation. PIK3R1 is suggested to be the potential gene for OA in overcoming tamoxifen resistance in breast cancer therapy.

Conclusion

The predicted molecular mechanism of OA in overcoming tamoxifen resistance involves inhibiting p85 activation, leading to the inhibition of the downstream activity of the PI3K signaling pathway, causing breast cancer to respond to tamoxifen therapy once again. Results of this study need to be validated by further studies, including in vitro and in vivo.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00341-4.

Epithelial to mesenchymal transition during mammalian neural crest cell delamination

Epithelial to mesenchymal transition (EMT) is a well-defined cellular process that was discovered in chicken embryos and described as "epithelial to mesenchymal transformation" [1]. During EMT, epithelial cells lose their epithelial features and acquire mesenchymal character with migratory potential. EMT has subsequently been shown to be essential for both developmental and pathological processes including embryo morphogenesis, wound healing, tissue fibrosis and cancer [2]. During the past 5 years, interest and study of EMT especially in cancer biology have increased exponentially due to the implied role of EMT in multiple aspects of malignancy such as cell invasion, survival, stemness, metastasis, therapeutic resistance and tumor heterogeneity [3]. Since the process of EMT in embryogenesis and cancer progression shares similar phenotypic changes, core transcription factors and molecular mechanisms, it has been proposed that the initiation and development of carcinoma could be attributed to abnormal activation of EMT factors usually required for normal embryo development. Therefore, developmental EMT mechanisms, whose timing, location, and tissue origin are strictly regulated, could prove useful for uncovering new insights into the phenotypic changes and corresponding gene regulatory control of EMT under pathological conditions. In this review, we initially provide an overview of the phenotypic and molecular mechanisms involved in EMT and discuss the newly emerging concept of epithelial to mesenchymal plasticity (EMP). Then we focus on our current knowledge of a classic developmental EMT event, neural crest cell (NCC) delamination, highlighting key differences in our understanding of NCC EMT between mammalian and non-mammalian species. Lastly, we highlight available tools and future directions to advance our understanding of mammalian NCC EMT.

Integrated Bioinformatics Analysis for Identification of the Hub Genes Linked with Prognosis of Ovarian Cancer Patients

BACKGROUND

One of the most usual gynecological state of tumor is ovarian cancer and is a major reason of gynecological tumor-related global mortality rate. There have been multiple risk elements related to ovarian cancer like the background of past cases associated with breast cancer or ovarian cancer, or excessive body weight issues, case history of smoking, and untimely menstruation or menopause. Because of unclear expressions, more than 70% of the ovarian cancer patient cases are determined during the early stage. Material and Methods. GSE38666, GSE40595, and GSE66957 were the three microarray datasets which were analyzed using GEO2R for screening the differentially expressed genes. GO, Kyoto Encyclopedia of Genes, and protein expression studies were performed for analysis of hub genes. Then, survival analysis was performed for all the hub genes.

RESULTS

From the dataset, a total of 199 differentially expressed genes (DEGs) were identified. Through the KEGG pathway study, it was noted that the DEGs are mainly linked with the AGE-RAGE signaling pathway, central carbon metabolism, and human papillomavirus infection. The survival analysis showed 4 highly expressed hub genes COL4A1, SDC1, CDKN2A, and TOP2A which correlated with overall survival in ovarian cancer patients. Moreover, the expression of the 4 hub genes was validated by the GEPIA database and the Human Protein Atlas.

CONCLUSION

The results have shown that all 4 hub genes were found to be upregulated in ovarian cancer tissues which predict poor prognosis in patients with ovarian cancer.

Development and Validation of Ten-RNA Binding Protein Signature Predicts Overall Survival in Osteosarcoma

Osteosarcoma is a malignant tumor that originates in the bones with the characteristics of high malignancy, predisposition to metastasis, and poor prognosis. RNA binding proteins (RBPs) are closely related to various tumors, but their relationship with osteosarcoma remains unclear. Based on GTEx and TARGET RNA sequencing data, we applied differential analysis to obtain RBP genes that are differentially expressed in osteosarcoma, and analyzed the functions of these RBPs. After applying univariate and LASSO Cox regression analysis, 10 key prognostic RBPs (TDRD6, TLR8, NXT2, EIF4E3, RPS27L, CPEB3, RBM34, TERT, RPS29, and ZC3HAV1) were screened, and an RBP prognostic risk assessment model for patients with osteosarcoma was established. The independent cohort GSE21257 was used for external verification, and the results showed that the signature has an excellent ability to predict prognosis. In addition, a nomogram that can be used for clinical evaluation was constructed. Finally, the expression levels of 10 prognostic RBPs in osteosarcoma cells and tissues were confirmed through experiments. Our study identified a ten-gene prognostic marker related to RBP, which is of great significance for adjusting the treatment strategy of patients with osteosarcoma and exploring prognostic markers.

Understanding Gene Expression and Transcriptome Profiling of COVID-19: An Initiative Towards the Mapping of Protective Immunity Genes Against SARS-CoV-2 Infection

The COVID-19 pandemic has created an urgent situation throughout the globe. Therefore, it is necessary to identify the differentially expressed genes (DEGs) in COVID-19 patients to understand disease pathogenesis and the genetic factor(s) responsible for inter-individual variability. The DEGs will help understand the disease's potential underlying molecular mechanisms and genetic characteristics, including the regulatory genes associated with immune response elements and protective immunity. This study aimed to determine the DEGs in mild and severe COVID-19 patients versus healthy controls. The Agilent-085982 Arraystar human lncRNA V5 microarray GEO dataset (GSE164805 dataset) was used for this study. We used statistical tools to identify the DEGs. Our 15 human samples dataset was divided into three groups: mild, severe COVID-19 patients and healthy control volunteers. We compared our result with three other published gene expression studies of COVID-19 patients. Along with significant DEGs, we developed an interactome map, a protein-protein interaction (PPI) pattern, a cluster analysis of the PPI network, and pathway enrichment analysis. We also performed the same analyses with the top-ranked genes from the three other COVID-19 gene expression studies. We also identified differentially expressed lncRNA genes and constructed protein-coding DEG-lncRNA co-expression networks. We attempted to identify the regulatory genes related to immune response elements and protective immunity. We prioritized the most significant 29 protein-coding DEGs. Our analyses showed that several DEGs were involved in forming interactome maps, PPI networks, and cluster formation, similar to the results obtained using data from the protein-coding genes from other investigations. Interestingly we found six lncRNAs (TALAM1, DLEU2, and UICLM CASC18, SNHG20, and GNAS) involved in the protein-coding DEG-lncRNA network; which might be served as potential biomarkers for COVID-19 patients. We also identified three regulatory genes from our study and 44 regulatory genes from the other investigations related to immune response elements and protective immunity. We were able to map the regulatory genes associated with immune elements and identify the virogenomic responses involved in protective immunity against SARS-CoV-2 infection during COVID-19 development.

Screening differentially expressed genes between endometriosis and ovarian cancer to find new biomarkers for endometriosis

AIM

Endometriosis is one of the most common reproductive system diseases, but the mechanisms of disease progression are still unclear. Due to its high recurrence rate, searching for potential therapeutic biomarkers involved in the pathogenesis of endometriosis is an urgent issue.

METHODS

Due to the similarities between endometriosis and ovarian cancer, four endometriosis datasets and one ovarian cancer dataset were downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified, followed by gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and protein-protein interaction (PPI) analyses. Then, we validated gene expression and performed survival analysis with ovarian serous cystadenocarcinoma (OV) datasets in TCGA/GTEx database, and searched for potential drugs in the Drug-Gene Interaction Database. Finally, we explored the miRNAs of key genes to find biomarkers associated with the recurrence of endometriosis.

RESULTS

In total, 104 DEGs were identified in the endometriosis datasets, and the main enriched GO functions included cell adhesion, extracellular exosome and actin binding. Fifty DEGs were identified between endometriosis and ovarian cancer datasets including 11 consistently regulated genes, and nine DEGs with significant expression in TCGA/GTEx. Only IGHM had both significant expression and an association with survival, three module DEGs and two significantly expressed DEGs had drug associations, and 10 DEGs had druggability.

CONCLUSIONS

ITGA7, ITGBL1 and SORBS1 may help us understand the invasive nature of endometriosis, and IGHM might be related to recurrence; moreover, these genes all may be potential therapeutic targets.KEY MESSAGEThis manuscript used a bioinformatics approach to find target genes for the treatment of endometriosis.This manuscript used a new approach to find target genes by drawing on common characteristics between ovarian cancer and endometriosis.We screened relevant therapeutic agents for target genes in the drug database, and performed histological validation of target genes with both expression and survival analysis difference in cancer databases.

Identification of key molecular markers of acute coronary syndrome using peripheral blood transcriptome sequencing analysis and mRNA-lncRNA co-expression network construction

Acute coronary syndrome (ACS) is a term used to describe major cardiovascular diseases, and treatment of in-stent restenosis in patients with ACS remains a major clinical challenge. Further investigation into molecular markers of ACS may aid early diagnosis, and the treatment of ACS and post-treatment recurrence. In the present study, total RNA was extracted from the peripheral blood samples of 3 patients with ACS, 3 patients with percutaneous coronary intervention (PCI)_non-restenosis, 3 patients with PCI_restenosis and 3 healthy controls. Subsequently, RNA library construction and high-throughput sequencing were performed. DESeq2 package in R was used to screen genes that were differentially expressed between the different samples. Moreover, the intersection of the differentially expressed mRNAs (DEmRNAs) and differentially expressed long noncoding RNAs (DElncRNAs) obtained. GeneCodis4.0 was used to perform function enrichment for DEmRNAs, and lncRNA-mRNA co-expression network was constructed. The GSE60993 dataset was utilized for diagnostic analysis, and the aforementioned investigations were verified using in vitro studies. Results of the present study revealed a large number of DEmRNAs and DElncRNAs in the different groups. We selected genes in the top 10 of differential expression and also involved in the co-expression of lncRNA-mRNA for diagnostic analysis in the GSE60993 dataset. The area under curve (AUC) of PDZK1IP1 (0.747), PROK2 (0.769) and LAMP3 (0.725) were all >0.7. These results indicated that the identified mRNAs and lncRNAs may act as potential clinical biomarkers, and more specifically, PDZK1IP1, PROK2 and LAMP3 may act as potential biomarkers for the diagnosis of ACS.

Expression of Microtubule-Associated Proteins in Relation to Prognosis and Efficacy of Immunotherapy in Non-Small Cell Lung Cancer

Background

Microtubule-associated proteins (MAPs) have been considered to play significant roles in the tumor evolution of non-small cell lung cancer (NSCLC). Nevertheless, mRNA transcription levels and prognostic value of distinct MAPs in patients with NSCLC remain to be clarified.

Methods

In this study, the Oncomine database, Gene Expression Profiling Interactive Analysis (GEPIA) database, and Human Protein Atlas were utilized to analyze the relationship between mRNA/protein expression of different MAPs and clinical characteristics in NSCLC patients, including tumor type and pathological stage. The correlation between the transcription level of MAPs and overall survival (OS) of NSCLC patients was analyzed by Kaplan–Meier plotter. Besides, 50 frequently altered neighbor genes of the MAPs were screened out, and a network has been constructed via the cBioPortal and Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) dataset. Meanwhile, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis on the expression data of MAPs and their 50 frequently altered neighbor genes in NSCLC tissues. Furthermore, The Cancer Immunome Atlas (TCIA) was utilized to analyze the relationship between MAP expression and the response to immunotherapy. Finally, we used reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to verify the expression of MAPs in 20 patients with NSCLC.

Results

The present study discovered that the mRNA transcription levels of MAP7/7D2 were enriched in NSCLC tissues, while those of the MAP2/4/6/7D3 were lower in NSCLC specimens than those in control specimens. The mRNA transcription level of MAP6 was significantly associated with the advanced stage of NSCLC. Besides, survival analysis indicated that higher mRNA expressions of MAP2/4/6/7/7D3 were correlated considerably with favorable OS of NSCLC patients, whereas increased mRNA expression levels of MAP1A/1S were associated with poor OS. Moreover, the expression of MAP1A/1B/1S/4/6/7D1/7D3 was significantly correlated with immunophenoscore (IPS) in NSCLC patients.

Conclusions

Our analysis indicated that MAP1A/1S could serve as potential personalized therapeutic targets for patients with NSCLC, and the enriched MAP2/4/6/7/7D3 expression could serve as a biomarker for favorable prognosis in NSCLC. Besides, the expression of MAP1A/1B/1S/4/6/7D1/7D3 was closely related to the response to immunotherapy. Taken together, MAP expression has potential application value in the clinical treatment and prognosis assessment of NSCLC patients, and further verifiable experiments can be conducted to verify our results.

Development and Verification of an Immune-Related Gene Pairs Prognostic Signature in Hepatocellular Carcinoma

With the increasing prevalence of Hepatocellular carcinoma (HCC) and the poor prognosis of immunotherapy, reliable immune-related gene pairs (IRGPs) prognostic signature is required for personalized management and treatment of patients. Gene expression profiles and clinical information of HCC patients were obtained from the TCGA and ICGC databases. The IRGPs are constructed using immune-related genes (IRGs) with large variations. The least absolute shrinkage and selection operator (LASSO) regression analysis was used to construct IRGPs signature. The IRGPs signature was verified through the ICGC cohort. 1,309 IRGPs were constructed from 90 IRGs with high variability. We obtained 50 IRGPs that were significantly connected to the prognosis and constructed a signature that included 17 IRGPs. In the TCGA and ICGC cohorts, patients were divided into high and low-risk patients by the IRGPs signature. The overall survival time of low-risk patients is longer than that of high-risk patients. After adjustment for clinical and pathological factors, multivariate analysis showed that the IRGPs signature is an independent prognostic factor. The Receiver operating characteristic (ROC) curve confirmed the accuracy of the signature. Besides, gene set enrichment analysis (GSEA) revealed that the signature is related to immune biological processes, and the immune microenvironment status is distinct in different risk patients. The proposed IRGPs signature can effectively assess the overall survival of HCC, and provide the relationship between the signature and the reactivity of immune checkpoint therapy and the sensitivity of targeted drugs, thereby providing new ideas for the diagnosis and treatment of the disease.

Comparison of critical biomarkers in 2 erectile dysfunction models based on GEO and NOS-cGMP-PDE5 pathway

BACKGROUND

Erectile dysfunction is a disease commonly caused by diabetes mellitus (DMED) and cavernous nerve injury (CNIED). Bioinformatics analyses including differentially expressed genes (DEGs), enriched functions and pathways (EFPs), and protein-protein interaction (PPI) networks were carried out in DMED and CNIED rats in this study. The critical biomarkers that may intervene in nitric oxide synthase (NOS, predominantly nNOS, ancillary eNOS, and iNOS)-cyclic guanosine monophosphate (cGMP)-phosphodiesterase 5 enzyme (PDE5) pathway, an important mechanism in erectile dysfunction treatment, were then explored for potential clinical applications.

METHODS

GSE2457 and GSE31247 were downloaded. Their DEGs with a |logFC (fold change)| > 0 were screened out. Database for Annotation, Visualization and Integrated Discovery (DAVID) online database was used to analyze the EFPs in Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes networks based on down-regulated and up-regulated DEGs respectively. PPI analysis of 2 datasets was performed in Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape. Interactions with an average score greater than 0.9 were chosen as the cutoff for statistical significance.

RESULTS

From a total of 1710 DEGs in GSE2457, 772 were down-regulated and 938 were up-regulated, in contrast to the 836 DEGs in GSE31247, from which 508 were down-regulated and 328 were up-regulated. The 25 common EFPs such as aging and response to hormone were identified in both models. PPI results showed that the first 10 hub genes in DMED were all different from those in CNIED.

CONCLUSIONS

The intervention of iNOS with the hub gene complement component 3 in DMED and the aging process in both DMED and CNIED deserves attention.

Biomarkers of Blood from Patients with Atherosclerosis Based on Bioinformatics Analysis

Atherosclerosis is a multifaceted disease characterized by the formation and accumulation of plaques that attach to arteries and cause cardiovascular disease and vascular embolism. A range of diagnostic techniques, including selective coronary angiography, stress tests, computerized tomography, and nuclear scans, assess cardiovascular disease risk and treatment targets. However, there is currently no simple blood biochemical index or biological target for the diagnosis of atherosclerosis. Therefore, it is of interest to find a biochemical blood marker for atherosclerosis. Three datasets from the Gene Expression Omnibus (GEO) database were analyzed to obtain differentially expressed genes (DEG) and the results were integrated using the Robustrankaggreg algorithm. The genes considered more critical by the Robustrankaggreg algorithm were put into their own data set and the data set system with cell classification information for verification. Twenty-one possible genes were screened out. Interestingly, we found a good correlation between RPS4Y1, EIF1AY, and XIST. In addition, we know the general expression of these genes in different cell types and whole blood cells. In this study, we identified BTNL8 and BLNK as having good clinical significance. These results will contribute to the analysis of the underlying genes involved in the progression of atherosclerosis and provide insights for the discovery of new diagnostic and evaluation methods.

Identification of Vital Hub Genes and Potential Molecular Pathways of Dermatomyositis by Bioinformatics Analysis

Dermatomyositis is an autoimmune disease characterized by severe symmetrical muscle dysfunction and pain. This study was aimed at discovering vital hub genes and potential molecular pathways of DM through bioinformatics analysis, which contributes to identifying potential diagnostic or therapeutic biomarkers and targets. In this study, a total of 915 DEGs in DM samples including 167 upregulated genes and 748 downregulated genes were screened out by the limma package based on the GSE142807 dataset from the Gene Expression Omnibus (GEO) database. Furthermore, the results of Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated that these downregulated genes were highly associated with the immune-related biological processes and pathways. Therefore, 41 genes closely related to DM were extracted for further study based on the subcluster analysis through the Molecular Complex Detection (MCODE) software plugin in Cytoscape. Ultimately, 10 hub genes (including ISG15, DDX58, IFIT3, CXCL10, and STAT1) were identified as the potential candidate biomarkers and targets. Besides, we found that the identified hub genes directly or indirectly communicated with each other via molecular signaling pathways on the protein and transcription level. In general, under the guidance of bioinformatics analysis, 10 vital hub genes and molecular mechanisms in DM were identified and the expression of proinflammatory factors and interferon family proteins and genes showed high association with DM, which might help provide a theoretical foundation for the development of point-to-point targeted therapy in the future treatment of DM.

Epigenetic Biomarkers of Transition from Metabolically Healthy Obesity to Metabolically Unhealthy Obesity Phenotype: A Prospective Study

Background: Identifying those parameters that could potentially predict the deterioration of metabolically healthy phenotype is a matter of debate. In this field, epigenetics, in particular DNA methylation deserves special attention. Results: The aim of the present study was to analyze the long-term evolution of methylation patterns in a subset of metabolically healthy subjects in order to search for epigenetic markers that could predict the progression to an unhealthy state. Twenty-six CpG sites were significantly differentially methylated, both at baseline and 11-year follow-up. These sites were related to 19 genes or pseudogenes; a more in-depth analysis of the methylation sites of these genes showed that CYP2E1 had 50% of the collected CpG sites differently methylated between stable metabolically healthy obesity (MHO) and unstable MHO, followed by HLA-DRB1 (33%), ZBTB45 (16%), HOOK3 (14%), PLCZ1 (14%), SLC1A1 (12%), MUC2 (12%), ZFPM2 (12.5%) and HLA-DQB2 (8%). Pathway analysis of the selected 26 CpG sites showed enrichment in pathways linked to th1 and th2 activation, antigen presentation, allograft rejection signals and metabolic processes. Higher methylation levels in the cg20707527 (ZFPM2) could have a protective effect against the progression to unstable MHO (OR: 0.21, 95%CI (0.067–0.667), p < 0.0001), whilst higher methylation levels in cg11445109 (CYP2E1) would increase the progression to MUO; OR: 2.72, 95%CI (1.094–6.796), p < 0.0014; respectively). Conclusions: DNA methylation status is associated with the stability/worsening of MHO phenotype. Two potential biomarkers of the transition to an unhealthy state were identified and deserve further investigation (cg20707527 and cg11445109). Moreover, the described differences in methylation could alter immune system-related pathways, highlighting these pathways as therapeutic targets to prevent metabolic deterioration in MHO patients.

Bio-informatics and in Vitro Experiments Reveal the Mechanism of Schisandrin A Against MDA-MB-231 cells

Schisandrin A (SchA) has been reported to have good anti-cancer effects. However, its anti-cancer mechanism in breast cancer remains unknown. This study aimed to explore the mechanism of SchA in breast cancer treatment using bio-informatics analysis and in vitro experiments. The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Gene Cards, and PharmMapper databases were used to screen the candidate targets of SchA against MDA-MB-231 cells selected as the tested cell line through MTT analysis. The functions and pathways of the targets were identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and further analyzed using DAVID 6.8.1 database. Network pharmacology analysis revealed 77 candidate targets, 31 signal pathways, and 208 GO entries (P < 0.05). The targets regulated serine-type endopeptidase and protein tyrosine kinase activities, thereby promoting the migration and inhibiting the apoptosis of MDA-MB-231 cells. Comprehensive analysis of the ‘Protein–Protein Interaction’ (PPI) and ‘Component-Targets-Pathways’ (C-T-P) networks constructed using Cytoscape 3.7.1 software revealed four core targets: EGFR, PIK3R1, MMP9 and Caspase 3. Their docking scores with SchA were subsequently investigated through molecular docking. The wound healing, Hoechst 33342/PI, and western blot assays confirmed that SchA significantly down-regulated EGFR, PIK3R1, and MMP9, but up-regulated cleaved-caspase 3, thus inhibiting the migration and promoting the apoptosis of MDA-MB-231 cells. Reckoning the findings of the study, SchA is a potential adjuvant treatment for breast cancer.

Computational Analysis of Gly482Ser Single-Nucleotide Polymorphism in PPARGC1A Gene Associated with CAD, NAFLD, T2DM, Obesity, Hypertension, and Metabolic Diseases

Peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PPARGC1A) regulates the expression of energy metabolism's genes and mitochondrial biogenesis. The essential roles of PPARGC1A encouraged the researchers to assess the relation between metabolism-related diseases and its variants. To study Gly482Ser (+1564G/A) single-nucleotide polymorphism (SNP) after PPARGC1A modeling, we substitute Gly482 for Ser482. Stability prediction tools showed that this substitution decreases the stability of PPARGC1A or has a destabilizing effect on this protein. We then utilized molecular dynamics simulation of both the Gly482Ser variant and wild type of the PPARGC1A protein to analyze the structural changes and to reveal the conformational flexibility of the PPARGC1A protein. We observed loss flexibility in the RMSD plot of the Gly482Ser variant, which was further supported by a decrease in the SASA value in the Gly482Ser variant structure of PPARGC1A and an increase of H-bond with the increase of β-sheet and coil and decrease of turn in the DSSP plot of the Gly482Ser variant. Such alterations may significantly impact the structural conformation of the PPARGC1A protein, and it might also affect its function. It showed that the Gly482Ser variant affects the PPARGC1A structure and makes the backbone less flexible to move. In general, molecular dynamics simulation (MDS) showed more flexibility in the native PPARGC1A structure. Essential dynamics (ED) also revealed that the range of eigenvectors in the conformational space has lower extension of motion in the Gly482Ser variant compared with WT. The Gly482Ser variant also disrupts PPARGC1A interaction. Due to this single-nucleotide polymorphism in PPARGC1A, it became more rigid and might disarray the structural conformation and catalytic function of the protein and might also induce type 2 diabetes mellitus (T2DM), coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). The results obtained from this study will assist wet lab research in expanding potent treatment on T2DM.

Meta-analysis identifying epithelial-derived transcriptomes predicts poor clinical outcome and immune infiltrations in ovarian cancer

BACKGROUND

Previous studies revealed that the epithelial component is associated with the modulation of the ovarian tumor microenvironment (TME). However, the identification of key transcriptional signatures of laser capture microdissected human ovarian cancer epithelia remains lacking.

METHODS

We identified the differentially expressed transcriptional signatures of human ovarian cancer epithelia by meta-analysis of GSE14407, GSE2765, GSE38666, GSE40595, and GSE54388. Then we investigated the enrichment of KEGG pathways that are associated with epithelia-derived transcriptomes. Finally, we investigated the correlation of key epithelia-hub genes with the survival prognosis and immune infiltrations. Finally, we investigated the genetic alterations of key prognostic hub genes and their diagnostic efficacy in ovarian cancer epithelia.

RESULTS

We identified 1339 differentially expressed genes (DEGs) in ovarian cancer epithelia including 541upregulated and 798 downregulated genes. We identified 21 (such as E2F4, FOXM1, TFDP1, E2F1, and SIN3A) and 11 (such as JUN, DDX4, FOSL1, NOC2L, and HMGA1) master transcriptional regulators (MTRs) that are interacted with upregulated and the downregulated genes in ovarian tumor epithelium, respectively. The STRING-based analysis identified hub genes (such as CDK1, CCNB1, AURKA, CDC20, and CCNA2) in ovarian cancer epithelia. The significant clusters of identified hub genes are associated with the enrichment of KEGG pathways including cell cycle, DNA replication, cytokine-cytokine receptor interaction, pathways in cancer, and focal adhesion. The upregulation of SCNN1A and CDCA3 and the downregulation of SOX6 are correlated with a shorter survival prognosis in ovarian cancer (OV). The expression level of SOX6 is negatively correlated with immune score and positively correlated with tumor purity in OV. Moreover, SOX6 is negatively correlated with the infiltration of TILs, CD8+ T cells, CD4+ Regulatory T cells, cytolytic activity, T cell activation, pDC, neutrophils, and macrophages in OV. Also, SOX6 is negatively correlated with various immune markers including CD8A, PRF1, GZMA, GZMB, NKG7, CCL3, and CCL4, indicating the immune regulatory efficiency of SOX6 in the TME of OV. Furthermore, SCNN1A, CDCA3, and SOX6 genes are genetically altered in OV and the expression levels of SCNN1A and SOX6 genes showed diagnostic efficacy in ovarian cancer epithelia.

CONCLUSIONS

The identified ovarian cancer epithelial-derived key transcriptional signatures are significantly correlated with survival prognosis and immune infiltrations, and may provide new insight into the diagnosis and treatment of epithelial ovarian cancer.

Network Pharmacology Combined with Bioinformatics to Investigate the Mechanisms and Molecular Targets of Astragalus Radix-Panax notoginseng Herb Pair on Treating Diabetic Nephropathy

Background

Astragalus Radix (AR)-Panax notoginseng (PN), a classical herb pair, has shown significant effects in treating diabetic nephropathy (DN). However, the intrinsic mechanism of AR-PN treating DN is still unclear. This study aims to illustrate the mechanism and molecular targets of AR-PN treating DN based on network pharmacology combined with bioinformatics.

Materials and Methods

The Traditional Chinese Medicine Systems Pharmacology database was used to screen bioactive ingredients of AR-PN. Subsequently, putative targets of bioactive ingredients were predicted utilizing the DrugBank database and converted into genes on UniProtKB database. DN-related targets were retrieved via analyzing published microarray data (GSE30528) from the Gene Expression Omnibus database. Protein-protein interaction networks of AR-PN putative targets and DN-related targets were established to identify candidate targets using Cytoscape 3.8.0. GO and KEGG enrichment analyses of candidate targets were reflected using a plugin ClueGO of Cytoscape. Molecular docking was performed using AutoDock Vina software, and the results were visualized by Pymol software. The diagnostic capacity of hub genes was verified by receiver operating characteristic (ROC) curves.

Results

Twenty-two bioactive ingredients and 189 putative targets of AR-PN were obtained. Eight hundred and fifty differently expressed genes related to DN were screened. The PPI network showed that 115 candidate targets of AR-PN against DN were identified. GO and KEGG analyses revealed that candidate targets of AR-PN against DN were mainly involved in the apoptosis, oxidative stress, cell cycle, and inflammation response, regulating the PI3K-Akt signaling pathway, cell cycle, and MAPK signaling pathway. Moreover, MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR were considered as the core potential therapeutic targets. Molecular docking demonstrated that these core targets had a great binding affinity with quercetin, kaempferol, isorhamnetin, and formononetin components. ROC curve analysis showed that AKT1, TP53, RELA, JUN, CDKN1A, and EGFR are effective in discriminating DN from controls.

Conclusions

AR-PN against DN may exert its renoprotective effects via various bioactive chemicals and the related pharmacological pathways, involving multiple molecular targets, which may be a promising herb pair treating DN. Nevertheless, these results should be further validated by experimental evidence.

Identification of Multiple Hub Genes and Pathways in Hepatocellular Carcinoma: A Bioinformatics Analysis

Hepatocellular carcinoma (HCC) is a common malignant tumor of the digestive system, and its early asymptomatic characteristic increases the difficulty of diagnosis and treatment. This study is aimed at obtaining some novel biomarkers with diagnostic and prognostic meaning and may find out potential therapeutic targets for HCC. We screen differentially expressed genes (DEGs) from the HCC gene expression profile GSE14520 using GEO2R. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted by using the clusterProfiler software while a protein-protein interaction (PPI) network was performed based on the STRING database. Then, prognosis analysis of hub genes was conducted using The Cancer Genome Atlas (TCGA) database. Quantitative real-time polymerase chain reaction (qRT-PCR) was utilized to further verify the expression of hub genes and explore the correlation between gene expression and clinicopathological parameters. A total of 1053 DEGs were captured, containing 497 upregulated genes and 556 downregulated genes. GO and KEGG analysis indicated that the downregulated DEGs were mainly enriched in the fatty acid catabolic process while upregulated DEGs were primarily enriched in the cell cycle. Simultaneously, ten hub genes (CYP3A4, UGT1A6, AOX1, UGT1A4, UGT2B15, CDK1, CCNB1, MAD2L1, CCNB2, and CDC20) were identified by the PPI network. Five prognosis-related hub genes (CYP3A4, CDK1, CCNB1, MAD2L1, and CDC20) were uncovered by the survival analysis based on TCGA database. The ten hub genes were further validated by qRT-PCR using samples obtained from our hospital. The prognosis-related hub genes such as CYP3A4, CDK1, CCNB1, MAD2L1, and CDC20 could be considered potential diagnosis biomarkers and prognosis targets for HCC. We also use Oncomine for further verification, and we found CCNB1, CCNB2, CDK1, and CYP3A4 which were highly expressed in HCC. Meanwhile, CCNB1, CCNB2, and CDK1 are highly expressed in almost all cancer types, which may play an important role in cancer. Still, further functional study should be conducted to explore the underlying mechanism and biological effect in the near future.

Identification of the correlations between interleukin-27 (IL-27) and immune-inflammatory imbalance in preterm birth

Preterm birth (PTB) is an immune-inflammatory disease that needs to be resolved. This study aimed to identify the role of interleukin-27 (IL-27), an immunomodulatory factor, in PTB and its associated mechanisms. Here, we analyzed the high-throughput of samples data from the maternal-fetal interface to the peripheral circulation obtained from public databases and reported that the elevated IL-27 was involved with the onset of PTB. Further bioinformatics analyses (e.g. GeneMANIA and GSEA) revealed that IL-27 overexpression in the peripheral circulation as well as maternal-fetal interface is related to the activation of the immune-inflammatory process represented by IFN-γ signaling, etc. In addition, IL-27 and immune infiltration correlation analysis demonstrated that IL-27 mediates this immune-inflammatory imbalance, plausibly mainly through monocyte-macrophage and neutrophils. This finding was further validated by analyzing additional datasets. Overall, this is the first study to elaborate on the role of IL-27-mediated immuno-inflammation in PTB from the perspective of bioinformatics, which may provide a novel strategy for the prevention and treatment of PTB.

Multi-omics landscape of circadian rhythm pathway alterations in Glioma

Circadian rhythm pathway was demonstrated pathological functions in glioma on single-gene level. We aim to depict the multi-omics landscape of circadian rhythm pathway alteration in glioma using bioinformatic analyses. Multi-omics data were obtained from “cBioPortal” database. Comparisons were done regarding clinical parameters, differential-expressed genes and functional annotations. A pathway index was generated using the expression data from TCGA and GTEx to quantify the general alteration level of the pathway with clinical association of circadian rhythm pathway index explored. A total of 30 genes were mapped on the circadian rhythm pathway. Genomic profile ofcircadian rhythm pathway genes exhibited distinct characteristics on multiple levels between lower grade glioma (LGG) and glioblastoma multiforme (GBM) patients. LGG patients presented significantly higher frequencies of multi-omics mutations, as well as significant clinical relevance, on single-gene level. Differential-expressed genes between LGG and GBM patients revealed different functions between subtypes that related to the alteration of circadian rhythm pathway. LGG have significantly higher pathway index than normal brain tissue, while GBM significantly lower than normal tissue (P < 0.01), indicating distinctly altered circadian pathway in LGG. Circadian rhythm pathway index correlated with the prognosis of LGG, but not GBM, patients, with higher score indicating better survival outcome (LGG: HR = 0.39, 95% CI: 0.26 − 0.59, P < 0.001). In conclusion, LGG have more multi-omics alterations of circadian rhythm pathway than GBM. Quantification of circadian rhythm pathway using pathway index demonstrated hyperactivated pathway status in LGG and correlated with the prognosis of LGG patients.

A Comprehensive Analysis of the Downregulation of miRNA-1827 and Its Prognostic Significance by Targeting SPTBN2 and BCL2L1 in Ovarian Cancer

Background: Previous studies demonstrated that miRNA-1827 could repress various cancers on proliferation, angiogenesis, and metastasis. However, little attention has been paid to its role in ovarian cancer as a novel biomarker or intervention target, especially its clinical significance and underlying regulatory network.

Methods: A meta-analysis of six microarrays was adopted here to determine the expression trend of miRNA-1827, and was further validated by gene expression profile data and cellular experiments. We explored the functional annotations through enrichment analysis for the differentially expressed genes targeted by miRNA-1827. Subsequently, we identified two hub genes, SPTBN2 and BCL2L1, based on interaction analysis using two online archive tools, miRWALK (it consolidates the resources of 12 miRNA-focused servers) and Gene Expression Profiling Interactive Analysis (GEPIA). Finally, we validated their characteristics and clinical significance in ovarian cancer.

Results: The comprehensive meta-analysis revealed that miRNA-1827 was markedly downregulated in clinical and cellular specimens. Transfection of the miRNA-1827 mimic could significantly inhibit cellular proliferation. Concerning its target genes, they were involved in diverse biological processes related to tumorigenesis, such as cell proliferation, migration, and the apoptosis signaling pathway. Moreover, interaction analysis proved that two hub genes, SPTBN2 and BCL2L1, were highly associated with poor prognosis in ovarian cancer.

Conclusion: These integrated bioinformatic analyses indicated that miRNA-1827 was dramatically downregulated in ovarian cancer as a tumor suppressor. The upregulation of its downstream modulators, SPTBN2 and BCL2L1, was associated with an unfavorable prognosis. Thus, the present study has identified miRNA-1827 as a potential intervention target for ovarian cancer based on our bioinformatic analysis processes.

Genes That Predict Poor Prognosis in Breast Cancer via Bioinformatical Analysis

Background

Breast cancer is one of the most commonly diagnosed cancers all over the world, and it is now the leading cause of cancer death among females. The aim of this study was to find DEGs (differentially expressed genes) which can predict poor prognosis in breast cancer and be effective targets for breast cancer patients via bioinformatical analysis.

Methods

GSE86374, GSE5364, and GSE70947 were chosen from the GEO database. DEGs between breast cancer tissues and normal breast tissues were picked out by GEO2R and Venn diagram software. Then, DAVID (Database for Annotation, Visualization, and Integrated Discovery) was used to analyze these DEGs in gene ontology (GO) including molecular function (MF), cellular component (CC), and biological process (BP) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway. Next, STRING (Search Tool for the Retrieval of Interacting Genes) was used to investigate potential protein-protein interaction (PPI) relationships among DEGs and these DEGs were analyzed by Molecular Complex Detection (MCODE) in Cytoscape. After that, UALCAN, GEPIA (gene expression profiling interactive analysis), and KM (Kaplan-Meier plotter) were used for the prognostic information and core genes were qualified.

Results

There were 96 upregulated genes and 98 downregulated genes in this study. 55 upregulated genes were selected as hub genes in the PPI network. For validation in UALCAN, GEPIA, and KM, 5 core genes (KIF4A, RACGAP1, CKS2, SHCBP1, and HMMR) were found to highly expressed in breast cancer tissues with poor prognosis. They differentially expressed between different subclasses of breast cancer.

Conclusion

These five genes (KIF4A, RACGAP1, CKS2, SHCBP1, and HMMR) could be potential targets for therapy in breast cancer and prediction of prognosis on the basis of bioinformatical analysis.

CXCL1 Clone Evolution Induced by the HDAC Inhibitor Belinostat Might Be a Favorable Prognostic Indicator in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is the most lethal subtype of breast cancer due to its lack of treatment options. Patients with TNBC frequently develop resistance to chemotherapy. As epigenetic-based antineoplastic drugs, histone deacetylase inhibitors (HDACis) have achieved particular efficacy in lymphoma but are less efficacious in solid tumors, and the resistance mechanism remains poorly understood. In this study, the GSE129944 microarray dataset from the Gene Expression Omnibus database was downloaded, and fold changes at the transcriptome level of a TNBC line (MDA-MB-231) after treatment with belinostat were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to identify the critical biological processes. Construction and analysis of the protein-protein interaction (PPI) network were performed to screen candidate genes related to cancer prognosis. A total of 465 DEGs were identified, including 240 downregulated and 225 upregulated genes. The cytokine-cytokine receptor pathway was identified as being significantly changed. Furthermore, the expression of CXCL1 was implicated as a favorable factor in the overall survival of breast cancer patients. With in vitro approaches, we also showed that belinostat could induce the expression of CXCL1 in another 2 TNBC cell lines (BT-549 and HCC-1937). We speculate that belinostat-induced CXCL1 expression could be one of the results of the stress clone evolution of cells after HDACi treatment. These findings provide new insights into clone evolution during HDACi treatment, which might guide us to a novel perspective that various mutation-targeted treatments should be implemented during the whole treatment cycle.

An integrative analysis to distinguish between emphysema (EML) and alpha-1 antitrypsin deficiency-related emphysema (ADL)-A systems biology approach

Lung Emphysema is an abnormal enlargement of the air sacs followed by the destruction of alveolar walls without any prominent fibrosis. This study primarily identifies the differentially expressed genes (DEGs), interactions between them, and their significant involvement in the activated signaling cascades. The dataset with ID GSE1122 (five normal lung tissue samples, five of usual emphysema, and five of alpha-1 antitrypsin deficiency-related emphysema) from the gene expression omnibus (GEO) was analyzed using the GEO2R tool. The physical association between the DEGs were mapped using the STRING tool and was visualized in the Cytoscape software. The enriched functional processes were identified with the ClueGO plugin's help from Cytoscape. Further integrative functional annotation was performed by implying the GeneGo Metacore™ to distinguish the enriched pathway maps, process networks, and GO processes. The results from this analysis revealed the critical signaling cascades that have been either activated or inhibited due to identified DEGs. We found the activated pathways such as immune response IL-1 signaling pathway, positive regulation of smooth muscle migration, BMP signaling pathway, positive regulation of leukocyte migration, NIK/NF-kappB signaling, and cytochrome-c oxidase activity. Finally, we mapped four crucial genes (CCL5, ALK, TAC1, CD74, and HLA-DOA) by comparing the functional annotations that could be significantly influential in emphysema molecular pathogenesis. Our study provides insights into the pathogenesis of emphysema and helps in developing potential drug targets against emphysema.

Cancer stem cells targets and combined therapies to prevent cancer recurrence

Cancer stem cells (CSCs) control the dynamics of tumorigenesis by self-renewal ability and differentiation potential. These properties contribute towards tumor malignancy, metastasis, cellular heterogeneity, and immune escape, which are regulated by multiple signaling pathways. The CSCs are chemoresistant and cause cancer recurrence, generally recognized as a small side-population that eventually leads to tumor relapse. Despite many treatment options available, none can be considered entirely efficient due to a lack of specificity and dose limitation. This review primarily highlights the processes involved in CSCs development and maintenance. Secondly, the current effective therapies based on stem cells, cell-free therapies that involve exosomes and miRNAs, and photodynamic therapy have been discussed. Also, the inhibitors that specifically target various signaling pathways, which can be used in combination to control CSCs kinetics have been highlighted. Conclusively, this comprehensive review is a detailed study of recently developed novel treatment strategies that will facilitate in coming up with better-targeted approaches against CSCs.

A systemic approach to explore the mechanisms of drug resistance and altered signaling cascades in extensively drug-resistant tuberculosis

BACKGROUND AND AIM

The persistence of extensively drug-resistant (XDR) strains of Mycobacterium tuberculosis (MTB) continue to pose a significant challenge to the treatment and control of tuberculosis infections worldwide. XDR-MTB strains exhibit resistance against first-line anti-TB drugs, fluoroquinolones, and second-line injectable drugs. The mechanisms of drug resistance of MTB remains poorly understood. Our study aims at identifying the differentially expressed genes (DEGs), associated gene networks, and signaling cascades involved in rendering this pathogen resistant to multiple drugs, namely, isoniazid, rifampicin, and capreomycin.

METHODS

We used the microarray dataset GSE53843. The GEO2R tool was used to prioritize the most significant DEGs (top 250) of each drug exposure sample between XDR strains and non-resistant strains. The validation of the 250 DEGs was performed using volcano plots. Protein-protein interaction networks of the DEGs were created using STRING and Cytoscape tools, which helped decipher the relationship between these genes. The significant DEGs were functionally annotated using DAVID and ClueGO. The concomitant biological processes (BP) and molecular functions (MF) were represented as dot plots.

RESULTS AND CONCLUSION

We identified relevant molecular pathways and biological processes, such as cell wall biogenesis, lipid metabolic process, ion transport, phosphopantetheine binding, and triglyceride lipase activity. These processes indicated the involvement of multiple interconnected mechanisms in drug resistance. Our study highlighted the impact of cell wall permeability, with the dysregulation of the mur family of proteins, as essential factors in the inference of resistance. Additionally, upregulation of genes responsible for ion transport such as ctpF, arsC, and nark3, emphasizes the importance of transport channels and efflux pumps in potentially driving out stress-inducing compounds. This study investigated the upregulation of the Lip family of proteins, which play a crucial role in triglyceride lipase activity. Thereby illuminating the potential role of drug-induced dormancy and subsequent resistance in the mycobacterial strains. Multiple mechanisms such as carboxylic acid metabolic process, NAD biosynthetic process, triglyceride lipase activity, phosphopantetheine binding, organic acid biosynthetic process, and growth of symbiont in host cell were observed to partake in resistance of XDR-MTB. This study ultimately provides a platform for important mapping targets for potential therapeutics against XDR-MTB.

Construction and Analysis of Survival-Associated Competing Endogenous RNA Network in Lung Adenocarcinoma

Increasing evidence has shown that noncoding RNAs play significant roles in the initiation, progression, and metastasis of tumours via participating in competing endogenous RNA (ceRNA) networks. However, the survival-associated ceRNA in lung adenocarcinoma (LUAD) remains poorly understood. In this study, we aimed to investigate the regulatory mechanisms underlying ceRNA in LUAD to identify novel prognostic factors. mRNA, lncRNA, and miRNA sequencing data obtained from the GDC data portal were utilized to identify differentially expressed (DE) RNAs. Survival-related RNAs were recognized using univariate Kaplan-Meier survival analysis. We performed functional enrichment analysis of survival-related mRNAs using the clusterProfiler package of R and STRING. lncRNA-miRNA and miRNA-mRNA interactions were predicted based on miRcode, Starbase, and miRanda. Subsequently, the survival-associated ceRNA network was constructed for LUAD. Multivariate Cox regression analysis was used to identify prognostic factors. Finally, we acquired 15 DE miRNAs, 49 DE lncRNAs, and 843 DE mRNAs associated with significant overall survival. Functional enrichment analysis indicated that survival-related DE mRNAs were enriched in cell cycle. The survival-associated lncRNA-miRNA-mRNA ceRNA network was constructed using five miRNAs, 49 mRNAs, and 21 lncRNAs. Furthermore, seven hub RNAs (LINC01936, miR-20a-5p, miR-31-5p, TNS1, TGFBR2, SMAD7, and NEDD4L) were identified based on the ceRNA network. LINC01936 and miR-31-5p were found to be significant using the multifactorial Cox regression model. In conclusion, we successfully constructed a survival-related lncRNA-miRNA-mRNA ceRNA regulatory network in LUAD and identified seven hub RNAs, which provide novel insights into the regulatory molecular mechanisms associated with survival of LUAD, and identified two independent prognostic predictors for LUAD.

Investigation of the Mechanism of Complement System in Diabetic Nephropathy via Bioinformatics Analysis

OBJECTIVES

Diabetic nephropathy (DN) is a major cause of end-stage renal disease (ESRD) throughout the world, and the identification of novel biomarkers via bioinformatics analysis could provide research foundation for future experimental verification and large-group cohort in DN models and patients.

METHODS

GSE30528, GSE47183, and GSE104948 were downloaded from Gene Expression Omnibus (GEO) database to find differentially expressed genes (DEGs). The difference of gene expression between normal renal tissues and DN renal tissues was firstly screened by GEO2R. Then, the protein-protein interactions (PPIs) of DEGs were performed by STRING database, the result was integrated and visualized via applying Cytoscape software, and the hub genes in this PPI network were selected by MCODE and topological analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out to determine the molecular mechanisms of DEGs involved in the progression of DN. Finally, the Nephroseq v5 online platform was used to explore the correlation between hub genes and clinical features of DN.

RESULTS

There were 64 DEGs, and 32 hub genes were identified, enriched pathways of hub genes involved in several functions and expression pathways, such as complement binding, extracellular matrix structural constituent, complement cascade related pathways, and ECM proteoglycans. The correlation analysis and subgroup analysis of 7 complement cascade-related hub genes and the clinical characteristics of DN showed that C1QA, C1QB, C3, CFB, ITGB2, VSIG4, and CLU may participate in the development of DN.

CONCLUSIONS

We confirmed that the complement cascade-related hub genes may be the novel biomarkers for DN early diagnosis and targeted treatment.

Network analysis of transcriptomics data for the prediction and prioritization of membrane-associated biomarkers for idiopathic pulmonary fibrosis (IPF) by bioinformatics approach

Idiopathic pulmonary fibrosis (IPF) is a rare yet crucial persistent lung disorder that actuates scarring of lung tissues, which makes breathing difficult. Smoking, environmental pollution, and certain viral infections could initiate lung scarring. However, the molecular mechanism involved in IPF remains elusive. To develop an efficient therapeutic arsenal against IPF, it is vital to understand the pathology and deviations in biochemical pathways that lead to disorder. In this study, we availed network analysis and other computational pipelines to delineate the prominent membrane proteins as diagnostic biomarkers and therapeutic targets for IPF. This study yielded a significant role of glycosaminoglycan binding, endothelin, and GABA-B receptor signaling pathway in IPF pathogenesis. Furthermore, ADCY8, CRH, FGB, GPR17, MCHR1, NMUR1, and SAA1 genes were found to be immensely involved with IPF, and the enrichment pathway analysis suggests that most of the pathways were corresponding to membrane transport and signal transduction functionalities. This analysis could help in better understanding the molecular mechanism behind IPF to develop an efficient therapeutic target or biomarkers for IPF.