Gene network in pulmonary tuberculosis based on bioinformatic analysis

Network pharmacological prediction of the mechanism of action of Shen-Zhu-Lian-Bai Decoction in the treatment of ulcerative colitis

The incidence of ulcerative colitis (UC) is on the rise globally. Shen-Zhu-Lian-Bai decoction (SZLBD) can relieve the clinical symptoms of UC. This study aimed to investigate the underlying molecular mechanism of SZLBD in the treatment of UC. The key treatment targets of SZLBD for UC were obtained based on the online database, and combined with the STRING database and Cytoscape 3.7.2 software, PPI network was constructed and visualized. The GEO database was utilized to validate the expression levels of core targets in UC. Metascape database GO functional annotation and KEGG pathway enrichment analysis. Molecular docking technology was used to verify the docking of core compounds with key targets. RT-qPCR and Western Blot were used to detect the expression of key targets in HCoEpiC cells for verification. After screening, 67 targets shared by SZLBD and UC were obtained. It is predicted that IL-6, IL-1B, and AKT1 might be the key targets of SZLBD in the treatment of UC. Quercetin was the main active ingredient. GEO results showed that the expression levels of IL-6, IL-1B and AKT1 were higher in the UC group compared to the control group. GO and KEGG analyses showed that these targets were related to apoptosis and inflammation. The results of molecular docking demonstrated that the AKT1 gene, a key target of quercetin, had the highest affinity of -9.2 kcal/mol. Cell experiments found that quercetin could affect the expression of IL-6, IL-1B, and AKT1. This study preliminarily explored and verified the mechanism of action of SZLBD in the treatment of UC, which provides a theoretical basis for subsequent in vivo mechanism studies.

Identification of Potential Differentially-Methylated/Expressed Genes in Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) is a chronic inflammatory lung disease that causes obstructed airflow from the lungs. DNA methylation can regulate gene expression. Understanding the potential molecular mechanism of COPD is of great importance. The aim of this study was to find differentially methylated/expressed genes in COPD. DNA methylation and gene expression profiles in COPD were downloaded from the dataset, followed by functional analysis of differentially-methylated/expressed genes. The potential diagnostic value of these differentially-methylated/expressed genes was determined by receiver operating characteristic (ROC) analysis. Expression validation of differentially-methylated/expressed genes was performed by in vitro experiment and extra online datasets. Totally, 81 hypermethylated-low expression genes and 121 hypomethylated-high expression genes were found in COPD. Among which, 9 core hypermethylated-low expression genes (CD247, CCR7, CD5, IKZF1, SLAMF1, IL2RB, CD3E, CD7 and IL7R) and 8 core hypomethylated-high expression genes (TREM1, AQP9, CD300LF, CLEC12A, NOD2, IRAK3, NLRP3 and LYZ) were identified in the protein-protein interaction (PPI) network. Moreover, these genes had a potential diagnostic utility for COPD. Some signaling pathways were identified in COPD, including T cell receptor signaling pathway, cytokine-cytokine receptor interaction, hematopoietic cell lineage, HTLV-I infection, endocytosis and Jak-STAT signaling pathway. In conclusion, differentially-methylated/expressed genes and involved signaling pathways are likely to be associated with the process of COPD.

Integrated bioinformatic analyses investigate macrophage-M1-related biomarkers and tuberculosis therapeutic drugs

Tuberculosis (TB) is a common infectious disease linked to host genetics and the innate immune response. It is vital to investigate new molecular mechanisms and efficient biomarkers for Tuberculosis because the pathophysiology of the disease is still unclear, and there aren't any precise diagnostic tools. This study downloaded three blood datasets from the GEO database, two of which (GSE19435 and 83456) were used to build a weighted gene co-expression network for searching hub genes associated with macrophage M1 by the CIBERSORT and WGCNA algorithms. Furthermore, 994 differentially expressed genes (DEGs) were extracted from healthy and TB samples, four of which were associated with macrophage M1, naming RTP4, CXCL10, CD38, and IFI44. They were confirmed as upregulation in TB samples by external dataset validation (GSE34608) and quantitative real-time PCR analysis (qRT-PCR). CMap was used to predict potential therapeutic compounds for tuberculosis using 300 differentially expressed genes (150 downregulated and 150 upregulated genes), and six small molecules (RWJ-21757, phenamil, benzanthrone, TG-101348, metyrapone, and WT-161) with a higher confidence value were extracted. We used in-depth bioinformatics analysis to investigate significant macrophage M1-related genes and promising anti-Tuberculosis therapeutic compounds. However, more clinical trials were necessary to determine their effect on Tuberculosis.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

Objective

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

Methods

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

Results

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

Conclusion

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Prognostic value of SOX9 in cervical cancer: Bioinformatics and experimental approaches

Among gynecological cancers, cervical cancer is a common malignancy and remains the leading cause of cancer-related death for women. However, the exact molecular pathogenesis of cervical cancer is not known. Hence, understanding the molecular mechanisms underlying cervical cancer pathogenesis will aid in the development of effective treatment modalities. In this research, we attempted to discern candidate biomarkers for cervical cancer by using multiple bioinformatics approaches. First, we performed differential expression analysis based on cervical squamous cell carcinoma and endocervical adenocarcinoma data from The Cancer Genome Atlas database, then used differentially expressed genes for weighted gene co-expression network construction to find the most relevant gene module for cervical cancer. Next, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed on the module genes, followed by using protein–protein interaction network analysis and Cytoscape to find the key gene. Finally, we validated the key gene by using multiple online sites and experimental methods. Through weighted gene co-expression network analysis, we found the turquoise module was the highest correlated module with cervical cancer diagnosis. The biological process of the module genes focused on cell proliferation, cell adhesion, and protein binding processes, while the Kyoto Encyclopedia of Genes and Genomes pathway of the module significantly enriched pathways related to cancer and cell circle. Among the module genes, SOX9 was identified as the hub gene, and its expression was associated with cervical cancer prognosis. We found the expression of SOX9 correlates with cancer-associated fibroblast immune infiltration in immune cells by Timer2.0. Furthermore, cancer-associated fibroblast infiltration is linked to cervical cancer patients’ prognosis. Compared to those in normal adjacent, immunohistochemical and real-time quantitative polymerase chain reaction (qPCR) showed that the protein and mRNA expression of SOX9 in cervical cancer were higher. Therefore, the SOX9 gene acts as an oncogene in cervical cancer, interactive with immune infiltration of cancer-associated fibroblasts, thereby affecting the prognosis of patients with cervical cancer.

An Epistatic Network Describes oppA and glgB as Relevant Genes for Mycobacterium tuberculosis

Mycobacterium tuberculosis is an acid-fast bacterium that causes tuberculosis worldwide. The role of epistatic interactions among different loci of the M. tuberculosis genome under selective pressure may be crucial for understanding the disease and the molecular basis of antibiotic resistance acquisition. Here, we analyzed polymorphic loci interactions by applying a model-free method for epistasis detection, SpydrPick, on a pan–genome-wide alignment created from a set of 254 complete reference genomes. By means of the analysis of an epistatic network created with the detected epistatic interactions, we found that glgB (α-1,4-glucan branching enzyme) and oppA (oligopeptide-binding protein) are putative targets of co-selection in M. tuberculosis as they were associated in the network with M. tuberculosis genes related to virulence, pathogenesis, transport system modulators of the immune response, and antibiotic resistance. In addition, our work unveiled potential pharmacological applications for genotypic antibiotic resistance inherent to the mutations of glgB and oppA as they epistatically interact with fprA and embC, two genes recently included as antibiotic-resistant genes in the catalog of the World Health Organization. Our findings showed that this approach allows the identification of relevant epistatic interactions that may lead to a better understanding of M. tuberculosis by deciphering the complex interactions of molecules involved in its metabolism, virulence, and pathogenesis and that may be applied to different bacterial populations.

An integrated analysis of prognostic and immune infiltrates for hub genes as potential survival indicators in patients with lung adenocarcinoma

Objective

Lung adenocarcinoma (LUAD) is one of the major subtypes of lung cancer that is associated with poor prognosis. The aim of this study was to identify useful biomarkers to enhance the treatment and diagnosis of LUAD.

Methods

GEO2R was used to identify common up-regulated differentially expressed genes (DEGs) in the GSE32863, GSE40791, and GSE75037 datasets. The DEGs were submitted to Metascape for gene ontology and pathway enrichment analysis as well as construction of the protein-protein interaction (PPI) network, while the molecular complex detection (MCODE) plug-in was employed to filter important subnetworks. The expression levels of the hub genes and their prognostic values were evaluated using the UALCAN, GEPIA2, and Kaplan-Meier plotter databases. The timer algorithm was utilized to determine the correlation between immune cell infiltration and the expression levels of hub genes in LUAD tissues. In addition, the hub gene mutation landscape and the correlation analysis with tumor mutational burden (TMB) score were evaluated using maftools package and ggstatsplot package in R software, respectively.

Results

We identified 156 common up-regulated DEGs, with gene ontology and pathway enrichment analysis indicating that they were mostly enriched in mitotic cell cycle process and cell cycle pathway. DEGs in the subnetwork with the largest number of genes were AURKB, CCNB2, CDC20, CDCA5, CDCA8, CENPF, and KNTC1. The seven hub genes were highly expressed in LUAD tissues and were associated with poor prognosis. These hub genes were negatively correlated with most immune cells. The somatic mutation landscape showed that AURKB, CDC20, CENPF, and KNTC1 had mutations and were positively correlated with TMB scores.

Conclusions

Our findings demonstrate that increased expression of seven hub genes is associated with poor prognosis for LUAD patients. Additionally, the TMB score indicates that the high expression of hub gene increases immune cell infiltration in patients with lung adenocarcinoma which may significantly improve response to immunotherapy.

ORN: Inferring patient-specific dysregulation status of pathway modules in cancer with OR-gate Network

Pathway level understanding of cancer plays a key role in precision oncology. However, the current amount of high-throughput data cannot support the elucidation of full pathway topology. In this study, instead of directly learning the pathway network, we adapted the probabilistic OR gate to model the modular structure of pathways and regulon. The resulting model, OR-gate Network (ORN), can simultaneously infer pathway modules of somatic alterations, patient-specific pathway dysregulation status, and downstream regulon. In a trained ORN, the differentially expressed genes (DEGs) in each tumour can be explained by somatic mutations perturbing a pathway module. Furthermore, the ORN handles one of the most important properties of pathway perturbation in tumours, the mutual exclusivity. We have applied the ORN to lower-grade glioma (LGG) samples and liver hepatocellular carcinoma (LIHC) samples in TCGA and breast cancer samples from METABRIC. Both datasets have shown abnormal pathway activities related to immune response and cell cycles. In LGG samples, ORN identified pathway modules closely related to glioma development and revealed two pathways closely related to patient survival. We had similar results with LIHC samples. Additional results from the METABRIC datasets showed that ORN could characterize critical mechanisms of cancer and connect them to less studied somatic mutations (e.g., BAP1, MIR604, MICAL3, and telomere activities), which may generate novel hypothesis for targeted therapy.