Hepatocellular carcinoma (HCC) is a devastating disease worldwide. Though many efforts have been made to elucidate the process of HCC, its molecular mechanisms of development remain elusive due to its complexity. To explore the stepwise carcinogenic process from pre-neoplastic lesions to the end stage of HCC, we employed weighted gene co-expression network analysis (WGCNA) which has been proved to be an effective method in many diseases to detect co-expressed modules and hub genes using eight pathological stages including normal, cirrhosis without HCC, cirrhosis, low-grade dysplastic, high-grade dysplastic, very early and early, advanced HCC and very advanced HCC. Among the eight consecutive pathological stages, five representative modules are selected to perform canonical pathway enrichment and upstream regulator analysis by using ingenuity pathway analysis (IPA) software. We found that cell cycle related biological processes were activated at four neoplastic stages, and the degree of activation of the cell cycle corresponded to the deterioration degree of HCC. The orange and yellow modules enriched in energy metabolism, especially oxidative metabolism, and the expression value of the genes decreased only at four neoplastic stages. The brown module, enriched in protein ubiquitination and ephrin receptor signaling pathways, correlated mainly with the very early stage of HCC. The darkred module, enriched in hepatic fibrosis/hepatic stellate cell activation, correlated with the cirrhotic stage only. The high degree hub genes were identified based on the protein-protein interaction (PPI) network and were verified by Kaplan-Meier survival analysis. The novel five high degree hub genes signature that was identified in our study may shed light on future prognostic and therapeutic approaches. Our study brings a new perspective to the understanding of the key pathways and genes in the dynamic changes of HCC progression. These findings shed light on further investigations.

BackgroundTumor microenvironment, in particular the stroma, plays an important role in breast cancer cell invasion and metastasis. Investigation of the molecular characteristics of breast cancer stroma may reveal targets for future study.MethodsThe transcriptome profiles of breast cancer stroma and normal breast stroma were compared to identify differentially expressed genes (DEGs). The method was analysis of GSE26910 and GSE10797 datasets. Common DEGs were identified and then analyses of enriched pathways and hub genes were performed.ResultsA total of 146 DEGs were common to GSE26910 and GSE10797. The enriched pathways were associated with “extracellular matrix (ECM) organization”, “ECM-receptor interaction” and “focal adhesion”. Network analysis identified six key genes, including JUN, FOS, ATF3, STAT1, COL1A1 and FN1. Notably, COL1A1 and FN1 were identified for the first time as cancer stromal key genes associated with breast cancer invasion and metastasis. Oncome analysis showed that the high expression levels of COL1A1 and FN1 correlated to an advanced stage of breast cancer and poor clinical outcomes.ConclusionsWe found that several conserved tumor stromal genes might regulate breast cancer invasion through ECM remodeling. The clinical outcome analyses of COL1A1 and FN1 suggest these two genes are promising targets for future studies.

ObjectivesThe aim of this study was to identify key pathological genes in osteoarthritis (OA).MethodsWe searched and downloaded mRNA expression data from the Gene Expression Omnibus database to identify differentially expressed genes (DEGs) of joint synovial tissues from OA and normal individuals. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analyses were used to assess the function of identified DEGs. The protein-protein interaction (PPI) network and transcriptional factors (TFs) regulatory network were used to further explore the function of identified DEGs. The quantitative real-time polymerase chain reaction (qRT-PCR) was applied to validate the result of bioinformatics analysis. Electronic validation was performed to verify the expression of selected DEGs. The diagnosis value of identified DEGs was accessed by receiver operating characteristic (ROC) analysis.ResultsA total of 1085 DEGs were identified. KEGG pathway analysis displayed that Wnt was a significantly enriched signalling pathway. Some hub genes with high interactions such as USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, and ZNF423 were identified in the PPI and TFs network. The results of qRT-PCR showed that GADD45B, ADAMTS1, and TFAM were down-regulated in joint synovial tissues of OA, which was consistent with the bioinformatics analysis. The expression levels of USP46, CPVL, FOSL2, and PTGS1 in electronic validation were compatible with the bio-informatics result. CPVL and TFAM had a potential diagnostic value for OA based on the ROC analysis.ConclusionThe deregulated genes including USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, ZNF423, ADAMTS1, and TFAM might be involved in the pathology of OA.Cite this article: X. Zhang, Y. Bu, B. Zhu, Q. Zhao, Z. Lv, B. Li, J. Liu. Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis. Bone Joint Res 2018;7:298–307. DOI: 10.1302/2046-3758.74.BJR-2017-0245.R1.