Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.     

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.     

SLC25A38 as a novel biomarker for metastasis and clinical outcome in uveal melanoma

Risk of metastasis is increased by the presence of chromosome 3 monosomy in uveal melanoma (UM). This study aimed to identify more accurate biomarker for risk of metastasis in UM. A total of 80 patients with UM from TCGA were assigned to two groups based on the metastatic status, and bioinformatic analyses were performed to search for critical genes for risk of metastasis. SLC25A38, located on chromosome 3, was the dominant downregulated gene in metastatic UM patients. Low expression of SLC25A38 was an independent predictive and prognostic factor in UM. The predictive potential of SLC25A38 expression was superior to that of pervious reported biomarkers in both TCGA cohort and {"type":"entrez-geo","attrs":{"text":"GSE22138","term_id":"22138"}}GSE22138 cohort. Subsequently, its role in promoting metastasis was explored in vitro and in vivo. Knock-out of SLC25A38 could enhance the migration ability of UM cells, and promote distant metastasis in mice models. Through the inhibition of CBP/HIF-mediated pathway followed by the suppression of pro-angiogenic factors, SLC25A38 was situated upstream of metastasis-related pathways, especially angiogenesis. Low expression of SLC25A38 promotes angiogenesis and metastasis, and identifies increased metastatic risk and worse survival in UM patients. This finding may further improve the accuracy of prognostic prediction for UM.Subject terms: Eye cancer, Prognostic markers     

Identification of core genes and outcomes in hepatocellular carcinoma by bioinformatics analysis

Hepatocellular carcinoma (HCC) is the most common malignant liver disease in the world. However, the mechanistic relationships among various genes and signaling pathways are still largely unclear. In this study, we aimed to elucidate potential core candidate genes and pathways in HCC. The expression profiles GSE14520, GSE25097, GSE29721, and GSE62232, which cover 606 tumor and 550 nontumour samples, were downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, HCC RNA-seq datasets were also downloaded from the Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were filtered using R software, and we performed gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis using the online databases DAVID 6.8 and KOBAS 3.0. Furthermore, the protein-protein interaction (PPI) network complex of these DEGs was constructed by Cytoscape software, the molecular complex detection (MCODE) plug-in and the online database STRING. First, a total of 173 DEGs (41 upregulated and 132 downregulated) were identified that were aberrantly expressed in both the GEO and TCGA datasets. Second, GO analysis revealed that most of the DEGs were significantly enriched in extracellular exosomes, cytosol, extracellular region, and extracellular space. Signaling pathway analysis indicated that the DEGs had common pathways in metabolism-related pathways, cell cycle, and biological oxidations. Third, 146 nodes were identified from the DEG PPI network complex, and two important modules with a high degree were detected using the MCODE plug-in. In addition, 10 core genes were identified, TOP2A, NDC80, FOXM1, HMMR, KNTC1, PTTG1, FEN1, RFC4, SMC4, and PRC1. Finally, Kaplan-Meier analysis of overall survival and correlation analysis were applied to these genes. The abovementioned findings indicate that the identified core genes and pathways in this bioinformatics analysis could significantly enrich our understanding of the development and recurrence of HCC; furthermore, these candidate genes and pathways could be therapeutic targets for HCC treatment.     

DNA methylation biomarkers for head and neck squamous cell carcinoma

DNA methylation plays an important role in the etiology and pathogenesis of head and neck squamous cell carcinoma (HNSCC). The current study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) by a comprehensive bioinformatics analysis. In addition, we screened for DEGs affected by DNA methylation modification and further investigated their prognostic values for HNSCC. We included microarray data of DNA methylation (GSE25093 and GSE33202) and gene expression (GSE23036 and GSE58911) from Gene Expression Omnibus. Aberrantly methylated-DEGs were analyzed with R software. The Cancer Genome Atlas (TCGA) RNA sequencing and DNA methylation (Illumina HumanMethylation450) databases were utilized for validation. In total, 27 aberrantly methylated genes accompanied by altered expression were identified. After confirmation by The Cancer Genome Atlas (TCGA) database, 2 hypermethylated-low-expression genes (FAM135B and ZNF610) and 2 hypomethylated-high-expression genes (HOXA9 and DCC) were identified. A receiver operating characteristic (ROC) curve confirmed the diagnostic value of these four methylated genes for HNSCC. Multivariate Cox proportional hazards analysis showed that FAM135B methylation was a favorable independent prognostic biomarker for overall survival of HNSCC patients.     

Gene regulation network analysis reveals core genes associated with survival in glioblastoma multiforme

Glioblastoma multiforme (GBM) is a very serious mortality of central nervous system cancer. The microarray data from GSE2223 , GSE4058 , GSE4290 , GSE13276 , GSE68848 and GSE70231 (389 GBM tumour and 67 normal tissues) and the RNA‐seq data from TCGA‐GBM dataset (169 GBM and five normal samples) were chosen to find differentially expressed genes (DEGs). RRA (Robust rank aggregation) method was used to integrate seven datasets and calculate 133 DEGs (82 up‐regulated and 51 down‐regulated genes). Subsequently, through the PPI (protein‐protein interaction) network and MCODE/ cytoHubba methods, we finally filtered out ten hub genes, including FOXM1, CDK4, TOP2A, RRM2, MYBL2, MCM2, CDC20, CCNB2, MYC and EZH2, from the whole network. Functional enrichment analyses of DEGs were conducted to show that these hub genes were enriched in various cancer‐related functions and pathways significantly. We also selected CCNB2, CDC20 and MYBL2 as core biomarkers, and further validated them in CGGA, HPA and CCLE database, suggesting that these three core hub genes may be involved in the origin of GBM. All these potential biomarkers for GBM might be helpful for illustrating the important role of molecular mechanisms of tumorigenesis in the diagnosis, prognosis and targeted therapy of GBM cancer.     

人肝细胞癌预后不良相关基因的生物信息学分析及其临床意义

目的:筛选肝细胞癌(HCC)预后不良相关基因,并探讨其临床意义。方法:在基因表达综合数据库(GEO)中获取符合分析条件的肝细胞癌全基因组表达谱数据并分析得到差异表达基因(DEGs),再运用生物学信息注释及可视化数据库 (DAVID) 和蛋白相互作用数据库 (String) 分别进行功能富集分析和蛋白质互作用网络的构建。利用癌症基因组图谱数据库(TCGA)和Cox比例风险回归模型对相关差异基因进行预后分析。结果:找到一个符合条件的人类HCC数据库 (GSE84402),共筛选出1141个差异表达基因(DEGs),其中上调基因720个,下调基因421个。基因功能富集分析和蛋白质互作用分析结果显示CDK1、CDC6、CCNA2、CHEK1、CENPE 、PIK3R1、RACGAP1、BIRC5、KIF11和CYP2B6为HCC预后的关键基因。TCGA数据库和Cox回归模型分析显示CDC6、PIK3R1、RACGAP1和KIF11的表达升高,CENPE的表达降低与HCC预后不良密切相关。结论:CDC6、CENPE、PIK3R1、RACGAP1和KIF11可能和HCC的预后不良相关,可作为未来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.     

The AP2M1 gene expression is a promising biomarker for predicting survival of patients with hepatocellular carcinoma

There is a growing need for the discovery of new prognostic factors for cases where the scoring and staging system of hepatocellular carcinoma (HCC) does not result in a clear definition. We analyzed whether AP-2 complex subunit mu (AP2M1) expression could be a new prognostic marker for HCC based on the roles of AP2M1 in influencing hepatocyte growth factor (HGF) promoter regulation and hepatitis C virus (HCV) assembly. Patient data were extracted from cohorts of the Gene Expression Omnibus (GSE10186), International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Differential expression value between matched cancer and normal liver was identified using ICGC cohort. Subsequently, we compared AP2M1 expression as a prognostic gene with other well-known prognostic genes for HCC, using the time-dependent area under the curve (AUC) of the Uno's C-index, the AUC value of the receiver operating characteristics at 5 years, Kaplan-Meier survival curve, and multivariate analysis. Particularly, TCGA and GSE10186 patients were divided into subgroups based on alcohol intake, hepatitis B, and C viral infections, and analyzed in the same methods. The AP2M1 expression values in patients with cancer were much higher than matched normal liver. The AP2M1 level showed excellent prognosis predictions in comparison with existing markers in the three independent cohorts (n = 647). In particular, it was more predictive of prognosis than other markers in alcohol intake and HCV infections. In conclusion, we were confident that AP2M1 provides sufficient value as a new prognostic marker for HCC especially patients with HCV infection and/or alcohol intake.     

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.     

MiR-139-5p/SLC7A11 inhibits the proliferation,invasion and metastasis of pancreatic carcinoma via PI3K/Akt signaling pathway

ObjectivePancreatic carcinoma (PANC) is one of the important aggressive cancers, with deficiency in effective therapeutics. The study aimed to investigate the effects and molecular mechanism of miR-139-5p/SLC7A11 on the proliferation and metastasis of PANC.MethodsBioinformatics was used to analyze the differentially expressed genes in the TCGA database. PANC cell lines with overexpressed miR-139-5p and Solute Carrier Family 7, Member 11 (SLC7A11) was established, and have been used to detect cell proliferation, invasion and metastasis of PANC Subsequently, bioinformatic analysis and dual luciferase reporter assay were performed to confirm that SLC7A11 was a target gene of miR-139-5p. Xenograft mice model was used to explore the functions of miR-139-5p in PANC tumorigenicity.ResultsMiR-139-5p could regulate and affect the protein expression of P13K and Akt associated with phosphatidylinositol signaling pathway by inhibiting SLC7A11. MiR-139-5p was found to be lowly expressed in PANC tissues, while SLC7A11 was highly expressed. Low expression of miR-139-5p and high expression of SLC7A11 were positively associated with poor clinical outcomes. PANC cell proliferation, invasion and metastasis could be inhibited by miR-139-5p overexpression and be promoted by SLC7A11 overexpression. MiR-139-5p overexpression could suppress PANC tumor growth and the expressions of SLC7A11, p-PI3K, p-Akt in tumor tissues. Therefore, the inhibitory of miR-139-5p to PANC cell proliferation, invasion and metastasis was partly due to its inhibiting effect on SLC7A11 expression.ConclusionOur study proves that miR-139-5p/SLC7A11 has important functions on PANC, suggesting that miR-139-5p can be used as a biomarker for PANC patients.     

A ferroptosis-related gene signature for overall survival prediction and immune infiltration in lung squamous cell carcinoma

Background: Ferroptosis is associated with cancer initiation and progression. However, the molecular mechanism and prognostic value of ferroptosis-related genes in lung squamous cell carcinoma (LUSC) are poorly understood.Methods: The mRNA expression profiles, methylation data, and clinical information of patients with LUSC were downloaded from TCGA and GEO database. Ferroptosis-related differentially expressed genes (DEGs) were identified between cancerous and non-cancerous tissues, and their prognostic value was systemically investigated by bioinformatic analyses.Results: A ferroptosis-related gene signature (ALOX5, TFRC, PHKG2, FADS2, NOX1) was constructed using multivariate Cox regression analysis and represented as a risk score. Overall survival (OS) probability was significantly lower in the high-risk group than in the low-risk group (P<0.001), and receiver operating characteristic curve showed a good predictive capacity (AUC = 0.739). The risk score was an independent prognostic factor for LUSC. A nomogram was constructed to predict the OS probabilities at 1, 3, and 5 years. High-risk score was associated with increased immune infiltration, lower methylation levels, higher immune checkpoint genes expression levels, and better chemotherapy response. Cell adhesion molecules, focal adhesion, and extracellular matrix receptor interaction were the main pathways in the high-risk group. The signature was validated using the TCGA test cohort, entire TCGA cohort, {"type":"entrez-geo","attrs":{"text":"GSE30219","term_id":"30219"}}GSE30219, {"type":"entrez-geo","attrs":{"text":"GSE157010","term_id":"157010"}}GSE157010, {"type":"entrez-geo","attrs":{"text":"GSE73403","term_id":"73403"}}GSE73403, and {"type":"entrez-geo","attrs":{"text":"GSE4573","term_id":"4573"}}GSE4573 datasets. The gene disorders in patients with LUSC were validated using real-time PCR and single-cell RNA sequencing analysis.Conclusions: A ferroptosis-related gene signature was constructed to predict OS probability in LUSC. This could facilitate novel therapeutic methods and guide individualized therapy.     

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.     

Beta-adrenergic pathway activation enhances aggressiveness and inhibits stemness in head and neck cancer

Chronic stress leads to the activation of the beta-adrenergic pathway. Its activation has been implicated in the progression of different types of cancer but its role on head and neck squamous cell carcinomas (HNSCCs) remains undefined. The aim of this study was to investigate the influence of the beta-adrenergic pathway activation in the progression of HNSCCs and offer a panel of potential treatments for patients with the active beta-adrenergic pathway. Five hundred and twenty TCGA patients with primary HNSCCs were divided in two groups: ADRB2^low / SLC6A2^low and ADRB2^high / SLC6A2^high. Differentially expressed genes (DEGs) were identified through differential expression analysis. The association of clinicopathological and genomic features between the groups was analyzed using a bioinformatic approach. Potential drugs for treatment of HNSCC were identified based on the DEGs. There was association between ADRB2 and SLC6A2 expressions with age, race, tumor site, histologic grade, perineural invasion, and HPV p16 status. It was identified 898 DEGs between the groups. High ADRB2/SLC6A2 expression stimulated HNSCC proliferation, adhesion, invasion, and angiogenesis. On the other hand, genes related to cell stemness were downregulated in patients with activation of the beta- adrenergic pathway. Finally, 56 FDA-approved antineoplastic and immunotherapeutic drugs were identified as potential targets for the personalized treatment.     

氨基酸转运载体SLC1A5在胃癌组织中的表达及其临床病理学意义

目的:研究氨基酸转运载体溶质载体家族1成员5(Solute Carrier Family 1 Member 5, SLC1A5)蛋白在胃癌组织中的表达情况,并探讨其与胃癌临床病理特征及预后的相关性。方法:收集进展期胃癌组织及对应癌旁组织90例,应用免疫组化技术检测SLC1A5在上述组织中的表达情况,并统计分析其表达与胃癌临床病理特征及预后的关系。同时通过基因数据库分析SLC1A5在胃癌组织和癌旁组织中表达情况及其对胃癌患者预后的影响。结果:与癌旁组织相比,胃癌组织中SLC1A5表达明显上调(P0.0001)。数据库研究也显示SLC1A5在胃癌组织中表达明显上调(GSE 65801,P=0.0046;GSE 63809,P0.0001;GSE 27342,P=0.0147)。胃癌组织中SLC1A5高表达与肿瘤大小(P0.05)、肿瘤浸润深度(P0.01)、淋巴结转移(P0.05)、TNM分期(P0.05)和Ki-67(P0.01)相关,而与年龄、性别、肿瘤位置及分化程度均无显著相关性(P0.05)。胃癌组织中SLC1A5表达强度与患者预后相关,表达越高,患者预后越差(总体生存率,P=0.0131;无复发生存率,P=0.0293)。数据库分析也显示SLC1A5高表达可明显缩短患者的总体生存期(GSE 14210,P=0.011;GSE 22377,P=0.0015)和无进展生存期(GSE 14210,P=0.0095;GSE 22377,P=0.0012)。结论:SLC1A5蛋白表达在胃癌组织中明显上调,且与肿瘤大小、肿瘤浸润深度、淋巴结转移及TNM分期有关。SLC1A5高表达与胃癌患者预后不良密切相关。     

Identification of hub genes and pathways in adrenocortical carcinoma by integrated bioinformatic analysis

Adrenocortical carcinoma (ACC), a rare malignant neoplasm originating from adrenal cortical cells, has high malignancy and few treatments. Therefore, it is necessary to explore the molecular mechanism of tumorigenesis, screen and verify potential biomarkers, which will provide new clues for the treatment and diagnosis of ACC. In this paper, three gene expression profiles (GSE10927, GSE12368 and GSE90713) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained using the Limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID. Protein‐protein interaction (PPI) network was evaluated by STRING database, and PPI network was constructed by Cytoscape. Finally, GEPIA was used to validate hub genes’ expression. Compared with normal adrenal tissues, 74 up‐regulated DEGs and 126 down‐regulated DEGs were found in ACC samples; GO analysis showed that up‐regulated DEGs were enriched in organelle fission, nuclear division, spindle, et al, while down‐regulated DEGs were enriched in angiogenesis, proteinaceous extracellular matrix and growth factor activity; KEGG pathway analysis showed that up‐regulated DEGs were significantly enriched in cell cycle, cellular senescence and progesterone‐mediated oocyte maturation; Nine hub genes (CCNB1, CDK1, TOP2A, CCNA2, CDKN3, MAD2L1, RACGAP1, BUB1 and CCNB2) were identified by PPI network; ACC patients with high expression of 9 hub genes were all associated with worse overall survival (OS). These hub genes and pathways might be involved in the tumorigenesis, which will offer the opportunities to develop the new therapeutic targets of ACC.     

Prognostic value of tumour microenvironment-related genes by TCGA database in rectal cancer

Rectal cancer is a common malignant tumour and the progression is highly affected by the tumour microenvironment (TME). This study intended to assess the relationship between TME and prognosis, and explore prognostic genes of rectal cancer. The gene expression profile of rectal cancer was obtained from TCGA and immune/stromal scores were calculated by Estimation of Stromal and Immune cells in Malignant Tumors using Expression data (ESTIMATE) algorithm. The correlation between immune/stromal scores and survival time as well as clinical characteristics were evaluated. Differentially expressed genes (DEGs) were identified according to the stromal/immune scores, and the functional enrichment analyses were conducted to explore functions and pathways of DEGs. The survival analyses were conducted to clarify the DEGs with prognostic value, and the protein-protein interaction (PPI) network was performed to explore the interrelation of prognostic DEGs. Finally, we validated prognostic DEGs using data from the Gene Expression Omnibus (GEO) database by PrognoScan, and we verified these genes at the protein levels using the Human Protein Atlas (HPA) databases. We downloaded gene expression profiles of 83 rectal cancer patients from The Cancer Genome Atlas (TCGA) database. The Kaplan-Meier plot demonstrated that low-immune score was associated with worse clinical outcome (P = .034), metastasis (M1 vs. M0, P = .031) and lymphatic invasion (+ vs. -, P < .001). A total of 540 genes were screened as DEGs with 539 up-regulated genes and 1 down-regulated gene. In addition, 60 DEGs were identified associated with overall survival. Functional enrichment analyses and PPI networks showed that the DEGs are mainly participated in immune process, and cytokine-cytokine receptor interaction. Finally, 19 prognostic genes were verified by GSE17536 and GSE17537 from GEO, and five genes (ADAM23, ARHGAP20, ICOS, IRF4, MMRN1) were significantly different in tumour tissues compared with normal tissues at the protein level. In summary, our study demonstrated the associations between TME and prognosis as well as clinical characteristics of rectal cancer. Moreover, we explored and verified microenvironment-related genes, which may be the potential key prognostic genes of rectal cancer. Further clinical samples and functional studies are needed to validate this finding.     

Identification of featured biomarkers in different types of lung cancer with DNA microarray

Lung cancer is a worldwide leading cause of cancer-related death. The aim of this study was to identify target genes and specific biomarkers for identification and treatment of different types of lung cancer with DNA microarray. Gene expression profile GSE6044 and miRNA microarray profile GSE17681 were downloaded from Gene Expression Omnibus database. The differentially expressed genes (DEGs) and miRNAs were screened with multtest package in R language. Then, functional enrichment analysis of identified DEGs was performed. Furthermore, the verified target genes based on screened miRNAs were selected from miRTarBase and miRecords databases. Then miRNA-target gene regulation network was constructed. APOE, CDC6 and ATP2B1were involved in most of the functions obtained for adenocarcinomas, small cell lung cancer and squamous cell carcinomas, respectively. The target DEGs of differentially expressed hsa-miR-29a included FGG in adenocarcinoma, RAN and COL4A1 in small cell lung cancer, GLUL in squamous cell carcinoma. The target DEGs of has-miR-7 were SNCA and SLC7A5 in adenocarcinoma and small cell lung cancer, respectively. ICAM1 and KIT were the target DEGs of hsa-miR-222 in adenocarcinoma and squamous cell carcinoma. The miRNAs and their differentially expressed target genes have the potential to be used in clinic for diagnosis and treatment of different kinds of lung cancer in the future.     

COPD 差异表达基因的生物信息学分析及在LUSC样本中的表达

为确定慢性阻塞性肺病(COPD)的分子标记物及COPD与肺鳞状细胞癌(LUSC)共存的差异表达基因,探寻COPD合并肺癌的预测因子,发现新的治疗靶点。本研究采用生物信息学方法,从GEO数据库中筛选3套基因芯片数据集,挖掘COPD患者小气道上皮细胞(SAEC)的差异表达基因(DEG)以及潜在的生物标记物,并通过基因本体(GO)、京都基因与基因组百科全书(KEGG)富集分析预测DEGs的功能及参与的代谢途径。继而对DEGs构建PPI网络,使用Cytoscape软件筛选子模块和Hub基因,并将Hub基因通过TCGA数据库分析其在LUSC中的差异表达情况及差异基因间的相关性。结果共获得52个上调基因和24个下调基因,代谢通路主要集中在细胞色素P450对外源物质的代谢、化学致癌、花生四烯酸代谢及甲状腺激素合成四条途径上,通过Cytoscape软件从PPI网络中筛选得到2个功能模块和10个Hub基因,进一步验证发现其中5个基因在TCGA数据库中的LUSC样本中同样差异表达。由此推测SPP1、ALDH3A1、SPRR3、KRT6A和SPRR1B 可能为COPD 分子标记物及COPD与LUSC共存的DEGs,从而为研究COPD和LUSC的发病机制及二者潜在关系奠定良好的基础。