A novel computational approach for predicting complex phenotypes in Drosophila (starvation-sensitive and sterile) by deriving their gene expression signatures from public data.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:33104720", "kind": "science-record", "revision": 7, "content_hash": "sha256:f5d24f4043b30a2bca00c5f39337746c534a73fad0de627081bdffe4c450ff0f", "created_at": "2026-11-01T15:41:14Z", "updated_at": "2026-11-01T22:02:24Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "A novel computational approach for predicting complex phenotypes in Drosophila (starvation-sensitive and sterile) by deriving their gene expression signatures from public data.", "doi": "10.1371/journal.pone.0240824", "pmid": "33104720", "authors": ["Dobril K Ivanov", "Gerrit Bostelmann", "Benoit Lan-Leung", "Julie Williams", "Linda Partridge", "Valentina Escott-Price", "Janet M Thornton"], "journal": "PloS one", "year": 2020, "abstract": "Many research teams perform numerous genetic, transcriptomic, proteomic and other types of omic experiments to understand molecular, cellular and physiological mechanisms of disease and health. Often (but not always), the results of these experiments are deposited in publicly available repository databases. These data records often include phenotypic characteristics following genetic and environmental perturbations, with the aim of discovering underlying molecular mechanisms leading to the phenotypic responses. A constrained set of phenotypic characteristics is usually recorded and these are mostly hypothesis driven of possible to record within financial or practical constraints. We present a novel proof-of-principal computational approach for combining publicly available gene-expression data from control/mutant animal experiments that exhibit a particular phenotype, and we use this approach to predict unobserved phenotypic characteristics in new experiments (data derived from EBI's ArrayExpress and ExpressionAtlas respectively). We utilised available microarray gene-expression data for two phenotypes (starvation-sensitive and sterile) in Drosophila. The data were combined using a linear-mixed effects model with the inclusion of consecutive principal components to account for variability between experiments in conjunction with Gene Ontology enrichment analysis. We present how available data can be ranked in accordance to a phenotypic likelihood of exhibiting these two phenotypes using random forest. The results from our study show that it is possible to integrate seemingly different gene-expression microarray data and predict a potential phenotypic manifestation with a relatively high degree of confidence (>80% AUC). This provides thus far unexplored opportunities for inferring unknown and unbiased phenotypic characteristics from already performed experiments, in order to identify studies for future analyses. Molecular mechanisms associated with gene and environment perturbations are intrinsically linked and give rise to a variety of phenotypic manifestations. Therefore, unravelling the phenotypic spectrum can help to gain insights into disease mechanisms associated with gene and environmental perturbations. Our approach uses public data that are set to increase in volume, thus providing value for money.", "claims": [{"id": "claim-method", "statement": "A random intercept linear mixed-effects model (LMEM) was used to generate normalised residuals for each gene within the Affymetrix Genome 2.0 microarray, accounting for a number of consecutive principal components.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33104720", "locator": "sections[8] \"Generation of the molecular signatures (linear-mixed effects model)\" paragraph 2", "quote": "A random intercept linear mixed-effects model (LMEM) was used to generate normalised residuals for each gene within the Affymetrix Genome 2.0 microarray, accounting for a number of consecutive principal components."}], "state": "complete"}, {"id": "claim-data-access", "statement": "All relevant data are within the manuscript and its Supporting Information files.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:33104720", "locator": "sections[23] \"Data Availability\" paragraph 2", "quote": "All relevant data are within the manuscript and its Supporting Information files."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "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:33104720"], "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:E-MTAB-3546"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "blocked", "blocked_on": "method:regression_modelling_mixed_effects", "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:E-MTAB-3546"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:run_analysis blocked on method:regression_modelling_mixed_effects"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:33104720", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1371/journal.pone.0240824", "sha256": "sha256:07ba292470dd04387831a46924c63fba435de662620e1ff3e5063f7a0a53ff5d", "bytes": 37860, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1371/journal.pone.0240824", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC7588067", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:E-MTAB-3546", "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=E-MTAB-3546", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "E-MTAB-3546", "version": null, "note": "All relevant data are within the manuscript and its Supporting Information files.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "All relevant data are within the manuscript and its Supporting Information files.", "access_statement_kind": "data_availability_statement", "sample_accessions": [], "sample_count": null, "group_sizes": null, "group_sizes_basis": null, "platform": null, "taxon": null, "title": null, "resolution_checked_at": "2026-10-31T18:28:59Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "asset:table:pmid:33104720:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:535f4136bfd080ee6b631aec6a3f7ccb3846f0ee7efb291d98a64402810dbd03", "bytes": 433, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[16]/table-wrap[1]#pone.0240824.t001", "version": null, "note": "One sample t-test for class probabilities (controls/mutants) in the two phenotypes following LOOCV.", "structure": {"pmcid": "PMC7588067", "section_index": 16, "section_title": "Determining the number of PCs for unwanted variation", "row_count": 4, "column_count": 3}, "table": {"label": "Table 1", "caption": "One sample t-test for class probabilities (controls/mutants) in the two phenotypes following LOOCV.", "rows": [["Class", "Phenotype"], ["starvation-sensitive p-value (μ = 0.5)", "Sterile p-value (μ = 0.5)"], ["Controls", "5.72x10-03", "3.87x10-03"], ["Mutants", "5.84x10-03", "3.10x10-03"]], "row_count": 4, "column_count": 3, "numeric_cells": 6, "header_rows": 2, "locator": "section[16]/table-wrap[1]#pone.0240824.t001"}}], "gaps": ["asset:data:geo:E-MTAB-3546 not downloaded", "asset:data:geo:E-MTAB-3546 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "regression_modelling", "specific": "mixed_effects", "label": "Regression modelling", "dsl_operations": [], "missing_operations": ["regression_modelling_mixed_effects"], "expressible": false, "parameters": {}, "data_shapes": ["predictors_and_outcome"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "A random intercept linear mixed-effects model (LMEM) was used to generate normalised residuals for each gene within the Affymetrix Genome 2.0 microarray, accounting for a number of consecutive principal components."}}, {"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"}}], "targets": [], "gaps": ["method.reader_differs_from_extractor", "method.regression_modelling_mixed_effects"]}, "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": false, "required_operations": [], "missing_operations": ["regression_modelling_mixed_effects"], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: mixed_effects."}, "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": ["E-MTAB-3546"], "resolved_accessions": ["E-MTAB-3546"], "statements": ["All relevant data are within the manuscript and its Supporting Information files."], "basis": "The paper names E-MTAB-3546 in EBI’s ExpressionAtlas and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:33104720", "data:geo:E-MTAB-3546", "asset:table:pmid:33104720:1"], "steps": []}, "must_build": {"dsl_operations": [{"operation": "regression_modelling_mixed_effects", "unlocks_papers": null, "priority": null}], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:E-MTAB-3546"]}, "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:23:44Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 3, "usage": {"prompt_tokens": 26350, "completion_tokens": 4804, "cached_tokens": 792}}, "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:33104720", "type": "paper", "label": "pmid:33104720"}, {"id": "source:pmid:33104720", "type": "manuscript", "label": "PMID 33104720 full text"}, {"id": "data:geo:E-MTAB-3546", "type": "dataset", "label": "GEO E-MTAB-3546"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:33104720", "type": "paper", "label": "pmid:33104720"}, {"id": "source:pmid:33104720", "type": "manuscript", "label": "PMID 33104720 full text"}, {"id": "data:geo:E-MTAB-3546", "type": "dataset", "label": "GEO E-MTAB-3546"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-33104720", "type": "extraction", "used": ["source:pmid:33104720"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-E-MTAB-3546", "type": "accession_resolution", "used": ["data:geo:E-MTAB-3546"], "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:33104720", "parent_hash": "sha256:17448f6b521c63fb421047960a26dbd80eb177a9075d8332c9af0c42cc6eb381", "fork_of": null, "generation": 6}, "gaps": ["asset:data:geo:E-MTAB-3546 licence not read", "asset:data:geo:E-MTAB-3546 not downloaded", "execution.not_run", "method.reader_differs_from_extractor", "method.regression_modelling_mixed_effects", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:run_analysis blocked on method:regression_modelling_mixed_effects"], "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": ["onAtlas respectively). We utilised available microarray gene-expression data for two phenotypes (sta", "ntegrate seemingly different gene-expression microarray data and predict a potential phenotypic mani"], "order_basis": "paper_text_order:1200", "id": "assay:1"}, {"role": "analysis", "definition_id": "regression_modelling", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["A random intercept linear mixed-effects model (LMEM) was used to generate normalised residuals for each gene within the Affymetrix Genome 2.0 microarray, accounting for a number of consecutive principal components."], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}