Identification potential biomarkers in pulmonary tuberculosis and latent infection based on bioinformatics analysis

Computational reassessment of RNA-seq data reveals key genes in active tuberculosis

BACKGROUND

Tuberculosis is a serious life-threatening disease among the top global health challenges and rapid and effective diagnostic biomarkers are vital for early diagnosis especially given the increasing prevalence of multidrug resistance.

METHODS

Two human whole blood microarray datasets, GSE42826 and GSE42830 were retrieved from publicly available gene expression omnibus (GEO) database. Deregulated genes (DEGs) were identified using GEO2R online tool and Gene Ontology (GO), protein-protein interaction (PPI) network analysis was performed using Metascape and STRING databases. Significant genes (n = 8) were identified using T-test/ANOVA and Molecular Complex Detection (MCODE) score ≥10, which was validated in GSE34608 dataset. The diagnostic potential of three biomarkers was assessed using Area Under Curve (AUC) of Receiver Operating Characteristic (ROC) plot. The transcriptional levels of these genes were also examined in a separate dataset GSE31348, to monitor the patterns of variation during tuberculosis treatment.

RESULTS

A total of 62 common DEGs (57 upregulated, 7 downregulated genes) were identified in two discovery datasets. GO functions and pathway enrichment analysis shed light on the functional roles of these DEGs in immune response and type-II interferon signaling. The genes in Module-1 (n = 18) were linked to innate immune response, interferon-gamma signaling. The common genes (n = 8) were validated in GSE34608 dataset, that corroborates the results obtained from discovery sets. The gene expression levels demonstrated responsiveness to Mtb infection during anti-TB therapy in GSE31348 dataset. In GSE34608 dataset, the expression levels of three specific genes, GBP5, IFITM3, and EPSTI1, emerged as potential diagnostic makers. In combination, these genes scored remarkable diagnostic performance with 100% sensitivity and 89% specificity, resulting in an impressive Area Under Curve (AUC) of 0.958. However, GBP5 alone showed the highest AUC of 0.986 with 100% sensitivity and 89% specificity.

CONCLUSIONS

The study presents valuable insights into the critical gene network perturbed during tuberculosis. These genes are determinants for assessing the effectiveness of an anti-TB response and distinguishing between active TB and healthy individuals. GBP5, IFITM3 and EPSTI1 emerged as candidate core genes in TB and holds potential as novel molecular targets for the development of interventions in the treatment of TB.

MICAL1 facilitates pancreatic cancer proliferation, migration, and invasion by activating WNT/β-catenin pathway

Background

MICAL1 is involved in the malignant processes of several types of cancer; however, the role of MICAL1 in pancreatic cancer (PC) has not been well-characterized. This study aimed to investigate the expression and function of MICAL1 in PC.

Methods

RT-qPCR and immunohistochemistry were used to detect MICAL1 expression in PC and adjacent nontumor tissues. Cell Counting Kit-8, EdU, clone formation, wound healing, and Transwell assays as well as animal models were used to investigate the effects of overexpression or inhibition of MICAL1 expression on the proliferation, invasion, and metastasis of PC cells. RNA-seq was used to explore the main pathway underlying the functions of MICAL1. Proteomics, mass spectrometry, and co-immunoprecipitation assays were used to investigate the interaction of proteins with MICAL1. Rescue experiments were conducted to validate these findings.

Results

Both MICAL1 mRNA and protein levels were upregulated in PC tissues compared with matched adjacent nontumor tissues. The expression level of MICAL1 was associated with the proliferative and metastatic status of PC. Repression of MICAL1 significantly inhibited PC cell growth, migration, and invasion in vitro and in vivo. RNA sequencing analysis indicated that MICAL1 was closely correlated with the WNT pathway. Overexpression of MICAL1 (1) promoted the phosphorylation of TBC1D1 at the Ser660 site, (2) facilitated the distribution of FZD7 on the cytomembrane, (3) inhibited the degradation of FZD7 in the lysosome, and (4) activated the WNT pathway.

Conclusions

MICAL1 was upregulated in PC and involved in stimulating the progression of PC cells by activating the WNT/β-catenin signaling pathway. Therefore, MICAL1 is a potential therapeutic target for PC.

Uncovering the mechanism of Huanglian-Wuzhuyu herb pair in treating nonalcoholic steatohepatitis based on network pharmacology and experimental validation

ETHNOPHARMACOLOGICAL RELEVANCE

The Huanglian-Wuzhuyu herb pair (HWHP) is a classic Chinese herbal formula consisting of the root of Coptis chinensis Franch and dried, nearly mature scented fruit of Tetradium ruticarpum (A.Juss.) T.G.Hartley. It is widely utilized to treat gastrointestinal and liver diseases such as diarrhea, dysentery, cholestasis, hepatocellular carcinoma, and nonalcoholic steatohepatitis (NASH). However, the mechanism of HWHP in treating NASH remains poorly understood.

AIM OF THE STUDY

This study investigated the mechanisms of HWHP in NASH treatment via network pharmacology and validated the results through in vivo experiment using mouse models.

MATERIALS AND METHODS

The compounds and targets corresponding to the active ingredients of HWHP were obtained from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) database. The genes associated with NASH were obtained from the DisGeNET database. Cytoscape software was employed to construct a "drug-ingredient-target-disease" network. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were employed to analyze the related signaling pathways affected by HWHP. Moreover, AutoDock software was used to assess the potential binding affinity between the key targets of the hub pathway and the bioactive compounds. Subsequently, in vivo experiment was conducted to verify the findings of network pharmacology.

RESULTS

A total of 41 active compounds and 198 targets of HWHP were screened, of which 51 common targets were related to NASH. GO functional enrichment analysis revealed that HWHP may affect NASH by modulating inflammatory response. KEGG pathway enrichment suggested that the NOD-like receptor (NLR) signaling pathway may play an important role in treating NASH. Molecular docking results demonstrated that most HWHP components were successfully docked to NLRP3 with good binding energy. In vivo experiments revealed that HWHP alleviated liver inflammation, improved liver steatosis, reduced TC, TG, LDL-C, ALT, and AST, decreased mRNA expressions of IL-6, IL-18, and TNF-α in the liver, and lowered the expressions of NLRP3, pro-IL-1β, and ASC protein. Also, immunohistochemical findings presented downregulation of caspase-1 and IL-1β by HWHP.

CONCLUSIONS

The results disclosed that HWHP ameliorates NASH in mice by reducing inflammation and liver steatosis via inhibition of NLRP3 inflammasome. This study revealed the mechanism of HWHP in treating NASH through experiments.

Comprehensive Analysis Identifies COL1A1, COL3A1, and POSTN as Key Genes Associated with Brain Metastasis in Patients with Breast Cancer

Objective

Brain metastasis (BM) is associated with a high mortality in patients with breast cancer (BC). Nevertheless, the molecular mechanisms of BM in BM remain uncertain. The study aims to identify the key genes in BC in relation to BM and to assess their prognostic value.

Methods

Two microarray datasets GSE125989 and GSE100534 were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) between primary BC and BM samples. The function enrichment analysis and protein-protein interaction (PPI) network were performed. We mapped hub genes into the Kaplan–Meier database for their correlations with BC survival.

Results

Venn diagram analysis showed an overlapped upregulated DEG and 18 overlapped downregulated ones between primary BC and BM samples. We constructed the PPI network, and top 5 hub genes were sorted out according to the node degree, including type I collagen α₁ chain (COL1A1), lumican (LUM), type III collagen α₁ chain (COL3A1), type V collagen α₂ chain (COL5A2), and periosteal protein (POSTN). The Kaplan–Meier database analysis found that COL1A1, COL3A1, and POSTN were significantly correlated with overall survival of BC patients.

Conclusion

The study suggests that COL1A1, COL3A1, and POSTN may be key genes associated with BM in patients with BC.

Pharmacological Mechanism of Ganlu Powder in the Treatment of NASH Based on Network Pharmacology and Molecular Docking

Nonalcoholic steatohepatitis (NASH), a progression of nonalcoholic fatty liver disease (NAFLD), is a clinical syndrome characterized by liver steatosis, inflammation, and hepatocellular damage. Ganlu powder (GLP) is a classic traditional Chinese medicine prescription that has shown favorable treatment effects on NASH. However, the underlying therapeutic mechanisms are still poorly understood. This study is aimed at exploring the potential mechanism of GLP in the treatment of NASH via network pharmacology and molecular docking. PubMed and CNKI databases were used to identify the components of GLP. Swiss and STITCH databases were employed to obtain corresponding drug targets. NASH targets were adopted from the Therapeutic Target Database (TTD), DisGeNET, DrugBank, GeneCards, and MalaCards databases. Cytoscape software was utilized to construct “drug-ingredient-target-disease” networks and the protein-protein interaction (PPI) network of GLP in NASH. AKT1 was identified as the key target. The GO functional enrichment analysis revealed that GLP might treat NASH by modulating the inflammatory response and regulating phosphatidylinositol 3-kinase signaling. The KEGG analysis showed that GLP might treat NASH by regulating the tumor necrosis factor (TNF) signal pathway by affecting the role of AKT1. According to the network pharmacology results, a virtual docking of active compounds with AKT1 was carried out, and the results indicated that the 7 components, berberine, epiberberine, jatrorrhizine, coptisine, palmatine, evodiamine, and rutecarpine, can bind stably with AKT1 and have higher binding energy than AKT1 inhibitors. The overall study findings suggest that GLP may treat NASH by regulating AKT1.

Oxidation and reduction of actin: Origin, impact in vitro and functional consequences in vivo

Actin is among the most abundant proteins in eukaryotic cells and assembles into dynamic filamentous networks regulated by many actin binding proteins. The actin cytoskeleton must be finely tuned, both in space and time, to fulfill key cellular functions such as cell division, cell shape changes, phagocytosis and cell migration. While actin oxidation by reactive oxygen species (ROS) at non-physiological levels are known for long to impact on actin polymerization and on the cellular actin cytoskeleton, growing evidence shows that direct and reversible oxidation/reduction of specific actin amino acids plays an important and physiological role in regulating the actin cytoskeleton. In this review, we describe which actin amino acid residues can be selectively oxidized and reduced in many different ways (e.g. disulfide bond formation, glutathionylation, carbonylation, nitration, nitrosylation and other oxidations), the cellular enzymes at the origin of these post-translational modifications, and the impact of actin redox modifications both in vitro and in vivo. We show that the regulated balance of oxidation and reduction of key actin amino acid residues contributes to the control of actin filament polymerization and disassembly at the subcellular scale and highlight how improper redox modifications of actin can lead to pathological conditions.

Identification of potential genes related to breast cancer brain metastasis in breast cancer patients

Brain metastases (BMs) usually develop in breast cancer (BC) patients. Thus, the molecular mechanisms of breast cancer brain metastasis (BCBM) are of great importance in designing therapeutic strategies to treat or prevent BCBM. The present study attempted to identify novel diagnostic and prognostic biomarkers of BCBM. Two datasets (GSE125989 and GSE100534) were obtained from the Gene Expression Omnibus (GEO) database to find differentially expressed genes (DEGs) in cases of BC with and without brain metastasis (BM). A total of 146 overlapping DEGs, including 103 up-regulated and 43 down-regulated genes, were identified. Functional enrichment analysis showed that these DEGs were mainly enriched for functions including extracellular matrix (ECM) organization and collagen catabolic fibril organization. Using protein-protein interaction (PPI) and principal component analysis (PCA) analysis, we identified ten key genes, including LAMA4, COL1A1, COL5A2, COL3A1, COL4A1, COL5A1, COL5A3, COL6A3, COL6A2, and COL6A1. Additionally, COL5A1, COL4A1, COL1A1, COL6A1, COL6A2, and COL6A3 were significantly associated with the overall survival of BC patients. Furthermore, COL6A3, COL5A1, and COL4A1 were potentially correlated with BCBM in human epidermal growth factor 2 (HER2) expression. Additionally, the miR-29 family might participate in the process of metastasis by modulating the cancer microenvironment. Based on datasets in the GEO database, several DEGs have been identified as playing potentially important roles in BCBM in BC patients.

Identification of Potential Biomarkers and Related Transcription Factors in Peripheral Blood of Tuberculosis Patients

Tuberculosis (TB), one major threat to humans, can infect one third of the worldwide population, and cause more than one million deaths each year. This study aimed to identify the effective diagnosis and therapy biomarkers of TB. Hence, we analyzed two microarray datasets (GSE54992 and GSE62525) derived from the Gene Expression Omnibus (GEO) database to find the differentially expressed genes (DEGs) of peripheral blood mononuclear cell (PBMC) between TB patients and healthy specimens. Functional and pathway enrichment of the DEGs were analyzed by Metascape database. Protein-protein interaction (PPI) network among the DEGs were constructed by STRING databases and visualized in Cytoscape software. The related transcription factors regulatory network of the DEGs was also constructed. A total of 190 DEGs including 36 up-regulated genes and 154 down-regulated genes were obtained in TB samples. Gene functional enrichment analysis showed that these DEGs were enriched in T cell activation, chemotaxis, leukocyte activation involved in immune response, cytokine secretion, head development, etc. The top six hub genes (namely, LRRK2, FYN, GART, CCR7, CXCR5, and FASLG) and two significant modules were got from PPI network of DEGs. Vital transcriptional factors, such as FoxC1 and GATA2, were discovered with close interaction with these six hub DEGs. By systemic bioinformatic analysis, many DEGs associated with TB were screened, and these identified hub DEGs may be potential biomarkers for diagnosis and treatment of TB in the future.

Gene network in pulmonary tuberculosis based on bioinformatic analysis

BACKGROUND

Pulmonary tuberculosis (PTB) is one of the serious infectious diseases worldwide; however, the gene network involved in the host response remain largely unclear.

METHODS

This study integrated two cohorts profile datasets GSE34608 and GSE83456 to elucidate the potential gene network and signaling pathways in PTB. Differentially expressed genes (DEGs) were obtained for Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis using Metascape database. Protein-Protein Interaction (PPI) network of DEGs was constructed by the online database the Search Tool for the Retrieval of Interacting Genes (STRING). Modules were identified by the plug-in APP Molecular Complex Detection (MCODE) in Cytoscape. GO and KEGG pathway of Module 1 were further analyzed by STRING. Hub genes were selected for further expression validation in dataset GSE19439. The gene expression level was also investigated in the dataset GSE31348 to display the change pattern during the PTB treatment.

RESULTS

Totally, 180 shared DEGs were identified from two datasets. Gene function and KEGG pathway enrichment revealed that DEGs mainly enriched in defense response to other organism, response to bacterium, myeloid leukocyte activation, cytokine production, etc. Seven modules were clustered based on PPI network. Module 1 contained 35 genes related to cytokine associated functions, among which 14 genes, including chemokine receptors, interferon-induced proteins and Toll-like receptors, were identified as hub genes. Expression levels of the hub genes were validated with a third dataset GSE19439. The signature of this core gene network showed significant response to Mycobacterium tuberculosis (Mtb) infection, and correlated with the gene network pattern during anti-PTB therapy.

CONCLUSIONS

Our study unveils the coordination of causal genes during PTB infection, and provides a promising gene panel for PTB diagnosis. As major regulators of the host immune response to Mtb infection, the 14 hub genes are also potential molecular targets for developing PTB drugs.

BATF Potentially Mediates Negative Regulation of PD-1/PD-Ls Pathway on T Cell Functions in Mycobacterium tuberculosis Infection

Background: Previously, we have found that blockade of PD-1/PD-Ls pathway could enhance CD4⁺ T cells-mediated protective immunity in patients with active tuberculosis (ATB). However, the mechanism of PD-1/PD-Ls pathway involved in negative regulation of anti-TB immunity has been still unclear. Recently, the study of human immunodeficiency virus (HIV) infection demonstrated that PD-1 could induce the expression of basic leucine zipper ATF-like transcription factor (BATF) to inhibit CD8⁺ T cell function. While the mechanism of immune regulation of BATF in Mycobacterium tuberculosis (M. tb) infection has not yet been elucidated.

Methods: We enrolled 104 participants including ATB patients (n = 66), latent tuberculosis infection (LTBI) (n = 16) and healthy control (HC) (n = 22). The expressions of BATF in peripheral blood CD4⁺ and CD8⁺ T cells from enrolled subjects were determined using flow cytometry. Intervention with PD-1/PD-Ls pathway was performed by using blocking antibodies or human PD-L1 fusion protein. Silencing BATF in peripheral blood mononuclear cells (PBMCs) by electroporation with siRNA. Real-time quantitative PCR, CFSE dilution assay and enzyme linked immunosorbent assay (ELISA) were employed to test T cell functions after BATF knockdown.

Results: The percentages of BATF⁺CD4⁺ (P = 0.0003 and P < 0.0001, respectively) and BATF⁺CD8⁺ (P = 0.0003 and P = 0.0003, respectively) cells were significantly increased in ATB patients compared with LTBI and HC. BATF-expressing PD-1⁺ T cells in CD4⁺ and CD8⁺ T cells were much higher in ATB group than those in LTBI group (P = 0.0426 and 0.0104, respectively) and HC group (P = 0.0133 and 0.0340, respectively). There was a positive correlation between BATF expression and PD-1 expression in ATB patients (for CD4⁺ T cells, r = 0.6761, P = 0.0158; for CD8⁺ T cells, r = 0.6104, P = 0.0350). BATF knockdown could enhance IL-2 and IFN-γ secretions (P = 0.0485 and 0.0473, respectively) and CD4⁺ T cells proliferation (P = 0.0041) in vitro.

Conclusions: In the context of tuberculosis, BATF mediates negative regulation of PD-1/PD-Ls pathway on T cell functions. BATF knockdown can improve cytokine secretion and cells proliferation in vitro.

The MICALs are a Family of F-actin Dismantling Oxidoreductases Conserved from Drosophila to Humans

Cellular form and function – and thus normal development and physiology – are specified via proteins that control the organization and dynamic properties of the actin cytoskeleton. Using the Drosophila model, we have recently identified an unusual actin regulatory enzyme, Mical, which is directly activated by F-actin to selectively post-translationally oxidize and destabilize filaments – regulating numerous cellular behaviors. Mical proteins are also present in mammals, but their actin regulatory properties, including comparisons among different family members, remain poorly defined. We now find that each human MICAL family member, MICAL-1, MICAL-2, and MICAL-3, directly induces F-actin dismantling and controls F-actin-mediated cellular remodeling. Specifically, each human MICAL selectively associates with F-actin, which directly induces MICALs catalytic activity. We also find that each human MICAL uses an NADPH-dependent Redox activity to post-translationally oxidize actin’s methionine (M) M44/M47 residues, directly dismantling filaments and limiting new polymerization. Genetic experiments also demonstrate that each human MICAL drives F-actin disassembly in vivo, reshaping cells and their membranous extensions. Our results go on to reveal that MsrB/SelR reductase enzymes counteract each MICAL’s effect on F-actin in vitro and in vivo. Collectively, our results therefore define the MICALs as an important phylogenetically-conserved family of catalytically-acting F-actin disassembly factors.