Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Cost and cost-effectiveness of a universal HIV testing and treatment intervention in Zambia and South Africa: evidence and projections from the HPTN 071 (PopART) trial.

BACKGROUND: The HPTN 071 (PopART) trial showed that a combination HIV prevention package including universal HIV testing and treatment (UTT) reduced population-level incidence of HIV compared with standard care. However, evidence is scarce on the costs and cost-effectiveness of such an intervention. METHODS: Using an individual-based model, we simulated the PopART intervention and standard care with antiretroviral therapy (ART) provided according to national guidelines for the 21 trial communities in Zambia and South Africa (for all individuals aged >14 years), with model parameters and primary cost data collected during the PopART trial and from published sources. Two intervention scenarios were modelled: annual rounds of PopART from 2014 to 2030 (PopART 2014-30; as the UNAIDS Fast-Track target year) and three rounds of PopART throughout the trial intervention period (PopART 2014-17). For each country, we calculated incremental cost-effectiveness ratios (ICERs) as the cost per disability-adjusted life-year (DALY) and cost per HIV infection averted. Cost-effectiveness acceptability curves were used to indicate the probability of PopART being cost-effective compared with standard care at different thresholds of cost per DALY averted. We also assessed budget impact by projecting undiscounted costs of the intervention compared with standard care up to 2030. FINDINGS: During 2014-17, the mean cost per person per year of delivering home-based HIV counselling and testing, linkage to care, promotion of ART adherence, and voluntary medical male circumcision via community HIV care providers for the simulated population was US$6·53 (SD 0·29) in Zambia and US$7·93 (0·16) in South Africa. In the PopART 2014-30 scenario, median ICERs for PopART delivered annually until 2030 were $2111 (95% credible interval [CrI] 1827-2462) per HIV infection averted in Zambia and $3248 (2472-3963) per HIV infection averted in South Africa; and $593 (95% CrI 526-674) per DALY averted in Zambia and $645 (538-757) per DALY averted in South Africa. In the PopART 2014-17 scenario, PopART averted one infection at a cost of $1318 (1098-1591) in Zambia and $2236 (1601-2916) in South Africa, and averted one DALY at $258 (225-298) in Zambia and $326 (266-391) in South Africa, when outcomes were projected until 2030. The intervention had almost 100% probability of being cost-effective at thresholds greater than $700 per DALY averted in Zambia, and greater than $800 per DALY averted in South Africa, in the PopART 2014-30 scenario. Incremental programme costs for annual rounds until 2030 were $46·12 million (for a mean of 341 323 people) in Zambia and $30·24 million (for a mean of 165 852 people) in South Africa. INTERPRETATION: Combination prevention with universal home-based testing can be delivered at low annual cost per person but accumulates to a considerable amount when scaled for a growing population. Combination prevention including UTT is cost-effective at thresholds greater than $800 per DALY averted and can be an efficient strategy to reduce HIV incidence in high-prevalence settings. FUNDING: US National Institutes of Health, President's Emergency Plan for AIDS Relief, International Initiative for Impact Evaluation, Bill & Melinda Gates Foundation.

Quantifying Transmission Heterogeneity Using Both Pathogen Phylogenies and Incidence Time Series.

Heterogeneity in individual-level transmissibility can be quantified by the dispersion parameter k of the offspring distribution. Quantifying heterogeneity is important as it affects other parameter estimates, it modulates the degree of unpredictability of an epidemic, and it needs to be accounted for in models of infection control. Aggregated data such as incidence time series are often not sufficiently informative to estimate k. Incorporating phylogenetic analysis can help to estimate k concurrently with other epidemiological parameters. We have developed an inference framework that uses particle Markov Chain Monte Carlo to estimate k and other epidemiological parameters using both incidence time series and the pathogen phylogeny. Using the framework to fit a modified compartmental transmission model that includes the parameter k to simulated data, we found that more accurate and less biased estimates of the reproductive number were obtained by combining epidemiological and phylogenetic analyses. However, k was most accurately estimated using pathogen phylogeny alone. Accurately estimating k was necessary for unbiased estimates of the reproductive number, but it did not affect the accuracy of reporting probability and epidemic start date estimates. We further demonstrated that inference was possible in the presence of phylogenetic uncertainty by sampling from the posterior distribution of phylogenies. Finally, we used the inference framework to estimate transmission parameters from epidemiological and genetic data collected during a poliovirus outbreak. Despite the large degree of phylogenetic uncertainty, we demonstrated that incorporating phylogenetic data in parameter inference improved the accuracy and precision of estimates.

Genomic Infectious Disease Epidemiology in Partially Sampled and Ongoing Outbreaks.

Genomic data are increasingly being used to understand infectious disease epidemiology. Isolates from a given outbreak are sequenced, and the patterns of shared variation are used to infer which isolates within the outbreak are most closely related to each other. Unfortunately, the phylogenetic trees typically used to represent this variation are not directly informative about who infected whom-a phylogenetic tree is not a transmission tree. However, a transmission tree can be inferred from a phylogeny while accounting for within-host genetic diversity by coloring the branches of a phylogeny according to which host those branches were in. Here we extend this approach and show that it can be applied to partially sampled and ongoing outbreaks. This requires computing the correct probability of an observed transmission tree and we herein demonstrate how to do this for a large class of epidemiological models. We also demonstrate how the branch coloring approach can incorporate a variable number of unique colors to represent unsampled intermediates in transmission chains. The resulting algorithm is a reversible jump Monte-Carlo Markov Chain, which we apply to both simulated data and real data from an outbreak of tuberculosis. By accounting for unsampled cases and an outbreak which may not have reached its end, our method is uniquely suited to use in a public health environment during real-time outbreak investigations. We implemented this transmission tree inference methodology in an R package called TransPhylo, which is freely available from https://github.com/xavierdidelot/TransPhylo.

Impact and Cost-Effectiveness of Point-Of-Care CD4 Testing on the HIV Epidemic in South Africa.

Rapid diagnostic tools have been shown to improve linkage of patients to care. In the context of infectious diseases, assessing the impact and cost-effectiveness of such tools at the population level, accounting for both direct and indirect effects, is key to informing adoption of these tools. Point-of-care (POC) CD4 testing has been shown to be highly effective in increasing the proportion of HIV positive patients who initiate ART. We assess the impact and cost-effectiveness of introducing POC CD4 testing at the population level in South Africa in a range of care contexts, using a dynamic compartmental model of HIV transmission, calibrated to the South African HIV epidemic. We performed a meta-analysis to quantify the differences between POC and laboratory CD4 testing on the proportion linking to care following CD4 testing. Cumulative infections averted and incremental cost-effectiveness ratios (ICERs) were estimated over one and three years. We estimated that POC CD4 testing introduced in the current South African care context can prevent 1.7% (95% CI: 0.4% - 4.3%) of new HIV infections over 1 year. In that context, POC CD4 testing was cost-effective 99.8% of the time after 1 year with a median estimated ICER of US$4,468/DALY averted. In healthcare contexts with expanded HIV testing and improved retention in care, POC CD4 testing only became cost-effective after 3 years. The results were similar when, in addition, ART was offered irrespective of CD4 count, and CD4 testing was used for clinical assessment. Our findings suggest that even if ART is expanded to all HIV positive individuals and HIV testing efforts are increased in the near future, POC CD4 testing is a cost-effective tool, even within a short time horizon. Our study also illustrates the importance of evaluating the potential impact of such diagnostic technologies at the population level, so that indirect benefits and costs can be incorporated into estimations of cost-effectiveness.

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.

Assessment of epidemic projections using recent HIV survey data in South Africa: a validation analysis of ten mathematical models of HIV epidemiology in the antiretroviral therapy era.

BACKGROUND: Mathematical models are widely used to simulate the effects of interventions to control HIV and to project future epidemiological trends and resource needs. We aimed to validate past model projections against data from a large household survey done in South Africa in 2012. METHODS: We compared ten model projections of HIV prevalence, HIV incidence, and antiretroviral therapy (ART) coverage for South Africa with estimates from national household survey data from 2012. Model projections for 2012 were made before the publication of the 2012 household survey. We compared adult (age 15-49 years) HIV prevalence in 2012, the change in prevalence between 2008 and 2012, and prevalence, incidence, and ART coverage by sex and by age groups between model projections and the 2012 household survey. FINDINGS: All models projected lower prevalence estimates for 2012 than the survey estimate (18·8%), with eight models' central projections being below the survey 95% CI (17·5-20·3). Eight models projected that HIV prevalence would remain unchanged (n=5) or decline (n=3) between 2008 and 2012, whereas prevalence estimates from the household surveys increased from 16·9% in 2008 to 18·8% in 2012 (difference 1·9, 95% CI -0·1 to 3·9). Model projections accurately predicted the 1·6 percentage point prevalence decline (95% CI -0·3 to 3·5) in young adults aged 15-24 years, and the 2·2 percentage point (0·5 to 3·9) increase in those aged 50 years and older. Models accurately represented the number of adults on ART in 2012; six of ten models were within the survey 95% CI of 1·54-2·12 million. However, the differential ART coverage between women and men was not fully captured; all model projections of the sex ratio of women to men on ART were lower than the survey estimate of 2·22 (95% CI 1·73-2·71). INTERPRETATION: Projections for overall declines in HIV epidemics during the ART era might have been optimistic. Future treatment and HIV prevention needs might be greater than previously forecasted. Additional data about service provision for HIV care could help inform more accurate projections. FUNDING: Bill & Melinda Gates Foundation.

Health benefits, costs, and cost-effectiveness of earlier eligibility for adult antiretroviral therapy and expanded treatment coverage: a combined analysis of 12 mathematical models.

BACKGROUND: New WHO guidelines recommend initiation of antiretroviral therapy for HIV-positive adults with CD4 counts of 500 cells per µL or less, a higher threshold than was previously recommended. Country decision makers have to decide whether to further expand eligibility for antiretroviral therapy accordingly. We aimed to assess the potential health benefits, costs, and cost-effectiveness of various eligibility criteria for adult antiretroviral therapy and expanded treatment coverage. METHODS: We used several independent mathematical models in four settings-South Africa (generalised epidemic, moderate antiretroviral therapy coverage), Zambia (generalised epidemic, high antiretroviral therapy coverage), India (concentrated epidemic, moderate antiretroviral therapy coverage), and Vietnam (concentrated epidemic, low antiretroviral therapy coverage)-to assess the potential health benefits, costs, and cost-effectiveness of various eligibility criteria for adult antiretroviral therapy under scenarios of existing and expanded treatment coverage, with results projected over 20 years. Analyses assessed the extension of eligibility to include individuals with CD4 counts of 500 cells per µL or less, or all HIV-positive adults, compared with the previous (2010) recommendation of initiation with CD4 counts of 350 cells per µL or less. We assessed costs from a health-system perspective, and calculated the incremental cost (in US$) per disability-adjusted life-year (DALY) averted to compare competing strategies. Strategies were regarded very cost effective if the cost per DALY averted was less than the country's 2012 per-head gross domestic product (GDP; South Africa: $8040; Zambia: $1425; India: $1489; Vietnam: $1407) and cost effective if the cost per DALY averted was less than three times the per-head GDP. FINDINGS: In South Africa, the cost per DALY averted of extending eligibility for antiretroviral therapy to adult patients with CD4 counts of 500 cells per µL or less ranged from $237 to $1691 per DALY averted compared with 2010 guidelines. In Zambia, expansion of eligibility to adults with a CD4 count threshold of 500 cells per µL ranged from improving health outcomes while reducing costs (ie, dominating the previous guidelines) to $749 per DALY averted. In both countries results were similar for expansion of eligibility to all HIV-positive adults, and when substantially expanded treatment coverage was assumed. Expansion of treatment coverage in the general population was also cost effective. In India, the cost for extending eligibility to all HIV-positive adults ranged from $131 to $241 per DALY averted, and in Vietnam extending eligibility to patients with CD4 counts of 500 cells per µL or less cost $290 per DALY averted. In concentrated epidemics, expanded access for key populations was also cost effective. INTERPRETATION: Our estimates suggest that earlier eligibility for antiretroviral therapy is very cost effective in low-income and middle-income settings, although these estimates should be revisited when more data become available. Scaling up antiretroviral therapy through earlier eligibility and expanded coverage should be considered alongside other high-priority health interventions competing for health budgets. FUNDING: Bill & Melinda Gates Foundation, WHO.

HPTN 071 (PopART): rationale and design of a cluster-randomised trial of the population impact of an HIV combination prevention intervention including universal testing and treatment - a study protocol for a cluster randomised trial.

BACKGROUND: Effective interventions to reduce HIV incidence in sub-Saharan Africa are urgently needed. Mathematical modelling and the HIV Prevention Trials Network (HPTN) 052 trial results suggest that universal HIV testing combined with immediate antiretroviral treatment (ART) should substantially reduce incidence and may eliminate HIV as a public health problem. We describe the rationale and design of a trial to evaluate this hypothesis. METHODS/DESIGN: A rigorously-designed trial of universal testing and treatment (UTT) interventions is needed because: i) it is unknown whether these interventions can be delivered to scale with adequate uptake; ii) there are many uncertainties in the models such that the population-level impact of these interventions is unknown; and ii) there are potential adverse effects including sexual risk disinhibition, HIV-related stigma, over-burdening of health systems, poor adherence, toxicity, and drug resistance.In the HPTN 071 (PopART) trial, 21 communities in Zambia and South Africa (total population 1.2 m) will be randomly allocated to three arms. Arm A will receive the full PopART combination HIV prevention package including annual home-based HIV testing, promotion of medical male circumcision for HIV-negative men, and offer of immediate ART for those testing HIV-positive; Arm B will receive the full package except that ART initiation will follow current national guidelines; Arm C will receive standard of care. A Population Cohort of 2,500 adults will be randomly selected in each community and followed for 3 years to measure the primary outcome of HIV incidence. Based on model projections, the trial will be well-powered to detect predicted effects on HIV incidence and secondary outcomes. DISCUSSION: Trial results, combined with modelling and cost data, will provide short-term and long-term estimates of cost-effectiveness of UTT interventions. Importantly, the three-arm design will enable assessment of how much could be achieved by optimal delivery of current policies and the costs and benefits of extending this to UTT. TRIAL REGISTRATION: ClinicalTrials.gov NCT01900977.

Health benefits, costs, and cost-effectiveness of earlier eligibility for adult antiretroviral therapy and expanded treatment coverage: a combined analysis of 12 mathematical models.

BACKGROUND: New WHO guidelines recommend ART initiation for HIV-positive persons with CD4 cell counts ≤500 cells/µL, a higher threshold than was previously recommended. Country decision makers must consider whether to further expand ART eligibility accordingly. METHODS: We used multiple independent mathematical models in four settings-South Africa, Zambia, India, and Vietnam-to evaluate the potential health impact, costs, and cost-effectiveness of different adult ART eligibility criteria under scenarios of current and expanded treatment coverage, with results projected over 20 years. Analyses considered extending eligibility to include individuals with CD4 ≤500 cells/µL or all HIV-positive adults, compared to the previous recommendation of initiation with CD4 ≤350 cells/µL. We assessed costs from a health system perspective, and calculated the incremental cost per DALY averted ($/DALY) to compare competing strategies. Strategies were considered 'very cost-effective' if the $/DALY was less than the country's per capita gross domestic product (GDP; South Africa: $8040, Zambia: $1425, India: $1489, Vietnam: $1407) and 'cost-effective' if $/DALY was less than three times per capita GDP. FINDINGS: In South Africa, the cost per DALY averted of extending ART eligibility to CD4 ≤500 cells/µL ranged from $237 to $1691/DALY compared to 2010 guidelines; in Zambia, expanded eligibility ranged from improving health outcomes while reducing costs (i.e. dominating current guidelines) to $749/DALY. Results were similar in scenarios with substantially expanded treatment access and for expanding eligibility to all HIV-positive adults. Expanding treatment coverage in the general population was therefore found to be cost-effective. In India, eligibility for all HIV-positive persons ranged from $131 to $241/DALY and in Vietnam eligibility for CD4 ≤500 cells/µL cost $290/DALY. In concentrated epidemics, expanded access among key populations was also cost-effective. INTERPRETATION: Earlier ART eligibility is estimated to be very cost-effective in low- and middle-income settings, although these questions should be revisited as further information becomes available. Scaling-up ART should be considered among other high-priority health interventions competing for health budgets. FUNDING: The Bill and Melinda Gates Foundation and World Health Organization.

HIV treatment as prevention: optimising the impact of expanded HIV treatment programmes.

Until now, decisions about how to allocate ART have largely been based on maximising the therapeutic benefit of ART for patients. Since the results of the HPTN 052 study showed efficacy of antiretroviral therapy (ART) in preventing HIV transmission, there has been increased interest in the benefits of ART not only as treatment, but also in prevention. Resources for expanding ART in the short term may be limited, so the question is how to generate the most prevention benefit from realistic potential increases in the availability of ART. Although not a formal systematic review, here we review different ways in which access to ART could be expanded by prioritising access to particular groups based on clinical or behavioural factors. For each group we consider (i) the clinical and epidemiological benefits, (ii) the potential feasibility, acceptability, and equity, and (iii) the affordability and cost-effectiveness of prioritising ART access for that group. In re-evaluating the allocation of ART in light of the new data about ART preventing transmission, the goal should be to create policies that maximise epidemiological and clinical benefit while still being feasible, affordable, acceptable, and equitable.

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.

Epidemic growth rate and household reproduction number in communities of households, schools and workplaces.

In this paper we present a novel and coherent modelling framework for the characterisation of the real-time growth rate in SIR models of epidemic spread in populations with social structures of increasing complexity. Known results about homogeneous mixing and multitype models are included in the framework, which is then extended to models with households and models with households and schools/workplaces. Efficient methods for the exact computation of the real-time growth rate are presented for the standard SIR model with constant infection and recovery rates (Markovian case). Approximate methods are described for a large class of models with time-varying infection rates (non-Markovian case). The quality of the approximation is assessed via comparison with results from individual-based stochastic simulations. The methodology is then applied to the case of influenza in models with households and schools/workplaces, to provide an estimate of a household-to-household reproduction number and thus asses the effort required to prevent an outbreak by targeting control policies at the level of households. The results highlight the risk of underestimating such effort when the additional presence of schools/workplaces is neglected. Our framework increases the applicability of models of epidemic spread in socially structured population by linking earlier theoretical results, mainly focused on time-independent key epidemiological parameters (e.g. reproduction numbers, critical vaccination coverage, epidemic final size) to new results on the epidemic dynamics.

Strategies for mitigating an influenza pandemic.

Development of strategies for mitigating the severity of a new influenza pandemic is now a top global public health priority. Influenza prevention and containment strategies can be considered under the broad categories of antiviral, vaccine and non-pharmaceutical (case isolation, household quarantine, school or workplace closure, restrictions on travel) measures. Mathematical models are powerful tools for exploring this complex landscape of intervention strategies and quantifying the potential costs and benefits of different options. Here we use a large-scale epidemic simulation to examine intervention options should initial containment of a novel influenza outbreak fail, using Great Britain and the United States as examples. We find that border restrictions and/or internal travel restrictions are unlikely to delay spread by more than 2-3 weeks unless more than 99% effective. School closure during the peak of a pandemic can reduce peak attack rates by up to 40%, but has little impact on overall attack rates, whereas case isolation or household quarantine could have a significant impact, if feasible. Treatment of clinical cases can reduce transmission, but only if antivirals are given within a day of symptoms starting. Given enough drugs for 50% of the population, household-based prophylaxis coupled with reactive school closure could reduce clinical attack rates by 40-50%. More widespread prophylaxis would be even more logistically challenging but might reduce attack rates by over 75%. Vaccine stockpiled in advance of a pandemic could significantly reduce attack rates even if of low efficacy. Estimates of policy effectiveness will change if the characteristics of a future pandemic strain differ substantially from those seen in past pandemics.