Estimating the reproduction number and transmission heterogeneity from the size distribution of clusters of identical pathogen sequences

Quantifying transmission intensity and heterogeneity is crucial to ascertain the threat posed by infectious diseases and inform the design of interventions. Methods that jointly estimate the reproduction number R and the dispersion parameter k have however mainly remained limited to the analysis of epidemiological clusters or contact tracing data, whose collection often proves difficult. Here, we show that clusters of identical sequences are imprinted by the pathogen offspring distribution, and we derive an analytical formula for the distribution of the size of these clusters. We develop and evaluate an inference framework to jointly estimate the reproduction number and the dispersion parameter from the size distribution of clusters of identical sequences. We then illustrate its application across a range of epidemiological situations. Finally, we develop a hypothesis testing framework relying on clusters of identical sequences to determine whether a given pathogen genetic subpopulation is associated with increased or reduced transmissibility. Our work provides tools to estimate the reproduction number and transmission heterogeneity from pathogen sequences without building a phylogenetic tree, thus making it easily scalable to large pathogen genome datasets.

Supercritical and homogenous transmission of monkeypox in the capital of China

Starting from May 31, 2023, the local transmission of monkeypox (Mpox) in mainland China began in Beijing. Till now, the transmission characteristics have not been explored. Based on the daily Mpox incidence data in the first 3 weeks of Beijing (from May 31 to June 21, 2023), we employed the instant-individual heterogeneity transmission model to simultaneously calculate the effective reproduction number (Re ) and the degree of heterogeneity (k) of the Beijing epidemic. We additionally simulated the monthly infection size in Beijing from July to November and compared with the reported data to project subsequent transmission dynamics. We estimated Re to be 1.68 (95% highest posterior density [HPD]: 1.12-2.41), and k to be 2.57 [95% HPD: 0.54-83.88], suggesting the transmission of Mpox in Beijing was supercritical and didn't have considerable transmission heterogeneity. We projected that Re fell in the range of 0.95-1.0 from July to November, highlighting more efforts needed to further reduce the Mpox transmissibility. Our findings revealed supercritical and homogeneous transmission of the Mpox epidemic in Beijing. Our results could serve as a reference for understanding and predicting the ongoing Mpox transmission in other regions of China and evaluating the effect of control measures.

Inapparent infections shape the transmission heterogeneity of dengue

Transmission heterogeneity, whereby a disproportionate fraction of pathogen transmission events result from a small number of individuals or geographic locations, is an inherent property of many, if not most, infectious disease systems. For vector-borne diseases, transmission heterogeneity is inferred from the distribution of the number of vectors per host, which could lead to significant bias in situations where vector abundance and transmission risk at the household do not correlate, as is the case with dengue virus (DENV). We used data from a contact tracing study to quantify the distribution of DENV acute infections within human activity spaces (AS), the collection of residential locations an individual routinely visits, and quantified measures of virus transmission heterogeneity from two consecutive dengue outbreaks (DENV-4 and DENV-2) that occurred in the city of Iquitos, Peru. Negative-binomial distributions and Pareto fractions showed evidence of strong overdispersion in the number of DENV infections by AS and identified super-spreading units (SSUs): i.e. AS where most infections occurred. Approximately 8% of AS were identified as SSUs, contributing to more than 50% of DENV infections. SSU occurrence was associated more with DENV-2 infection than with DENV-4, a predominance of inapparent infections (74% of all infections), households with high Aedes aegypti mosquito abundance, and high host susceptibility to the circulating DENV serotype. Marked heterogeneity in dengue case distribution, and the role of inapparent infections in defining it, highlight major challenges faced by reactive interventions if those transmission units contributing the most to transmission are not identified, prioritized, and effectively treated.

Factors related to human-vector contact that modify the likelihood of malaria transmission during a contained Plasmodium falciparum outbreak in Praia, Cabo Verde

Background

Determining the reproductive rate and how it varies over time and space (RT) provides important insight to understand transmission of a given disease and inform optimal strategies for controlling or eliminating it. Estimating RT for malaria is difficult partly due to the widespread use of interventions and immunity to disease masking incident infections. A malaria outbreak in Praia, Cabo Verde in 2017 provided a unique opportunity to estimate RT directly, providing a proxy for the intensity of vector-human contact and measure the impact of vector control measures.

Methods

Out of 442 confirmed malaria cases reported in 2017 in Praia, 321 (73%) were geolocated and informed this analysis. RT was calculated using the joint likelihood of transmission between two cases, based on the time (serial interval) and physical distance (spatial interval) between them. Log-linear regression was used to estimate factors associated with changes in RT, including the impact of vector control interventions. A geostatistical model was developed to highlight areas receptive to transmission where vector control activities could be focused in future to prevent or interrupt transmission.

Results

The RT from individual cases ranged between 0 and 11 with a median serial- and spatial-interval of 34 days [interquartile range (IQR): 17-52] and 1,347 m (IQR: 832-1,985 m), respectively. The number of households receiving indoor residual spraying (IRS) 4 weeks prior was associated with a reduction in RT by 0.84 [95% confidence interval (CI) 0.80-0.89; p-value <0.001] in the peak-and post-epidemic compared to the pre-epidemic period.

Conclusions

Identifying the effect of reduced human-vector contact through IRS is essential to determining optimal intervention strategies that modify the likelihood of malaria transmission and can inform optimal intervention strategies to accelerate time to elimination. The distance within which two cases are plausibly linked is important for the potential scale of any reactive interventions as well as classifying infections as imported or introduced and confirming malaria elimination.

Model-based Analysis of Tuberculosis Genotype Clusters in the United States Reveals High Degree of Heterogeneity in Transmission and State-level Differences Across California, Florida, New York, and Texas

BACKGROUND

Reductions in tuberculosis (TB) transmission have been instrumental in lowering TB incidence in the United States. Sustaining and augmenting these reductions are key public health priorities.

METHODS

We fit mechanistic transmission models to distributions of genotype clusters of TB cases reported to the Centers for Disease Control and Prevention during 2012-2016 in the United States and separately in California, Florida, New York, and Texas. We estimated the mean number of secondary cases generated per infectious case (R0) and individual-level heterogeneity in R0 at state and national levels and assessed how different definitions of clustering affected these estimates.

RESULTS

In clusters of genotypically linked TB cases that occurred within a state over a 5-year period (reference scenario), the estimated R0 was 0.29 (95% confidence interval [CI], .28-.31) in the United States. Transmission was highly heterogeneous; 0.24% of simulated cases with individual R0 >10 generated 19% of all recent secondary transmissions. R0 estimate was 0.16 (95% CI, .15-.17) when a cluster was defined as cases occurring within the same county over a 3-year period. Transmission varied across states: estimated R0s were 0.34 (95% CI, .3-.4) in California, 0.28 (95% CI, .24-.36) in Florida, 0.19 (95% CI, .15-.27) in New York, and 0.38 (95% CI, .33-.46) in Texas.

CONCLUSIONS

TB transmission in the United States is characterized by pronounced heterogeneity at the individual and state levels. Improving detection of transmission clusters through incorporation of whole-genome sequencing and identifying the drivers of this heterogeneity will be essential to reducing TB transmission.

Superspreading potential of COVID-19 outbreak seeded by Omicron variants of SARS-CoV-2 in Hong Kong

Using two early transmission chains in Hong Kong, the estimated R and k were 1.34 (95%CrI: 0.94-2.19) and 0.33 (95%CrI: 0.17-0.62) respectively, inferring 20.3% (95%CrI: 12.7%-29.6%) cases were responsible for 80% of the transmissions of the Omicron epidemic. Compared with Omicron BA.1, Omicron BA.2 had a greater superspreading potential.

Apathogenic proxies for transmission dynamics of a fatal virus

Identifying drivers of transmission-especially of emerging pathogens-is a formidable challenge for proactive disease management efforts. While close social interactions can be associated with microbial sharing between individuals, and thereby imply dynamics important for transmission, such associations can be obscured by the influences of factors such as shared diets or environments. Directly-transmitted viral agents, specifically those that are rapidly evolving such as many RNA viruses, can allow for high-resolution inference of transmission, and therefore hold promise for elucidating not only which individuals transmit to each other, but also drivers of those transmission events. Here, we tested a novel approach in the Florida panther, which is affected by several directly-transmitted feline retroviruses. We first inferred the transmission network for an apathogenic, directly-transmitted retrovirus, feline immunodeficiency virus (FIV), and then used exponential random graph models to determine drivers structuring this network. We then evaluated the utility of these drivers in predicting transmission of the analogously transmitted, pathogenic agent, feline leukemia virus (FeLV), and compared FIV-based predictions of outbreak dynamics against empirical FeLV outbreak data. FIV transmission was primarily driven by panther age class and distances between panther home range centroids. FIV-based modeling predicted FeLV dynamics similarly to common modeling approaches, but with evidence that FIV-based predictions captured the spatial structuring of the observed FeLV outbreak. While FIV-based predictions of FeLV transmission performed only marginally better than standard approaches, our results highlight the value of proactively identifying drivers of transmission-even based on analogously-transmitted, apathogenic agents-in order to predict transmission of emerging infectious agents. The identification of underlying drivers of transmission, such as through our workflow here, therefore holds promise for improving predictions of pathogen transmission in novel host populations, and could provide new strategies for proactive pathogen management in human and animal systems.

SARS-CoV-2 Droplet and Airborne Transmission Heterogeneity

The spread dynamics of the SARS-CoV-2 virus have not yet been fully understood after two years of the pandemic. The virus's global spread represented a unique scenario for advancing infectious disease research. Consequently, mechanistic epidemiological theories were quickly dismissed, and more attention was paid to other approaches that considered heterogeneity in the spread. One of the most critical advances in aerial pathogens transmission was the global acceptance of the airborne model, where the airway is presented as the epicenter of the spread of the disease. Although the aerodynamics and persistence of the SARS-CoV-2 virus in the air have been extensively studied, the actual probability of contagion is still unknown. In this work, the individual heterogeneity in the transmission of 22 patients infected with COVID-19 was analyzed by close contact (cough samples) and air (environmental samples). Viral RNA was detected in 2/19 cough samples from patient subgroups, with a mean Ct (Cycle Threshold in Quantitative Polymerase Chain Reaction analysis) of 25.7 ± 7.0. Nevertheless, viral RNA was only detected in air samples from 1/8 patients, with an average Ct of 25.0 ± 4.0. Viral load in cough samples ranged from 7.3 × 105 to 8.7 × 108 copies/mL among patients, while concentrations between 1.1-4.8 copies/m3 were found in air, consistent with other reports in the literature. In patients undergoing follow-up, no viral load was found (neither in coughs nor in the air) after the third day of symptoms, which could help define quarantine periods in infected individuals. In addition, it was found that the patient's Ct should not be considered an indicator of infectiousness, since it could not be correlated with the viral load disseminated. The results of this work are in line with proposed hypotheses of superspreaders, which can attribute part of the heterogeneity of the spread to the oversized emission of a small percentage of infected people.

Genotype and sex-based host variation in behaviour and susceptibility drives population disease dynamics

Host heterogeneity in pathogen transmission is widespread and presents a major hurdle to predicting and minimizing disease outbreaks. Using Drosophila melanogaster infected with Drosophila C virus as a model system, we integrated experimental measurements of social aggregation, virus shedding, and disease-induced mortality from different genetic lines and sexes into a disease modelling framework. The experimentally measured host heterogeneity produced substantial differences in simulated disease outbreaks, providing evidence for genetic and sex-specific effects on disease dynamics at a population level. While this was true for homogeneous populations of single sex/genetic line, the genetic background or sex of the index case did not alter outbreak dynamics in simulated, heterogeneous populations. Finally, to explore the relative effects of social aggregation, viral shedding and mortality, we compared simulations where we allowed these traits to vary, as measured experimentally, to simulations where we constrained variation in these traits to the population mean. In this context, variation in infectiousness, followed by social aggregation, was the most influential component of transmission. Overall, we show that host heterogeneity in three host traits dramatically affects population-level transmission, but the relative impact of this variation depends on both the susceptible population diversity and the distribution of population-level variation.

Evaluating Transmission Heterogeneity and Super-Spreading Event of COVID-19 in a Metropolis of China

COVID-19 caused rapid mass infection worldwide. Understanding its transmission characteristics, including heterogeneity and the emergence of super spreading events (SSEs) where certain individuals infect large numbers of secondary cases, is of vital importance for prediction and intervention of future epidemics. Here, we collected information of all infected cases (135 cases) between 21 January and 26 February 2020 from official public sources in Tianjin, a metropolis of China, and grouped them into 43 transmission chains with the largest chain of 45 cases and the longest chain of four generations. Utilizing a heterogeneous transmission model based on branching process along with a negative binomial offspring distribution, we estimated the reproductive number R and the dispersion parameter k (lower value indicating higher heterogeneity) to be 0.67 (95% CI: 0.54-0.84) and 0.25 (95% CI: 0.13-0.88), respectively. A super-spreader causing six infections was identified in Tianjin. In addition, our simulation allowing for heterogeneity showed that the outbreak in Tianjin would have caused 165 infections and sustained for 7.56 generations on average if no control measures had been taken by local government since 28 January. Our results highlighted more efforts are needed to verify the transmission heterogeneity of COVID-19 in other populations and its contributing factors.

The primary case is not enough: Variation among individuals, groups and social networks modify bacterial transmission dynamics

The traits of the primary case of an infectious disease outbreak, and the circumstances for their aetiology, potentially influence the trajectory of transmission dynamics. However, these dynamics likely also depend on the traits of the individuals with whom the primary case interacts. We used the social spider Stegodyphus dumicola to test how the traits of the primary case, group phenotypic composition and group size interact to facilitate the transmission of a GFP-labelled cuticular bacterium. We also compared bacterial transmission across experimentally generated "daisy-chain" vs. "star" networks of social interactions. Finally, we compared social network structure across groups of different sizes. Groups of 10 spiders experienced more bacterial transmission events compared to groups of 30 spiders, regardless of groups' behavioural composition. Groups containing only one bold spider experienced the lowest levels of bacterial transmission regardless of group size. We found no evidence for the traits of the primary case influencing any transmission dynamics. In a second experiment, bacteria were transmitted to more individuals in experimentally induced star networks than in daisy-chains, on which transmission never exceeded three steps. In both experimental network types, transmission success depended jointly on the behavioural traits of the interacting individuals; however, the behavioural traits of the primary case were only important for transmission on star networks. Larger social groups exhibited lower interaction density (i.e. had a low ratio of observed to possible connections) and were more modular, i.e. they had more connections between nodes within a subgroup and fewer connections across subgroups. Thus, larger groups may restrict transmission by forming fewer interactions and by isolating subgroups that interacted with the primary case. These findings suggest that accounting for the traits of single exposed hosts has less power in predicting transmission dynamics compared to the larger scale factors of the social groups in which they reside. Factors like group size and phenotypic composition appear to alter social interaction patterns, which leads to differential transmission of microbes.

Combining contact tracing with targeted indoor residual spraying significantly reduces dengue transmission

The widespread transmission of dengue viruses (DENV), coupled with the alarming increase of birth defects and neurological disorders associated with Zika virus, has put the world in dire need of more efficacious tools for Aedes aegypti-borne disease mitigation. We quantitatively investigated the epidemiological value of location-based contact tracing (identifying potential out-of-home exposure locations by phone interviews) to infer transmission foci where high-quality insecticide applications can be targeted. Space-time statistical modeling of data from a large epidemic affecting Cairns, Australia, in 2008-2009 revealed a complex pattern of transmission driven primarily by human mobility (Cairns accounted for ~60% of virus transmission to and from residents of satellite towns, and 57% of all potential exposure locations were nonresidential). Targeted indoor residual spraying with insecticides in potential exposure locations reduced the probability of future DENV transmission by 86 to 96%, compared to unsprayed premises. Our findings provide strong evidence for the effectiveness of combining contact tracing with residual spraying within a developed urban center, and should be directly applicable to areas with similar characteristics (for example, southern USA, Europe, or Caribbean countries) that need to control localized Aedes-borne virus transmission or to protect pregnant women's homes in areas with active Zika transmission. Future theoretical and empirical research should focus on evaluation of the applicability and scalability of this approach to endemic areas with variable population size and force of DENV infection.

Coupled Heterogeneities and Their Impact on Parasite Transmission and Control

Most host-parasite systems exhibit remarkable heterogeneity in the contribution to transmission of certain individuals, locations, host infectious states, or parasite strains. While significant advancements have been made in the understanding of the impact of transmission heterogeneity in epidemic dynamics and parasite persistence and evolution, the knowledge base of the factors contributing to transmission heterogeneity is limited. We argue that research efforts should move beyond considering the impact of single sources of heterogeneity and account for complex couplings between conditions with potential synergistic impacts on parasite transmission. Using theoretical approaches and empirical evidence from various host-parasite systems, we investigate the ecological and epidemiological significance of couplings between heterogeneities and discuss their potential role in transmission dynamics and the impact of control.

Individual differences in boldness influence patterns of social interactions and the transmission of cuticular bacteria among group-mates

Despite the importance of host attributes for the likelihood of associated microbial transmission, individual variation is seldom considered in studies of wildlife disease. Here, we test the influence of host phenotypes on social network structure and the likelihood of cuticular bacterial transmission from exposed individuals to susceptible group-mates using female social spiders (Stegodyphus dumicola). Based on the interactions of resting individuals of known behavioural types, we assessed whether individuals assorted according to their behavioural traits. We found that individuals preferentially interacted with individuals of unlike behavioural phenotypes. We next applied a green fluorescent protein-transformed cuticular bacterium,Pantoeasp., to individuals and allowed them to interact with an unexposed colony-mate for 24 h. We found evidence for transmission of bacteria in 55% of cases. The likelihood of transmission was influenced jointly by the behavioural phenotypes of both the exposed and susceptible individuals: transmission was more likely when exposed spiders exhibited higher 'boldness' relative to their colony-mate, and when unexposed individuals were in better body condition. Indirect transmission via shared silk took place in only 15% of cases. Thus, bodily contact appears key to transmission in this system. These data represent a fundamental step towards understanding how individual traits influence larger-scale social and epidemiological dynamics.