Screening and identification of genes related to ferroptosis in keratoconus.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:37626095", "kind": "science-record", "revision": 5, "content_hash": "sha256:e1f16ace478e43a769d0c2d9e38982cba02b651aded9589fcf449ab836391aa9", "created_at": "2026-11-01T19:02:25Z", "updated_at": "2026-11-01T22:52:35Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Screening and identification of genes related to ferroptosis in keratoconus.", "doi": "10.1038/s41598-023-41194-2", "pmid": "37626095", "authors": ["Xiaojun Wu", "Qing Deng", "Zhe Han", "Feixue Ni", "Daxi Sun", "Yuxue Xu"], "journal": "Scientific reports", "year": 2023, "abstract": "Corneal keratoconus (KC) is a dilated (ectatic) corneal disease characterized by a central thinning of the cornea, which causes protrusion into a conical shape that seriously affects vision. However, due to the complex etiology of keratoconus, its entire mechanism remains unclear and there is no mechanism-directed treatment method. Ferroptosis is a novel programmed cell death mechanism related to lipid peroxidation, stress, and amino acid metabolism, which plays a crucial role in various diseases. This study aimed to explore the relationship between keratoconus and ferroptosis, to provide new insights into the mechanism of keratoconus development, and potential treatment options based on further elucidation of this mechanism. The corresponding mRNA microarray expression matrix data of KC patients were obtained from GEO database (GSE204791). Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb. GO and GSEA were performed to analyze differential pathways, ssGSEA was used to determine immune status, and then, feasible drugs were predicted by gene-drug network. Additionally, we predicted the miRNA and IncRNA of hub genes to identify the underlying mechanism of disease so as to predict treatment for the disease. The epithelial transcriptome from keratoconus tissue mRNA microarray data (GSE204791) was extracted for the main analysis, including eight epithelial cells and eight epithelial control cells. The differential genes that were overlapped by WGCAN, SVM-RFE and FRGs were mainly related to oxidative stress, immune regulation, cellular inflammation, and metal ion transport. Through further analysis, aldo-keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages. Then, gene-drug interaction network analysis and miRNA prediction were performed through the website. It was concluded that four immune-related drugs (INDOMETHACIN, DAUNORUBICIN, DOXORUBICIN, DOCETAXEL) and a miRNA (has-miR-184) were screened to predict potential drugs and targets for disease treatment. To our knowledge, this was the first report of KC being associated with ferroptosis and prompted search for differential genes to predict drug targets of gene immunotherapy. Our findings provided insight and a solid basis for the analysis and treatment of KC.", "claims": [{"id": "claim-method", "statement": "Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37626095", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets generated or analyzed during the current study are freely available for the research community to access from the GEO (GSE204791) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE204791).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37626095", "locator": "sections[28] \"notes\" paragraph 1", "quote": "The datasets generated or analyzed during the current study are freely available for the research community to access from the GEO (GSE204791) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE204791)."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Through further analysis, aldo–keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:37626095", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Through further analysis, aldo–keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "Through further analysis, aldo–keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages.", "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:37626095"], "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:GSE204791"], "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:GSE204791"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:37626095", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1038/s41598-023-41194-2", "sha256": "sha256:5c01d70a5ff07e30f92c5a3fd7cf37fd4e99592fd92cb155c035efc0d1dfec4a", "bytes": 30846, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1038/s41598-023-41194-2", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10457308", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE204791", "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=GSE204791", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE204791", "version": null, "note": "The datasets generated or analyzed during the current study are freely available for the research community to access from the GEO (GSE204791) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE204791).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets generated or analyzed during the current study are freely available for the research community to access from the GEO (GSE204791) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE204791).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM6193900", "GSM6193903", "GSM6193906", "GSM6193908", "GSM6193911", "GSM6193913", "GSM6193916", "GSM6193918", "GSM6193920", "GSM6193923", "GSM6193925", "GSM6193928", "GSM6193930", "GSM6193933", "GSM6193935", "GSM6193938", "GSM6193940", "GSM6193943", "GSM6193945", "GSM6193948", "GSM6193950", "GSM6193953", "GSM6193955", "GSM6193957", "GSM6193960", "GSM6193962", "GSM6193965", "GSM6193967", "GSM6193970", "GSM6193972", "GSM6193975"], "sample_count": 31, "group_sizes": null, "group_sizes_basis": null, "platform": "21185", "taxon": "Homo sapiens", "title": "Altered regulation of mRNA and miRNA expression in epithelial and stromal tissue of keratoconus corneas [RNA]", "resolution_checked_at": "2026-11-01T19:09:46Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE204791 not downloaded", "asset:data:geo:GSE204791 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "machine_learning", "specific": "support_vector_machine", "label": "Machine learning", "dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb."}}, {"id": "method-2", "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": 0.8, "band": "high", "quote": "limma"}}, {"id": "method-3", "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": ["method.reader_differs_from_extractor"]}, "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": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: support_vector_machine."}, "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": ["GSE204791"], "resolved_accessions": ["GSE204791"], "statements": ["The datasets generated or analyzed during the current study are freely available for the research community to access from the GEO (GSE204791) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE204791)."], "basis": "The paper names GSE204791 in GEO and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "assets": ["source:pmid:37626095", "data:geo:GSE204791"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE204791"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:51:26Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 9199, "completion_tokens": 2133, "cached_tokens": 264}}, "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:37626095", "type": "paper", "label": "pmid:37626095"}, {"id": "source:pmid:37626095", "type": "manuscript", "label": "PMID 37626095 full text"}, {"id": "data:geo:GSE204791", "type": "dataset", "label": "GEO GSE204791"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:37626095", "type": "paper", "label": "pmid:37626095"}, {"id": "source:pmid:37626095", "type": "manuscript", "label": "PMID 37626095 full text"}, {"id": "data:geo:GSE204791", "type": "dataset", "label": "GEO GSE204791"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-37626095", "type": "extraction", "used": ["source:pmid:37626095"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE204791", "type": "accession_resolution", "used": ["data:geo:GSE204791"], "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:37626095", "parent_hash": "sha256:d5fc182337871f5165ed3d30401f4f74a94aaebcf7de3e7b261eac9efa519f3c", "fork_of": null, "generation": 4}, "gaps": ["asset:data:geo:GSE204791 licence not read", "asset:data:geo:GSE204791 not downloaded", "execution.not_run", "method.reader_differs_from_extractor", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "microarray", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["on of this mechanism. The corresponding mRNA microarray expression matrix data of KC patients were o", "l transcriptome from keratoconus tissue mRNA microarray data (GSE204791) was extracted for the main"], "order_basis": "paper_text_order:791", "id": "assay:1"}, {"role": "analysis", "definition_id": "machine_learning", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["train_classifier", "cross_validate", "roc_auc"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb."], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}