Weighted gene correlation network analysis identifies RSAD2, HERC5, and CCL8 as prognostic candidates for breast cancer

As the most commonly diagnosed malignant tumor in female population, the prognosis of breast cancer is affected by complex gene interaction networks. In this research weighted gene co-expression network analysis (WGCNA) would be utilized to build a gene co-expression network to identify potential biomarkers for prediction the prognosis of patients with breast cancer. We downloaded GSE25065 from Gene Expression Omnibus database as the test set. GSE25055 and GSE42568 were utilized to validate findings in the research. Seven modules were established in the GSE25065 by utilizing average link hierarchical clustering. Three hub genes, RSAD2, HERC5, and CCL8 were screened out from the significant module (R ² = 0.44), which were considerably interrelated to worse prognosis. Within test dataset GSE25065, RSAD2, and CCL8 were correlated with tumor stage, grade, and lymph node metastases, whereas HERC5 was correlated with lymph node metastases and tumor grade. In the validation dataset GSE25055 and RSAD2 expression was correlated with tumor grade, stage, and size, whereas HERC5 was related to tumor stage and tumor grade, and CCL8 was associated with tumor size and tumor grade. Multivariable survival analysis demonstrated that RSAD2, HERC5, and CCL8 were independent risk factors. In conclusion, the WGCNA analysis conducted in this study screened out novel prognostic biomarkers of breast cancer. Meanwhile, further in vivo and in vitro studies are required to make the clear molecular mechanisms.     

Identification of gene expression profiles and key genes in subchondral bone of osteoarthritis using weighted gene coexpression network analysis

Osteoarthritis (OA) significantly influences the quality life of people around the world. It is urgent to find an effective way to understand the genetic etiology of OA. We used weighted gene coexpression network analysis (WGCNA) to explore the key genes involved in the subchondral bone pathological process of OA. Fifty gene expression profiles of GSE51588 were downloaded from the Gene Expression Omnibus database. The OA‐associated genes and gene ontologies were acquired from JuniorDoc. Weighted gene coexpression network analysis was used to find disease‐related networks based on 21756 gene expression correlation coefficients, hub‐genes with the highest connectivity in each module were selected, and the correlation between module eigengene and clinical traits was calculated. The genes in the traits‐related gene coexpression modules were subject to functional annotation and pathway enrichment analysis using ClusterProfiler. A total of 73 gene modules were identified, of which, 12 modules were found with high connectivity with clinical traits. Five modules were found with enriched OA‐associated genes. Moreover, 310 OA‐associated genes were found, and 34 of them were among hub‐genes in each module. Consequently, enrichment results indicated some key metabolic pathways, such as extracellular matrix (ECM)‐receptor interaction (hsa04512), focal adhesion (hsa04510), the phosphatidylinositol 3'‐kinase (PI3K)‐Akt signaling pathway (PI3K‐AKT) (hsa04151), transforming growth factor beta pathway, and Wnt pathway. We intended to identify some core genes, collagen (COL)6A3, COL6A1, ITGA11, BAMBI, and HCK, which could influence downstream signaling pathways once they were activated. In this study, we identified important genes within key coexpression modules, which associate with a pathological process of subchondral bone in OA. Functional analysis results could provide important information to understand the mechanism of OA.     

Identification of SLITRK6 as a Novel Biomarker in hepatocellular carcinoma by comprehensive bioinformatic analysis

Hepatocellular carcinoma (HCC) is the most common primary malignancy of the adult liver and morbidity are increasing in recent years, however, there is still no effective strategy to prevent and diagnose HCC. Therefore, it is urgent to research the effective biomarker to predict clinical outcomes of HCC tumorigenesis. In the current study, differentially expressed genes in HCC and normal tissues were investigated using the Gene Expression Omnibus (GEO) dataset GSE144269 and The Cancer Genome Atlas (TCGA). Gene differential expression analysis and weighted correlation network analysis (WGCNA) methods were used to identify nine and 16 key gene modules from the GEO dataset and TCGA dataset, respectively, in which the green module in the GEO dataset and magenta module in TCGA were significantly correlated with HCC occurrence. Third, the enrichment score of gene function annotation results showed that these two key modules focus on the positive regulation of inflammatory response and cell differentiation, etc. Besides, PPI network analysis, mutation analysis, and survival analysis found that SLITRK6 had high connectivity, and its mutation significantly impacted overall survival. In addition, SLITRK6 was found to be low expressed in tumor cells. To summarize, SLITRK6 mutation was found to significantly affect the occurrence and prognosis of HCC. SLITRK6 was confirmed as a new potential gene target for HCC, which may provide a new theoretical basis for personalized diagnosis and chemotherapy of HCC in the future.     

伴HBV感染性肝癌调控枢纽基因筛选与初步鉴定

慢性乙型肝炎病毒(Hepatitis B virus,HBV)感染引起的原发性肝癌涉及多种基因、转录本和蛋白质的相互作用及调控。从单个基因的角度来看,某个基因的表达量的改变只能对肝癌发生发展的局部作出解释而无法从整体行为进行深入和全面的探索,无法满足高度复杂性的调控研究需要。筛选乙肝相关性肝癌的基因芯片数据获取差异表达基因后,应用加权基因共表达网络分析算法构建基因共表达网络,识别与肝癌发生相关的模块,利用可视化筛选枢纽基因,并针对枢纽基因进行基因本体富集分析和初步验证。富集分析和文献挖掘一致发现,某些枢纽基因确实与多种癌症的发生与发展存在显著的关联。权重基因共表达网络分析方法被证明是一个高效的系统生物学方法,应用该方法发现了新的HBV相关性肝癌枢纽基因。经实验验证,发现枢纽基因SHARPIN促进细胞迁移。该研究对肝癌发生的调控机制以及发现HBV慢性感染导致肝癌的新型诊断标志物和(或)药物作用靶点提供了新的视野。     

Identification of hub genes in colorectal cancer based on weighted gene co-expression network analysis and clinical data from The Cancer Genome Atlas

Colorectal cancer (CRC) is one of the most common tumors worldwide and is associated with high mortality. Here we performed bioinformatics analysis, which we validated using immunohistochemistry in order to search for hub genes that might serve as biomarkers or therapeutic targets in CRC. Based on data from The Cancer Genome Atlas (TCGA), we identified 4832 genes differentially expressed between CRC and normal samples (1562 up-regulated and 3270 down-regulated in CRC). Gene ontology (GO) analysis showed that up-regulated genes were enriched mainly in organelle fission, cell cycle regulation, and DNA replication; down-regulated genes were enriched primarily in the regulation of ion transmembrane transport and ion homeostasis. Weighted gene co-expression network analysis (WGCNA) identified eight gene modules that were associated with clinical characteristics of CRC patients, including brown and blue modules that were associated with cancer onset. Analysis of the latter two hub modules revealed the following six hub genes: adhesion G protein-coupled receptor B3 (BAI3, also known as ADGRB3), cyclin F (CCNF), cytoskeleton-associated protein 2 like (CKAP2L), diaphanous-related formin 3 (DIAPH3), oxysterol binding protein-like 3 (OSBPL3), and RERG-like protein (RERGL). Expression levels of these hub genes were associated with prognosis, based on Kaplan–Meier survival analysis of data from the Gene Expression Profiling Interactive Analysis database. Immunohistochemistry of CRC tumor tissues confirmed that OSBPL3 is up-regulated in CRC. Our findings suggest that CCNF, DIAPH3, OSBPL3, and RERGL may be useful as therapeutic targets against CRC. BAI3 and CKAP2L may be novel biomarkers of the disease.     

Systematic prediction of key genes for ovarian cancer by co-expression network analysis

Ovarian cancer (OC) is the most lethal gynaecological malignancy, characterized by high recurrence and mortality. However, the mechanisms of its pathogenesis remain largely unknown, hindering the investigation of the functional roles. This study sought to identify key hub genes that may serve as biomarkers correlated with prognosis. Here, we conduct an integrated analysis using the weighted gene co-expression network analysis (WGCNA) to explore the clinically significant gene sets and identify candidate hub genes associated with OC clinical phenotypes. The gene expression profiles were obtained from the MERAV database. Validations of candidate hub genes were performed with RNASeqV2 data and the corresponding clinical information available from The Cancer Genome Atlas (TCGA) database. In addition, we examined the candidate genes in ovarian cancer cells. Totally, 19 modules were identified and 26 hub genes were extracted from the most significant module (R² = .53) in clinical stages. Through the validation of TCGA data, we found that five hub genes (COL1A1, DCN, LUM, POSTN and THBS2) predicted poor prognosis. Receiver operating characteristic (ROC) curves demonstrated that these five genes exhibited diagnostic efficiency for early-stage and advanced-stage cancer. The protein expression of these five genes in tumour tissues was significantly higher than that in normal tissues. Besides, the expression of COL1A1 was associated with the TAX resistance of tumours and could be affected by the autophagy level in OC cell line. In conclusion, our findings identified five genes could serve as biomarkers related to the prognosis of OC and may be helpful for revealing pathogenic mechanism and developing further research.     

Fifteen hub genes associated with progression and prognosis of clear cell renal cell carcinoma identified by coexpression analysis

Renal cell carcinoma (RCC) is the most common type of renal tumor, and the clear cell renal cell carcinoma (ccRCC) is the most frequent subtype. In this study, our aim is to identify potential biomarkers that could effectively predict the prognosis and progression of ccRCC. First, we used The Cancer Genome Atlas (TCGA) RNA-sequencing (RNA-seq) data of ccRCC to identify 2370 differentially expressed genes (DEGs). Second, the DEGs were used to construct a coexpression network by weighted gene coexpression network analysis (WGCNA). Moreover, we identified the yellow module, which was strongly related to the histologic grade and pathological stage of ccRCC. Then, the functional annotation of the yellow module and single-samples gene-set enrichment analysis of DEGs were performed and mainly enriched in cell cycle. Subsequently, 18 candidate hub genes were screened through WGCNA and protein–protein interaction (PPI) network analysis. After verification of TCGA’s ccRCC data set, Gene Expression Omnibus (GEO) data set (GSE73731) and tissue validation, we finally identified 15 hub genes that can actually predict the progression of ccRCC. In addition, by using survival analysis, we found that patients of ccRCC with high expression of each hub gene were more likely to have poor prognosis than those with low expression. The receiver operating characteristic curve showed that each hub gene could effectively distinguish between localized and advanced ccRCC. In summary, our study indicates that 15 hub genes have great predictive value for the prognosis and progression of ccRCC, and may contribute to the exploration of the pathogenesis of ccRCC.     

基于SNP共表达网络肝癌分子分型及预后分析

基于SNP突变数据与mRNA表达谱关联分析,构建一种肝癌分子分型方法并对比不同分型预后的差异,并对不同分型肝癌的发生发展机制进一步研究。首先通过TCGA数据库收集359例肝细胞癌患者的SNP突变数据和mRNA表达数据,采用Wilcoxon秩和检验,筛选突变后差异表达基因,并通过生物信息学工具String和Cytoscape 构建差异表达基因的蛋白互作网络,筛选连接度最高的10个Hub基因。利用Consensus Cluster Plus软件包,基于Hub基因mRNA表达水平构建NMF分子分型模型,再结合生存数据评估各分型患者的预后。最后利用加权基因共表达网络分析(WGCNA),识别与肝癌分子分型相关的模块,并针对关键模块的基因进行通路富集,从而对不同分型肝癌的基因表达谱进行比较。结果:NMF模型将肝癌分为高危、低危2个分型,其中CDKN2A和FOXO1基因对分型贡献度高。生存分析显示低危组患者的生存情况显著优于高危组,高危组富集多个与肿瘤细胞侵蚀、转移、复发过程相关的信号通路,低危组则与细胞周期和胰液分泌相关。本研究在无先验性信息的前提下,基于突变后显著差异表达的Hub基因表达水平构建的肝癌分子分型对肝癌患者预后评估具有一定的指导意义,其中CDKN2A和FOXO1突变是肝癌患者的不良预后因素,针对二者的靶向药研发,可能为肝癌患者提供新的治疗策略。     

环氧化物水解酶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在肝细胞癌中呈显著低表达,提示其可能在肝细胞癌发生发展过程中发挥抑癌基因作用,但具体作用机制还需进一步验证。     

Prognostic genes of hepatocellular carcinoma based on gene coexpression network analysis

Hepatocellular carcinoma (HCC) is the most common subtype in liver cancer whose prognosis is affected by malignant progression associated with complex gene interactions. However, there is currently no available biomarkers associated with HCC progression in clinical application. In our study, RNA sequencing expression data of 50 normal samples and 374 tumor samples was analyzed and 9225 differentially expressed genes were screened. Weighted gene coexpression network analysis was then conducted and the blue module we were interested was identified by calculating the correlations between 17 gene modules and clinical features. In the blue module, the calculation of topological overlap was applied to select the top 30 genes and these 30 genes were divided into the green group (11 genes) and the yellow group (19 genes) through searching whether these genes were validated by in vitro or in vivo experiments. The genes in the green group which had never been validated by any experiments were recognized as hub genes. These hub genes were subsequently validated by a new data set GSE76427 and KM Plotter Online Tool, and the results indicated that 10 genes (FBXO43, ARHGEF39, MXD3, VIPR1, DNASE1L3, PHLDA1, CSRNP1, ADR2B, C1RL, and CDC37L1) could act as prognosis and progression biomarkers of HCC. In summary, 10 genes who have never been mentioned in HCC were identified to be associated with malignant progression and prognosis of patients. These findings may contribute to the improvement of the therapeutic decision, risk stratification, and prognosis prediction for HCC patients.     

Identification of gene coexpression modules,hub genes,and pathways related to spinal cord injury using integrated bioinformatics methods

Spinal cord injury (SCI) is characterized by dramatic neurons loss and axonal regeneration suppression. The underlying mechanism associated with SCI-induced immune suppression is still unclear. Weighted gene coexpression network analysis (WGCNA) is now widely applied for the identification of the coexpressed modules, hub genes, and pathways associated with clinic traits of diseases. We performed this study to identify hub genes associated with SCI development. Gene Expression Omnibus (GEO) data sets GSE45006 and GSE20907 were downloaded and the significant correlativity and connectivity between them were detected using WGCNA. Three significant consensus modules, including 567 eigengenes, were identified from the master GSE45006 data following the preconditions of approximate scale-free topology for WGCNA. Further bioinformatics analysis showed these eigengenes were involved in inflammatory and immune responses in SCI. Three hub genes Rac2, Itgb2, and Tyrobp and one pathway “natural killer cell-mediated cytotoxicity” were identified following short time-series expression miner, protein-protein interaction network, and functional enrichment analysis. Gradually upregulated expression patterns of Rac2, Itgb2, and Tyrobp genes at 0, 3, 7, and 14 days after SCI were confirmed based on GSE45006 and GSE20907 data set. Finally, we found that Rac2, Itgb2, and Tyrobp genes might take crucial roles in SCI development through the “natural killer cell–mediated cytotoxicity” pathway.     

Identification of biomarkers correlated with hypertrophic cardiomyopathy with co-expression analysis

Hypertrophic cardiomyopathy (HCM) is reported to be the most common genetic heart disease. To identify key module and candidate biomarkers correlated with clinical prognosis of patients with HCM, we carried out this study with co-expression analysis. To construct a co-expression network of hub genes correlated with HCM, the Weighted Gene Co-expression Network Analysis (WGCNA) was performed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by Database for Annotation, Visualization and Integrated Discovery (DAVID). The protein-protein interaction network analysis of central genes was performed to recognize the interactions of central genes. Gene set enrichment analyses were carried out to discover the possible mechanisms involved in the pathways promoted by hub genes. To validate the hub genes, quantitative real-time polymerase chain reaction (RT-PCR) was performed. Based on the results of topological overlap measure based clustering, 2,351 differentially expressed genes (DEGs) were identified. Those genes were included in six different modules. Of these modules, the yellow and the blue modules showed a pivotal correlation with HCM. DEGs were enriched in immune system procedure associated GO terms and KEGG pathways. We identified nine hub genes (TYROBP, STAT3, CSF1R, ITGAM, SYK, ITGB2, LILRB2, LYN, and HCK) affected the immune system significantly. Among the genes we validated with RT-PCR, TYROBP, CSF1R, and SYK showed significant increasing expression levels in model HCM rats. In conclusion, we identified two modules and nine hub genes, which were prominently associated with HCM. We found that immune system may play a crucial role in the HCM. Accordingly, those genes and pathways might become therapeutic targets with clinical usefulness in the future.     

Multi-omic analyses of hepatocellular carcinoma to determine immunological characteristics and key nodes in gene-expression network

Hepatocellular carcinoma (HCC) is a common malignant tumor worldwide, but effective immunotherapy is still limited for those affected. Therefore, there is an urgent need to explore the specific mechanisms governing tumor immunity to improve the survival rate for those diagnosed with HCC. In the present study, we performed a new immune stratification of HCC samples into two subclasses (A and B) from The Cancer Genome Atlas and the International Cancer Genome Consortium databases, and comprehensive multi-omic analyses of major histocompatibility complex genes, gene copy-number variations, somatic mutations, DNA methylation, and non-coding RNAs. Subclass A was found to have a higher survival rate compared with subclass B, and there were significant immunological differences between the two clusters. Based on these differences, we identified DRD1 and MYCN as key hub genes in the immune-phenotype gene expression regulatory network. These results provide novel ideas and evidence for HCC regulatory mechanisms that may improve immunotherapy for this cancer.     

Development and validation of a 14-gene signature for prognosis prediction in hepatocellular carcinoma

Worldwide, hepatocellular carcinoma (HCC) remains a crucial medical problem. Precise and concise prognostic models are urgently needed because of the intricate gene variations among liver cancer cells. We conducted this study to identify a prognostic gene signature with biological significance. We applied two algorithms to generate differentially expressed genes (DEGs) between HCC and normal specimens in The Cancer Genome Atlas cohort (training set included) and performed enrichment analyses to expound on their biological significance. A protein-protein interactions network was established based on the STRING online tool. We then used Cytoscape to screen hub genes in crucial modules. A multigene signature was constructed by Cox regression analysis of hub genes to stratify the prognoses of HCC patients in the training set. The prognostic value of the multigene signature was externally validated in two other sets from Gene Expression Omnibus (GSE14520 and GSE76427), and its role in recurrence prediction was also investigated. A total of 2000 DEGs were obtained, including 1542 upregulated genes and 458 downregulated genes. Subsequently, we constructed a 14-gene signature on the basis of 56 hub genes, which was a good predictor of overall survival. The prognostic signature could be replicated in GSE14520 and GSE76427. Moreover, the 14-gene signature could be applied for recurrence prediction in the training set and GSE14520. In summary, the 14-gene signature extracted from hub genes was involved in some of the HCC-related signalling pathways; it not only served as a predictive signature for HCC outcome but could also be used to predict HCC recurrence.     

Six genes as potential diagnosis and prognosis biomarkers for hepatocellular carcinoma through data mining

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors and the third of cancer mortality worldwide. Although the study of HCC has made great progress, the molecular mechanism and signal pathways of HCC are not yet clear. Therefore, it is necessary to investigate the early diagnosis and prognosis biomarkers for HCC. The aim of this study is to screen the relevant genes and study the association of gene expression with the survival status of HCC patients using bioinformatics approaches, in the hope of establishing marker genes for diagnosis and prognosis of HCC. The gene expression data and corresponding clinical information of HCC samples were downloaded from the The Cancer Genome Atlas database. We performed to study the relationship between gene expression and prognosis of HCC and screen significantly relevant genes associated with prognosis of HCC by analyzing survival and function enrichment of genes. In this study, we collected 421 samples with gene expression data, including 371 tumor samples and 50 normal samples. By using single factor Cox regression analysis, we screened 1,197 genes significantly associated with survival time in the modeling data containing 117 samples and also searched six genes as the best markers to predict living status of HCC patients. Besides, we established score system of survival risk of HCC. Our study recognized six genes (PGBD3, PGM5P3-AS1, RNF5, UTP11, BAG6, and KCND2) to be significantly associated with diagnosis and prognosis of HCC, providing novel targets for studying potential mechanism about the progression of HCC.     

lncRNA A1BG-AS1 suppresses proliferation and invasion of hepatocellular carcinoma cells by targeting miR-216a-5p

Extensive evidence indicate that long noncoding RNAs (lncRNAs) regulates the tumorigenesis and progression of hepatocellular carcinoma (HCC). However, the expression and biological function of lncRNA A1BG antisense RNA 1 (A1BG-AS1) were poorly known in HCC. Here, we found the underexpression of A1BG-AS1 in HCC via analysis of The Cancer Genome Atlas database. Further analyses confirmed that A1BG-AS1 expression in HCC was markedly lower than that in noncancerous tissues based on our HCC cohort. Clinical association analysis revealed that low A1BG-AS1 expression correlated with poor prognostic features, such as microvascular invasion, high tumor grade, and advanced tumor stage. Follow-up data indicated that low A1BG-AS1 level evidently correlated with poor clinical outcomes of HCC patients. Moreover, forced expression of A1BG-AS1 repressed proliferation, migration, and invasion of HCC cells in vitro. Conversely, A1BG-AS1 knockdown promoted these malignant behaviors in HepG2 cells. Mechanistically, A1BG-AS1 functioned as a competing endogenous RNA by directly sponging miR-216a-5p in HCC cells. Notably, miR-216a-5p restoration rescued A1BG-AS1 attenuated proliferation, migration and invasion of HCCLM3 cells. A1BG-AS1 positively regulated the levels of phosphatase and tensin homolog and SMAD family member 7, which were reduced by miR-216a-5p in HCC cells. Altogether, we conclude that A1BG-AS1 exerts a tumor suppressive role in HCC progression.     

Rat Hepatocytes Weighted Gene Co-Expression Network Analysis Identifies Specific Modules and Hub Genes Related to Liver Regeneration after Partial Hepatectomy

The recovery of liver mass is mainly mediated by proliferation of hepatocytes after 2/3 partial hepatectomy (PH) in rats. Studying the gene expression profiles of hepatocytes after 2/3 PH will be helpful to investigate the molecular mechanisms of liver regeneration (LR). We report here the first application of weighted gene co-expression network analysis (WGCNA) to analyze the biological implications of gene expression changes associated with LR. WGCNA identifies 12 specific gene modules and some hub genes from hepatocytes genome-scale microarray data in rat LR. The results suggest that upregulated MCM5 may promote hepatocytes proliferation during LR; BCL3 may play an important role by activating or inhibiting NF-kB pathway; MAPK9 may play a permissible role in DNA replication by p38 MAPK inactivation in hepatocytes proliferation stage. Thus, WGCNA can provide novel insight into understanding the molecular mechanisms of LR.     

Key modules and hub genes identified by coexpression network analysis for revealing novel biomarkers for larynx squamous cell carcinoma

Larynx squamous cell carcinoma (LSCC) is the second most aggressive head and neck squamous cell carcinoma. Numerous genes have been identified to be aberrantly expressed during the development of LSCC. However, currently, researchers focus more on the individual molecule and downstream genes, leaving the coexpression among genes and key upstream disease driver genes unexploited. In this study, we applied weighted gene coexpression analysis (WGCNA) to decipher potential hub genes driving the development of LSCC. After downloading of LSCC microarray profile from gene expression omnibus, different expression analysis was performed, which was used to conduct functional enrichment analysis. Then, we applied WGCNA to highlight the hub genes which were relevant to the carcinogenesis and progression. A total of 2858 differentially expressed genes were identified in LSCC samples compared with adjacent non-neoplastic tissues. WGCNA revealed three LSCC set-specific modules having significant Kyoto Encyclopedia of Genes and Genomes enrichment effect, including pink, cyan, and black module. Nine hub genes were identified to be crucial in LSCC onset and progression, which may assist clinical decisions and serve as potential targets for LSCC treatment.     

埃及伊蚊不同组织的基因共表达模式分析

【目的】利用权重基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)探索埃及伊蚊Aedes aegypti不同组织基因共表达模式。【方法】从NCBI SRA数据库中选择埃及伊蚊不同组织的转录组数据中具代表性的9种组织(雌雄成蚊的触角和脑,雌蚊的喙、下颚须和卵巢,雄成蚊的前足、中足、后足和腹部末端)的双端测序数据;经过缺失值移除以及方差计算后,筛选出方差最大的5 000个基因,利用R软件中WGCNA包建立埃及伊蚊成蚊不同组织的基因共表达网络并划分模块;然后利用clusterProfiler包对组织特异性模块内的基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopediaof Genes and Genomes)富集分析,并用Cytoscape软件中的CytoHubba插件筛选共表达模块内的hub基因。【结果】从埃及伊蚊成蚊不同组织中共鉴定出11个基因共表达模块,在雌蚊触角、喙、卵巢、下颚须以及雄蚊脑、腹部末端组织中各鉴定出1个特异性表达模块,雄蚊前足、中足和后足组织中无特异性表达模块。6个组织特异性表达模块内基因功能注释到组织生物学功能;其中,雌蚊触角特异性green模块内基因具有气味结合和嗅觉受体活性等功能;雌蚊喙特异性purple模块内基因具有丝氨酸型肽链内切酶活性和丝氨酸水解酶活性等功能;雄蚊脑特异性blue模块内基因在生物学过程调节、信号转导和神经系统过程等生物学过程中发挥主要作用。利用CytoHubba进一步鉴定出所选组织特异性共表达模块中具有高连通性的hub基因,包括AAEL010426, AAEL002896, AAEL002600, AAEL000961, AAEL007784和AAEL006429。【结论】本研究依据埃及伊蚊不同组织转录组数据,利用WGCNA方法发现了许多重要的基因共表达模块。本研究的结果为蚊虫基因共表达模式分析提供新思路和方法基础,对探究蚊虫不同组织特有的基因资源信息以及功能基因生物信息学研究有参考价值。     

Identification of castration‐resistant prostate cancer‐related hub genes using weighted gene co‐expression network analysis

Prostate cancer is the most common malignancy in urinary system and brings heavy burdens in men. We downloaded gene expression profile of mRNA and related clinical data of GSE70768 data set from public database. Weighted gene co‐expression network analysis (WGCNA) was used to identify the relationships between gene modules and clinical features, as well as the candidate genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were developed to investigate the potential functions of related hub genes. Importantly, basic experiments were performed to verify the relationship between hub genes and the phenotype previously identified. Lastly, copy number variation (CNV) analysis was conducted to explore the genetical alteration. WGCNA identified that black module was the most relevant module which was tightly related to castration‐resistant prostate cancer (CRPC) phenotype. KEGG and GO analysis results revealed genes in black module were mainly related to RNA splicing. Additionally, 9 genes were chosen as hub genes and heterogeneous nuclear ribonucleoprotein A2/B1 (HNRNPA2B1), golgin A8 family member B (GOLGA8B) and mitogen‐activated protein kinase 8 interacting protein 3 (MAPK8IP3) were identified to be associated with PCa progression and prognosis. Moreover, all above three genes were highly expressed in CRPC‐like cells and their suppression led to hindered cell proliferation in vitro. Finally, CNV analysis found that amplification was the main type of alteration of the 3 hub genes. Our study found that HNRNPA2B1, GOLGA8B and MAPK8IP3 were identified to be tightly associated with tumour progression and prognosis, and further researches are needed before clinical application.