{"schema":"https://automatedassays.local/schemas/real-data-evidence-bundle/1.0.0","schema_version":"1.0.0","bundle_public_id":"976195ee-994c-4c75-8e3d-aab80fdd3bab","generated_at":"2026-11-03T17:43:58.000000Z","mode":"real_data","label":"Differential expression (edgeR quasi-likelihood)","record":{"record_id":"pmid:27508061","science_record_content_hash":null,"science_record_source":null},"operation":{"name":"rna_seq_differential_expression","version":"0.1.0","parameters":{"paper":"pmid:27508061","paper_version":"2","stage":"differential_expression","source_files":{"targets.txt":"sha256:88ae42c0539403073f0bbbd53287dd9c9cf61ddd5b1a42856243d7d99113cb8b","edgeRQL.R":"sha256:3b11478ac314c4456af8b4ba7af7992e938488ffc7bec88d3f73f42aed4c0bba","edgeRQL.Rnw":"sha256:defd3acb73dcbfe885ca02209ac9437b785b4ccf12c7ba09a8a8bb311173902d"},"software":{"edgeR":"3.14.0","limma":"3.28.21","org.Mm.eg.db":"3.3.0","GO.db":"3.3.0","gplots":"3.0.1"},"deviations":{"kegg_transport":"KEGG's REST endpoint now redirects http to https, and the pinned R's url() reader cannot follow a redirect, so the gene-to-pathway table kegga reads is fetched over https and passed to limma::kegga as the gene.pathway and pathway.names arguments that limma documents for supplying those tables directly. The species argument, the gene sets and the analysis are the authors' own.","r_patch_release":"The paper ran on R 3.3.1 and this platform runs on R 3.3.2.","limma_patch_release":"The paper ran on limma 3.28.6 and this platform runs on limma 3.28.21, the last patch of the same Bioconductor 3.3 series.","annotation_dbi_patch_release":"The paper ran on AnnotationDbi 1.34.3 and this platform runs on AnnotationDbi 1.38.0, the last patch of the same Bioconductor 3.3 series.","operating_system":"The paper ran on Windows 7 x64 and this platform runs on Linux."}}},"dataset":{"asset_id":"GEO:GSE60450:Lactation-GenewiseCounts","source_url":"https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE60450","licence":{"state":"open"}},"result":{"stage":"differential_expression","label":"Differential expression (edgeR quasi-likelihood)","status":"complete","numbers":{"stage":"differential_expression","package":"edgeR 3.14.0","limma":"limma 3.28.21","design":"model.matrix(~0+group), one column per cell type and status","contrasts":{"primary":"B.lactating - B.pregnant","luminal":["L.pregnant - L.lactating","L.virgin - L.lactating","L.virgin - L.pregnant"],"interaction":"(L.lactating - L.pregnant) - (B.lactating - B.pregnant)"},"dispersion":"estimateDisp(robust = TRUE)","fit":"glmQLFit(robust = TRUE)","test":"glmQLFTest","treat_lfc":0.5849625007,"qlf_fdr_0.05":{"down":2757,"not_significant":10408,"up":2488,"total":5245},"treat_fdr_0.05":{"down":723,"not_significant":14530,"up":400,"total":1123},"top_gene":{"entrez_id":"12992","symbol":"Csn1s2b","logFC":6.084638853,"logCPM":10.18750496,"fdr":7.479370032e-7},"top_table":[{"entrez_id":"12992","symbol":"Csn1s2b","logFC":6.084638853,"logCPM":10.18750496,"fdr":7.479370032e-7},{"entrez_id":"211577","symbol":"Mrgprf","logFC":5.145799852,"logCPM":2.749212987,"fdr":7.990433387e-7},{"entrez_id":"226101","symbol":"Myof","logFC":2.322206353,"logCPM":6.445570291,"fdr":7.990433387e-7},{"entrez_id":"381290","symbol":"Atp2b4","logFC":2.143691947,"logCPM":6.150633497,"fdr":7.990433387e-7},{"entrez_id":"140474","symbol":"Muc4","logFC":-7.174703375,"logCPM":6.056779838,"fdr":8.145409077e-7},{"entrez_id":"231830","symbol":"Micall2","logFC":-2.250634292,"logCPM":5.188178793,"fdr":1.149727669e-6},{"entrez_id":"24117","symbol":"Wif1","logFC":-1.820902768,"logCPM":6.770643875,"fdr":1.595267172e-6},{"entrez_id":"12740","symbol":"Cldn4","logFC":-5.321129504,"logCPM":9.872591634,"fdr":1.744562078e-6},{"entrez_id":"21953","symbol":"Tnni2","logFC":5.751919396,"logCPM":3.863649801,"fdr":1.744562078e-6},{"entrez_id":"231991","symbol":"Creb5","logFC":2.573613406,"logCPM":4.871705678,"fdr":1.831355452e-6}],"df_prior":{"Min.":3,"1st Qu.":6.766,"Median":6.766,"Mean":6.626,"3rd Qu.":6.766,"Max.":6.766},"common_dispersion":0.01344354583,"trended_dispersion":[]},"artefact_sha256":{"top_tags_qlf.tsv":"sha256:874f770e0b4f171c2697c7eb5f19edf635c63d2476985722090b93af93887246","top_tags_treat.tsv":"sha256:d3a7184113f19753d7ff2dcad36e43b107ebb068ca9baca3ec63544320a6b9d4","de_counts.tsv":"sha256:1b06df4231d37a618f7510a9ab512b77193264f35984e464dd0c03951685348d","bcv.png":"sha256:0689a5070d150d7c4930ac409de42249112b7c7bbcc92b3a935394d38da16776","qldisp.png":"sha256:c8df377af2b98d69a85cf652101a2fc75402d9dddde596cdae06ecaf0389d1b2","md_fit.png":"sha256:9be8dea8f650d766f82b0d0ccede3337330fdae847e530bc25644fa6641d0f73","md_treat.png":"sha256:1f96236b44b61bf326898c2cfed8550a22cf6e0869beae6d54cf4480304a5c52"},"comparison":{"stage":"differential_expression","label":"Differential expression (edgeR quasi-likelihood)","source":"The paper prints the leading table of the basal lactating against basal pregnant comparison, the counts of differentially expressed genes from both the quasi-likelihood test and the TREAT test, and the leading table of the TREAT test.","criteria":[{"name":"leading_gene","statement":"The leading gene of the basal lactating against basal pregnant comparison is Csn1s2b, with a large positive log fold change.","comparison":"exact on the symbol and the sign of the fold change","published":{"symbol":"Csn1s2b","direction":"positive"},"verdict":"reproduced","agreements":["the leading gene"],"differences":[],"compared":1,"detail":"1 of 1 agree with the paper"},{"name":"qlf_top_table","statement":"The leading table of the quasi-likelihood test carries the same ten genes, in the same order, with the same log fold changes and the same false discovery rates as the paper prints.","comparison":"the log fold change to the two decimals the paper prints and the false discovery rate to the three significant figures the paper prints","published":[{"symbol":"Csn1s2b","logFC":6.08,"FDR":7.48e-7},{"symbol":"Mrgprf","logFC":5.15,"FDR":7.99e-7},{"symbol":"Myof","logFC":2.32,"FDR":7.99e-7},{"symbol":"Atp2b4","logFC":2.14,"FDR":7.99e-7},{"symbol":"Muc4","logFC":-7.17,"FDR":8.15e-7},{"symbol":"Micall2","logFC":-2.25,"FDR":1.15e-6},{"symbol":"Wif1","logFC":-1.82,"FDR":1.6e-6},{"symbol":"Cldn4","logFC":-5.32,"FDR":1.74e-6},{"symbol":"Tnni2","logFC":5.75,"FDR":1.74e-6},{"symbol":"Creb5","logFC":2.57,"FDR":1.83e-6}],"verdict":"reproduced","agreements":["Csn1s2b","Mrgprf","Myof","Atp2b4","Muc4","Micall2","Wif1","Cldn4","Tnni2","Creb5"],"differences":[],"compared":10,"detail":"10 of 10 agree with the paper"},{"name":"qlf_de_gene_counts","statement":"The quasi-likelihood test finds the same number of differentially expressed genes at a false discovery rate of 5%, in each direction.","comparison":"exact","published":{"down":2757,"not_significant":10408,"up":2488,"total":5245},"verdict":"reproduced","agreements":["down","not_significant","up","total"],"differences":[],"compared":4,"detail":"4 of 4 agree with the paper"},{"name":"treat_de_gene_counts","statement":"The TREAT test at a fold change of 1.5 finds the same number of differentially expressed genes at a false discovery rate of 5%, in each direction.","comparison":"exact","published":{"down":723,"not_significant":14530,"up":400,"total":1123},"verdict":"reproduced","agreements":["down","not_significant","up","total"],"differences":[],"compared":4,"detail":"4 of 4 agree with the paper"},{"name":"de_gene_count","statement":"The comparison identifies more than 5000 differentially expressed genes at a false discovery rate of 5%.","comparison":"the paper's own statement, checked against the exact counts above","published":{"more_than":5000},"verdict":"reproduced","agreements":["the number of genes"],"differences":[],"compared":1,"detail":"1 of 1 agree with the paper"},{"name":"test_used","statement":"The comparison is a quasi-likelihood F-test, not a likelihood ratio test.","comparison":"exact","published":"glmQLFTest with robust dispersions","verdict":"reproduced","agreements":["the method"],"differences":[],"compared":1,"detail":"1 of 1 agree with the paper"}],"summary":{"reproduced":6,"not_reproduced":0,"not_scored":0}},"log_sha256":"sha256:a7bf7a9203e798e569fcb8ebc783935cd4d2276869892a6ed00c9a7af29b6719","started_at":"2026-11-03T17:43:09.972641Z","completed_at":"2026-11-03T17:43:58.598189Z"},"integrity":{"dataset_sha256":"sha256:46cd415ac9b82853efb608bddc98cddbf8c7a908807877d448d9c73419b615cc","result_sha256":"sha256:194ba83e9468a04462de791795cd0bdf4bbdc992bb5cea0b420fa659fe8ecc1f","bundle_hash":"ffd49280b7bcfa17e938f077569600a4d38d3dce873176d5f4f8f5aa230659da"},"replay_contract":{"description":"Re-running the operation on the dataset identified above must reproduce the recorded result hash exactly.","dataset_asset_id":"GEO:GSE60450:Lactation-GenewiseCounts","dataset_source_url":"https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE60450","expected_dataset_sha256":"sha256:46cd415ac9b82853efb608bddc98cddbf8c7a908807877d448d9c73419b615cc","expected_result_sha256":"sha256:194ba83e9468a04462de791795cd0bdf4bbdc992bb5cea0b420fa659fe8ecc1f","operation":"rna_seq_differential_expression","operation_version":"0.1.0","design":{"paper":"pmid:27508061","paper_version":"2","stage":"differential_expression","source_files":{"targets.txt":"sha256:88ae42c0539403073f0bbbd53287dd9c9cf61ddd5b1a42856243d7d99113cb8b","edgeRQL.R":"sha256:3b11478ac314c4456af8b4ba7af7992e938488ffc7bec88d3f73f42aed4c0bba","edgeRQL.Rnw":"sha256:defd3acb73dcbfe885ca02209ac9437b785b4ccf12c7ba09a8a8bb311173902d"},"software":{"edgeR":"3.14.0","limma":"3.28.21","org.Mm.eg.db":"3.3.0","GO.db":"3.3.0","gplots":"3.0.1"},"deviations":{"kegg_transport":"KEGG's REST endpoint now redirects http to https, and the pinned R's url() reader cannot follow a redirect, so the gene-to-pathway table kegga reads is fetched over https and passed to limma::kegga as the gene.pathway and pathway.names arguments that limma documents for supplying those tables directly. The species argument, the gene sets and the analysis are the authors' own.","r_patch_release":"The paper ran on R 3.3.1 and this platform runs on R 3.3.2.","limma_patch_release":"The paper ran on limma 3.28.6 and this platform runs on limma 3.28.21, the last patch of the same Bioconductor 3.3 series.","annotation_dbi_patch_release":"The paper ran on AnnotationDbi 1.34.3 and this platform runs on AnnotationDbi 1.38.0, the last patch of the same Bioconductor 3.3 series.","operating_system":"The paper ran on Windows 7 x64 and this platform runs on Linux."}},"platform":{"accession":null,"source_url":null,"sha256":null}}}