Liberty of your own prognostic gene trademark off their systematic variables from inside the TCGA
Study population
In expose data, we checked and you will installed mRNA phrase chip investigation from HCC structures about vederlo GEO databases utilising the keywords away from “hepatocellular carcinoma” and you can “Homo sapiens”. Half a dozen microarray datasets (GSE121248, GSE84402, GSE65372, GSE51401, GSE45267 and you will GSE14520 (according to the GPL571 program) had been obtained getting DEGs research. Details of brand new GEO datasets used in this study are provided from inside the Dining table step 1. RNA-sequencing analysis of 371 HCC frameworks and 50 normal frameworks normalized of the log2 sales was received on Cancer tumors Genome Atlas (TCGA) getting evaluating the fresh incorporated DEGs regarding half dozen GEO datasets and strengthening gene prognostic habits. GSE14520 datasets (in accordance with the GPL3921 system) provided 216 HCC structures having done health-related suggestions and you will mRNA phrase analysis to have outside recognition of the prognostic gene signature. After leaving out TCGA cases with incomplete logical advice, 233 HCC people making use of their done many years, intercourse, gender, tumor amount, Western Combined Panel for the Cancer (AJCC) pathologic tumor phase, vascular invasion, Operating system status and you will date information had been integrated to possess univariable and multivariable Cox regression studies. Mutation research have been obtained from the newest cBioPortal for Malignant tumors Genomics .
Operating regarding gene expression data
To integrated gene expression chip data downloaded from the GEO datasets, we firstly conducted background correction, quartile normalization for the raw data followed by log2 transformation to obtain normally distributed expression values. The DEGs between HCC tissues and non-tumor tissues were identified using the “Limma” package in R . The thresholds of absolute value of the log2 fold change (logFC) > 1 and adjusted P value < 0.05 were adopted. Mean expression values were applied for genes with multiprobes. Then, we used the robust rank aggregation (RRA) method to finally identify overlapping DEGs (P < 0.05) from the six GEO datasets.
Design out of a potential prognostic trademark
To identify the prognostic genes, we firstly sifted 341 patients from the TCGA Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort with follow-up times of more than 30 days. Then, univariable Cox regression survival analysis was performed based on the overlapping DEGs. A value of P < 0.01 in the univariable Cox regression analysis was considered statistically significant. Subsequently, the prognostic gene signature was constructed by Lasso?penalized Cox regression analysis , and the optimal values of the penalty parameter alpha were determined through 10-times cross-validations by using R package “glmnet” . Based on the optimal alpha value, a twelve-gene prognostic signature with corresponding coefficients was selected, and a risk score was calculated for each TCGA-LIHC patient. Next, the HCC patients were divided into two or three groups based on the optimal cutoff of the risk score determined by “survminer” package in R and X-Tile software. To assess the performance of the twelve-gene prognostic signature, the Kaplan–Meier estimator curves and the C-index comparing the predicted and observed OS were calculated using package “survival” in R. Time-dependent receiver operating characteristic (ROC) curve analysis was also conducted by using the R packages “pROC” and “survivalROC” . Then, the GSE14520 datasets with complete clinical information was used to validate the prognostic performance of twelve-gene signature. The GSE14520 external validation datasets was based on the GPL3921 platform of the Affymetrix HT Human Genome U133A Array Plate Set (HT_HG-U133A, Affymetrix, Santa Clara, CA, United States).
The risk score and other clinical variants, including age, body mass index (BMI), sex, tumor grade, the AJCC pathologic tumor stage, vascular invasion, residual tumor status and AFP value, were analyzed by univariable Cox regression analysis. Next, we conducted a multivariable Cox regression model that combined the risk score and the above clinical indicators (P value < 0.2) to assess the predictive performance. The univariable and multivariable Cox regression analysis were performed with TCGA-LIHC patients (n = 234) that had complete clinical information.
0 دیدگاه