EndoTime: non-categorical timing estimates for luteal endometrium.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:35092277", "kind": "science-record", "revision": 5, "content_hash": "sha256:b56484576ccce3c51df7c1e19a72c7b90fca69d5bf4844e1a9c58b96537b8665", "created_at": "2026-11-01T17:36:49Z", "updated_at": "2026-11-01T22:05:41Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "EndoTime: non-categorical timing estimates for luteal endometrium.", "doi": "10.1093/humrep/deac006", "pmid": "35092277", "authors": ["Julia Lipecki", "Andrew E Mitchell", "Joanne Muter", "Emma S Lucas", "Komal Makwana", "Katherine Fishwick", "Joshua Odendaal", "Amelia Hawkes", "Pavle Vrljicak", "Jan J Brosens", "Sascha Ott"], "journal": "Human reproduction (Oxford, England)", "year": 2022, "abstract": "STUDY QUESTION: Can the accuracy of timing of luteal phase endometrial biopsies based on urinary ovulation testing be improved by measuring the expression of a small number of genes and a continuous, non-categorical modelling approach?\nSUMMARY ANSWER: Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample.\nWHAT IS KNOWN ALREADY: Commercially available endometrial timing approaches based on gene expression require large gene sets and use a categorical approach that classifies samples as pre-receptive, receptive or post-receptive.\nSTUDY DESIGN, SIZE, DURATION: Gene expression was measured by RTq-PCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4-12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RTq-PCR as well as RNA-sequencing.\nPARTICIPANTS/MATERIALS, SETTING, METHODS: A computational procedure, named 'EndoTime', was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points.\nMAIN RESULTS AND THE ROLE OF CHANCE: The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e-6; DPP4: P < 3.7e-13). Pearson correlation between EndoTime estimates for the same sample set but based on RTq-PCR or RNA-sequencing data showed a high degree of congruency between the two (P = 8.6e-10, R2 = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (P > 0.05).\nLARGE SCALE DATA: The RTq-PCR data files are available via the GitHub repository for the EndoTime software at https://github.com/AE-Mitchell/EndoTime, as is the code used for pre-processing of RTq-PCR data. The RNA-sequencing data are available on GEO (accession GSE180485).\nLIMITATIONS, REASONS FOR CAUTION: Timing estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Our approach, in its current iteration, performs batch correction such that larger sample batches impart greater accuracy to timing estimations. In theory, our method requires endometrial samples obtained at different days in the luteal phase. In practice, however, this is not a concern as timings based on urinary ovulation testing are associated with a sufficient level of noise to ensure that a variety of time points will be sampled.\nWIDER IMPLICATIONS OF THE FINDINGS: Our method is the first to assay the temporal state of luteal-phase endometrial samples on a continuous domain. It is freely available with fully shared data and open-source software. EndoTime enables accurate temporal profiling of any gene in luteal endometrial samples for a wide range of research applications and, potentially, clinical use.\nSTUDY FUNDING/COMPETING INTEREST(S): This study was supported by a Wellcome Trust Investigator Award (Grant/Award Number: 212233/Z/18/Z) and the Tommy's National Miscarriage Research Centre. None of the authors have any competing interests. J.L. was funded by the Biotechnology and Biological Sciences Research Council (UK) through the Midlands Integrative Biology Training Partnership (MIBTP, BB/M01116X/1).", "claims": [{"id": "claim-method", "statement": "The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:35092277", "locator": "sections[0] \"Abstract\" paragraph 6", "quote": "The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13)."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The RNA-seq data are available in the Gene Expression Omnibus under accession GSE180485. The RTq-PCR data are available in the GitHub repository for the EndoTime software at https://www.github.com/AE-Mitchell/EndoTime.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:35092277", "locator": "sections[21] \"Data availability\" paragraph 2", "quote": "The RNA-seq data are available in the Gene Expression Omnibus under accession GSE180485. The RTq-PCR data are available in the GitHub repository for the EndoTime software at https://www.github.com/AE-Mitchell/EndoTime."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:35092277", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample.", "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:35092277"], "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:GSE180485"], "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:group_comparison_wilcoxon", "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:GSE180485"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:run_analysis blocked on method:group_comparison_wilcoxon"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:35092277", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1093/humrep/deac006", "sha256": "sha256:63a811306e03dba7ab3cf56ddda4575179e6d873d102ebcbe28a3a0ab8140160", "bytes": 50451, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1093/humrep/deac006", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC8971653", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE180485", "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=GSE180485", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE180485", "version": null, "note": "The RNA-seq data are available in the Gene Expression Omnibus under accession GSE180485. The RTq-PCR data are available in the GitHub repository for the EndoTime software at https://www.github.com/AE-Mitchell/EndoTime.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The RNA-seq data are available in the Gene Expression Omnibus under accession GSE180485. The RTq-PCR data are available in the GitHub repository for the EndoTime software at https://www.github.com/AE-Mitchell/EndoTime.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM5463770", "GSM5463771", "GSM5463772", "GSM5463773", "GSM5463774", "GSM5463775", "GSM5463776", "GSM5463777", "GSM5463778", "GSM5463779", "GSM5463780", "GSM5463781", "GSM5463782", "GSM5463783", "GSM5463784", "GSM5463785", "GSM5463786", "GSM5463787", "GSM5463788", "GSM5463789", "GSM5463790", "GSM5463791", "GSM5463792", "GSM5463793", "GSM5463794", "GSM5463795", "GSM5463796", "GSM5463797", "GSM5463798", "GSM5463799", "GSM5463800", "GSM5463801", "GSM5463802", "GSM5463803", "GSM5463804", "GSM5463805", "GSM5463806", "GSM5463807", "GSM5463808", "GSM5463809", "GSM5463810", "GSM5463811", "GSM5463812", "GSM5463813", "GSM5463814", "GSM5463815", "GSM5463816", "GSM5463817", "GSM5463818", "GSM5463819", "GSM5463820", "GSM5463821", "GSM5463822", "GSM5463823", "GSM5463824", "GSM5463825", "GSM5463826", "GSM5463827", "GSM5463828", "GSM5463829", "GSM5463830", "GSM5463831", "GSM5463832", "GSM5463833", "GSM5463834", "GSM5463835", "GSM5463836", "GSM5463837", "GSM5463838", "GSM5463839", "GSM5463840", "GSM5463841", "GSM5463842", "GSM5463843", "GSM5463844", "GSM5463845", "GSM5463846", "GSM5463847", "GSM5463848", "GSM5463849", "GSM5463850", "GSM5463851", "GSM5463852", "GSM5463853", "GSM5463854", "GSM5463855", "GSM5463856", "GSM5463857", "GSM5463858", "GSM5463859", "GSM5463860", "GSM5463861", "GSM5463862", "GSM5463863", "GSM5463864", "GSM5463865", "GSM5463866", "GSM5463867", "GSM5463868", "GSM5463869", "GSM5463870", "GSM5463871", "GSM5463872", "GSM5463873", "GSM5463874", "GSM5463875", "GSM5463876", "GSM5463877"], "sample_count": 108, "group_sizes": null, "group_sizes_basis": null, "platform": "20301", "taxon": "Homo sapiens", "title": "EndoTime: Non-categorical timing estimates for endometrial tissue during the secretory phase", "resolution_checked_at": "2026-11-01T18:06:43Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE180485 not downloaded", "asset:data:geo:GSE180485 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "group_comparison", "specific": "wilcoxon", "label": "Group comparison", "dsl_operations": [], "missing_operations": ["group_comparison_wilcoxon"], "expressible": false, "parameters": {}, "data_shapes": ["two_or_more_groups_continuous"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13)."}}, {"id": "method-2", "family": "correlation_analysis", "specific": "pearson", "label": "Correlation analysis", "dsl_operations": ["pearson_r"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 0.7, "band": "medium", "quote": "Pearson correlation"}}, {"id": "method-3", "family": "correlation_analysis", "specific": "spearman", "label": "Correlation analysis", "dsl_operations": ["spearman_rho"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 0.8, "band": "high", "quote": "Spearman correlation"}}], "targets": [], "gaps": ["method.group_comparison_wilcoxon"]}, "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": ["group_comparison_wilcoxon"], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: wilcoxon."}, "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": ["GSE180485"], "resolved_accessions": ["GSE180485"], "statements": ["The RNA-seq data are available in the Gene Expression Omnibus under accession GSE180485. The RTq-PCR data are available in the GitHub repository for the EndoTime software at https://www.github.com/AE-Mitchell/EndoTime."], "basis": "The paper names GSE180485 in Gene Expression Omnibus and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:35092277", "data:geo:GSE180485"], "steps": []}, "must_build": {"dsl_operations": [{"operation": "group_comparison_wilcoxon", "unlocks_papers": null, "priority": null}], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE180485"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T17:31:34Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 12354, "completion_tokens": 1754, "cached_tokens": 256}}, "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:35092277", "type": "paper", "label": "pmid:35092277"}, {"id": "source:pmid:35092277", "type": "manuscript", "label": "PMID 35092277 full text"}, {"id": "data:geo:GSE180485", "type": "dataset", "label": "GEO GSE180485"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:35092277", "type": "paper", "label": "pmid:35092277"}, {"id": "source:pmid:35092277", "type": "manuscript", "label": "PMID 35092277 full text"}, {"id": "data:geo:GSE180485", "type": "dataset", "label": "GEO GSE180485"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-35092277", "type": "extraction", "used": ["source:pmid:35092277"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE180485", "type": "accession_resolution", "used": ["data:geo:GSE180485"], "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:35092277", "parent_hash": "sha256:fd032ebc0ef1016e8325b8121d9741a47f8003989c2aa7ce2243845f8c95a220", "fork_of": null, "generation": 4}, "gaps": ["asset:data:geo:GSE180485 licence not read", "asset:data:geo:GSE180485 not downloaded", "execution.not_run", "method.group_comparison_wilcoxon", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:run_analysis blocked on method:group_comparison_wilcoxon"], "assays": [{"role": "data_generation", "definition_id": "rna_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["samples were profiled by RTq-PCR as well as RNA-sequencing. PARTICIPANTS/MATERIALS, SETTING, METHODS: A", "the same sample set but based on RTq-PCR or RNA-sequencing data showed a high degree of congruency betw", "used for pre-processing of RTq-PCR data. The RNA-sequencing data are available on GEO (accession GSE1804"], "order_basis": "paper_text_order:1014", "id": "assay:1"}, {"role": "data_generation", "definition_id": "ihc", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["ogen, and one part fixed in 10% formalin for immunohistochemistry. Two sets of RTq-PCR data were utilized in t"], "order_basis": "paper_text_order:5174", "id": "assay:2"}, {"role": "data_generation", "definition_id": "qpcr", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["ize and apply minor batch correction to –ΔCT qPCR gene expression data and transform reported"], "order_basis": "paper_text_order:13521", "id": "assay:3"}, {"role": "analysis", "definition_id": "group_comparison", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13)."], "order_basis": "consumes_earlier_output", "id": "assay:4"}]}