Inferring transmission trees of infectious disease outbreaks with multiple introductions and transmission routes
Bibliographic record
Abstract
Like a detective in a murder mystery seeks the culprit, researchers studying infectious disease transmission aim to uncover who infected whom. Tracking links between infected individuals and identifying transmission events is essential for understanding many aspects of an outbreak. To solve this question, epidemiologists use a toolkit that includes two classes of methods: identifying transmission events through contact information and identifying them by tracing the evolution of the pathogen within the population. In densely sampled outbreaks, where genetic sequences are available for almost all cases, methods have been developed that infer transmission events from genetic and epidemiological information. Yet the assumptions behind these methods can limit their use across different outbreak settings. This thesis focuses on how infection of cases is represented in transmission models. One part examines the role of multiple introductions into a population. Another part studies the role of contact in transmission. Many current models assume a single introduction and treat all contact in the same way, which reduces their value for real outbreak analysis. To address these issues, this work expands the phybreak model to allow multiple introductions and to incorporate contact information, improving both practical use and accuracy. Multiple introductions are common in infectious disease outbreaks, yet many inference methods assume only one index case. To resolve this, I extended the transmission inference model in the R package phybreak with a “history host”, an artificial source for the index cases. This allows the entire outbreak to be represented in one phylogenetic tree. This approach successfully identified repeated introductions of SARS CoV 2 into mink farms (Chapter 2) and showed that ongoing introductions were a key driver of the 2022 Mpox outbreak in Slovenia (Chapter 3). The method also performed well during real time analysis, where strongly supported introductions matched those later confirmed. Recognizing repeated introductions is important for effective public health action, since outside sources often sustain an outbreak. Contact information is another essential source for understanding transmission, especially when used together with whole genome sequencing. In studies of Mycobacterium tuberculosis, cases are commonly grouped using SNP distance cutoffs, but these depend in part on the quality of contact tracing. In Chapter 4, phylodynamic methods were used to infer transmission events in the Netherlands. This supported the use of common SNP thresholds as a first step in identifying possible transmission pairs. Chapter 5 introduced another extension of phybreak that estimates the contribution of different types of contact to transmission. When applied to SARS CoV 2 in Dutch mink farms, the model showed that shared personnel accounted for most transmission whenever present, while other contact types were less influential. The phybreak model is best suited for outbreaks with relatively few cases, nearly complete observation, and genetic sequences for most individuals. It captures transmission, observation, within host processes, and mutation, and is particularly useful when multiple introductions occur. Although challenges remain, the extensions presented here strengthen the use of phybreak for reconstruction of transmission events and support informed public health decisions.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.003 | 0.014 |
| Meta-epidemiology (narrow) | 0.001 | 0.001 |
| Meta-epidemiology (broad) | 0.001 | 0.002 |
| Bibliometrics | 0.003 | 0.002 |
| Science and technology studies | 0.001 | 0.001 |
| Scholarly communication | 0.002 | 0.003 |
| Open science | 0.001 | 0.001 |
| Research integrity | 0.001 | 0.002 |
| Insufficient payload (model declined to judge) | 0.003 | 0.001 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".