Quantifying the heterogeneity and determinants of Ebola Zaire transmission during the 10th outbreak in DRC, 2018-2020.

Background. The 2018-2020 Ebola virus disease outbreak in the Democratic Republic of the Congo (DRC) was the country's largest, and the second largest globally, amid armed conflict and community mistrust. Transmission heterogeneity (superspreading) is recognised in Ebola epidemics, but empirical estimates of its extent and determinants remain scarce for DRC outbreaks. We quantified transmission heterogeneity and its determinants during this outbreak. Methods. In this retrospective observational study, we reconstructed transmission chains for confirmed and probable cases (Aug 1, 2018, to June 25, 2020) using routinely collected Ministry of Health and Medecins Sans Frontieres surveillance data. We modelled the offspring distribution with a Bayesian negative binomial framework, correcting for incomplete contact tracing, to estimate the effective reproduction number (Reff), dispersion parameter (k), and proportion of cases responsible for 80% of transmission (prop80), overall, by subgroup, and over time. Individual-level determinants were assessed with a regression extension, adjusting for covariates. Findings. Among 3481 cases, 2008 transmission events linked 2402 (69%) individuals into 415 chains (median size 3, range 2-102). Overall Reff was 1.00 (95% CI 0.92-1.08) with k 0.29 (0.26-0.32); 17.8% of cases generated 80% of transmission. Overdispersion stayed stable despite fluctuating Reff. Non-isolation (IRR 1.79), death outside a treatment centre (IRR 4.34), and unfollowed contact status (IRR up to 4.48) predicted more secondary cases; vaccination cut transmission by about 60%. Interpretation. Epidemiological investigations linked 67.7% (2356/3481) of cases into 415 transmission chains (median size 3, range 2-102); linkage to a known infector fell to 10% during the November 2018-February 2019 period of peak insecurity. Transmission was heterogeneous overall, with dispersion parameter k of 0.29 (0.26-0.32), such that 17.8% (16.8-18.9) of cases generated 80% of onward transmission confirming superspreading as a stable, structural feature of Ebola dynamics. Critically, k remained stable throughout the outbreak, including during periods of elevated Reff, indicating that transmission surges reflected intensification of the same underlying process rather than new superspreading contexts, and that Reff alone is an insufficient summary of epidemic potential. Regression analyses identified predominantly modifiable determinants: cases not isolated in an Ebola treatment centre (IRR 1.79 [1.54-2.06]) or who died outside one (IRR 4.34 [3.47-5.31]) generated substantially more secondary cases, as did those registered as contacts but not followed up (IRR 3.01 [2.39-3.70]) or unregistered altogether (IRR 4.48 [3.67-5.38]) relative to actively followed-up contacts. Vaccination reduced onward transmission by 60-64% (IRR 0.36-0.40). These findings indicate that transmission was shaped less by gaps in epidemiological knowledge than by the operational reach of contact tracing, isolation, and vaccination delivery, particularly during periods of insecurity.