Search

Background: Mixed infection with multiple strains of the same pathogen in a single host can present clinical and analytical challenges. Whole genome sequence (WGS) data can identify signals of multiple strains in samples, though the precision of previous methods can be improved. Here, we present MixInfect2, a new tool to accurately detect mixed samples from Mycobacterium tuberculosis short-read WGS data. We then evaluate three approaches for reconstructing the underlying mixed constituent strain sequences. This allows these samples to be included in downstream analysis to gain insights into the epidemiology and transmission of mixed infections., Methods: We employed a Gaussian mixture model to cluster allele frequencies at mixed sites (hSNPs) in each sample to identify signals of multiple strains. Building upon our previous tool, MixInfect, we increased the accuracy of classifying in vitro mixed samples through multiple improvements to the bioinformatic pipeline. Major and minor proportion constituent strains were reconstructed using three approaches and assessed by comparing the estimated sequence to the known constituent strain sequence. Lastly, mixed infections in a real-world Mycobacterium tuberculosis population from Moldova were detected with MixInfect2 and clusters of recent transmission that included major and minor constituent strains were built., Results: All 36/36 in vitro mixed and 12/12 non-mixed samples were correctly classified with MixInfect2, and major strain proportions were estimated with high accuracy (within 3% of the true strain proportion), outperforming previous tools. Reconstructed major strain sequences closely matched the true constituent sequence by taking the allele at the highest frequency at hSNPs, while the best-performing approach to reconstruct the minor proportion strain sequence was identifying the closest non-mixed isolate in the same population, though no approach was effective when the minor strain proportion was at 5%. Finally, fewer mixed infections were identified in Moldova than previous estimates (6.6% vs 17.4%) and we found multiple instances where the constituent strains of mixed samples were present in transmission clusters., Conclusions: MixInfect2 accurately detects samples with evidence of mixed infection from short-read WGS data and provides an excellent estimate of the mixture proportions. While there are limitations in reconstructing the constituent strain sequences of mixed samples, we present recommendations for the best approach to include these isolates in further analyses., Competing Interests: Declarations. Ethics approval and consent to participate: This study was approved by the Ethics Committee of Research of the Phthisiopneumology Institute in Moldova and the Yale University Human Investigation Committee (no. 2000023071). Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Background: Mycobacterium tuberculosis complex (MTBC) species evolve slowly, so isolates from individuals linked in transmission often have identical or nearly identical genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared MTBC variation could help overcome this problem. Previous studies have reported MTBC diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear. Here, we aimed to quantify within-host MTBC variation and assess whether such information improves transmission inferences., Methods: We conducted a retrospective genomic epidemiology study in which we reanalysed publicly available sequence data from household transmission studies published in PubMed from database inception until Jan 31, 2024, for which both genomic and epidemiological contact data were available, using household membership as a proxy for transmission linkage. We quantified minority variants (ie, positions with two or more alleles each supported by at least five-fold coverage and with a minor allele frequency of 1% or more) outside of PE and PPE genes, within individual samples and shared across samples. We used receiver operator characteristic (ROC) curves to compare the performance of a general linear model for household membership that included shared minority variants and one that included only fixed genetic differences., Findings: We identified three MTBC household transmission studies with publicly available whole-genome sequencing data and epidemiological linkages: a household transmission study in Vitória, Brazil (Colangeli et al), a retrospective population-based study of paediatric tuberculosis in British Columbia, Canada (Guthrie et al), and a retrospective population-based study in Oxfordshire, England (Walker et al). We found moderate levels of minority variation present in MTBC sequence data from cultured isolates that varied significantly across studies: mean 168·6 minority variants (95% CI 151·4-185·9) for the Colangeli et al dataset, 5·8 (1·5-10·2) for Guthrie et al (p<0·0001, Wilcoxon rank sum test, vs Colangeli et al), and 7·1 (2·4-11·9) for Walker et al (p<0·0001, Wilcoxon rank sum test, vs Colangeli et al). Isolates from household pairs shared more minority variants than did randomly selected pairs of isolates: mean 97·7 shared minority variants (79·1-116·3) versus 9·8 (8·6-11·0) in Colangeli et al, 0·8 (0·1-1·5) versus 0·2 (0·1-0·2) in Guthrie et al, and 0·7 (0·1-1·3) versus 0·2 (0·2-0·2) in Walker et al (all p<0·0001, Wilcoxon rank sum test). Shared within-host variation was significantly associated with household membership (odds ratio 1·51 [95% CI 1·30-1·71], p<0·0001), for one standard deviation increase in shared minority variants. Models that included shared within-host variation versus models without within-host variation improved the accuracy of predicting household membership in all three studies: area under the ROC curve 0·95 versus 0·92 for the Colangeli et al study, 0·99 versus 0·95 for the Guthrie et al study, and 0·93 versus 0·91 for the Walker et al study., Interpretation: Within-host MTBC variation persists through culture of sputum and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlight the need to optimise approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference., Funding: National Institutes of Health., Competing Interests: Declaration of interests JRA reports support from the US National Institutes of Health (NIH); receiving donated Cepheid disposable open cartridges from Cepheid to Stanford University for NIH-funded tuberculosis research; participation in two Data Safety Monitoring Boards for NIH-funded tuberculosis trials; and participation in the Science Advisory Board for an NIH-funded tuberculosis network. JC reports grants from the NIH, Valneva/Butantan Institution, Coalition for Epidemic Preparedness Innovations–Sabin Institute, MSD, and Sanofi Pasteur for clinical and epidemiological studies; honoraria from Pfizer; participation in the Brazil Advisory Board for the mRNA-1273 vaccine Moderna–Zodiac, the Latin America and Brazilian Advisory Board for paxlovid (nirmatrelvir–ritonavir; Pfizer), the Brazil Advisory Board for the Qdenga vaccine (Takeda), and the Takeda Global Dengue Steering Committee. All other authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Diversity plays an important role in various domains, including conservation, whether it describes diversity within a population or diversity over a set of species. While various strategies for measuring among-species diversity have emerged (e.g. Phylogenetic Diversity (PD), Split System Diversity (SSD) and entropy-based methods), extensions to populations are rare. An understudied problem is how to assess the diversity of a collection of populations where each has its own internal diversity. Relying solely on measures that treat each population as a monomorphic lineage (like a species) can be misleading. To address this problem, we present four population-level diversity assessment approaches: Pooling, Averaging, Pairwise Differencing, and Fixing. These approaches can be used to extend any diversity measure that is primarily defined for a group of individuals to a collection of populations. We then apply the approaches to two measures of diversity that have been used in conservation-Heterozygosity (Het) and Split System Diversity (SSD)-across a dataset comprising SNP data for 50 anadromous Atlantic salmon populations. We investigate agreement and disagreement between these measures of diversity when used to identify optimal sets of populations for conservation, on both the observed data, and randomized and simulated datasets. The similarity and differences of the maximum-diversity sets as well as the pairwise correlations among our proposed measures emphasize the need to clearly define what aspects of biodiversity we aim to both measure and optimize, to ensure meaningful and effective conservation decisions., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Abhari et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Plasmids play a major role in the transfer of antimicrobial resistance (AMR) genes among bacteria via horizontal gene transfer. The identification of plasmids in short-read assemblies is a challenging problem and a very active research area. Plasmid binning aims at detecting, in a draft genome assembly, groups (bins) of contigs likely to originate from the same plasmid. Several methods for plasmid binning have been developed recently, such as PlasBin-flow, HyAsP, gplas, MOB-suite, and plasmidSPAdes. This motivates the problem of evaluating the performances of plasmid binning methods, either against a given ground truth or between them., Results: We describe PlasEval, a novel method aimed at comparing the results of plasmid binning tools. PlasEval computes a dissimilarity measure between two sets of plasmid bins, that can originate either from two plasmid binning tools, or from a plasmid binning tool and a ground truth set of plasmid bins. The PlasEval dissimilarity accounts for the contig content of plasmid bins, the length of contigs and is repeat-aware. Moreover, the dissimilarity score computed by PlasEval is broken down into several parts, that allows to understand qualitative differences between the compared sets of plasmid bins. We illustrate the use of PlasEval by benchmarking four recently developed plasmid binning tools-PlasBin-flow, HyAsP, gplas, and MOB-recon-on a data set of 53 E. coli bacterial genomes., Conclusion: Analysis of the results of plasmid binning methods using PlasEval shows that their behaviour varies significantly. PlasEval can be used to decide which specific plasmid binning method should be used for a specific dataset. The disagreement between different methods also suggests that the problem of plasmid binning on short-read contigs requires further research. We believe that PlasEval can prove to be an effective tool in this regard. PlasEval is publicly available at https://github.com/acme92/PlasEval., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare that they have no competing interests., (© 2024. The Author(s).)

Identifying individuals with tuberculosis (TB) with a high risk of onward transmission can guide disease prevention and public health strategies. Here, we train classification models to predict the first sampled isolates in Mycobacterium tuberculosis transmission clusters from demographic and disease data. We find that supervised learning, in particular balanced random forests, can be used to develop predictive models to identify people with TB that are more likely associated with TB cluster growth, with good model performance and AUCs of ≥ 0.75. We also identified the most important patient and disease characteristics in the best performing classification model, including host demographics, site of infection, TB lineage, and age at diagnosis. This framework can be used to develop predictive tools for the early assessment of potential cluster growth to prioritise individuals for enhanced follow-up with the aim of reducing transmission chains., Competing Interests: Competing interests The authors declare no competing interests. Ethics declarations Ethics were obtained from the University of British Columbia (certificate H12-00910) and informed consent for participation in the study was not required, as determined by institutional REB review., (© 2024. The Author(s).)

Hypermutated proviruses, which arise in a single Human Immunodeficiency Virus (HIV) replication cycle when host antiviral APOBEC3 proteins introduce extensive guanine to adenine mutations throughout the viral genome, persist in all people living with HIV receiving antiretroviral therapy (ART). However, hypermutated sequences are routinely excluded from phylogenetic trees because their extensive mutations complicate phylogenetic inference, and as a result, we know relatively little about their within-host evolutionary origins and dynamics. Using >1400 longitudinal single-genome-amplified HIV env-gp120 sequences isolated from six women over a median of 18 years of follow-up-including plasma HIV RNA sequences collected over a median of 9 years between seroconversion and ART initiation, and >500 proviruses isolated over a median of 9 years on ART-we evaluated three approaches for masking hypermutation in nucleotide alignments. Our goals were to (i) reconstruct phylogenies that can be used for molecular dating and (ii) phylogenetically infer the integration dates of hypermutated proviruses persisting during ART. Two of the approaches (stripping all positions containing putative APOBEC3 mutations from the alignment or replacing individual putative APOBEC3 mutations in hypermutated sequences with the ambiguous base R) consistently normalized tree topologies, eliminated erroneous clustering of hypermutated proviruses, and brought env -intact and hypermutated proviruses into comparable ranges with respect to multiple tree-based metrics. Importantly, these corrected trees produced integration date estimates for env -intact proviruses that were highly concordant with those from benchmark trees that excluded hypermutated sequences, supporting the use of these corrected trees for molecular dating. Subsequent molecular dating of hypermutated proviruses revealed that these sequences spanned a wide within-host age range, with the oldest ones dating to shortly after infection. This indicates that hypermutated proviruses, like other provirus types, begin to be seeded into the proviral pool immediately following infection and can persist for decades. In two of the six participants, hypermutated proviruses differed from env -intact ones in terms of their age distributions, suggesting that different provirus types decay at heterogeneous rates in some hosts. These simple approaches to reconstruct hypermutated provirus' evolutionary histories reveal insights into their in vivo origins and longevity toward a more comprehensive understanding of HIV persistence during ART., Competing Interests: None declared., (© The Author(s) 2025. Published by Oxford University Press.)

We unified the recent literature with the goal to contribute to the discussion on how genetic diversity might best be conserved. We argue that this decision will be guided by how genomic variation is distributed among manageable populations (i.e., its spatial structure), the degree to which adaptive potential is best predicted by variation across the entire genome or the subset of that variation that is identified as putatively adaptive (i.e., its genomic structure), and whether we are managing species as single entities or as collections of diversifying lineages. The distribution of genetic variation and our ultimate goal will have practical implications for on-the-ground management. If adaptive variation is largely polygenic or responsive to change, its spatial structure might be broadly governed by the forces determining genome-wide variation (linked selection, drift, and gene flow), making measurement and prioritization straightforward. If we are managing species as single entities, then population-level prioritization schemes are possible so as to maximize future pooled genetic variation. We outline one such scheme based on the popular Shapley value from cooperative game theory that considers the relative genetic contribution of a population to an unknown future collection of populations., Competing Interests: The authors declare there are no competing interests.

Background: Understanding how infectious disease transmission varies from person to person, including associations with age and contact behavior, can help design effective control strategies. Within households, transmission may be highly variable because of differing transmission risks by age, household size, and individual contagiousness. Our aim was to disentangle those factors by fitting mathematical models to SARS-CoV-2 household survey and serologic data., Methods: We surveyed members of 3,381 Utah households from January-April 2021 and performed SARS-CoV-2 antibody testing on all available members. We paired these data with a probabilistic model of household importation and transmission composed of a novel combination of transmission variability and age- and size-structured heterogeneity. We calculated maximum likelihood estimates of mean and variability of household transmission probability between household members in different age groups and different household sizes, simultaneously with importation probability and probabilities of false negative and false positive test results., Results: 12.8% of individual participants, residing in 17.4% of the participating households, showed serologic evidence of prior infection or reported a prior positive test on the survey. Serologically positive individuals in younger age groups were less likely than older adults to have tested positive during their infection according to our survey results. Our model results suggested that adolescents and young adults (ages 13-24) acquired SARS-CoV-2 infection outside the household at a rate substantially higher than younger children and older adults. Our estimate of the household secondary attack rate (HSAR) among adults aged 45 and older exceeded HSARs to and/or from younger age groups at a given household size. We found lower HSAR in households with more members, independent of age differences. The age-specific HSAR patterns we found could not be explained by age-dependent biological susceptibility and transmissibility alone, suggesting that age groups contacted each other at different rates within households., Conclusions: We disentangled several factors contributing to age-specific infection risk, including non-household exposure, within-household exposure to specific age groups, and household size. Within-household contact rate differences played a significant role in driving household transmission epidemiology. These findings provide nuanced insights for understanding community outbreak patterns and mechanisms of differential infection risk., (© 2024. The Author(s).)

Background: The COVID-19 pandemic underlined the need for pandemic planning but also brought into focus the use of mathematical modelling to support public health decisions. The types of models needed (compartment, agent-based, importation) are described. Best practices regarding biological realism (including the need for multidisciplinary expert advisors to modellers), model complexity, consideration of uncertainty and communications to decision-makers and the public are outlined., Methods: A narrative review was developed from the experiences of COVID-19 by members of the Public Health Agency of Canada External Modelling Network for Infectious Diseases (PHAC EMN-ID), a national community of practice on mathematical modelling of infectious diseases for public health., Results: Modelling can best support pandemic preparedness in two ways: 1) by modelling to support decisions on resource needs for likely future pandemics by estimating numbers of infections, hospitalized cases and cases needing intensive care, associated with epidemics of "hypothetical-yet-plausible" pandemic pathogens in Canada; and 2) by having ready-to-go modelling methods that can be readily adapted to the features of an emerging pandemic pathogen and used for long-range forecasting of the epidemic in Canada, as well as to explore scenarios to support public health decisions on the use of interventions., Conclusion: There is a need for modelling expertise within public health organizations in Canada, linked to modellers in academia in a community of practice, within which relationships built outside of times of crisis can be applied to enhance modelling during public health emergencies. Key challenges to modelling for pandemic preparedness include the availability of linked public health, hospital and genomic data in Canada.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform the public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This portal has been coupled with other resources, such as Viral AI, and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this portal (https://virusseq-dataportal.ca/), including its contextual data not available elsewhere, and the Duotang (https://covarr-net.github.io/duotang/duotang.html), a web platform that presents key genomic epidemiology and modelling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the portal (COVID-MVP, CoVizu), are all open source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Clostridioides difficile infection (CDI) is a significant public health threat, associated with antibiotic-induced disruption of the normally protective gastrointestinal microbiota. CDI is thought to occur in two stages: acquisition of asymptomatic colonization from ingesting C. difficile bacteria followed by progression to symptomatic CDI caused by toxins produced during C. difficile overgrowth. The degree to which disruptive antibiotic exposure increases susceptibility at each stage is uncertain, which might contribute to divergent published projections of the impact of hospital antibiotic stewardship interventions on CDI. Here, we model C. difficile transmission and CDI among hospital inpatients, including exposure to high-CDI-risk antibiotics and their effects on each stage of CDI epidemiology. We derive the mathematical relationship, using a deterministic model, between those parameters and observed equilibrium levels of colonization, CDI, and risk ratio of CDI among certain antibiotic-exposed patients relative to patients with no recent antibiotic exposure. We then quantify the sensitivity of projected antibiotic stewardship intervention impacts to alternate assumptions. We find that two key parameters, the antibiotic effects on susceptibility to colonization and to CDI progression, are not identifiable given the data frequently available. Furthermore, the effects of antibiotic stewardship interventions are sensitive to their assumed values. Thus, discrepancies between different projections of antibiotic stewardship interventions may be largely due to model assumptions. Data supporting improved quantification of mechanistic antibiotic effects on CDI epidemiology are needed to understand stewardship effects better., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Mitchell et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Infectious disease dynamics are driven by the complex interplay of epidemiological, ecological, and evolutionary processes. Accurately modeling these interactions is crucial for understanding pathogen spread and informing public health strategies. However, existing simulators often fail to capture the dynamic interplay between these processes, resulting in oversimplified models that do not fully reflect real-world complexities in which the pathogen's genetic evolution dynamically influences disease transmission. We introduce the epidemiological-ecological-evolutionary simulator (e3SIM), an open-source framework that concurrently models the transmission dynamics and molecular evolution of pathogens within a host population while integrating environmental factors. Using an agent-based, discrete-generation, forward-in-time approach, e3SIM incorporates compartmental models, host-population contact networks, and quantitative-trait models for pathogens. This integration allows for realistic simulations of disease spread and pathogen evolution. Key features include a modular and scalable design, flexibility in modeling various epidemiological and population-genetic complexities, incorporation of time-varying environmental factors, and a user-friendly graphical interface. We demonstrate e3SIM's capabilities through simulations of realistic outbreak scenarios with SARS-CoV-2 and Mycobacterium tuberculosis , illustrating its flexibility for studying the genomic epidemiology of diverse pathogen types.

COVID-19 has become endemic, with dynamics that reflect the waning of immunity and re-exposure, by contrast to the epidemic phase driven by exposure in immunologically naïve populations. Endemic does not, however, mean constant. Further evolution of SARS-CoV-2, as well as changes in behavior and public health policy, continue to play a major role in the endemic load of disease and mortality. In this article, we analyze evolutionary models to explore the impact that a newly arising variant can have on the short-term and longer-term endemic load, characterizing how these impacts depend on the transmission and immunological properties of the variants. We describe how evolutionary changes in the virus will increase the endemic load most for a persistently immune-escape variant, by an intermediate amount for a more transmissible variant, and least for a transiently immune-escape variant. Balancing the tendency for evolution to favor variants that increase the endemic load, we explore the impact of vaccination strategies and non-pharmaceutical interventions that can counter these increases in the impact of disease. We end with some open questions about the future of COVID-19 as an endemic disease., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses., Competing Interests: J.T.S. receives research funding from Oxford Nanopore Technologies (ONT) and has received travel support to attend and speak at meetings organized by ONT, and is on the Scientific Advisory Board of Day Zero Diagnostics.

Background: Spatial analysis at different levels can help understand spatial variation of human immunodeficiency virus (HIV) infection, disease drivers, and targeted interventions. Combining spatial analysis and the evaluation of the determinants of the HIV burden in Southern African countries is essential for a better understanding of the disease dynamics in high-burden settings., Methods: The study countries were selected based on the availability of demographic and health surveys (DHS) and corresponding geographic coordinates. We used multivariable regression to evaluate the determinants of HIV burden and assessed the presence and nature of HIV spatial autocorrelation in six Southern African countries., Results: The overall prevalence of HIV for each country varied between 11.3% in Zambia and 22.4% in South Africa. The HIV prevalence rate was higher among female respondents in all six countries. There were reductions in prevalence estimates in most countries yearly from 2011 to 2020. The hotspot cluster findings show that the major cities in each country are the key sites of high HIV burden. Compared with female respondents, the odds of being HIV positive were lesser among the male respondents. The probability of HIV infection was higher among those who had sexually transmitted infections (STI) in the last 12 months, divorced and widowed individuals, and women aged 25 years and older., Conclusions: Our research findings show that analysis of survey data could provide reasonable estimates of the wide-ranging spatial structure of the HIV epidemic in Southern African countries. Key determinants such as individuals who are divorced, middle-aged women, and people who recently treated STIs, should be the focus of HIV prevention and control interventions. The spatial distribution of high-burden areas for HIV in the selected countries was more pronounced in the major cities. Interventions should also be focused on locations identified as hotspot clusters., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Adetokunboh, Are. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The projected trajectory of multidrug resistant tuberculosis (MDR-TB) epidemics depends on the reproductive fitness of circulating strains of MDR M. tuberculosis (Mtb). Previous efforts to characterize the fitness of MDR Mtb have found that Mtb strains of the Beijing sublineage (Lineage 2.2.1) may be more prone to develop resistance and retain fitness in the presence of resistance-conferring mutations than other lineages. Using Mtb genome sequences from all culture-positive cases collected over two years in Moldova, we estimate the fitness of Ural (Lineage 4.2) and Beijing strains, the two lineages in which MDR is concentrated in the country. We estimate that the fitness of MDR Ural strains substantially exceeds that of other susceptible and MDR strains, and we identify several mutations specific to these MDR Ural strains. Our findings suggest that MDR Ural Mtb has been transmitting efficiently in Moldova and poses a substantial risk of spreading further in the region., (© 2024. The Author(s).)

Celebrities can generate substantial attention and influence public interest in species. Using a large-scale examination of publicly available data, we assessed whether species across 6 taxonomic groups received more page views on Wikipedia when the species was named after a celebrity than when it was not. We conducted our analysis for 4 increasingly strict thresholds of how many average daily Wikipedia page views a celebrity had (1, 10, 100, or 1000 views). Overall, we found a high probability (0.96-0.98) that species named after celebrities had more page views than their closest relatives that were not named after celebrities, irrespective of the celebrity threshold. The multiplicative effect on species' page views was larger but more uncertain as celebrity page-view thresholds increased. The range for thresholds of 1 to 1000 was 1.08 (95% credible interval [CI] 1.00-1.18) to 1.76 (95% CI 0.96-2.80), respectively. The hierarchical estimates for the taxa tended to be positive. The strongest effects were for invertebrates, followed by amphibians, reptiles, fish, and mammals, whereas the weakest effect was for birds at lower page-view thresholds. Our results suggest that naming species after celebrities could be particularly significant for those belonging to taxonomic groups that are generally less popular than others (e.g., invertebrates). Celebrities may further influence the effectiveness of this marketing strategy, depending on their likability and connection to the species named after them. Eponyms may serve as a reminder of the disproportionate power dynamics between populations and some namesakes' untenable actions. However, they also provide an opportunity to recognize remarkable individuals and promote equity, inclusivity, and diversity in taxonomic practice. We encourage taxonomists to examine whether naming threatened species after celebrities could affect conservation support, especially for species that are otherwise typically overlooked by the public., (© 2023 The Authors. Conservation Biology published by Wiley Periodicals LLC on behalf of Society for Conservation Biology.)

The recovery of top predators is thought to have cascading effects on vegetated ecosystems and their geomorphology 1,2 , but the evidence for this remains correlational and intensely debated 3,4 . Here we combine observational and experimental data to reveal that recolonization of sea otters in a US estuary generates a trophic cascade that facilitates coastal wetland plant biomass and suppresses the erosion of marsh edges-a process that otherwise leads to the severe loss of habitats and ecosystem services 5,6 . Monitoring of the Elkhorn Slough estuary over several decades suggested top-down control in the system, because the erosion of salt marsh edges has generally slowed with increasing sea otter abundance, despite the consistently increasing physical stress in the system (that is, nutrient loading, sea-level rise and tidal scour 7-9 ). Predator-exclusion experiments in five marsh creeks revealed that sea otters suppress the abundance of burrowing crabs, a top-down effect that cascades to both increase marsh edge strength and reduce marsh erosion. Multi-creek surveys comparing marsh creeks pre- and post-sea otter colonization confirmed the presence of an interaction between the keystone sea otter, burrowing crabs and marsh creeks, demonstrating the spatial generality of predator control of ecosystem edge processes: densities of burrowing crabs and edge erosion have declined markedly in creeks that have high levels of sea otter recolonization. These results show that trophic downgrading could be a strong but underappreciated contributor to the loss of coastal wetlands, and suggest that restoring top predators can help to re-establish geomorphic stability., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Motivation: Transcriptomic long-read (LR) sequencing is an increasingly cost-effective technology for probing various RNA features. Numerous tools have been developed to tackle various transcriptomic sequencing tasks (e.g. isoform and gene fusion detection). However, the lack of abundant gold-standard datasets hinders the benchmarking of such tools. Therefore, the simulation of LR sequencing is an important and practical alternative. While the existing LR simulators aim to imitate the sequencing machine noise and to target specific library protocols, they lack some important library preparation steps (e.g. PCR) and are difficult to modify to new and changing library preparation techniques (e.g. single-cell LRs)., Results: We present TKSM, a modular and scalable LR simulator, designed so that each RNA modification step is targeted explicitly by a specific module. This allows the user to assemble a simulation pipeline as a combination of TKSM modules to emulate a specific sequencing design. Additionally, the input/output of all the core modules of TKSM follows the same simple format (Molecule Description Format) allowing the user to easily extend TKSM with new modules targeting new library preparation steps., Availability and Implementation: TKSM is available as an open source software at https://github.com/vpc-ccg/tksm., (© The Author(s) 2024. Published by Oxford University Press.)

The mechanisms behind vaccine-induced strain replacement in the pneumococcus remain poorly understood. There is emerging evidence that distinct pneumococcal lineages can co-colonise for significant time periods, and that novel recombinants can readily emerge during natural colonisation. Despite this, patterns of post-vaccine replacement are indicative of competition between specific lineages. Here, we develop a multiscale transmission model to investigate explicitly how within host dynamics shape observed ecological patterns, both pre- and post-vaccination. Our model framework explores competition between and within strains defined by distinct antigenic, metabolic and resistance profiles. We allow for strains to freely co-colonise and recombine within hosts, and consider how each of these types may contribute to a strain's overall fitness. Our results suggest that antigenic and resistance profiles are key drivers of post-vaccine success., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Mulberry et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Hyperparasites (species which parasitize other parasites) are common in natural populations, affecting many parasitic taxa, including: eukaryotic parasites; bacterial and fungal pathogens. Hyperparasitism is therefore likely to shape the ecology and evolution of many host-parasite systems, representing a promising method for biocontrol (e.g., treating antimicrobial resistant infections). However, the eco-evolutionary consequences of hyperparasitism have received little attention. We use a host-parasite-hyperparasite model to explore how introducing a hyperparasite drives the evolution of parasite virulence, and what impact this has on the host population. We show when the introduction of a hyperparasite selects for higher or lower parasite virulence, and the changes in virulence experienced by the host population. Crucially, we show that variation in the direct effects of hyperparasites on virulence and transmission, and the probability of cotransmission, can lead to a previously unseen hysteresis effect, whereby small shifts in hyperparasite characteristics can lead to sudden shifts in parasite virulence. We also show that hyperparasites can induce diversification in parasite virulence, leading to the coexistence of high and low virulence strains. Our results show hyperparasites can have dramatic effects on the evolution of parasite virulence, and that the use of hyperparasites in biocontrol should be approached with caution., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

Background: Outbreak response modelling often involves collaboration among academics, and experts from governmental and non-governmental organizations. We conducted a systematic review of modelling studies on human vaccine-preventable disease (VPD) outbreaks to identify patterns in modelling practices between two collaboration types. We complemented this with a mini comparison of foot-and-mouth disease (FMD), a veterinary disease that is controllable by vaccination., Methods: We searched three databases for modelling studies that assessed the impact of an outbreak response. We extracted data on author affiliation type (academic institution, governmental, and non-governmental organizations), location studied, and whether at least one author was affiliated to the studied location. We also extracted the outcomes and interventions studied, and model characteristics. Included studies were grouped into two collaboration types: purely academic (papers with only academic affiliations), and mixed (all other combinations) to help investigate differences in modelling patterns between collaboration types in the human disease literature and overall differences with FMD collaboration practices., Results: Human VPDs formed 227 of 252 included studies. Purely academic collaborations dominated the human disease studies (56%). Notably, mixed collaborations increased in the last seven years (2013-2019). Most studies had an author affiliated to an institution in the country studied (75.2%) but this was more likely among the mixed collaborations. Contrasted to the human VPDs, mixed collaborations dominated the FMD literature (56%). Furthermore, FMD studies more often had an author with an affiliation to the country studied (92%) and used complex model design, including stochasticity, and model parametrization and validation., Conclusion: The increase in mixed collaboration studies over the past seven years could suggest an increase in the uptake of modelling for outbreak response decision-making. We encourage more mixed collaborations between academic and non-academic institutions and the involvement of locally affiliated authors to help ensure that the studies suit local contexts., Competing Interests: Declaration of competing interest The authors report no declarations of interest., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

The serial interval of an infectious disease is an important variable in epidemiology. It is defined as the period of time between the symptom onset times of the infector and infectee in a direct transmission pair. Under partially sampled data, purported infector-infectee pairs may actually be separated by one or more unsampled cases in between. Misunderstanding such pairs as direct transmissions will result in overestimating the length of serial intervals. On the other hand, two cases that are infected by an unseen third case (known as coprimary transmission) may be classified as a direct transmission pair, leading to an underestimation of the serial interval. Here, we introduce a method to jointly estimate the distribution of serial intervals factoring in these two sources of error. We simultaneously estimate the distribution of the number of unsampled intermediate cases between purported infector-infectee pairs, as well as the fraction of such pairs that are coprimary. We also extend our method to situations where each infectee has multiple possible infectors, and show how to factor this additional source of uncertainty into our estimates. We assess our method's performance on simulated data sets and find that our method provides consistent and robust estimates. We also apply our method to data from real-life outbreaks of four infectious diseases and compare our results with published results. With similar accuracy, our method of estimating serial interval distribution provides unique advantages, allowing its application in settings of low sampling rates and large population sizes, such as widespread community transmission tracked by routine public health surveillance., Competing Interests: Declaration of competing interest The authors have declared no competing interests., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

This paper deals with the second quantization of interacting relativistic Fermionic and Bosonic fields in the arena of discrete phase space and continuous time. The mathematical formulation involves partial difference equations. The corresponding Feynman diagrams and a new [Formula: see text]-matrix theory is developed. In the special case of proton-proton Møller scattering via an exchange of a neutral meson, the explicit second order element [Formula: see text] is deduced. In the approximation of very low external three-momenta, a new Yukawa potential is explicitly derived from [Formula: see text]. Moreover, it is rigorously proved that this new Yukawa potential is divergence-free. The mass parameter of the exchanged meson may be set to zero to obtain a type of scalar Boson exchange between hypothetical Fermions. This provides a limiting case of a new Coulomb type potential directly from the new singularity free Yukawa potential. A divergence-free Coulomb potential between two Fermions at two discrete points is shown to be proportional to the Euler beta function. Within this relativistic discrete phase space continuous time, a single quanta is shown to occupy the hyper-tori [Formula: see text] where [Formula: see text] is a circle of radius [Formula: see text]., (© 2023. The Author(s).)

The seasonal influenza (flu) vaccine is designed to protect against those influenza viruses predicted to circulate during the upcoming flu season, but identifying which viruses are likely to circulate is challenging. We use features from phylogenetic trees reconstructed from hemagglutinin (HA) and neuraminidase (NA) sequences, together with a support vector machine, to predict future circulation. We obtain accuracies of 0.75 to 0.89 (AUC 0.83 to 0.91) over 2016-2020. We explore ways to select potential candidates for a seasonal vaccine and find that the machine learning model has a moderate ability to select strains that are close to future populations. However, consensus sequences among the most recent 3 years also do well at this task. We identify similar candidate strains to those proposed by the World Health Organization, suggesting that this approach can help inform vaccine strain selection.

Identification of plasmids from sequencing data is an important and challenging problem related to antimicrobial resistance spread and other One-Health issues. We provide a new architecture for identifying plasmid contigs in fragmented genome assemblies built from short-read data. We employ graph neural networks (GNNs) and the assembly graph to propagate the information from nearby nodes, which leads to more accurate classification, especially for short contigs that are difficult to classify based on sequence features or database searches alone. We trained plASgraph2 on a data set of samples from the ESKAPEE group of pathogens. plASgraph2 either outperforms or performs on par with a wide range of state-of-the-art methods on testing sets of independent ESKAPEE samples and samples from related pathogens. On one hand, our study provides a new accurate and easy to use tool for contig classification in bacterial isolates; on the other hand, it serves as a proof-of-concept for the use of GNNs in genomics. Our software is available at https://github.com/cchauve/plasgraph2 and the training and testing data sets are available at https://github.com/fmfi-compbio/plasgraph2-datasets., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Sielemann, Sielemann, Brejová, Vinař and Chauve.)

Objectives: Clustering pathogen sequence data is a common practice in epidemiology to gain insights into the genetic diversity and evolutionary relationships among pathogens. We can find groups of cases with a shared transmission history and common origin, as well as identifying transmission hotspots. Motivated by the experience of clustering SARS-CoV-2 cases using whole genome sequence data during the COVID-19 pandemic to aid with public health investigation, we investigated how differences in epidemiology and sampling can influence the composition of clusters that are identified., Methods: We performed genomic clustering on simulated SARS-CoV-2 outbreaks produced with different transmission rates and levels of genomic diversity, along with varying the proportion of cases sampled., Results: In single outbreaks with a low transmission rate, decreasing the sampling fraction resulted in multiple, separate clusters being identified where intermediate cases in transmission chains are missed. Outbreaks simulated with a high transmission rate were more robust to changes in the sampling fraction and largely resulted in a single cluster that included all sampled outbreak cases. When considering multiple outbreaks in a sampled jurisdiction seeded by different introductions, low genomic diversity between introduced cases caused outbreaks to be merged into large clusters. If the transmission and sampling fraction, and diversity between introductions was low, a combination of the spurious break-up of outbreaks and the linking of closely related cases in different outbreaks resulted in clusters that may appear informative, but these did not reflect the true underlying population structure. Conversely, genomic clusters matched the true population structure when there was relatively high diversity between introductions and a high transmission rate., Conclusion: Differences in epidemiology and sampling can impact our ability to identify genomic clusters that describe the underlying population structure. These findings can help to guide recommendations for the use of pathogen clustering in public health investigations., Competing Interests: Declaration of Competing Interest The authors declare they have no competing interests., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Serial intervals - the time between symptom onset in infector and infectee - are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals' exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2-3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities., (© 2023. Springer Nature Limited.)

The reconstruction of phylogenetic networks is an important but challenging problem in phylogenetics and genome evolution, as the space of phylogenetic networks is vast and cannot be sampled well. One approach to the problem is to solve the minimum phylogenetic network problem, in which phylogenetic trees are first inferred, and then the smallest phylogenetic network that displays all the trees is computed. The approach takes advantage of the fact that the theory of phylogenetic trees is mature, and there are excellent tools available for inferring phylogenetic trees from a large number of biomolecular sequences. A tree-child network is a phylogenetic network satisfying the condition that every nonleaf node has at least one child that is of indegree one. Here, we develop a new method that infers the minimum tree-child network by aligning lineage taxon strings in the phylogenetic trees. This algorithmic innovation enables us to get around the limitations of the existing programs for phylogenetic network inference. Our new program, named ALTS, is fast enough to infer a tree-child network with a large number of reticulations for a set of up to 50 phylogenetic trees with 50 taxa that have only trivial common clusters in about a quarter of an hour on average., (© 2023 Zhang et al.; Published by Cold Spring Harbor Laboratory Press.)

Motivation: The analysis of bacterial isolates to detect plasmids is important due to their role in the propagation of antimicrobial resistance. In short-read sequence assemblies, both plasmids and bacterial chromosomes are typically split into several contigs of various lengths, making identification of plasmids a challenging problem. In plasmid contig binning, the goal is to distinguish short-read assembly contigs based on their origin into plasmid and chromosomal contigs and subsequently sort plasmid contigs into bins, each bin corresponding to a single plasmid. Previous works on this problem consist of de novo approaches and reference-based approaches. De novo methods rely on contig features such as length, circularity, read coverage, or GC content. Reference-based approaches compare contigs to databases of known plasmids or plasmid markers from finished bacterial genomes., Results: Recent developments suggest that leveraging information contained in the assembly graph improves the accuracy of plasmid binning. We present PlasBin-flow, a hybrid method that defines contig bins as subgraphs of the assembly graph. PlasBin-flow identifies such plasmid subgraphs through a mixed integer linear programming model that relies on the concept of network flow to account for sequencing coverage, while also accounting for the presence of plasmid genes and the GC content that often distinguishes plasmids from chromosomes. We demonstrate the performance of PlasBin-flow on a real dataset of bacterial samples., Availability and Implementation: https://github.com/cchauve/PlasBin-flow., (© The Author(s) 2023. Published by Oxford University Press.)

The problem of computing the Elementary Flux Modes (EFMs) and Minimal Cut Sets (MCSs) of metabolic network is a fundamental one in metabolic networks. A key insight is that they can be understood as a dual pair of monotone Boolean functions (MBFs). Using this insight, this computation reduces to the question of generating from an oracle a dual pair of MBFs. If one of the two sets (functions) is known, then the other can be computed through a process known as dualization. Fredman and Khachiyan provided two algorithms, which they called simply A and B that can serve as an engine for oracle-based generation or dualization of MBFs. We look at efficiencies available in implementing their algorithm B, which we will refer to as FK -B. Like their algorithm A, FK -B certifies whether two given MBFs in the form of Conjunctive Normal Form and Disjunctive Normal Form are dual or not, and in case of not being dual it returns a conflicting assignment (CA), that is, an assignment that makes one of the given Boolean functions True and the other one False . The FK -B algorithm is a recursive algorithm that searches through the tree of assignments to find a CA. If it does not find any CA, it means that the given Boolean functions are dual. In this article, we propose six techniques applicable to the FK -B and hence to the dualization process. Although these techniques do not reduce the time complexity, they considerably reduce the running time in practice. We evaluate the proposed improvements by applying them to compute the MCSs from the EFMs in the 19 small- and medium-sized models from the BioModels database along with 4 models of biomass synthesis in Escherichia coli that were used in an earlier computational survey Haus et al. (2008).

Defensive symbionts in the host microbiome can confer protection from infection or reduce the harms of being infected by a parasite. Defensive symbionts are therefore promising agents of biocontrol that could be used to control or ameliorate the impact of infectious diseases. Previous theory has shown how symbionts can evolve along the parasitism-mutualism continuum to confer greater or lesser protection to their hosts and in turn how hosts may coevolve with their symbionts to potentially form a mutualistic relationship. However, the consequences of introducing a defensive symbiont for parasite evolution and how the symbiont may coevolve with the parasite have received relatively little theoretical attention. Here, we investigate the ecological and evolutionary implications of introducing a tolerance-conferring defensive symbiont into an established host-parasite system. We show that while the defensive symbiont may initially have a positive impact on the host population, parasite and symbiont evolution tend to have a net negative effect on the host population in the long term. This is because the introduction of the defensive symbiont always selects for an increase in parasite virulence and may cause diversification into high- and low-virulence strains. Even if the symbiont experiences selection for greater host protection, this simply increases selection for virulence in the parasite, resulting in a net negative effect on the host population. Our results therefore suggest that tolerance-conferring defensive symbionts may be poor biocontrol agents for population-level infectious disease control., Competing Interests: The authors declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE) and European Society for Evolutionary Biology (ESEN).)

Theoretical models of the evolution of parasites and their hosts have shaped our understanding of infectious disease dynamics for over 40 years. Many theoretical models assume that the underlying ecological dynamics are at equilibrium or constant, yet we know that in a great many systems there are fluctuations in the ecological dynamics owing to a variety of intrinsic or extrinsic factors. Here, we discuss the challenges presented when modelling evolution in systems with fluctuating ecological dynamics and summarize the main approaches that have been developed to study host-parasite evolution in such systems. We provide an in-depth guide to one of the methods by applying it to two worked examples of host evolution that have not previously been studied in the literature: when cycles occur owing to seasonal forcing in competition, and when the presence of a free-living parasite causes cycles, with accompanying interactive Python code provided. We review the findings of studies that have explored host-parasite evolution when ecological dynamics fluctuate, and point to areas of future research. Throughout we stress the importance of feedbacks between the ecological and evolutionary dynamics in driving the outcomes of infectious disease systems. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

In the management of infectious disease outbreaks, grouping cases into clusters and understanding their underlying epidemiology are fundamental tasks. In genomic epidemiology, clusters are typically identified either using pathogen sequences alone or with sequences in combination with epidemiological data such as location and time of collection. However, it may not be feasible to culture and sequence all pathogen isolates, so sequence data may not be available for all cases. This presents challenges for identifying clusters and understanding epidemiology, because these cases may be important for transmission. Demographic, clinical and location data are likely to be available for unsequenced cases, and comprise partial information about their clustering. Here, we use statistical modelling to assign unsequenced cases to clusters already identified by genomic methods, assuming that a more direct method of linking individuals, such as contact tracing, is not available. We build our model on pairwise similarity between cases to predict whether cases cluster together, in contrast to using individual case data to predict the cases' clusters. We then develop methods that allow us to determine whether a pair of unsequenced cases are likely to cluster together, to group them into their most probable clusters, to identify which are most likely to be members of a specific (known) cluster, and to estimate the true size of a known cluster given a set of unsequenced cases. We apply our method to tuberculosis data from Valencia, Spain. Among other applications, we find that clustering can be predicted successfully using spatial distance between cases and whether nationality is the same. We can identify the correct cluster for an unsequenced case, among 38 possible clusters, with an accuracy of approximately 35 %, higher than both direct multinomial regression (17 %) and random selection (< 5 %).

Alternative splicing (AS) is an important mechanism in the development of many cancers, as novel or aberrant AS patterns play an important role as an independent onco-driver. In addition, cancer-specific AS is potentially an effective target of personalized cancer therapeutics. However, detecting AS events remains a challenging task, especially if these AS events are novel. This is exacerbated by the fact that existing transcriptome annotation databases are far from being comprehensive, especially with regard to cancer-specific AS. Additionally, traditional sequencing technologies are severely limited by the short length of the generated reads, which rarely spans more than a single splice junction site. Given these challenges, transcriptomic long-read (LR) sequencing presents a promising potential for the detection and discovery of AS. We present Freddie, a computational annotation-independent isoform discovery and detection tool. Freddie takes as input transcriptomic LR sequencing of a sample alongside its genomic split alignment and computes a set of isoforms for the given sample. It then partitions the input reads into sets that can be processed independently and in parallel. For each partition, Freddie segments the genomic alignment of the reads into canonical exon segments. The goal of this segmentation is to be able to represent any potential isoform as a subset of these canonical exons. This segmentation is formulated as an optimization problem and is solved with a dynamic programming algorithm. Then, Freddie reconstructs the isoforms by jointly clustering and error-correcting the reads using the canonical segmentation as a succinct representation. The clustering and error-correcting step is formulated as an optimization problem-the Minimum Error Clustering into Isoforms (MErCi) problem-and is solved using integer linear programming (ILP). We compare the performance of Freddie on simulated datasets with other isoform detection tools with varying dependence on annotation databases. We show that Freddie outperforms the other tools in its accuracy, including those given the complete ground truth annotation. We also run Freddie on a transcriptomic LR dataset generated in-house from a prostate cancer cell line with a matched short-read RNA-seq dataset. Freddie results in isoforms with a higher short-read cross-validation rate than the other tested tools. Freddie is open source and available at https://github.com/vpc-ccg/freddie/., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Innate, infection-preventing resistance often varies between host life stages. Juveniles are more resistant than adults in some species, whereas the opposite pattern is true in others. This variation cannot always be explained by prior exposure or physiological constraints and so it has been hypothesized that trade-offs with other life-history traits may be involved. However, little is known about how trade-offs between various life-history traits and resistance at different life stages affect the evolution of age-specific resistance. Here, we use a mathematical model to explore how trade-offs with natural mortality, reproduction and maturation combine to affect the evolution of resistance at different life stages. Our results show that certain combinations of trade-offs have substantial effects on whether adults or juveniles are more resistant, with trade-offs between juvenile resistance and adult reproduction inherently more costly than trade-offs involving maturation or mortality (all else being equal), resulting in consistent evolution of lower resistance at the juvenile stage even when infection causes a lifelong fecundity reduction. Our model demonstrates how the differences between patterns of age-structured resistance seen in nature may be explained by variation in the trade-offs involved and our results suggest conditions under which trade-offs tend to select for lower resistance in juveniles than adults.

Non-pharmaceutical interventions (NPIs), such as social distancing and contact tracing, are important public health measures that can reduce pathogen transmission. In addition to playing a crucial role in suppressing transmission, NPIs influence pathogen evolution by mediating mutation supply, restricting the availability of susceptible hosts, and altering the strength of selection for novel variants. Yet it is unclear how NPIs might affect the emergence of novel variants that are able to escape pre-existing immunity (partially or fully), are more transmissible or cause greater mortality. We analyse a stochastic two-strain epidemiological model to determine how the strength and timing of NPIs affect the emergence of variants with similar or contrasting life-history characteristics to the wild type. We show that, while stronger and timelier NPIs generally reduce the likelihood of variant emergence, it is possible for more transmissible variants with high cross-immunity to have a greater probability of emerging at intermediate levels of NPIs. This is because intermediate levels of NPIs allow an epidemic of the wild type that is neither too small (facilitating high mutation supply), nor too large (leaving a large pool of susceptible hosts), to prevent a novel variant from becoming established in the host population. However, since one cannot predict the characteristics of a variant, the best strategy to prevent emergence is likely to be an implementation of strong, timely NPIs., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Foundation for Evolution, Medicine, and Public Health.)

Objectives/aims: Prolonged infections of immunocompromised individuals have been proposed as a crucial source of new variants of SARS-CoV-2 during the COVID-19 pandemic. In principle, sustained within-host antigenic evolution in immunocompromised hosts could allow novel immune escape variants to emerge more rapidly, but little is known about how and when immunocompromised hosts play a critical role in pathogen evolution., Materials and Methods: Here, we use a simple mathematical model to understand the effects of immunocompromised hosts on the emergence of immune escape variants in the presence and absence of epistasis., Conclusions: We show that when the pathogen does not have to cross a fitness valley for immune escape to occur (no epistasis), immunocompromised individuals have no qualitative effect on antigenic evolution (although they may accelerate immune escape if within-host evolutionary dynamics are faster in immunocompromised individuals). But if a fitness valley exists between immune escape variants at the between-host level (epistasis), then persistent infections of immunocompromised individuals allow mutations to accumulate, therefore, facilitating rather than simply speeding up antigenic evolution. Our results suggest that better genomic surveillance of infected immunocompromised individuals and better global health equality, including improving access to vaccines and treatments for individuals who are immunocompromised (especially in lower- and middle-income countries), may be crucial to preventing the emergence of future immune escape variants of SARS-CoV-2., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Foundation for Evolution, Medicine, and Public Health.)

The role of schools in the spread of SARS-CoV-2 is controversial, with some claiming they are an important driver of the pandemic and others arguing that transmission in schools is negligible. School cluster reports that have been collected in various jurisdictions are a source of data about transmission in schools. These reports consist of the name of a school, a date, and the number of students known to be infected. We provide a simple model for the frequency and size of clusters in this data, based on random arrivals of index cases at schools who then infect their classmates with a highly variable rate, fitting the overdispersion evident in the data. We fit our model to reports from four Canadian provinces, providing estimates of mean and dispersion for cluster size, as well as the distribution of the instantaneous transmission parameter β , whilst factoring in imperfect ascertainment. According to our model with parameters estimated from the data, in all four provinces (i) more than 65% of non-index cases occur in the 20% largest clusters, and (ii) reducing instantaneous transmission rate and the number of contacts a student has at any given time are effective in reducing the total number of cases, whereas strict bubbling (keeping contacts consistent over time) does not contribute much to reduce cluster sizes. We predict strict bubbling to be more valuable in scenarios with substantially higher transmission rates., Competing Interests: PT, SP No competing interests declared, CC Reviewing editor, eLife, (© 2022, Tupper et al.)

Background: The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada., Results: We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data., Conclusions: Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination., (© 2022. The Author(s).)

While the negative effects that pathogens have on their hosts are well-documented in humans and agricultural systems, direct evidence of pathogen-driven impacts in wild host populations is scarce and mixed. Here, to determine how the strength of pathogen-imposed selection depends on spatial structure, we analyze growth rates across approximately 4000 host populations of a perennial plant through time coupled with data on pathogen presence-absence. We find that infection decreases growth more in the isolated than well-connected host populations. Our inoculation study reveals isolated populations to be highly susceptible to disease while connected host populations support the highest levels of resistance diversity, regardless of their disease history. A spatial eco-evolutionary model predicts that non-linearity in the costs to resistance may be critical in determining this pattern. Overall, evolutionary feedbacks define the ecological impacts of disease in spatially structured systems with host gene flow being more important than disease history in determining the outcome., (© 2022. The Author(s).)

The field of genomic epidemiology is rapidly growing as many jurisdictions begin to deploy whole-genome sequencing (WGS) in their national or regional pathogen surveillance programmes. WGS data offer a rich view of the shared ancestry of a set of taxa, typically visualized with phylogenetic trees illustrating the clusters or subtypes present in a group of taxa, their relatedness and the extent of diversification within and between them. When methicillin-resistant Staphylococcus aureus (MRSA) arose and disseminated widely, phylogenetic trees of MRSA-containing types of S. aureus had a distinctive 'comet' shape, with a 'comet head' of recently adapted drug-resistant isolates in the context of a 'comet tail' that was predominantly drug-sensitive. Placing an S. aureus isolate in the context of such a 'comet' helped public health laboratories interpret local data within the broader setting of S. aureus evolution. In this work, we ask what other tree shapes, analogous to the MRSA comet, are present in bacterial WGS datasets. We extract trees from large bacterial genomic datasets, visualize them as images and cluster the images. We find nine major groups of tree images, including the 'comets', star-like phylogenies, 'barbell' phylogenies and other shapes, and comment on the evolutionary and epidemiological stories these shapes might illustrate. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Genomic epidemiology plays an ever-increasing role in our understanding of and response to the spread of infectious pathogens. Phylogeography, the reconstruction of the historical location and movement of pathogens from the evolutionary relationships among sampled pathogen sequences, can inform policy decisions related to viral movement among jurisdictions. However, phylogeographic reconstruction is impacted by the fact that the sampling and virus sequencing policies differ among jurisdictions, and these differences can cause bias in phylogeographic reconstructions. Here we assess the potential impacts of geographic-based sampling bias on estimated viral locations in the past, and on whether key viral movements can be detected. We quantify the effect of bias using simulated phylogenies with known geographic histories, and determine the impact of the biased sampling and of the underlying migration rate on the accuracy of estimated past viral locations. We find that overall, the accuracy of phylogeographic reconstruction is high, particularly when the migration rate is low. However, results depend on sampling, and sampling bias can have a large impact on the numbers and nature of estimated migration events. We apply these insights to the geographic spread of Ebolavirus in the 2014-2016 West Africa epidemic. This work highlights how sampling policy can both impact geographic inference and be optimized to best ensure the accuracy of specific features of geographic spread., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

An essential consideration in urban transit facility planning is service efficiency and accessibility. Previous research has shown that reducing the number of facilities along a route may increase efficiency but decrease accessibility. Striking a balance between these two is a critical consideration in transit planning. Transit facility consolidation is a cost-effective way to improve the quality of service by strategically determining the desirable allocation of a limited number of facilities. This paper develops an optimization framework that integrates Geographical Information systems (GIS), decision-making analysis, and quantum technologies for addressing the problem of facility consolidation. Our proposed framework includes a novel mathematical model that captures non-linear interactions between facilities and surrounding demand nodes, inter-facility competition, ridership demand and spatial coverage. The developed model can harness the power of quantum effects such as superposition and quantum tunnelling and enables transportation planners to utilize the most recent hardware solutions such as quantum and digital annealers, coherent Ising Machines and gate-based universal quantum computers. This study presents a real-world application of the framework to the public transit facility redundancy problem in the British Columbia Vancouver metropolitan area. We demonstrate the effectiveness of our framework by reducing the number of facilities by 40% while maintaining the same service accessibility. Additionally, we showcase the ability of the proposed mathematical model to take advantage of quantum annealing and classical optimization techniques., Competing Interests: The authors have declared that no competing interests exist.

Background: COVID-19 vaccination is a key public health measure in the pandemic response. The rapid evolution of SARS-CoV-2 variants introduce new groups of spike protein mutations. These new mutations are thought to aid in the evasion of vaccine-induced immunity and render vaccines less effective. However, not all spike mutations contribute equally to vaccine escape. Previous studies associate mutations with vaccine breakthrough infections (BTI), but information at the population level remains scarce. We aimed to identify spike mutations associated with SARS-CoV-2 vaccine BTI in a community setting during the emergence and predominance of the Delta-variant., Methods: This case-control study used both genomic, and epidemiological data from a provincial COVID-19 surveillance program. Analyses were stratified into two periods approximating the emergence and predominance of the Delta-variant, and restricted to primary SARS-CoV-2 infections from either unvaccinated individuals, or those infected ≥14 days after their second vaccination dose in a community setting. Each sample's spike mutations were concatenated into a unique spike mutation profile (SMP). Penalized logistic regression was used to identify spike mutations and SMPs associated with SARS-CoV-2 vaccine BTI in both time periods., Results and Discussion: This study reports population level relative risk estimates, between 2 and 4-folds, of spike mutation profiles associated with BTI during the emergence and predominance of the Delta-variant, which comprised 19,624 and 17,331 observations, respectively. The identified mutations cover multiple spike domains including the N-terminal domain (NTD), receptor binding domain (RBD), S1/S2 cleavage region, fusion peptide and heptad regions. Mutations in these different regions imply various mechanisms contribute to vaccine escape. Our profiling method identifies naturally occurring spike mutations associated with BTI, and can be applied to emerging SARS-CoV-2 variants with novel groups of spike mutations., Competing Interests: NJ participates in Abbvie Advisory Board meeting. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Fibke, Joffres, Tyson, Colijn, Janjua, Fjell, Prystajecky, Jassem and Sbihi.)

Single-cell RNA sequencing allows for characterizing the gene expression landscape at the cell type level. However, because of its use of short-reads, it is severely limited at detecting full-length features of transcripts such as alternative splicing. New library preparation techniques attempt to extend single-cell sequencing by utilizing both long-reads and short-reads. These techniques split the library material, after it is tagged with cellular barcodes, into two pools: one for short-read sequencing and one for long-read sequencing. However, the challenge of utilizing these techniques is that they require matching the cellular barcodes sequenced by the erroneous long-reads to the cellular barcodes detected by the short-reads. To overcome this challenge, we introduce scTagger, a computational method to match cellular barcodes data from long-reads and short-reads. We tested scTagger against another state-of-the-art tool on both real and simulated datasets, and we demonstrate that scTagger has both significantly better accuracy and time efficiency., Competing Interests: The authors declare no competing interests., (© 2022 The Authors.)