Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England.

The SARS-CoV-2 lineage B.1.1.7, designated variant of concern (VOC) 202012/01 by Public Health England1, was first identified in the UK in late summer to early autumn 20202. Whole-genome SARS-CoV-2 sequence data collected from community-based diagnostic testing for COVID-19 show an extremely rapid expansion of the B.1.1.7 lineage during autumn 2020, suggesting that it has a selective advantage. Here we show that changes in VOC frequency inferred from genetic data correspond closely to changes inferred by S gene target failures (SGTF) in community-based diagnostic PCR testing. Analysis of trends in SGTF and non-SGTF case numbers in local areas across England shows that B.1.1.7 has higher transmissibility than non-VOC lineages, even if it has a different latent period or generation time. The SGTF data indicate a transient shift in the age composition of reported cases, with cases of B.1.1.7 including a larger share of under 20-year-olds than non-VOC cases. We estimated time-varying reproduction numbers for B.1.1.7 and co-circulating lineages using SGTF and genomic data. The best-supported models did not indicate a substantial difference in VOC transmissibility among different age groups, but all analyses agreed that B.1.1.7 has a substantial transmission advantage over other lineages, with a 50% to 100% higher reproduction number.

Estimating fine age structure and time trends in human contact patterns from coarse contact data: The Bayesian rate consistency model.

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), large-scale social contact surveys are now longitudinally measuring the fundamental changes in human interactions in the face of the pandemic and non-pharmaceutical interventions. Here, we present a model-based Bayesian approach that can reconstruct contact patterns at 1-year resolution even when the age of the contacts is reported coarsely by 5 or 10-year age bands. This innovation is rooted in population-level consistency constraints in how contacts between groups must add up, which prompts us to call the approach presented here the Bayesian rate consistency model. The model can also quantify time trends and adjust for reporting fatigue emerging in longitudinal surveys through the use of computationally efficient Hilbert Space Gaussian process priors. We illustrate estimation accuracy on simulated data as well as social contact data from Europe and Africa for which the exact age of contacts is reported, and then apply the model to social contact data with coarse information on the age of contacts that were collected in Germany during the COVID-19 pandemic from April to June 2020 across five longitudinal survey waves. We estimate the fine age structure in social contacts during the early stages of the pandemic and demonstrate that social contact intensities rebounded in an age-structured, non-homogeneous manner. The Bayesian rate consistency model provides a model-based, non-parametric, computationally tractable approach for estimating the fine structure and longitudinal trends in social contacts and is applicable to contemporary survey data with coarsely reported age of contacts as long as the exact age of survey participants is reported.

Inferring HIV transmission patterns from viral deep-sequence data via latent typed point processes.

Viral deep-sequencing data play a crucial role toward understanding disease transmission network flows, providing higher resolution compared to standard Sanger sequencing. To more fully utilize these rich data and account for the uncertainties in outcomes from phylogenetic analyses, we propose a spatial Poisson process model to uncover human immunodeficiency virus (HIV) transmission flow patterns at the population level. We represent pairings of individuals with viral sequence data as typed points, with coordinates representing covariates such as gender and age and point types representing the unobserved transmission statuses (linkage and direction). Points are associated with observed scores on the strength of evidence for each transmission status that are obtained through standard deep-sequence phylogenetic analysis. Our method is able to jointly infer the latent transmission statuses for all pairings and the transmission flow surface on the source-recipient covariate space. In contrast to existing methods, our framework does not require preclassification of the transmission statuses of data points, and instead learns them probabilistically through a fully Bayesian inference scheme. By directly modeling continuous spatial processes with smooth densities, our method enjoys significant computational advantages compared to previous methods that rely on discretization of the covariate space. We demonstrate that our framework can capture age structures in HIV transmission at high resolution, bringing valuable insights in a case study on viral deep-sequencing data from Southern Uganda.

Global minimum estimates of children affected by COVID-19-associated orphanhood and deaths of caregivers: a modelling study.

BACKGROUND: The COVID-19 pandemic priorities have focused on prevention, detection, and response. Beyond morbidity and mortality, pandemics carry secondary impacts, such as children orphaned or bereft of their caregivers. Such children often face adverse consequences, including poverty, abuse, and institutionalisation. We provide estimates for the magnitude of this problem resulting from COVID-19 and describe the need for resource allocation. METHODS: We used mortality and fertility data to model minimum estimates and rates of COVID-19-associated deaths of primary or secondary caregivers for children younger than 18 years in 21 countries. We considered parents and custodial grandparents as primary caregivers, and co-residing grandparents or older kin (aged 60-84 years) as secondary caregivers. To avoid overcounting, we adjusted for possible clustering of deaths using an estimated secondary attack rate and age-specific infection-fatality ratios for SARS-CoV-2. We used these estimates to model global extrapolations for the number of children who have experienced COVID-19-associated deaths of primary and secondary caregivers. FINDINGS: Globally, from March 1, 2020, to April 30, 2021, we estimate 1 134 000 children (95% credible interval 884 000-1 185 000) experienced the death of primary caregivers, including at least one parent or custodial grandparent. 1 562 000 children (1 299 000-1 683 000) experienced the death of at least one primary or secondary caregiver. Countries in our study set with primary caregiver death rates of at least one per 1000 children included Peru (10·2 per 1000 children), South Africa (5·1), Mexico (3·5), Brazil (2·4), Colombia (2·3), Iran (1·7), the USA (1·5), Argentina (1·1), and Russia (1·0). Numbers of children orphaned exceeded numbers of deaths among those aged 15-50 years. Between two and five times more children had deceased fathers than deceased mothers. INTERPRETATION: Orphanhood and caregiver deaths are a hidden pandemic resulting from COVID-19-associated deaths. Accelerating equitable vaccine delivery is key to prevention. Psychosocial and economic support can help families to nurture children bereft of caregivers and help to ensure that institutionalisation is avoided. These data show the need for an additional pillar of our response: prevent, detect, respond, and care for children. FUNDING: UK Research and Innovation (Global Challenges Research Fund, Engineering and Physical Sciences Research Council, Medical Research Council), UK National Institute for Health Research, US National Institutes of Health, and Imperial College London.

Phylogenies from dynamic networks.

The relationship between the underlying contact network over which a pathogen spreads and the pathogen phylogenetic trees that are obtained presents an opportunity to use sequence data to learn about contact networks that are difficult to study empirically. However, this relationship is not explicitly known and is usually studied in simulations, often with the simplifying assumption that the contact network is static in time, though human contact networks are dynamic. We simulate pathogen phylogenetic trees on dynamic Erdos-Renyi random networks and on two dynamic networks with skewed degree distribution, of which one is additionally clustered. We use tree shape features to explore how adding dynamics changes the relationships between the overall network structure and phylogenies. Our tree features include the number of small substructures (cherries, pitchforks) in the trees, measures of tree imbalance (Sackin index, Colless index), features derived from network science (diameter, closeness), as well as features using the internal branch lengths from the tip to the root. Using principal component analysis we find that the network dynamics influence the shapes of phylogenies, as does the network type. We also compare dynamic and time-integrated static networks. We find, in particular, that static network models like the widely used Barabasi-Albert model can be poor approximations for dynamic networks. We explore the effects of mis-specifying the network on the performance of classifiers trained identify the transmission rate (using supervised learning methods). We find that both mis-specification of the underlying network and its parameters (mean degree, turnover rate) have a strong adverse effect on the ability to estimate the transmission parameter. We illustrate these results by classifying HIV trees with a classifier that we trained on simulated trees from different networks, infection rates and turnover rates. Our results point to the importance of correctly estimating and modelling contact networks with dynamics when using phylodynamic tools to estimate epidemiological parameters.

Improved Algorithmic Complexity for the 3SEQ Recombination Detection Algorithm.

Identifying recombinant sequences in an era of large genomic databases is challenging as it requires an efficient algorithm to identify candidate recombinants and parents, as well as appropriate statistical methods to correct for the large number of comparisons performed. In 2007, a computation was introduced for an exact nonparametric mosaicism statistic that gave high-precision P values for putative recombinants. This exact computation meant that multiple-comparisons corrected P values also had high precision, which is crucial when performing millions or billions of tests in large databases. Here, we introduce an improvement to the algorithmic complexity of this computation from O(mn3) to O(mn2), where m and n are the numbers of recombination-informative sites in the candidate recombinant. This new computation allows for recombination analysis to be performed in alignments with thousands of polymorphic sites. Benchmark runs are presented on viral genome sequence alignments, new features are introduced, and applications outside recombination analysis are discussed.

PHYLOSCANNER: Inferring Transmission from Within- and Between-Host Pathogen Genetic Diversity.

A central feature of pathogen genomics is that different infectious particles (virions and bacterial cells) within an infected individual may be genetically distinct, with patterns of relatedness among infectious particles being the result of both within-host evolution and transmission from one host to the next. Here, we present a new software tool, phyloscanner, which analyses pathogen diversity from multiple infected hosts. phyloscanner provides unprecedented resolution into the transmission process, allowing inference of the direction of transmission from sequence data alone. Multiply infected individuals are also identified, as they harbor subpopulations of infectious particles that are not connected by within-host evolution, except where recombinant types emerge. Low-level contamination is flagged and removed. We illustrate phyloscanner on both viral and bacterial pathogens, namely HIV-1 sequenced on Illumina and Roche 454 platforms, HCV sequenced with the Oxford Nanopore MinION platform, and Streptococcus pneumoniae with sequences from multiple colonies per individual. phyloscanner is available from https://github.com/BDI-pathogens/phyloscanner.

Molecular Epidemiology of HIV-1 Subtype B Reveals Heterogeneous Transmission Risk: Implications for Intervention and Control.

Background: The impact of HIV pre-exposure prophylaxis (PrEP) depends on infections averted by protecting vulnerable individuals as well as infections averted by preventing transmission by those who would have been infected if not receiving PrEP. Analysis of HIV phylogenies reveals risk factors for transmission, which we examine as potential criteria for allocating PrEP. Methods: We analyzed 6912 HIV-1 partial pol sequences from men who have sex with men (MSM) in the United Kingdom combined with global reference sequences and patient-level metadata. Population genetic models were developed that adjust for stage of infection, global migration of HIV lineages, and changing incidence of infection through time. Models were extended to simulate the effects of providing susceptible MSM with PrEP. Results: We found that young age <25 years confers higher risk of HIV transmission (relative risk = 2.52 [95% confidence interval, 2.32-2.73]) and that young MSM are more likely to transmit to one another than expected by chance. Simulated interventions indicate that 4-fold more infections can be averted over 5 years by focusing PrEP on young MSM. Conclusions: Concentrating PrEP doses on young individuals can avert more infections than random allocation.

Phylogenetic Tools for Generalized HIV-1 Epidemics: Findings from the PANGEA-HIV Methods Comparison.

Viral phylogenetic methods contribute to understanding how HIV spreads in populations, and thereby help guide the design of prevention interventions. So far, most analyses have been applied to well-sampled concentrated HIV-1 epidemics in wealthy countries. To direct the use of phylogenetic tools to where the impact of HIV-1 is greatest, the Phylogenetics And Networks for Generalized HIV Epidemics in Africa (PANGEA-HIV) consortium generates full-genome viral sequences from across sub-Saharan Africa. Analyzing these data presents new challenges, since epidemics are principally driven by heterosexual transmission and a smaller fraction of cases is sampled. Here, we show that viral phylogenetic tools can be adapted and used to estimate epidemiological quantities of central importance to HIV-1 prevention in sub-Saharan Africa. We used a community-wide methods comparison exercise on simulated data, where participants were blinded to the true dynamics they were inferring. Two distinct simulations captured generalized HIV-1 epidemics, before and after a large community-level intervention that reduced infection levels. Five research groups participated. Structured coalescent modeling approaches were most successful: phylogenetic estimates of HIV-1 incidence, incidence reductions, and the proportion of transmissions from individuals in their first 3 months of infection correlated with the true values (Pearson correlation > 90%), with small bias. However, on some simulations, true values were markedly outside reported confidence or credibility intervals. The blinded comparison revealed current limits and strengths in using HIV phylogenetics in challenging settings, provided benchmarks for future methods' development, and supports using the latest generation of phylogenetic tools to advance HIV surveillance and prevention.

A Method to Estimate the Size and Characteristics of HIV-positive Populations Using an Individual-based Stochastic Simulation Model.

It is important not only to collect epidemiologic data on HIV but to also fully utilize such information to understand the epidemic over time and to help inform and monitor the impact of policies and interventions. We describe and apply a novel method to estimate the size and characteristics of HIV-positive populations. The method was applied to data on men who have sex with men living in the UK and to a pseudo dataset to assess performance for different data availability. The individual-based simulation model was calibrated using an approximate Bayesian computation-based approach. In 2013, 48,310 (90% plausibility range: 39,900-45,560) men who have sex with men were estimated to be living with HIV in the UK, of whom 10,400 (6,160-17,350) were undiagnosed. There were an estimated 3,210 (1,730-5,350) infections per year on average between 2010 and 2013. Sixty-two percent of the total HIV-positive population are thought to have viral load <500 copies/ml. In the pseudo-epidemic example, HIV estimates have narrower plausibility ranges and are closer to the true number, the greater the data availability to calibrate the model. We demonstrate that our method can be applied to settings with less data, however plausibility ranges for estimates will be wider to reflect greater uncertainty of the data used to fit the model.

Dispersion of the HIV-1 Epidemic in Men Who Have Sex with Men in the Netherlands: A Combined Mathematical Model and Phylogenetic Analysis.

BACKGROUND: The HIV-1 subtype B epidemic amongst men who have sex with men (MSM) is resurgent in many countries despite the widespread use of effective combination antiretroviral therapy (cART). In this combined mathematical and phylogenetic study of observational data, we aimed to find out the extent to which the resurgent epidemic is the result of newly introduced strains or of growth of already circulating strains. METHODS AND FINDINGS: As of November 2011, the ATHENA observational HIV cohort of all patients in care in the Netherlands since 1996 included HIV-1 subtype B polymerase sequences from 5,852 patients. Patients who were diagnosed between 1981 and 1995 were included in the cohort if they were still alive in 1996. The ten most similar sequences to each ATHENA sequence were selected from the Los Alamos HIV Sequence Database, and a phylogenetic tree was created of a total of 8,320 sequences. Large transmission clusters that included ≥10 ATHENA sequences were selected, with a local support value ≥ 0.9 and median pairwise patristic distance below the fifth percentile of distances in the whole tree. Time-varying reproduction numbers of the large MSM-majority clusters were estimated through mathematical modeling. We identified 106 large transmission clusters, including 3,061 (52%) ATHENA and 652 Los Alamos sequences. Half of the HIV sequences from MSM registered in the cohort in the Netherlands (2,128 of 4,288) were included in 91 large MSM-majority clusters. Strikingly, at least 54 (59%) of these 91 MSM-majority clusters were already circulating before 1996, when cART was introduced, and have persisted to the present. Overall, 1,226 (35%) of the 3,460 diagnoses among MSM since 1996 were found in these 54 long-standing clusters. The reproduction numbers of all large MSM-majority clusters were around the epidemic threshold value of one over the whole study period. A tendency towards higher numbers was visible in recent years, especially in the more recently introduced clusters. The mean age of MSM at diagnosis increased by 0.45 years/year within clusters, but new clusters appeared with lower mean age. Major strengths of this study are the high proportion of HIV-positive MSM with a sequence in this study and the combined application of phylogenetic and modeling approaches. Main limitations are the assumption that the sampled population is representative of the overall HIV-positive population and the assumption that the diagnosis interval distribution is similar between clusters. CONCLUSIONS: The resurgent HIV epidemic amongst MSM in the Netherlands is driven by several large, persistent, self-sustaining, and, in many cases, growing sub-epidemics shifting towards new generations of MSM. Many of the sub-epidemics have been present since the early epidemic, to which new sub-epidemics are being added.

Impact of cross-protective vaccines on epidemiological and evolutionary dynamics of influenza.

Large-scale immunization has profoundly impacted control of many infectious diseases such as measles and smallpox because of the ability of vaccination campaigns to maintain long-term herd immunity and, hence, indirect protection of the unvaccinated. In the case of human influenza, such potential benefits of mass vaccination have so far proved elusive. The central difficulty is a considerable viral capacity for immune escape; new pandemic variants, as well as viral escape mutants in seasonal influenza, compromise the buildup of herd immunity from natural infection or deployment of current vaccines. Consequently, most current influenza vaccination programs focus mainly on protection of specific risk groups, rather than mass prophylactic protection. Here, we use epidemiological models to show that emerging vaccine technologies, aimed at broad-spectrum protection, could qualitatively alter this picture. We demonstrate that sustained immunization with such vaccines could--through potentially lowering transmission rates and improving herd immunity--significantly moderate both influenza pandemic and seasonal epidemics. More subtly, phylodynamic models indicate that widespread cross-protective immunization could slow the antigenic evolution of seasonal influenza; these effects have profound implications for a transition to mass vaccination strategies against human influenza, and for the management of antigenically variable viruses in general.

Increasing intra- and inter-subtype HIV diversity despite declining HIV incidence in Uganda.

HIV incidence has been declining in Africa with scale-up of HIV interventions. However, there is limited data on HIV evolutionary trends in African populations with waning epidemics. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a twenty-five-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from persons living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994 to 2018) and four hyperendemic Lake Victoria fishing communities (2011 to 2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was estimated using the Shannon diversity index and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Evolutionary dynamics were assessed among demographic and behavioral sub-groups, including by migration status. 9,931 HIV sequences were available from 4,999 PLHIV, including 3,060 and 1,939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (p<0.001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (p<0.001). The proportion of viruses classified as recombinants significantly increased by more than four-fold. Inter-subtype HIV diversity has generally increased. While p24 intra-subtype genetic diversity and divergence leveled off after 2014, diversity and divergence of gp41 increased through 2017. Inter- and intra-subtype viral diversity increased across all population sub-groups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics in declining African HIV epidemics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

Age and gender profiles of HIV infection burden and viraemia: novel metrics for HIV epidemic control in African populations with high antiretroviral therapy coverage.

Introduction: To prioritize and tailor interventions for ending AIDS by 2030 in Africa, it is important to characterize the population groups in which HIV viraemia is concentrating. Methods: We analysed HIV testing and viral load data collected between 2013-2019 from the open, population-based Rakai Community Cohort Study (RCCS) in Uganda, to estimate HIV seroprevalence and population viral suppression over time by gender, one-year age bands and residence in inland and fishing communities. All estimates were standardized to the underlying source population using census data. We then assessed 95-95-95 targets in their ability to identify the populations in which viraemia concentrates. Results: Following the implementation of Universal Test and Treat, the proportion of individuals with viraemia decreased from 4.9% (4.6%-5.3%) in 2013 to 1.9% (1.7%-2.2%) in 2019 in inland communities and from 19.1% (18.0%-20.4%) in 2013 to 4.7% (4.0%-5.5%) in 2019 in fishing communities. Viraemia did not concentrate in the age and gender groups furthest from achieving 95-95-95 targets. Instead, in both inland and fishing communities, women aged 25-29 and men aged 30-34 were the 5-year age groups that contributed most to population-level viraemia in 2019, despite these groups being close to or had already achieved 95-95-95 targets. Conclusions: The 95-95-95 targets provide a useful benchmark for monitoring progress towards HIV epidemic control, but do not contextualize underlying population structures and so may direct interventions towards groups that represent a marginal fraction of the population with viraemia.

Demographics of sources of HIV-1 transmission in Zambia: a molecular epidemiology analysis in the HPTN 071 PopART study.

BACKGROUND: In the last decade, universally available antiretroviral therapy (ART) has led to greatly improved health and survival of people living with HIV in sub-Saharan Africa, but new infections continue to appear. The design of effective prevention strategies requires the demographic characterisation of individuals acting as sources of infection, which is the aim of this study. METHODS: Between 2014 and 2018, the HPTN 071 PopART study was conducted to quantify the public health benefits of ART. Viral samples from 7124 study participants in Zambia were deep-sequenced as part of HPTN 071-02 PopART Phylogenetics, an ancillary study. We used these sequences to identify likely transmission pairs. After demographic weighting of the recipients in these pairs to match the overall HIV-positive population, we analysed the demographic characteristics of the sources to better understand transmission in the general population. FINDINGS: We identified a total of 300 likely transmission pairs. 178 (59·4%) were male to female, with 130 (95% CI 110-150; 43·3%) from males aged 25-40 years. Overall, men transmitted 2·09-fold (2·06-2·29) more infections per capita than women, a ratio peaking at 5·87 (2·78-15·8) in the 35-39 years source age group. 40 (26-57; 13·2%) transmissions linked individuals from different communities in the trial. Of 288 sources with recorded information on drug resistance mutations, 52 (38-69; 18·1%) carried viruses resistant to first-line ART. INTERPRETATION: HIV-1 transmission in the HPTN 071 study communities comes from a wide range of age and sex groups, and there is no outsized contribution to new infections from importation or drug resistance mutations. Men aged 25-39 years, underserved by current treatment and prevention services, should be prioritised for HIV testing and ART. FUNDING: National Institute of Allergy and Infectious Diseases, US President's Emergency Plan for AIDS Relief, International Initiative for Impact Evaluation, Bill & Melinda Gates Foundation, National Institute on Drug Abuse, and National Institute of Mental Health.

Temporal dynamics and drivers of durable HIV viral load suppression and persistent high- and low-level viraemia during Universal Test and Treat scale-up in Uganda: a population-based study.

INTRODUCTION: Population-level data on durable HIV viral load suppression (VLS) following the implementation of Universal Test and Treat (UTT) in Africa are limited. We assessed trends in durable VLS and viraemia among persons living with HIV in 40 Ugandan communities during the UTT scale-up. METHODS: In 2015-2020, we measured VLS (<200 RNA copies/ml) among participants in the Rakai Community Cohort Study, a longitudinal population-based HIV surveillance cohort in southern Uganda. Persons with unsuppressed viral loads were characterized as having low-level (200-999 copies/ml) or high-level (≥1000 copies/ml) viraemia. Individual virologic outcomes were assessed over two consecutive RCCS survey visits (i.e. visit-pairs; ∼18-month visit intervals) and classified as durable VLS (<200 copies/ml at both visits), new/renewed VLS (<200 copies/ml at follow-up only), viral rebound (<200 copies/ml at initial visit only) or persistent viraemia (≥200 copies/ml at both visits). Population prevalence of each outcome was assessed over calendar time. Community-level prevalence and individual-level predictors of persistent high-level viraemia were also assessed using multivariable Poisson regression with generalized estimating equations. RESULTS: Overall, 3080 participants contributed 4604 visit-pairs over three survey rounds. Most visit-pairs (72.4%) exhibited durable VLS, with few (2.5%) experiencing viral rebound. Among those with any viraemia at the initial visit (23.5%, n = 1083), 46.9% remained viraemic through follow-up, 91.3% of which was high-level viraemia. One-fifth (20.8%) of visit-pairs exhibiting persistent high-level viraemia self-reported antiretroviral therapy (ART) use for ≥12 months. Prevalence of persistent high-level viraemia varied substantially across communities and was significantly elevated among young persons aged 15-29 years (vs. 40- to 49-year-olds; adjusted risk ratio [adjRR] = 2.96; 95% confidence interval [95% CI]: 2.21-3.96), males (vs. females; adjRR = 2.40, 95% CI: 1.87-3.07), persons reporting inconsistent condom use with non-marital/casual partners (vs. persons with marital/permanent partners only; adjRR = 1.38, 95% CI: 1.10-1.74) and persons reporting hazardous alcohol use (adjRR = 1.09, 95% CI: 1.03-1.16). The prevalence of persistent high-level viraemia was highest among males <30 years (32.0%). CONCLUSIONS: Following universal ART provision, most persons living with HIV in south-central Uganda are durably suppressed. Among persons exhibiting any viraemia, nearly half exhibited high-level viraemia for ≥12 months and reported higher-risk behaviours associated with onward HIV transmission. Intensified efforts linking individuals to HIV treatment services could accelerate momentum towards HIV epidemic control.

Longitudinal population-level HIV epidemiologic and genomic surveillance highlights growing gender disparity of HIV transmission in Uganda.

HIV incidence in eastern and southern Africa has historically been concentrated among girls and women aged 15-24 years. As new cases decline with HIV interventions, population-level infection dynamics may shift by age and gender. Here, we integrated population-based surveillance of 38,749 participants in the Rakai Community Cohort Study and longitudinal deep-sequence viral phylogenetics to assess how HIV incidence and population groups driving transmission have changed from 2003 to 2018 in Uganda. We observed 1,117 individuals in the incidence cohort and 1,978 individuals in the transmission cohort. HIV viral suppression increased more rapidly in women than men, however incidence declined more slowly in women than men. We found that age-specific transmission flows shifted: whereas HIV transmission to girls and women (aged 15-24 years) from older men declined by about one-third, transmission to women (aged 25-34 years) from men that were 0-6 years older increased by half in 2003 to 2018. Based on changes in transmission flows, we estimated that closing the gender gap in viral suppression could have reduced HIV incidence in women by half in 2018. This study suggests that HIV programmes to increase HIV suppression in men are critical to reduce incidence in women, close gender gaps in infection burden and improve men's health in Africa.

Phylodynamic inference and model assessment with approximate bayesian computation: influenza as a case study.

A key priority in infectious disease research is to understand the ecological and evolutionary drivers of viral diseases from data on disease incidence as well as viral genetic and antigenic variation. We propose using a simulation-based, Bayesian method known as Approximate Bayesian Computation (ABC) to fit and assess phylodynamic models that simulate pathogen evolution and ecology against summaries of these data. We illustrate the versatility of the method by analyzing two spatial models describing the phylodynamics of interpandemic human influenza virus subtype A(H3N2). The first model captures antigenic drift phenomenologically with continuously waning immunity, and the second epochal evolution model describes the replacement of major, relatively long-lived antigenic clusters. Combining features of long-term surveillance data from The Netherlands with features of influenza A (H3N2) hemagglutinin gene sequences sampled in northern Europe, key phylodynamic parameters can be estimated with ABC. Goodness-of-fit analyses reveal that the irregularity in interannual incidence and H3N2's ladder-like hemagglutinin phylogeny are quantitatively only reproduced under the epochal evolution model within a spatial context. However, the concomitant incidence dynamics result in a very large reproductive number and are not consistent with empirical estimates of H3N2's population level attack rate. These results demonstrate that the interactions between the evolutionary and ecological processes impose multiple quantitative constraints on the phylodynamic trajectories of influenza A(H3N2), so that sequence and surveillance data can be used synergistically. ABC, one of several data synthesis approaches, can easily interface a broad class of phylodynamic models with various types of data but requires careful calibration of the summaries and tolerance parameters.

Intrinsic randomness in epidemic modelling beyond statistical uncertainty.

Uncertainty can be classified as either aleatoric (intrinsic randomness) or epistemic (imperfect knowledge of parameters). The majority of frameworks assessing infectious disease risk consider only epistemic uncertainty. We only ever observe a single epidemic, and therefore cannot empirically determine aleatoric uncertainty. Here, we characterise both epistemic and aleatoric uncertainty using a time-varying general branching process. Our framework explicitly decomposes aleatoric variance into mechanistic components, quantifying the contribution to uncertainty produced by each factor in the epidemic process, and how these contributions vary over time. The aleatoric variance of an outbreak is itself a renewal equation where past variance affects future variance. We find that, superspreading is not necessary for substantial uncertainty, and profound variation in outbreak size can occur even without overdispersion in the offspring distribution (i.e. the distribution of the number of secondary infections an infected person produces). Aleatoric forecasting uncertainty grows dynamically and rapidly, and so forecasting using only epistemic uncertainty is a significant underestimate. Therefore, failure to account for aleatoric uncertainty will ensure that policymakers are misled about the substantially higher true extent of potential risk. We demonstrate our method, and the extent to which potential risk is underestimated, using two historical examples.

Assessment of COVID-19 as the Underlying Cause of Death Among Children and Young People Aged 0 to 19 Years in the US.

Importance: COVID-19 was the underlying cause of death for more than 940 000 individuals in the US, including at least 1289 children and young people (CYP) aged 0 to 19 years, with at least 821 CYP deaths occurring in the 1-year period from August 1, 2021, to July 31, 2022. Because deaths among US CYP are rare, the mortality burden of COVID-19 in CYP is best understood in the context of all other causes of CYP death. Objective: To determine whether COVID-19 is a leading (top 10) cause of death in CYP in the US. Design, Setting, and Participants: This national population-level cross-sectional epidemiological analysis for the years 2019 to 2022 used data from the US Centers for Disease Control and Prevention Wide-Ranging Online Data for Epidemiologic Research (WONDER) database on underlying cause of death in the US to identify the ranking of COVID-19 relative to other causes of death among individuals aged 0 to 19 years. COVID-19 deaths were considered in 12-month periods between April 1, 2020, and August 31, 2022, compared with deaths from leading non-COVID-19 causes in 2019, 2020, and 2021. Main Outcomes and Measures: Cause of death rankings by total number of deaths, crude rates per 100 000 population, and percentage of all causes of death, using the National Center for Health Statistics 113 Selected Causes of Death, for ages 0 to 19 and by age groupings (<1 year, 1-4 years, 5-9 years, 10-14 years, 15-19 years). Results: There were 821 COVID-19 deaths among individuals aged 0 to 19 years during the study period, resulting in a crude death rate of 1.0 per 100 000 population overall; 4.3 per 100 000 for those younger than 1 year; 0.6 per 100 000 for those aged 1 to 4 years; 0.4 per 100 000 for those aged 5 to 9 years; 0.5 per 100 000 for those aged 10 to 14 years; and 1.8 per 100 000 for those aged 15 to 19 years. COVID-19 mortality in the time period of August 1, 2021, to July 31, 2022, was among the 10 leading causes of death in CYP aged 0 to 19 years in the US, ranking eighth among all causes of deaths, fifth in disease-related causes of deaths (excluding unintentional injuries, assault, and suicide), and first in deaths caused by infectious or respiratory diseases when compared with 2019. COVID-19 deaths constituted 2% of all causes of death in this age group. Conclusions and Relevance: The findings of this study suggest that COVID-19 was a leading cause of death in CYP. It caused substantially more deaths in CYP annually than any vaccine-preventable disease historically in the recent period before vaccines became available. Various factors, including underreporting and not accounting for COVID-19's role as a contributing cause of death from other diseases, mean that these estimates may understate the true mortality burden of COVID-19. The findings of this study underscore the public health relevance of COVID-19 to CYP. In the likely future context of sustained SARS-CoV-2 circulation, appropriate pharmaceutical and nonpharmaceutical interventions (eg, vaccines, ventilation, air cleaning) will continue to play an important role in limiting transmission of the virus and mitigating severe disease in CYP.