Identification of Ten Core Hub Genes as Potential Biomarkers and Treatment Target for Hepatoblastoma.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:33868991", "kind": "science-record", "revision": 8, "content_hash": "sha256:aa90cd6850872cbdb1ee528da7ad3766575dcea29a3c86669a2d2724b9add46a", "created_at": "2026-11-01T15:41:32Z", "updated_at": "2026-11-01T22:46:32Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Identification of Ten Core Hub Genes as Potential Biomarkers and Treatment Target for Hepatoblastoma.", "doi": "10.3389/fonc.2021.591507", "pmid": "33868991", "authors": ["Rui Sun", "Simin Li", "Ke Zhao", "Mei Diao", "Long Li"], "journal": "Frontiers in oncology", "year": 2021, "abstract": "BACKGROUND: This study aimed to systematically investigate gene signatures for hepatoblastoma (HB) and identify potential biomarkers for its diagnosis and treatment.\nMATERIALS AND METHODS: GSE131329 and GSE81928 were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between hepatoblastoma and normal samples were identified using the Limma package in R. Then, the similarity of network traits between two sets of genes was analyzed by weighted gene correlation network analysis (WGCNA). Cytoscape was used to visualize and select hub genes. PPI network of hub genes was construed by Cytoscape. GO enrichment and KEGG pathway analyses of hub genes were carried out using ClueGO. The random forest classifier was constructed based on the hub genes using the GSE131329 dataset as the training set, and its reliability was validated using the GSE81928 dataset. The resulting core hub genes were combined with the InnateDB database to identify the innate core genes.\nRESULTS: A total of 4244 DEGs in HB were identified. WGCNA identified four modules that were significantly correlated with the disease status. A total of 114 hub genes were obtained within the top 20 genes of each node rank. 6982 relation pairs and 3700 nodes were contained in the PPI network of 114 hub genes. GO enrichment and KEGG pathway analyses of hub genes were focused on MAPK, cell cycle, p53, and other crucial pathways involved in HB. A random forest classifier was constructed using the 114 hub genes as feature genes, resulting in a 95.5% true positive rate when classifying HB and normal samples. A total of 35 core hub genes were obtained through the mean decrease in accuracy and mean decrease Gini of the random forest model. The classification efficiency of the random forest model was 81.4%. Finally, CDK1, TOP2A, ADRA1A, FANCI, XRCC1, TPX2, CCNB2, CDK4, GLYATL1, and CFHR3 were identified by cross-comparison with the InnateDB database.\nCONCLUSION: Our study established a random forest classifier that identified 10 core genes in HB. These findings may be beneficial for the diagnosis, prediction, and targeted therapy of HB.", "claims": [{"id": "claim-method", "statement": "The Limma package in R was used to identify DEGs between HB and non-tumor samples.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33868991", "locator": "sections[4] \"Identification of Differentially Expressed Genes (DEGs)\" paragraph 2", "quote": "The Limma package in R was used to identify DEGs between HB and non-tumor samples."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material .", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33868991", "locator": "sections[20] \"Data Availability Statement\" paragraph 2", "quote": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material ."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": null, "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:33868991"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "complete", "blocked_on": null, "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE131329"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE131329"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:33868991", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.3389/fonc.2021.591507", "sha256": "sha256:5e1997238f69d61d80ea1762fabd38a25d0063ed92fd95219182682ad68bee97", "bytes": 30538, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.3389/fonc.2021.591507", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC8047669", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE131329", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE131329", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE131329", "version": null, "note": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material .", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material .", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM3770503", "GSM3770504", "GSM3770505", "GSM3770506", "GSM3770507", "GSM3770508", "GSM3770509", "GSM3770510", "GSM3770511", "GSM3770512", "GSM3770513", "GSM3770514", "GSM3770515", "GSM3770516", "GSM3770517", "GSM3770518", "GSM3770519", "GSM3770520", "GSM3770521", "GSM3770522", "GSM3770523", "GSM3770524", "GSM3770525", "GSM3770526", "GSM3770527", "GSM3770528", "GSM3770529", "GSM3770530", "GSM3770531", "GSM3770532", "GSM3770533", "GSM3770534", "GSM3770535", "GSM3770536", "GSM3770537", "GSM3770538", "GSM3770539", "GSM3770540", "GSM3770541", "GSM3770542", "GSM3770543", "GSM3770544", "GSM3770545", "GSM3770546", "GSM3770547", "GSM3770548", "GSM3770549", "GSM3770550", "GSM3770551", "GSM3770552", "GSM3770553", "GSM3770554", "GSM3770555", "GSM3770556", "GSM3770557", "GSM3770558", "GSM3770559", "GSM3770560", "GSM3770561", "GSM3770562", "GSM3770563", "GSM3770564", "GSM3770565", "GSM3770566", "GSM3770567", "GSM3770568", "GSM3770569"], "sample_count": 67, "group_sizes": null, "group_sizes_basis": null, "platform": "6244", "taxon": "Homo sapiens", "title": "Expression data from hepatoblastoma tissues", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE81928", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81928", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE81928", "version": null, "note": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material .", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material .", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM2177840", "GSM2177841", "GSM2177842", "GSM2177843", "GSM2177844", "GSM2177845", "GSM2177846", "GSM2177847", "GSM2177848", "GSM2177849", "GSM2177850", "GSM2177851", "GSM2177852", "GSM2177853", "GSM2177854", "GSM2177855", "GSM2177856", "GSM2177857", "GSM2177858", "GSM2177859", "GSM2177860", "GSM2177861", "GSM2177862", "GSM2177863", "GSM2177864", "GSM2177865", "GSM2177866", "GSM2177867", "GSM2177868", "GSM2177869", "GSM2177870", "GSM2177871"], "sample_count": 32, "group_sizes": null, "group_sizes_basis": null, "platform": "16791", "taxon": "Homo sapiens", "title": "Hepatic transcriptome of pediatric hepatoblastoma.", "resolution_checked_at": "2026-11-01T15:39:55Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "asset:table:pmid:33868991:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:99e221ce5b969fdaddd1969ba777980235f3782f8c821f084cbbe5225224df8a", "bytes": 635, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[3]/table-wrap[1]#T1", "version": null, "note": "The characteristics of the 2 datasets in the study.", "structure": {"pmcid": "PMC8047669", "section_index": 3, "section_title": "Acquisition of Microarray Profiles", "row_count": 3, "column_count": 7}, "table": {"label": "Table 1", "caption": "The characteristics of the 2 datasets in the study.", "rows": [["Datesets", "Country", "Researchers/ References", "Experiment type", "Tumor site", "Sample size (normal/tumor)", "Platform"], ["GSE131329", "Japan", "Contributed by Hiyama E, et al.", "Expression profiling by array", "hepatoblastoma", "67 (14/53)", "GPL6244 [HuGene-1_0-st] Affymetrix Human Gene 1.0 ST Array"], ["GSE81928", "USA", "(10)", "Expression profiling by high throughput sequencing", "hepatoblastoma", "26 (3/23)", "GPL16791 Illumina HiSeq 2500"]], "row_count": 3, "column_count": 7, "numeric_cells": 7, "header_rows": 1, "locator": "section[3]/table-wrap[1]#T1"}}, {"id": "asset:table:pmid:33868991:2", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:ff6f874c6c7cd749fc2ba563f625f8d61954e39f7d200be552b9d543513e2d51", "bytes": 2078, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[14]/table-wrap[1]#T2", "version": null, "note": "Network topological characteristic of top 20 nodes in PPI network.", "structure": {"pmcid": "PMC8047669", "section_index": 14, "section_title": "PPI Network Construction of Hub Genes", "row_count": 20, "column_count": 9}, "table": {"label": "Table 2", "caption": "Network topological characteristic of top 20 nodes in PPI network.", "rows": [["Gene", "label", "Degree", "Average ShortestPath Length", "BetweennessCentrality", "ClosenessCentrality", "ClusteringCoefficient", "Stress", "TopologicalCoefficient"], ["RPS2", "Turquoise", "304", "2.796578", "0.088088", "0.35758", "0.002084", "35651276", "0.011482"], ["PPP2R1A", "Turquoise", "285", "2.733026", "0.10859", "0.365895", "0.001112", "32731778", "0.008927"], ["CDK1", "Green", "277", "2.606464", "0.086532", "0.383662", "0.007953", "23033244", "0.008038"], ["FBL", "Turquoise", "272", "2.764259", "0.083119", "0.361761", "0.003337", "26783624", "0.009266"], ["PLK1", "Turquoise", "271", "2.749593", "0.087714", "0.36369", "0.002952", "24705604", "0.009257"], ["TRIM28", "Turquoise", "251", "2.847094", "0.07047", "0.351235", "0.000829", "35607200", "0.02181"], ["CDK4", "Turquoise", "215", "2.763987", "0.064648", "0.361796", "0.002478", "17600494", "0.009869"], ["PRMT1", "Turquoise", "174", "2.833786", "0.051374", "0.352885", "0.003654", "15009342", "0.013589"], ["SF3A2", "Brown", "156", "2.933188", "0.048426", "0.340926", "0", "11812770", "0.028122"], ["ITCH", "Brown", "155", "2.893808", "0.049782", "0.345565", "0.000922", "10898838", "0.016011"], ["ANLN", "Green", "149", "2.937534", "0.055011", "0.340422", "0.000998", "7721342", "0.020761"], ["USP15", "Brown", "146", "2.891092", "0.04818", "0.34589", "0.001606", "9324218", "0.015811"], ["CCNB1", "Green", "144", "2.747691", "0.02827", "0.363942", "0.019814", "5878754", "0.013833"], ["EHMT2", "Turquoise", "132", "2.939164", "0.038808", "0.340233", "0", "8115194", "0.02823"], ["CCNA2", "Green", "126", "2.771863", "0.026821", "0.360768", "0.01981", "4666808", "0.015196"], ["USP9X", "Brown", "125", "2.898968", "0.034266", "0.34495", "0.001419", "9220858", "0.020271"], ["HCFC1", "Brown", "120", "2.885117", "0.032388", "0.346606", "0.003641", "7730212", "0.018792"], ["KIF11", "Green", "118", "2.847366", "0.039087", "0.351202", "0.001883", "5902744", "0.013225"], ["TOP2A", "Green", "106", "2.892178", "0.020279", "0.34576", "0.004672", "7719634", "0.024587"]], "row_count": 20, "column_count": 9, "numeric_cells": 136, "header_rows": 1, "locator": "section[14]/table-wrap[1]#T2"}}], "gaps": ["asset:data:geo:GSE131329 not downloaded", "asset:data:geo:GSE131329 licence not read", "asset:data:geo:GSE81928 not downloaded", "asset:data:geo:GSE81928 licence not read"]}, "analysis": {"state": "complete", "methods": [{"id": "method-1", "family": "differential_expression", "specific": "limma_voom", "label": "Differential expression analysis", "dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["count_matrix_with_condition"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "The Limma package in R was used to identify DEGs between HB and non-tumor samples."}}, {"id": "method-2", "family": "machine_learning", "specific": "random_forest", "label": "Machine learning", "dsl_operations": [], "missing_operations": ["machine_learning_random_forest"], "expressible": false, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 0.7, "band": "medium", "quote": "random forest"}}, {"id": "method-3", "family": "network_analysis", "specific": "ppi_network", "label": "Network / pathway analysis", "dsl_operations": ["build_graph", "centrality_scores", "pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.7, "band": "medium", "quote": "PPI"}}, {"id": "method-4", "family": "network_analysis", "specific": "pathway_enrichment", "label": "Network / pathway analysis", "dsl_operations": ["pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.7, "band": "medium", "quote": "GO enrichment"}}, {"id": "method-5", "family": "diagnostic_evaluation", "specific": "roc", "label": "Diagnostic test evaluation", "dsl_operations": ["confusion_counts", "roc_curve", "sensitivity_specificity", "solve_confusion_matrix", "classification_metrics", "published_precision_envelope"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["reference_standard_and_test"], "evidence": {"score": 0.8, "band": "high", "quote": "ROC curve"}}], "targets": [], "gaps": []}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: limma_voom."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R1", "reachable": true, "accessions": ["GSE131329", "GSE81928"], "resolved_accessions": ["GSE131329", "GSE81928"], "statements": ["The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material ."], "basis": "The paper names GSE131329, GSE81928 in Gene Expression Omnibus (GEO) and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "assets": ["source:pmid:33868991", "data:geo:GSE131329", "asset:table:pmid:33868991:1", "asset:table:pmid:33868991:2"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE131329"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-10-31T18:29:08Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 22236, "completion_tokens": 1122, "cached_tokens": 784}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:33868991", "type": "paper", "label": "pmid:33868991"}, {"id": "source:pmid:33868991", "type": "manuscript", "label": "PMID 33868991 full text"}, {"id": "data:geo:GSE131329", "type": "dataset", "label": "GEO GSE131329"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:33868991", "type": "paper", "label": "pmid:33868991"}, {"id": "source:pmid:33868991", "type": "manuscript", "label": "PMID 33868991 full text"}, {"id": "data:geo:GSE131329", "type": "dataset", "label": "GEO GSE131329"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-33868991", "type": "extraction", "used": ["source:pmid:33868991"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE131329", "type": "accession_resolution", "used": ["data:geo:GSE131329"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:33868991", "parent_hash": "sha256:fd61c653790f11f3ed527c8625f13138f3dbef44ec2f8127cfae359aa51a617e", "fork_of": null, "generation": 7}, "gaps": ["asset:data:geo:GSE131329 licence not read", "asset:data:geo:GSE131329 not downloaded", "asset:data:geo:GSE81928 licence not read", "asset:data:geo:GSE81928 not downloaded", "execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "analysis", "definition_id": "differential_expression", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["The Limma package in R was used to identify DEGs between HB and non-tumor samples."], "order_basis": "analysis_only", "id": "assay:1"}]}