Corpus Engine Towards a reproducible .science file.
← All .science files

pmid:40489546

Detecting SARS-CoV-2 cryptic lineages using publicly available whole genome wastewater sequencing data.

Reinier Suarez · Devon A Gregory · David A Baker · Clayton A Rushford · Torin L Hunter · Nicholas R Minor · Clayton M Russ · Emma E Copen · David H O'Connor · Marc C Johnson

PLoS pathogens · 2025 doi:10.1371/journal.ppat.1012850

Paper

complete

Experiment

partial

Data

partial

Analysis

absent

Execution

absent

Execution evidence has not yet been added to this .science file. Explore its methods, data and open questions to plan the next step.
Record details sci_version · id · kind · revision · content_hash · schema · created_at · updated_at · generator
id
pmid:40489546
kind
science-record
content_hash
sha256:2871ccb2aa2b0b5cb69182a2c4b32cbe88ba16e0d76b681822523b0ef1564319
schema
https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json
created_at
2026-11-01T20:30:36Z
updated_at
2026-11-01T20:30:36Z
generator object · 2 fields
name
corpus-assay-engine
run_id
reading-sci-2026-11-01
paper object · 10 fields
state
complete
title
Detecting SARS-CoV-2 cryptic lineages using publicly available whole genome wastewater sequencing data.
doi
10.1371/journal.ppat.1012850
pmid
40489546
authors array · 10 items
[0]
Reinier Suarez
[1]
Devon A Gregory
[2]
David A Baker
[3]
Clayton A Rushford
[4]
Torin L Hunter
[5]
Nicholas R Minor
[6]
Clayton M Russ
[7]
Emma E Copen
[8]
David H O'Connor
[9]
Marc C Johnson
journal
PLoS pathogens
year
2025
abstract
Beginning in early 2021, unique and highly divergent lineages of SARS-CoV-2 were sporadically found in wastewater sewersheds using a sequencing strategy focused on amplifying the most rapidly evolving region of SARS-CoV-2, the receptor binding domain (RBD). Because these RBD sequences did not match known circulating strains and their source was not known, we termed them "cryptic lineages". To date, more than 20 cryptic lineages have been identified using the RBD-focused sequencing strategy. Here, we identified and characterized additional cryptic lineages from SARS-CoV-2 wastewater sequences submitted to NCBI's Sequence Read Archives (SRA). Wastewater sequence datasets were screened for individual sequence reads that contained combinations of mutations frequently found in cryptic lineages but not contemporary circulating lineages. Using this method, we identified 18 cryptic lineages that appeared in multiple (2-81) samples from the same sewershed, including 12 that were not previously reported. Partial consensus sequences were generated for each cryptic lineage by extracting and mapping sequences containing cryptic-specific mutations. Surprisingly, seven of the mutations that appeared convergently in cryptic lineages were reversions to sequences that were highly conserved in SARS-CoV-2-related enteric bat Sarbecoviruses. The apparent reversion to bat Sarbecovirus sequences is consistent with the notion that SARS-CoV-2 adaptation to replicate efficiently in respiratory tissues preceded the COVID-19 pandemic.
claims array · 2 items
[0] object · 6 fields
id
claim-data-access
statement
To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author.
kind
methodological
target
null
evidence array · 1 item
[0] object · 3 fields
asset
source:pmid:40489546
locator
sections[10] "Supporting information" paragraph 7
quote
To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author.
state
complete
[1] object · 6 fields
id
claim-target-result
statement
Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported.
kind
qualitative
target
null
evidence array · 1 item
[0] object · 3 fields
asset
source:pmid:40489546
locator
sections[0] "Abstract" paragraph 1
quote
Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported.
state
complete
gaps array · 0 items
experiment object · 8 fields
state
partial
objective
Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported.
steps array · 3 items
[0] object · 13 fields
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.
parameters object · 0 fields
inputs array · 1 item
[0]
source:pmid:40489546
outputs array · 1 item
[0]
the published estimate and its evidence
instruments array · 0 items
reagents array · 0 items
[1] object · 13 fields
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.
parameters object · 0 fields
inputs array · 0 items
outputs array · 1 item
[0]
the source data with its repository provenance
instruments array · 0 items
reagents array · 0 items
[2] object · 13 fields
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.
parameters object · 0 fields
inputs array · 0 items
outputs array · 1 item
[0]
the replicated result and its comparison
instruments array · 0 items
reagents array · 0 items
materials array · 0 items
instruments array · 0 items
reagents array · 0 items
conditions object · 0 fields
gaps array · 1 item
[0]
step:obtain_source_data blocked on data:accession_absent
data object · 3 fields
state
partial
assets array · 1 item
[0] object · 15 fields
id
source:pmid:40489546
role
manuscript
format
xml
media_type
application/xml
uri
https://doi.org/10.1371/journal.ppat.1012850
sha256
sha256:95f635c338268ead8ecb8d4ed1d055c3aeb93a5cf1925369bfaf3fb8a920c227
bytes
35122
access
open
state
resolved
reader
null
reader_status
reader_required
identifier
10.1371/journal.ppat.1012850
version
null
note
Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.
structure object · 4 fields
pmcid
PMC12176291
source
europe_pmc_fulltext
coverage
whole_body_without_references
text_hash_scope
rendered_section_body
gaps array · 0 items
analysis object · 4 fields
state
absent
methods array · 0 items
targets array · 0 items
gaps array · 2 items
[0]
method
[1]
operations
capability object · 3 fields
state
complete
requirements array · 1 item
[0] object · 6 fields
capability_id
capability:computational_analysis
label
Computational analysis platform
assay_families array · 9 items
[0]
meta_analysis
[1]
machine_learning
[2]
survival_analysis
[3]
statistical_regression
[4]
pca
[5]
molecular_docking
[6]
network_analysis
[7]
differential_expression
[8]
gwas
candidate
Automated Assays Laravel application and simulated adapters
gaps array · 0 items
readiness object · 8 fields
state
blocked
tier
R4
blocked_on
data:accession_absent
axes object · 3 fields
method object · 5 fields
tier
R4
expressible
false
required_operations array · 0 items
missing_operations array · 0 items
basis
The reader quoted the method the paper names (not named), but the corpus's method extractor could not map it onto a DSL operation, so the method is published but not expressible.
capability object · 4 fields
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 object · 6 fields
tier
R4
reachable
false
accessions array · 0 items
resolved_accessions array · 0 items
statements array · 1 item
[0]
To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author.
basis
The paper names no accession to its source data, so the data has to come from the authors.
can_run object · 3 fields
dsl_operations array · 0 items
assets array · 1 item
[0]
source:pmid:40489546
steps array · 1 item
[0]
run_analysis
must_build object · 5 fields
dsl_operations array · 0 items
converters array · 0 items
readers array · 0 items
capabilities array · 0 items
data_gaps array · 1 item
[0]
data:geo:None
gaps array · 1 item
[0]
data:accession_absent
review object · 5 fields
state
partial
machine_critique object · 5 fields
kind
extraction_review
reader object · 5 fields
model
MiniMax-M3.1-Flash
provider
minimax
provider_family
minimax
read_at
2026-11-01T19:54:03Z
independence
single_reader
evidence_check object · 2 fields
method
null
availability
verbatim
calls
5
usage object · 3 fields
prompt_tokens
9285
completion_tokens
1658
cached_tokens
256
stages array · 1 item
[0] object · 4 fields
id
extraction
state
complete
reader
MiniMax-M3.1-Flash
provider
minimax
gaps array · 3 items
[0]
review.independent_reader
[1]
review.machine_critique
[2]
review.human_critique
execution object · 5 fields
state
absent
status
not_run
mode
null
runs array · 0 items
gaps array · 1 item
[0]
execution.not_run
provenance object · 7 fields
entities array · 3 items
[0] object · 3 fields
id
pmid:40489546
type
paper
label
pmid:40489546
[1] object · 3 fields
id
source:pmid:40489546
type
manuscript
label
PMID 40489546 full text
[2] object · 3 fields
id
data:geo:None
type
dataset
label
GEO None
activities array · 3 items
[0] object · 4 fields
id
activity-reading-sci-2026-11-01
type
classification
used array · 3 items
[0] object · 3 fields
id
pmid:40489546
type
paper
label
pmid:40489546
[1] object · 3 fields
id
source:pmid:40489546
type
manuscript
label
PMID 40489546 full text
[2] object · 3 fields
id
data:geo:None
type
dataset
label
GEO None
agent
agent:corpus-assay-engine
[1] object · 4 fields
id
activity-read-40489546
type
extraction
used array · 1 item
[0]
source:pmid:40489546
agent
agent:minimax-reader
[2] object · 4 fields
id
activity-resolve-None
type
accession_resolution
used array · 1 item
[0]
data:geo:None
agent
agent:corpus-assay-engine
agents array · 2 items
[0] object · 3 fields
id
agent:corpus-assay-engine
type
software
label
Corpus assay engine
[1] object · 3 fields
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 object · 4 fields
gaps array · 7 items
[0]
execution.not_run
[1]
method
[2]
operations
[3]
review.human_critique
[4]
review.independent_reader
[5]
review.machine_critique
[6]
step:obtain_source_data blocked on data:accession_absent
Complete file
{"sci_version": "1.0", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:40489546", "kind": "science-record", "revision": 1, "content_hash": "sha256:2871ccb2aa2b0b5cb69182a2c4b32cbe88ba16e0d76b681822523b0ef1564319", "created_at": "2026-11-01T20:30:36Z", "updated_at": "2026-11-01T20:30:36Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Detecting SARS-CoV-2 cryptic lineages using publicly available whole genome wastewater sequencing data.", "doi": "10.1371/journal.ppat.1012850", "pmid": "40489546", "authors": ["Reinier Suarez", "Devon A Gregory", "David A Baker", "Clayton A Rushford", "Torin L Hunter", "Nicholas R Minor", "Clayton M Russ", "Emma E Copen", "David H O'Connor", "Marc C Johnson"], "journal": "PLoS pathogens", "year": 2025, "abstract": "Beginning in early 2021, unique and highly divergent lineages of SARS-CoV-2 were sporadically found in wastewater sewersheds using a sequencing strategy focused on amplifying the most rapidly evolving region of SARS-CoV-2, the receptor binding domain (RBD). Because these RBD sequences did not match known circulating strains and their source was not known, we termed them \"cryptic lineages\". To date, more than 20 cryptic lineages have been identified using the RBD-focused sequencing strategy. Here, we identified and characterized additional cryptic lineages from SARS-CoV-2 wastewater sequences submitted to NCBI's Sequence Read Archives (SRA). Wastewater sequence datasets were screened for individual sequence reads that contained combinations of mutations frequently found in cryptic lineages but not contemporary circulating lineages. Using this method, we identified 18 cryptic lineages that appeared in multiple (2-81) samples from the same sewershed, including 12 that were not previously reported. Partial consensus sequences were generated for each cryptic lineage by extracting and mapping sequences containing cryptic-specific mutations. Surprisingly, seven of the mutations that appeared convergently in cryptic lineages were reversions to sequences that were highly conserved in SARS-CoV-2-related enteric bat Sarbecoviruses. The apparent reversion to bat Sarbecovirus sequences is consistent with the notion that SARS-CoV-2 adaptation to replicate efficiently in respiratory tissues preceded the COVID-19 pandemic.", "claims": [{"id": "claim-data-access", "statement": "To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:40489546", "locator": "sections[10] \"Supporting information\" paragraph 7", "quote": "To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:40489546", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Using this method, we identified 18 cryptic lineages that appeared in multiple (2–81) samples from the same sewershed, including 12 that were not previously reported.", "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:40489546"], "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:40489546", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1371/journal.ppat.1012850", "sha256": "sha256:95f635c338268ead8ecb8d4ed1d055c3aeb93a5cf1925369bfaf3fb8a920c227", "bytes": 35122, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1371/journal.ppat.1012850", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12176291", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}], "gaps": []}, "analysis": {"state": "absent", "methods": [], "targets": [], "gaps": ["method", "operations"]}, "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": "R4", "expressible": false, "required_operations": [], "missing_operations": [], "basis": "The reader quoted the method the paper names (not named), but the corpus's method extractor could not map it onto a DSL operation, so the method is published but not expressible."}, "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": ["To comply with submission guidelines, many low-abundance polymorphisms were excluded from the analysis. A comprehensive report is available upon request from the corresponding author."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:40489546"], "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:54:03Z", "independence": "single_reader"}, "evidence_check": {"method": null, "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 9285, "completion_tokens": 1658, "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:40489546", "type": "paper", "label": "pmid:40489546"}, {"id": "source:pmid:40489546", "type": "manuscript", "label": "PMID 40489546 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:40489546", "type": "paper", "label": "pmid:40489546"}, {"id": "source:pmid:40489546", "type": "manuscript", "label": "PMID 40489546 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-40489546", "type": "extraction", "used": ["source:pmid:40489546"], "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": null, "parent_hash": null, "fork_of": null, "generation": 0}, "gaps": ["execution.not_run", "method", "operations", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent"]}