Statistical Learning of Large-Scale Genetic Data: How to Run a Genome-Wide Association Study of Gene-Expression Data Using the 1000 Genomes Project Data.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:38495080", "kind": "science-record", "revision": 2, "content_hash": "sha256:f4c80adb0a82a6b54b4ad193f159cfb36f01eb16b43ea5e4fd8ccc1768088ab2", "created_at": "2026-11-01T19:20:26Z", "updated_at": "2026-11-01T21:32:42Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Statistical Learning of Large-Scale Genetic Data: How to Run a Genome-Wide Association Study of Gene-Expression Data Using the 1000 Genomes Project Data.", "doi": "10.1007/s12561-023-09375-9", "pmid": "38495080", "authors": ["Anton Sugolov", "Eric Emmenegger", "Andrew D Paterson", "Lei Sun"], "journal": "Statistics in biosciences", "year": 2024, "abstract": "Teaching statistics through engaging applications to contemporary large-scale datasets is essential to attracting students to the field. To this end, we developed a hands-on, week-long workshop for senior high-school or junior undergraduate students, without prior knowledge in statistical genetics but with some basic knowledge in data science, to conduct their own genome-wide association study (GWAS). The GWAS was performed for open source gene expression data, using publicly available human genetics data. Assisted by a detailed instruction manual, students were able to obtain ∼1.4 million p-values from a real scientific study, within several days. This early motivation kept students engaged in learning the theories that support their results, including regression, data visualization, results interpretation, and large-scale multiple hypothesis testing. To further their learning motivation by emphasizing the personal connection to this type of data analysis, students were encouraged to make short presentations about how GWAS has provided insights into the genetic basis of diseases that are present in their friends or families. The appended open source, step-by-step instruction manual includes descriptions of the datasets used, the software needed, and results from the workshop. Additionally, scripts used in the workshop are archived on Github and Zenodo to further enhance reproducible research and training.", "claims": [{"id": "claim-method", "statement": "Since gene expression data is continuous, students should specify a linear regression, with PLINK v1.9 command –linear or PLINK v2.0 command –glm.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:38495080", "locator": "sections[5] \"Overview of the Workshop Content\" paragraph 4", "quote": "Since gene expression data is continuous, students should specify a linear regression, with PLINK v1.9 command –linear or PLINK v2.0 command –glm."}], "state": "complete"}, {"id": "claim-data-access", "statement": "Finally, the tested workshop datasets and other materials were also made publicly available on Zenodo.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:38495080", "locator": "sections[9] \"Discussion\" paragraph 3", "quote": "Finally, the tested workshop datasets and other materials were also made publicly available on Zenodo."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Students obtained approximately 1.4 million p-values from a real scientific study within several days.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:38495080", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Assisted by a detailed instruction manual, students were able to obtain ∼1.4 million p-values from a real scientific study, within several days."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Students obtained approximately 1.4 million p-values from a real scientific study within several days.", "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:38495080"], "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": "blocked", "blocked_on": "data:accession_absent", "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": [], "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": [], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:obtain_source_data blocked on data:accession_absent"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:38495080", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1007/s12561-023-09375-9", "sha256": "sha256:3d03e2190c95b44697c7766611d8c4dc4338a1f9ee9894a2ae8a620d4d355ec3", "bytes": 39474, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1007/s12561-023-09375-9", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10940486", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}], "gaps": []}, "analysis": {"state": "complete", "methods": [{"id": "method-1", "family": "regression_modelling", "specific": "ols", "label": "Regression modelling", "dsl_operations": ["ols_fit"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["predictors_and_outcome"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Since gene expression data is continuous, students should specify a linear regression, with PLINK v1.9 command –linear or PLINK v2.0 command –glm."}}, {"id": "method-2", "family": "genetic_association", "specific": "gwas", "label": "Genetic association (GWAS / linkage / MR)", "dsl_operations": ["gwas_locus_scan"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["genotype_phenotype_table"], "evidence": {"score": 0.7, "band": "medium", "quote": "Genome-Wide Association Study"}}], "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": "blocked", "tier": "R4", "blocking_axis": "data", "blocked_on": "data:accession_absent", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["ols_fit"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: ols."}, "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": "R4", "reachable": false, "accessions": [], "resolved_accessions": [], "statements": ["Finally, the tested workshop datasets and other materials were also made publicly available on Zenodo."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": ["ols_fit"], "assets": ["source:pmid:38495080"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:None"]}, "gaps": ["data:accession_absent"]}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T19:12:15Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 11480, "completion_tokens": 1834, "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:38495080", "type": "paper", "label": "pmid:38495080"}, {"id": "source:pmid:38495080", "type": "manuscript", "label": "PMID 38495080 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:38495080", "type": "paper", "label": "pmid:38495080"}, {"id": "source:pmid:38495080", "type": "manuscript", "label": "PMID 38495080 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-38495080", "type": "extraction", "used": ["source:pmid:38495080"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-None", "type": "accession_resolution", "used": ["data:geo:None"], "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:38495080", "parent_hash": "sha256:d63cb33ed86d425da895c69d9e658bf78315addab578a49d5ea72bd0657030af", "fork_of": null, "generation": 1}, "gaps": ["execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent"], "assays": [{"role": "analysis", "definition_id": "regression_modelling", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["ols_fit"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Since gene expression data is continuous, students should specify a linear regression, with PLINK v1.9 command –linear or PLINK v2.0 command –glm."], "order_basis": "analysis_only", "id": "assay:1"}]}