Classification and characterisation of extracellular vesicles-related tuberculosis subgroups and immune cell profiles

Around the world, tuberculosis (TB) remains one of the most common causes of morbidity and mortality. The molecular mechanism of Mycobacterium tuberculosis (Mtb) infection is still unclear. Extracellular vesicles (EVs) play a key role in the onset and progression of many disease states and can serve as effective biomarkers or therapeutic targets for the identification and treatment of TB patients. We analysed the expression profile to better clarify the EVs characteristics of TB and explored potential diagnostic markers to distinguish TB from healthy control (HC). Twenty EVs-related differentially expressed genes (DEGs) were identified, and 17 EVs-related DEGs were up-regulated and three DEGs were down-regulated in TB samples, which were related to immune cells. Using machine learning, a nine EVs-related gene signature was identified and two EVs-related subclusters were defined. The single-cell RNA sequence (scRNA-seq) analysis further confirmed that these hub genes might play important roles in TB pathogenesis. The nine EVs-related hub genes had excellent diagnostic values and accurately estimated TB progression. TB's high-risk group had significantly enriched immune-related pathways, and there were substantial variations in immunity across different groups. Furthermore, five potential drugs were predicted for TB using CMap database. Based on the EVs-related gene signature, the TB risk model was established through a comprehensive analysis of different EV patterns, which can accurately predict TB. These genes could be used as novel biomarkers to distinguish TB from HC. These findings lay the foundation for further research and design of new therapeutic interventions aimed at treating this deadly infectious disease.     

环氧化物水解酶2在肝细胞癌中的表达水平与功能作用

本研究通过公共数据和实验数据,全面分析环氧化物水解酶2(epoxide hydrolase 2, EPHX2)在肝细胞癌中的表达情况、功能作用以及预后意义。利用GEO和MitoCarta数据集,筛选肝细胞癌中呈差异表达的线粒体相关基因;利用TCGA数据库分析EPHX2及其相关基因在肝细胞癌中的表达水平;运行R包绘制Kaplan-Meier生存曲线和功能富集分析;基于STRING和GSEA构建蛋白质互作网络和基因集富集分析;荧光定量PCR和GEO数据集验证EPHX2在肝细胞癌中的表达水平。本研究共筛选得到15个在肝细胞癌中呈差异表达的线粒体相关基因。EPHX2在肝细胞癌组织中的表达水平显著降低(P<0.01)。EPHX2表达水平与肝癌患者性别、分期和级别有关,而与年龄、T分期等因素无关。与EPHX2低表达组肝癌患者相比,EPHX2高表达组肝癌患者预后较好。功能富集结果显示,EPHX2与补体途径、脂肪酸降解等信号通路有关。蛋白质互作网络结果显示,EPHX2与HAO1、AGXT、ACOX1、GSTκ1、SCP-2、CAT、CYP2C8,CYP2C9,CYP2B6,和CYP2J2等密切相关。GSEA结果显示,EPHX2低表达组与肝癌细胞增殖、肝癌复发等基因集正相关。荧光定量PCR和GEO数据集验证结果显示,EPHX2在肝细胞癌组织和肝癌细胞株中呈显著低表达。EPHX2在肝细胞癌中呈显著低表达,提示其可能在肝细胞癌发生发展过程中发挥抑癌基因作用,但具体作用机制还需进一步验证。     

Identification of a six-gene prognostic signature for oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignancies worldwide, and its morbidity and mortality have increased in the near term. Consequently, the purpose of the present study was to identify the notable differentially expressed genes (DEGs) involved in their pathogenesis to obtain new biomarkers or potential therapeutic targets for OSCC. The gene expression profiles of the microarray datasets GSE85195, GSE23558, and GSE10121 were obtained from the Gene Expression Omnibus (GEO) database. After screening the DEGs in each GEO dataset, 249 DEGs in OSCC tissues were obtained. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway enrichment analysis was employed to explore the biological functions and pathways of the above DEGs. A protein–protein interaction network was constructed to obtain a central gene. The corresponding total survival information was analyzed in patients with oral cancer from The Cancer Genome Atlas (TCGA). A total of six candidate genes (CXCL10, OAS2, IFIT1, CCL5, LRRK2, and PLAUR) closely related to the survival rate of patients with oral cancer were identified, and expression verification and overall survival analysis of six genes were performed based on TCGA database. Time-dependent receiver operating characteristic curve analysis yields predictive accuracy of the patient's overall survival. At the same time, the six genes were further verified by quantitative real-time polymerase chain reaction using samples obtained from the patients recruited to the present study. In conclusion, the present study identified the prognostic signature of six genes in OSCC for the first time via comprehensive bioinformatics analysis, which could become potential prognostic markers for OCSS and may provide potential therapeutic targets for tumors.     

CT23 knockdown attenuating malignant behaviors of hepatocellular carcinoma cell is associated with upregulation of metallothionein 1

The cancer-testis antigen 23 (CT23) gene has been reported in association with the pathogenesis and progress of hepatocellular carcinoma (HCC). However, the alterations of gene expression profiling induced by CT23 knockdown in HCC cells remains largely unknown. In this study, the RNA interfering (RNAi) method was used to silence CT23 expression in BEL-7404 cells. Microarray analysis was performed on mRNA extracted from the CT23 knockdown cells and the control cells to determine the alterations of gene expression profiles. The result showed a total of 1051 genes expressed differentially (two-fold change), including 470 genes upregulated and 581 gene downregulated in the CT23 knockdown cells. A bioinformatic analysis showed that the functional differentially expressed genes (DEGs) were linked to cell proliferation, migration, and apoptosis, and metallothionein 1 (MT1) attained the maximum enrichment scores in functional annotation, classification, and pathway analysis of DEGs. Furthermore, Western blot analysis and cell behaviors assays verified that CT23 modulates cell proliferation, migration, and apoptosis by regulating MT1 expression in HCC cells and non-neoplastic hepatocytes. In summary, downregulated CT23 gene in BEL-7404 cells might change the expressions of carcinogenesis and progression related genes in HCC by upregulating MT1 expression, which would provide insight into searching for a novel therapeutic target for HCC.     

Identification of Prognostic Biomarkers by Combined mRNA and miRNA Expression Microarray Analysis in Pancreatic Cancer

Pancreatic cancer is the fourth leading cause for cancer-related death, and early diagnosis is one key to improve the survival rate of this disease. Molecular biomarkers are an important method for diagnostic use in pancreatic cancer. We used data from three mRNA microarray datasets and a microRNA dataset (GSE16515, GSE15471, GSE28735, and GSE41372) to identify potential key genes. Differentially expressed genes (DEGs) and microRNAs (DEMs) were identified. Functional, pathway enrichment, and protein-protein interaction analyses were performed on common DEGs across all datasets. The target genes of the DEMs were identified. DEMs targets that were also DEGs were further scrutinized using overall survival analysis. A total of 236 DEGs and 21 DEMs were identified. There were a total of four DEGs (ECT2, NR5A2, NRP2, and TGFBI), which were also predicted target genes of DEMs. Overall survival analysis showed that high expression levels of three of these genes (ECT2, NRP2, and TGFBI) were associated with poor overall survival for pancreatic cancer patients. The basic expression of DEGs in pancreas stood lower level in various organ tissues. The expression of ECT2 and NRP2 was higher in different pancreatic cancer cell lines than normal pancreas cell line. Knockout of ECT2 by Crispr Cas9 gene editing system decreased proliferation and migration ability in pancreatic cancer cell line MiaPaCa2. In conclusion, we think that data mining method can do well in biomarker screening, and ECT2 and NRP2 can play as potential biomarker or therapy target by Crispr Cas9 in pancreatic cancer.     

基于TCGA数据库分析甲状腺癌基因表达谱

为分析甲状腺癌基因表达谱,筛选疾病相关的基因标志物。基于肿瘤基因组图谱(TCGA)数据库中的甲状腺癌基因表达数据,运用R/Bioconductor统计平台进行数据处理与统计学分析。分别应用edgeR算法和limma算法选取肿瘤组织与对照组间倍数改变 > 2,P< 0.05的基因作为差异基因;进一步运用Medcalc统计软件进行受试者工作特征曲线(ROC)分析,鉴定出有诊断标志物潜在应用价值的基因标志物。通过两种运算方法筛选出甲状腺癌组织中存在着1 945个差异基因(上调基因1 033个,下调基因912个);根据差异倍数进一步鉴定出11个基因在肿瘤组织中表达上调,且对鉴别肿瘤组与对照组有较好的应用价值。本研究分析了TCGA中的甲状腺癌表达谱数据,鉴定出了与疾病诊断显著相关的差异表达基因,能够为探索疾病发生发展机制及寻找新型分子标志物提供依据。     

Identification of differentially expressed genes and biological characteristics of colorectal cancer by integrated bioinformatics analysis

Colorectal cancer (CRC) ranks as one of the most common malignant tumors worldwide. Its mortality rate has remained high in recent years. Therefore, the aim of this study was to identify significant differentially expressed genes (DEGs) involved in its pathogenesis, which may be used as novel biomarkers or potential therapeutic targets for CRC. The gene expression profiles of GSE21510, GSE32323, GSE89076, and GSE113513 were downloaded from the Gene Expression Omnibus (GEO) database. After screening DEGs in each GEO data set, we further used the robust rank aggregation method to identify 494 significant DEGs including 212 upregulated and 282 downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by DAVID and the KOBAS online database, respectively. These DEGs were shown to be significantly enriched in different cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network. The module analysis was performed by the MCODE plug-in of Cytoscape based on the whole network. We finally filtered out seven hub genes by the cytoHubba plug-in, including PPBP, CCL28, CXCL12, INSL5, CXCL3, CXCL10, and CXCL11. The expression validation and survival analysis of these hub genes were analyzed based on The Cancer Genome Atlas database. In conclusion, the robust DEGs associated with the carcinogenesis of CRC were screened through the GEO database, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for CRC.     

Loss of cytoplasmic CDK1 predicts poor survival in human lung cancer and confers chemotherapeutic resistance

The dismal lethality of lung cancer is due to late stage at diagnosis and inherent therapeutic resistance. The incorporation of targeted therapies has modestly improved clinical outcomes, but the identification of new targets could further improve clinical outcomes by guiding stratification of poor-risk early stage patients and individualizing therapeutic choices. We hypothesized that a sequential, combined microarray approach would be valuable to identify and validate new targets in lung cancer. We profiled gene expression signatures during lung epithelial cell immortalization and transformation, and showed that genes involved in mitosis were progressively enhanced in carcinogenesis. 28 genes were validated by immunoblotting and 4 genes were further evaluated in non-small cell lung cancer tissue microarrays. Although CDK1 was highly expressed in tumor tissues, its loss from the cytoplasm unexpectedly predicted poor survival and conferred resistance to chemotherapy in multiple cell lines, especially microtubule-directed agents. An analysis of expression of CDK1 and CDK1-associated genes in the NCI60 cell line database confirmed the broad association of these genes with chemotherapeutic responsiveness. These results have implications for personalizing lung cancer therapy and highlight the potential of combined approaches for biomarker discovery.     

Comprehensive analysis of microRNA-messenger RNA regulatory network in gemcitabine-resistant bladder cancer cells

Chemotherapy is still a standard treatment of unresectable bladder cancer or distant metastases. The chemotherapy resistance always occurs after a period of treatment indicating poor prognosis. The current study aimed to explore the molecular mechanism of chemoresistance in bladder cancer cells. The gene expression profiles of GSE77883, including three untreated T24 cells samples and three gemcitabine-resistant T24 cells samples, was downloaded from Gene Expression Omnibus database. The screening of differentially expressed genes (DEGs), gene function analysis, and interaction prediction between microRNAs (miRNAs) and DEGs were performed by R software. The protein-protein interaction (PPI) and miRNA-DEGs networks were constructed and visualized by Cytoscape software. Then, the small molecules, with potential synergistic or antagonistic effects to gemcitabine resistance, were identified using the Connectivity Map database. Finally, gemcitabine-resistant T24 cell line was established and key genes were validated by quantitative real-time polymerase chain reaction (qRT-PCR). In total, 536 upregulated and 513 downregulated genes were screened and mainly enriched in oxidative stress response and signaling pathways related to extracellular matrix–receptor interaction and focal adhesion. PPI network showed interleukin 6, tumor necrosis factor, kinesin family member 11, and BUB1 mitotic checkpoint serine/threonine kinase B were key genes. The miRNA-DEGs regulatory networks included 18 miRNAs and 185 DEGs, including miR-182-5p, miR-590-3p, miR-320a and serum- and glucocorticoid-regulated kinase 1 (SGK1). Then, the related key genes and miRNAs were confirmed by qRT-PCR. Furthermore, 81 small molecules with antagonistic or synergistic effect to GEM were screened. We have investigated the molecular mechanisms driving GEM-resistance in bladder cancer cells that would contribute to the development of chemotherapy for advanced bladder cancer.     

Temporal expression and functional analysis of long non‐coding RNAs in colorectal cancer initiation

Long non‐coding RNAs (lncRNAs) have potential applications in clinical diagnosis and targeted cancer therapies. However, the expression profile of lncRNAs in colorectal cancer (CRC) initiation is still unclear. In this study, the expression profiles of lncRNAs and mRNAs were determined by microarray at specific tumour stages in an AOM/DSS‐induced primary colon cancer model. The temporal expression of lncRNAs was analysed by K‐means clustering. Additionally, weighted correlation network analysis (WGCNA) and gene ontology analysis were performed to construct co‐expression networks and establish functions of the identified lncRNAs and mRNAs. Our results suggested that 4307 lncRNAs and 5798 mRNAs are deregulated during CRC initiation. These differential expression genes (DEGs) exhibited a clear correlation with the differential stage of tumour initiation. WGCNA results suggested that a series of hub lncRNAs are involved in regulating cell stemness, colon inflammation, oxidative stress response and cell death at each stage. Among them, lncRNA H19 was up‐regulated in colon tumours and correlated with poor patient prognosis. Collectively, we have been the first to demonstrate the temporal expression and function of lncRNAs in CRC initiation. These results provide novel diagnosis and therapy targets for CRC.     

Identification of key genes and pathways of thyroid cancer by integrated bioinformatics analysis

Thyroid cancer is a common endocrine malignancy with a rapidly increasing incidence worldwide. Although its mortality is steady or declining because of earlier diagnoses, its survival rate varies because of different tumour types. Thus, the aim of this study was to identify key biomarkers and novel therapeutic targets in thyroid cancer. The expression profiles of GSE3467, GSE5364, GSE29265 and GSE53157 were downloaded from the Gene Expression Omnibus database, which included a total of 97 thyroid cancer and 48 normal samples. After screening significant differentially expressed genes (DEGs) in each data set, we used the robust rank aggregation method to identify 358 robust DEGs, including 135 upregulated and 224 downregulated genes, in four datasets. Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analyses of DEGs were performed by DAVID and the KOBAS online database, respectively. The results showed that these DEGs were significantly enriched in various cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network, and modules analysis was performed. Finally, we filtered out five hub genes, including LPAR5, NMU, FN1, NPY1R, and CXCL12, from the whole network. Expression validation and survival analysis of these hub genes based on the The Cancer Genome Atlas database suggested the robustness of the above results. In conclusion, these results provided novel and reliable biomarkers for thyroid cancer, which will be useful for further clinical applications in thyroid cancer diagnosis, prognosis and targeted therapy.     

Network cluster analysis of protein–protein interaction network identified biomarker for early onset colorectal cancer

Colorectal cancer (CRC) is a major cause of morbidity and mortality throughout the world. However, the genetic alterations and molecular mechanism of the early onset CRCs are not fully investigated. The present study aimed to characterize early onset CRC by analyzing its gene expression compared with normal controls and to identify network-based biomarkers of early onset CRC. The gene expression profiles of early onset CRC were downloaded from Gene Expression Omnibus and the differentially expressed genes (DEGs) in CRC patients were identified. Then, a protein–protein interaction (PPI) network was constructed and the clusters in PPI were analyzed by ClusterONE. Furthermore, the gene ontology functional analysis and pathway enrichment analysis were conducted to the modules in PPI network. A systems biology approach integrating microarray data and PPI was further applied to construct a PPI network in CRC. Total 631 DEGs were identified from the early onset CRC compared to healthy controls. These genes were found to be involved in several biological processes, including cell communication, cell proliferation, cell shape and apoptosis. Five functional modules which may play important roles in the initiation of early onset CRC were identified from the PPI network. Functional annotation revealed that these five modules were involved in the pathways of signal transduction, carcinogenesis and metastasis. The hub nodes of these five modules, CDC42, TEX11, QKI, CAV1 and FN1, may serve as the biomarkers of early onset CRC and have the potential to be targets for therapeutic intervention. However, further investigations are still needed to confirm our findings.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

Small interfering RNA–mediated gene suppression as a therapeutic intervention in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the lethal and difficult-to-cure cancers worldwide. Owing to the late diagnosis and drug resistance of malignant hepatocytes, treatment of this cancer by conventional chemotherapy agents is challenging, and researchers are seeking new alternative treatment options to overcome therapy resistance in this neoplasm. RNA interference (RNAi) is a potent and specific approach in targeting gene expression and has emerged as a novel therapeutic tool for many diseases, including cancers. Small interfering RNA (siRNA) is a type of RNAi that is produced intracellularly from exogenous synthetic oligonucleotides and can selectively knock down target gene expression in a sequence-specific manner. Various factors play roles in the initiation and progression of HCC and provide multiple candidate targets for siRNA intervention. In addition, due to the liver's unique architecture and availability of some hepatic siRNA delivery methods, this organ has received much more attention as a target tissue for such oligonucleotide action. Recent advances in designing nanoparticle systems for the in vivo delivery of siRNAs have markedly enhanced the potency of siRNA-mediated gene silencing under clinical development for HCC therapy. The utility of siRNAs as anti-HCC agents is the subject of the current review. siRNA-based gene therapies could be one of the main feasible approaches for HCC therapy in the future.