Lineage frequency time series reveal elevated levels of genetic drift in SARS-CoV-2 transmission in England

Genetic drift in infectious disease transmission results from randomness of transmission and host recovery or death. The strength of genetic drift for SARS-CoV-2 transmission is expected to be high due to high levels of superspreading, and this is expected to substantially impact disease epidemiology and evolution. However, we don't yet have an understanding of how genetic drift changes over time or across locations. Furthermore, noise that results from data collection can potentially confound estimates of genetic drift. To address this challenge, we develop and validate a method to jointly infer genetic drift and measurement noise from time-series lineage frequency data. Our method is highly scalable to increasingly large genomic datasets, which overcomes a limitation in commonly used phylogenetic methods. We apply this method to over 490,000 SARS-CoV-2 genomic sequences from England collected between March 2020 and December 2021 by the COVID-19 Genomics UK (COG-UK) consortium and separately infer the strength of genetic drift for pre-B.1.177, B.1.177, Alpha, and Delta. We find that even after correcting for measurement noise, the strength of genetic drift is consistently, throughout time, higher than that expected from the observed number of COVID-19 positive individuals in England by 1 to 3 orders of magnitude, which cannot be explained by literature values of superspreading. Our estimates of genetic drift suggest low and time-varying establishment probabilities for new mutations, inform the parametrization of SARS-CoV-2 evolutionary models, and motivate future studies of the potential mechanisms for increased stochasticity in this system.

Strategy to develop broadly effective multivalent COVID-19 vaccines against emerging variants based on Ad5/35 platform

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron strain has evolved into highly divergent variants with several sub-lineages. These newly emerging variants threaten the efficacy of available COVID-19 vaccines. To mitigate the occurrence of breakthrough infections and re-infections, and more importantly, to reduce the disease burden, it is essential to develop a strategy for producing updated multivalent vaccines that can provide broad neutralization against both currently circulating and emerging variants. We developed bivalent vaccine AdCLD-CoV19-1 BA.5/BA.2.75 and trivalent vaccines AdCLD-CoV19-1 XBB/BN.1/BQ.1.1 and AdCLD-CoV19-1 XBB.1.5/BN.1/BQ.1.1 using an Ad5/35 platform-based non-replicating recombinant adenoviral vector. We compared immune responses elicited by the monovalent and multivalent vaccines in mice and macaques. We found that the BA.5/BA.2.75 bivalent and the XBB/BN.1/BQ.1.1 and XBB.1.5/BN.1/BQ.1.1 trivalent vaccines exhibited improved cross-neutralization ability compared to their respective monovalent vaccines. These data suggest that the developed multivalent vaccines enhance immunity against circulating Omicron subvariants and effectively elicit neutralizing antibodies across a broad spectrum of SARS-CoV-2 variants.

Fitness, growth and transmissibility of SARS-CoV-2 genetic variants

The massive scale of the global SARS-CoV-2 sequencing effort created new opportunities and challenges for understanding SARS-CoV-2 evolution. Rapid detection and assessment of new variants has become one of the principal objectives of genomic surveillance of SARS-CoV-2. Because of the pace and scale of sequencing, new strategies have been developed for characterizing fitness and transmissibility of emerging variants. In this Review, I discuss a wide range of approaches that have been rapidly developed in response to the public health threat posed by emerging variants, ranging from new applications of classic population genetics models to contemporary synthesis of epidemiological models and phylodynamic analysis. Many of these approaches can be adapted to other pathogens and will have increasing relevance as large-scale pathogen sequencing becomes a regular feature of many public health systems.

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.

Epidemiological drivers of transmissibility and severity of SARS-CoV-2 in England

As the SARS-CoV-2 pandemic progressed, distinct variants emerged and dominated in England. These variants, Wildtype, Alpha, Delta, and Omicron were characterized by variations in transmissibility and severity. We used a robust mathematical model and Bayesian inference framework to analyse epidemiological surveillance data from England. We quantified the impact of non-pharmaceutical interventions (NPIs), therapeutics, and vaccination on virus transmission and severity. Each successive variant had a higher intrinsic transmissibility. Omicron (BA.1) had the highest basic reproduction number at 8.3 (95% credible interval (CrI) 7.7-8.8). Varying levels of NPIs were crucial in controlling virus transmission until population immunity accumulated. Immune escape properties of Omicron decreased effective levels of immunity in the population by a third. Furthermore, in contrast to previous studies, we found Alpha had the highest basic infection fatality ratio (2.9%, 95% CrI 2.7-3.2), followed by Delta (2.2%, 95% CrI 2.0-2.4), Wildtype (1.2%, 95% CrI 1.1-1.2), and Omicron (0.7%, 95% CrI 0.6-0.8). Our findings highlight the importance of continued surveillance. Long-term strategies for monitoring and maintaining effective immunity against SARS-CoV-2 are critical to inform the role of NPIs to effectively manage future variants with potentially higher intrinsic transmissibility and severe outcomes.

The use of representative community samples to assess SARS-CoV-2 lineage competition: Alpha outcompetes Beta and wild-type in England from January to March 2021

Genomic surveillance for SARS-CoV-2 lineages informs our understanding of possible future changes in transmissibility and vaccine efficacy and will be a high priority for public health for the foreseeable future. However, small changes in the frequency of one lineage over another are often difficult to interpret because surveillance samples are obtained using a variety of methods all of which are known to contain biases. As a case study, using an approach which is largely free of biases, we here describe lineage dynamics and phylogenetic relationships of the Alpha and Beta variant in England during the first 3 months of 2021 using sequences obtained from a random community sample who provided a throat and nose swab for rt-PCR as part of the REal-time Assessment of Community Transmission-1 (REACT-1) study. Overall, diversity decreased during the first quarter of 2021, with the Alpha variant (first identified in Kent) becoming predominant, driven by a reproduction number 0.3 higher than for the prior wild-type. During January, positive samples were more likely to be Alpha in those aged 18 to 54 years old. Although individuals infected with the Alpha variant were no more likely to report one or more classic COVID-19 symptoms compared to those infected with wild-type, they were more likely to be antibody-positive 6 weeks after infection. Further, viral load was higher in those infected with the Alpha variant as measured by cycle threshold (Ct) values. The presence of infections with non-imported Beta variant (first identified in South Africa) during January, but not during February or March, suggests initial establishment in the community followed by fade-out. However, this occurred during a period of stringent social distancing. These results highlight how sequence data from representative community surveys such as REACT-1 can augment routine genomic surveillance during periods of lineage diversity.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organization as Alpha. Originating in early autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK and the imposition of new restrictions, in particular, the English national lockdown in November 2020. While these interventions succeeded in reducing the absolute number of cases, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of the SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically infected individual. We conclude that the latter provides the best explanation of the observed behaviour and dynamics of the variant, although the individual need not be immunocompromised, as persistently infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs and find that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations and a lack of the rapid evolutionary rate on its ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms), it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Evaluating whole HIV-1 genome sequence for estimation of incidence and migration in a rural South African community

Background: South Africa has the largest number of people living with HIV (PLWHIV) in the world, with HIV prevalence and transmission patterns varying greatly between provinces. Transmission between regions is still poorly understood, but phylodynamics of HIV-1 evolution can reveal how many infections are attributable to contacts outside a given community. We analysed whole genome HIV-1 genetic sequences to estimate incidence and the proportion of transmissions between communities in Hlabisa, a rural South African community. Methods: We separately analysed HIV-1 for gag, pol, and env genes sampled from 2,503 PLWHIV. We estimated time-scaled phylogenies by maximum likelihood under a molecular clock model. Phylodynamic models were fitted to time-scaled trees to estimate transmission rates, effective number of infections, incidence through time, and the proportion of infections imported to Hlabisa. We also partitioned time-scaled phylogenies with significantly different distributions of coalescent times. Results: Phylodynamic analyses showed similar trends in epidemic growth rates between 1980 and 1990. Model-based estimates of incidence and effective number of infections were consistent across genes. Parameter estimates with gag were generally smaller than those estimated with pol and env. When estimating the proportions of new infections in Hlabisa from immigration or transmission from external sources, our posterior median estimates were 85% (95% credible interval (CI) = 78%–92%) for gag, 62% (CI = 40%–78%) for pol, and 77% (CI = 58%–90%) for env in 2015. Analysis of phylogenetic partitions by gene showed that most close global reference sequences clustered within a single partition. This suggests local evolving epidemics or potential unmeasured heterogeneity in the population. Conclusions: We estimated consistent epidemic dynamic trends for gag, pol and env genes using phylodynamic models. There was a high probability that new infections were not attributable to endogenous transmission within Hlabisa, suggesting high inter-connectedness between communities in rural South Africa.

Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study

BACKGROUND

The omicron variant (B.1.1.529) of SARS-CoV-2 has demonstrated partial vaccine escape and high transmissibility, with early studies indicating lower severity of infection than that of the delta variant (B.1.617.2). We aimed to better characterise omicron severity relative to delta by assessing the relative risk of hospital attendance, hospital admission, or death in a large national cohort.

METHODS

Individual-level data on laboratory-confirmed COVID-19 cases resident in England between Nov 29, 2021, and Jan 9, 2022, were linked to routine datasets on vaccination status, hospital attendance and admission, and mortality. The relative risk of hospital attendance or admission within 14 days, or death within 28 days after confirmed infection, was estimated using proportional hazards regression. Analyses were stratified by test date, 10-year age band, ethnicity, residential region, and vaccination status, and were further adjusted for sex, index of multiple deprivation decile, evidence of a previous infection, and year of age within each age band. A secondary analysis estimated variant-specific and vaccine-specific vaccine effectiveness and the intrinsic relative severity of omicron infection compared with delta (ie, the relative risk in unvaccinated cases).

FINDINGS

The adjusted hazard ratio (HR) of hospital attendance (not necessarily resulting in admission) with omicron compared with delta was 0·56 (95% CI 0·54-0·58); for hospital admission and death, HR estimates were 0·41 (0·39-0·43) and 0·31 (0·26-0·37), respectively. Omicron versus delta HR estimates varied with age for all endpoints examined. The adjusted HR for hospital admission was 1·10 (0·85-1·42) in those younger than 10 years, decreasing to 0·25 (0·21-0·30) in 60-69-year-olds, and then increasing to 0·47 (0·40-0·56) in those aged at least 80 years. For both variants, past infection gave some protection against death both in vaccinated (HR 0·47 [0·32-0·68]) and unvaccinated (0·18 [0·06-0·57]) cases. In vaccinated cases, past infection offered no additional protection against hospital admission beyond that provided by vaccination (HR 0·96 [0·88-1·04]); however, for unvaccinated cases, past infection gave moderate protection (HR 0·55 [0·48-0·63]). Omicron versus delta HR estimates were lower for hospital admission (0·30 [0·28-0·32]) in unvaccinated cases than the corresponding HR estimated for all cases in the primary analysis. Booster vaccination with an mRNA vaccine was highly protective against hospitalisation and death in omicron cases (HR for hospital admission 8-11 weeks post-booster vs unvaccinated: 0·22 [0·20-0·24]), with the protection afforded after a booster not being affected by the vaccine used for doses 1 and 2.

INTERPRETATION

The risk of severe outcomes following SARS-CoV-2 infection is substantially lower for omicron than for delta, with higher reductions for more severe endpoints and significant variation with age. Underlying the observed risks is a larger reduction in intrinsic severity (in unvaccinated individuals) counterbalanced by a reduction in vaccine effectiveness. Documented previous SARS-CoV-2 infection offered some protection against hospitalisation and high protection against death in unvaccinated individuals, but only offered additional protection in vaccinated individuals for the death endpoint. Booster vaccination with mRNA vaccines maintains over 70% protection against hospitalisation and death in breakthrough confirmed omicron infections.

FUNDING

Medical Research Council, UK Research and Innovation, Department of Health and Social Care, National Institute for Health Research, Community Jameel, and Engineering and Physical Sciences Research Council.

Diagnosis delays in the UK according to pre or postmigration acquisition of HIV

OBJECTIVES

The aim of this study was to evaluate whether infection occurred pre or postmigration and the associated diagnosis delay in migrants diagnosed with HIV in the UK.

DESIGN

We analyzed a cohort of individuals diagnosed with HIV in the UK in 2014-2016 born in Africa or elsewhere in Europe. Inclusion criteria were arrival within 15 years before diagnosis, availability of HIV pol sequence, and viral subtype shared by at least 10 individuals.

METHODS

We examined phylogenies for evidence of infection after entry into the UK and incorporated this information into a Bayesian analysis of timing of infection using biomarkers of CD4+ cell count, avidity assays, proportion of ambiguous nucleotides in viral sequences, and last negative test dates where available.

RESULTS

One thousand, two hundred and fifty-six individuals were included. The final model indicated that HIV was acquired postmigration for most MSM born in Europe (posterior expectation 65%, 95% credibility interval 64-67%) or Africa (65%, 62-69%), whereas a minority (20-30%) of men and women with heterosexual transmission acquired HIV postmigration. Estimated diagnosis delays were lower for MSM than for those with heterosexual transmission, and were lower for those with postmigration infection across all subgroups. For MSM acquiring HIV postmigration, the estimated mean time to diagnosis was less than one year, but for those who acquired HIV premigration, the mean time from infection to diagnosis was more than five years for all subgroups.

CONCLUSION

Acquisition of HIV postmigration is common, particularly among MSM, calling for prevention efforts aimed at migrant communities. Delays in diagnosis reinforce the need for targeted testing initiatives.

Genomic assessment of quarantine measures to prevent SARS-CoV-2 importation and transmission

Mitigation of SARS-CoV-2 transmission from international travel is a priority. We evaluated the effectiveness of travellers being required to quarantine for 14-days on return to England in Summer 2020. We identified 4,207 travel-related SARS-CoV-2 cases and their contacts, and identified 827 associated SARS-CoV-2 genomes. Overall, quarantine was associated with a lower rate of contacts, and the impact of quarantine was greatest in the 16-20 age-group. 186 SARS-CoV-2 genomes were sufficiently unique to identify travel-related clusters. Fewer genomically-linked cases were observed for index cases who returned from countries with quarantine requirement compared to countries with no quarantine requirement. This difference was explained by fewer importation events per identified genome for these cases, as opposed to fewer onward contacts per case. Overall, our study demonstrates that a 14-day quarantine period reduces, but does not completely eliminate, the onward transmission of imported cases, mainly by dissuading travel to countries with a quarantine requirement.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases 1-3 . The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions 4,5 . Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta's invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Bayesian inference of clonal expansions in a dated phylogeny

Microbial population genetics models often assume that all lineages are constrained by the same population size dynamics over time. However, many neutral and selective events can invalidate this assumption and can contribute to the clonal expansion of a specific lineage relative to the rest of the population. Such differential phylodynamic properties between lineages result in asymmetries and imbalances in phylogenetic trees that are sometimes described informally but which are difficult to analyze formally. To this end, we developed a model of how clonal expansions occur and affect the branching patterns of a phylogeny. We show how the parameters of this model can be inferred from a given dated phylogeny using Bayesian statistics, which allows us to assess the probability that one or more clonal expansion events occurred. For each putative clonal expansion event, we estimate its date of emergence and subsequent phylodynamic trajectory, including its long-term evolutionary potential which is important to determine how much effort should be placed on specific control measures. We demonstrate the applicability of our methodology on simulated and real data sets. Inference under our clonal expansion model can reveal important features in the evolution and epidemiology of infectious disease pathogens. [Clonal expansion; genomic epidemiology; microbial population genomics; phylodynamics.]

Genomic reconstruction of the SARS-CoV-2 epidemic in England

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus leads to new variants that warrant timely epidemiological characterization. Here we use the dense genomic surveillance data generated by the COVID-19 Genomics UK Consortium to reconstruct the dynamics of 71 different lineages in each of 315 English local authorities between September 2020 and June 2021. This analysis reveals a series of subepidemics that peaked in early autumn 2020, followed by a jump in transmissibility of the B.1.1.7/Alpha lineage. The Alpha variant grew when other lineages declined during the second national lockdown and regionally tiered restrictions between November and December 2020. A third more stringent national lockdown suppressed the Alpha variant and eliminated nearly all other lineages in early 2021. Yet a series of variants (most of which contained the spike E484K mutation) defied these trends and persisted at moderately increasing proportions. However, by accounting for sustained introductions, we found that the transmissibility of these variants is unlikely to have exceeded the transmissibility of the Alpha variant. Finally, B.1.617.2/Delta was repeatedly introduced in England and grew rapidly in early summer 2021, constituting approximately 98% of sampled SARS-CoV-2 genomes on 26 June 2021.

Changing composition of SARS-CoV-2 lineages and rise of Delta variant in England

BACKGROUND

Since its emergence in Autumn 2020, the SARS-CoV-2 Variant of Concern (VOC) B.1.1.7 (WHO label Alpha) rapidly became the dominant lineage across much of Europe. Simultaneously, several other VOCs were identified globally. Unlike B.1.1.7, some of these VOCs possess mutations thought to confer partial immune escape. Understanding when and how these additional VOCs pose a threat in settings where B.1.1.7 is currently dominant is vital.

METHODS

We examine trends in the prevalence of non-B.1.1.7 lineages in London and other English regions using passive-case detection PCR data, cross-sectional community infection surveys, genomic surveillance, and wastewater monitoring. The study period spans from 31st January 2021 to 15th May 2021.

FINDINGS

Across data sources, the percentage of non-B.1.1.7 variants has been increasing since late March 2021. This increase was initially driven by a variety of lineages with immune escape. From mid-April, B.1.617.2 (WHO label Delta) spread rapidly, becoming the dominant variant in England by late May.

INTERPRETATION

The outcome of competition between variants depends on a wide range of factors such as intrinsic transmissibility, evasion of prior immunity, demographic specificities and interactions with non-pharmaceutical interventions. The presence and rise of non-B.1.1.7 variants in March likely was driven by importations and some community transmission. There was competition between non-B.1.17 variants which resulted in B.1.617.2 becoming dominant in April and May with considerable community transmission. Our results underscore that early detection of new variants requires a diverse array of data sources in community surveillance. Continued real-time information on the highly dynamic composition and trajectory of different SARS-CoV-2 lineages is essential to future control efforts.

FUNDING

National Institute for Health Research, Medicines and Healthcare products Regulatory Agency, DeepMind, EPSRC, EA Funds programme, Open Philanthropy, Academy of Medical Sciences Bill,Melinda Gates Foundation, Imperial College Healthcare NHS Trust, The Novo Nordisk Foundation, MRC Centre for Global Infectious Disease Analysis, Community Jameel, Cancer Research UK, Imperial College COVID-19 Research Fund, Medical Research Council, Wellcome Sanger Institute.

Changing composition of SARS-CoV-2 lineages and rise of Delta variant in England

BACKGROUND

Since its emergence in Autumn 2020, the SARS-CoV-2 Variant of Concern (VOC) B.1.1.7 (WHO label Alpha) rapidly became the dominant lineage across much of Europe. Simultaneously, several other VOCs were identified globally. Unlike B.1.1.7, some of these VOCs possess mutations thought to confer partial immune escape. Understanding when and how these additional VOCs pose a threat in settings where B.1.1.7 is currently dominant is vital.

METHODS

We examine trends in the prevalence of non-B.1.1.7 lineages in London and other English regions using passive-case detection PCR data, cross-sectional community infection surveys, genomic surveillance, and wastewater monitoring. The study period spans from 31st January 2021 to 15th May 2021.

FINDINGS

Across data sources, the percentage of non-B.1.1.7 variants has been increasing since late March 2021. This increase was initially driven by a variety of lineages with immune escape. From mid-April, B.1.617.2 (WHO label Delta) spread rapidly, becoming the dominant variant in England by late May.

INTERPRETATION

The outcome of competition between variants depends on a wide range of factors such as intrinsic transmissibility, evasion of prior immunity, demographic specificities and interactions with non-pharmaceutical interventions. The presence and rise of non-B.1.1.7 variants in March likely was driven by importations and some community transmission. There was competition between non-B.1.17 variants which resulted in B.1.617.2 becoming dominant in April and May with considerable community transmission. Our results underscore that early detection of new variants requires a diverse array of data sources in community surveillance. Continued real-time information on the highly dynamic composition and trajectory of different SARS-CoV-2 lineages is essential to future control efforts.

FUNDING

National Institute for Health Research, Medicines and Healthcare products Regulatory Agency, DeepMind, EPSRC, EA Funds programme, Open Philanthropy, Academy of Medical Sciences Bill,Melinda Gates Foundation, Imperial College Healthcare NHS Trust, The Novo Nordisk Foundation, MRC Centre for Global Infectious Disease Analysis, Community Jameel, Cancer Research UK, Imperial College COVID-19 Research Fund, Medical Research Council, Wellcome Sanger Institute.

Spatiotemporal invasion dynamics of SARS-CoV-2 lineage B.1.1.7 emergence

Understanding the causes and consequences of the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern is crucial to pandemic control yet difficult to achieve because they arise in the context of variable human behavior and immunity. We investigated the spatial invasion dynamics of lineage B.1.1.7 by jointly analyzing UK human mobility, virus genomes, and community-based polymerase chain reaction data. We identified a multistage spatial invasion process in which early B.1.1.7 growth rates were associated with mobility and asymmetric lineage export from a dominant source location, enhancing the effects of B.1.1.7's increased intrinsic transmissibility. We further explored how B.1.1.7 spread was shaped by nonpharmaceutical interventions and spatial variation in previous attack rates. Our findings show that careful accounting of the behavioral and epidemiological context within which variants of concern emerge is necessary to interpret correctly their observed relative growth rates.

HIV-1 Transmission linkages among persons with incident infection to inform public health surveillance

BACKGROUND

We evaluated features of HIV transmission networks involving persons diagnosed during incident HIV infection (IHI) to assess network-based opportunities to curtail onward transmission.

METHODS

Transmission networks were constructed using partial pol sequences reported to North Carolina surveillance among persons with recent (2014-2018) and past (<2014) HIV diagnoses. IHI were defined as documented acute infections or seroconversion. Demographic and virologic features of HIV genetic clusters (<1.5% pairwise genetic distance) involving ≥ 1 IHI were assessed. Persons with viral genetic links and who had diagnoses >90 days prior to an IHI were further characterized. We assessed named partner outcomes among IHI index persons using contact tracing data.

FINDINGS

Of 4,405 HIV diagnoses 2014-2018 with sequences, there were 323 (7%) IHI index persons; most were male (88%), Black (65%), young (68% <30 years), and reported sex with men (MSM) risk (79%). Index persons were more likely to be cluster members compared to non-index persons diagnosed during the same period (72% vs. 49%). In total, 162 clusters were identified involving 233 IHI, 577 recent diagnoses, and 163 past diagnoses. Most IHI cases (53%) had viral linkages to ≥1 previously diagnosed person without evidence of HIV viral suppression in the year prior to the diagnosis of the IHI index. In contact tracing, only 53% IHI cases named an HIV-positive contact, resulting in 0.5 previously diagnosed persons detected per IHI investigated. When combined with viral analyses, the detection rate of viremic previously diagnosed persons increased to 1.3.

INTERPRETATION

Integrating public health with molecular epidemiology, revealed that more than half of IHI have viral links to persons with previously diagnosed unsuppressed HIV infection which was largely unrecognized by traditional contact tracing. Enhanced partner services to support engagement and retention in HIV care and improved case finding supported by rapid phylogenetic analysis are tools to substantially reduce onward HIV transmission.

HIV genetic diversity informs stage of HIV-1 infection among patients receiving antiretroviral therapy in Botswana

Background

Human immunodeficiency virus (HIV)-1 genetic diversity increases during infection and can help infer the time elapsed since infection. However, the effect of antiretroviral treatment (ART) on the inference remains unknown.

Methods

Participants with estimated duration of HIV-1 infection based on repeated testing were sourced from cohorts in Botswana (n = 1944). Full-length HIV genome sequencing was performed from proviral deoxyribonucleic acid. We optimized a machine learning model to classify infections as < or >1 year based on viral genetic diversity, demographic, and clinical data.

Results

The best predictive model included variables for genetic diversity of HIV-1 gag, pol, and env, viral load, age, sex, and ART status. Most participants were on ART. Balanced accuracy was 90.6% (95% confidence interval, 86.7%–94.1%). We tested the algorithm among newly diagnosed participants with or without documented negative HIV tests. Among those without records, those who self-reported a negative HIV test within <1 year were more frequently classified as recent than those who reported a test >1 year previously. There was no difference in classification between those self-reporting a negative HIV test <1 year, whether or not they had a record.

Conclusions

These results indicate that recency of HIV-1 infection can be inferred from viral sequence diversity even among patients on suppressive ART.

Reduction in mobility and COVID-19 transmission

In response to the COVID-19 pandemic, countries have sought to control SARS-CoV-2 transmission by restricting population movement through social distancing interventions, thus reducing the number of contacts. Mobility data represent an important proxy measure of social distancing, and here, we characterise the relationship between transmission and mobility for 52 countries around the world. Transmission significantly decreased with the initial reduction in mobility in 73% of the countries analysed, but we found evidence of decoupling of transmission and mobility following the relaxation of strict control measures for 80% of countries. For the majority of countries, mobility explained a substantial proportion of the variation in transmissibility (median adjusted R-squared: 48%, interquartile range - IQR - across countries [27-77%]). Where a change in the relationship occurred, predictive ability decreased after the relaxation; from a median adjusted R-squared of 74% (IQR across countries [49-91%]) pre-relaxation, to a median adjusted R-squared of 30% (IQR across countries [12-48%]) post-relaxation. In countries with a clear relationship between mobility and transmission both before and after strict control measures were relaxed, mobility was associated with lower transmission rates after control measures were relaxed indicating that the beneficial effects of ongoing social distancing behaviours were substantial.

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.

Database of epidemic trends and control measures during the first wave of COVID-19 in mainland China

OBJECTIVES

In this data collation study, we aimed to provide a comprehensive database describing the epidemic trends and responses during the first wave of coronavirus disease 2019 (COVID-19) throughout the main provinces in China.

METHODS

From mid-January to March 2020, we extracted publicly available data regarding the spread and control of COVID-19 from 31 provincial health authorities and major media outlets in mainland China. Based on these data, we conducted descriptive analyses of the epidemic in the six most-affected provinces.

RESULTS

School closures, travel restrictions, community-level lockdown, and contact tracing were introduced concurrently around late January but subsequent epidemic trends differed among provinces. Compared with Hubei, the other five most-affected provinces reported a lower crude case fatality ratio and proportion of critical and severe hospitalised cases. From March 2020, as the local transmission of COVID-19 declined, switching the focus of measures to the testing and quarantine of inbound travellers may have helped to sustain the control of the epidemic.

CONCLUSIONS

Aggregated indicators of case notifications and severity distributions are essential for monitoring an epidemic. A publicly available database containing these indicators and information regarding control measures is a useful resource for further research and policy planning in response to the COVID-19 epidemic.

Emerging phylogenetic structure of the SARS-CoV-2 pandemic

Since spilling over into humans, SARS-CoV-2 has rapidly spread across the globe, accumulating significant genetic diversity. The structure of this genetic diversity and whether it reveals epidemiological insights are fundamental questions for understanding the evolutionary trajectory of this virus. Here, we use a recently developed phylodynamic approach to uncover phylogenetic structures underlying the SARS-CoV-2 pandemic. We find support for three SARS-CoV-2 lineages co-circulating, each with significantly different demographic dynamics concordant with known epidemiological factors. For example, Lineage C emerged in Europe with a high growth rate in late February, just prior to the exponential increase in cases in several European countries. Non-synonymous mutations that characterize Lineage C occur in functionally important gene regions responsible for viral replication and cell entry. Even though Lineages A and B had distinct demographic patterns, they were much more difficult to distinguish. Continuous application of phylogenetic approaches to track the evolutionary epidemiology of SARS-CoV-2 lineages will be increasingly important to validate the efficacy of control efforts and monitor significant evolutionary events in the future.

Estimates of the severity of coronavirus disease 2019: a model-based analysis

BACKGROUND

In the face of rapidly changing data, a range of case fatality ratio estimates for coronavirus disease 2019 (COVID-19) have been produced that differ substantially in magnitude. We aimed to provide robust estimates, accounting for censoring and ascertainment biases.

METHODS

We collected individual-case data for patients who died from COVID-19 in Hubei, mainland China (reported by national and provincial health commissions to Feb 8, 2020), and for cases outside of mainland China (from government or ministry of health websites and media reports for 37 countries, as well as Hong Kong and Macau, until Feb 25, 2020). These individual-case data were used to estimate the time between onset of symptoms and outcome (death or discharge from hospital). We next obtained age-stratified estimates of the case fatality ratio by relating the aggregate distribution of cases to the observed cumulative deaths in China, assuming a constant attack rate by age and adjusting for demography and age-based and location-based under-ascertainment. We also estimated the case fatality ratio from individual line-list data on 1334 cases identified outside of mainland China. Using data on the prevalence of PCR-confirmed cases in international residents repatriated from China, we obtained age-stratified estimates of the infection fatality ratio. Furthermore, data on age-stratified severity in a subset of 3665 cases from China were used to estimate the proportion of infected individuals who are likely to require hospitalisation.

FINDINGS

Using data on 24 deaths that occurred in mainland China and 165 recoveries outside of China, we estimated the mean duration from onset of symptoms to death to be 17·8 days (95% credible interval [CrI] 16·9-19·2) and to hospital discharge to be 24·7 days (22·9-28·1). In all laboratory confirmed and clinically diagnosed cases from mainland China (n=70 117), we estimated a crude case fatality ratio (adjusted for censoring) of 3·67% (95% CrI 3·56-3·80). However, after further adjusting for demography and under-ascertainment, we obtained a best estimate of the case fatality ratio in China of 1·38% (1·23-1·53), with substantially higher ratios in older age groups (0·32% [0·27-0·38] in those aged <60 years vs 6·4% [5·7-7·2] in those aged ≥60 years), up to 13·4% (11·2-15·9) in those aged 80 years or older. Estimates of case fatality ratio from international cases stratified by age were consistent with those from China (parametric estimate 1·4% [0·4-3·5] in those aged <60 years [n=360] and 4·5% [1·8-11·1] in those aged ≥60 years [n=151]). Our estimated overall infection fatality ratio for China was 0·66% (0·39-1·33), with an increasing profile with age. Similarly, estimates of the proportion of infected individuals likely to be hospitalised increased with age up to a maximum of 18·4% (11·0-37·6) in those aged 80 years or older.

INTERPRETATION

These early estimates give an indication of the fatality ratio across the spectrum of COVID-19 disease and show a strong age gradient in risk of death.

FUNDING

UK Medical Research Council.

Estimates of the severity of coronavirus disease 2019: a model-based analysis

BACKGROUND

In the face of rapidly changing data, a range of case fatality ratio estimates for coronavirus disease 2019 (COVID-19) have been produced that differ substantially in magnitude. We aimed to provide robust estimates, accounting for censoring and ascertainment biases.

METHODS

We collected individual-case data for patients who died from COVID-19 in Hubei, mainland China (reported by national and provincial health commissions to Feb 8, 2020), and for cases outside of mainland China (from government or ministry of health websites and media reports for 37 countries, as well as Hong Kong and Macau, until Feb 25, 2020). These individual-case data were used to estimate the time between onset of symptoms and outcome (death or discharge from hospital). We next obtained age-stratified estimates of the case fatality ratio by relating the aggregate distribution of cases to the observed cumulative deaths in China, assuming a constant attack rate by age and adjusting for demography and age-based and location-based under-ascertainment. We also estimated the case fatality ratio from individual line-list data on 1334 cases identified outside of mainland China. Using data on the prevalence of PCR-confirmed cases in international residents repatriated from China, we obtained age-stratified estimates of the infection fatality ratio. Furthermore, data on age-stratified severity in a subset of 3665 cases from China were used to estimate the proportion of infected individuals who are likely to require hospitalisation.

FINDINGS

Using data on 24 deaths that occurred in mainland China and 165 recoveries outside of China, we estimated the mean duration from onset of symptoms to death to be 17·8 days (95% credible interval [CrI] 16·9-19·2) and to hospital discharge to be 24·7 days (22·9-28·1). In all laboratory confirmed and clinically diagnosed cases from mainland China (n=70 117), we estimated a crude case fatality ratio (adjusted for censoring) of 3·67% (95% CrI 3·56-3·80). However, after further adjusting for demography and under-ascertainment, we obtained a best estimate of the case fatality ratio in China of 1·38% (1·23-1·53), with substantially higher ratios in older age groups (0·32% [0·27-0·38] in those aged <60 years vs 6·4% [5·7-7·2] in those aged ≥60 years), up to 13·4% (11·2-15·9) in those aged 80 years or older. Estimates of case fatality ratio from international cases stratified by age were consistent with those from China (parametric estimate 1·4% [0·4-3·5] in those aged <60 years [n=360] and 4·5% [1·8-11·1] in those aged ≥60 years [n=151]). Our estimated overall infection fatality ratio for China was 0·66% (0·39-1·33), with an increasing profile with age. Similarly, estimates of the proportion of infected individuals likely to be hospitalised increased with age up to a maximum of 18·4% (11·0-37·6) in those aged 80 years or older.

INTERPRETATION

These early estimates give an indication of the fatality ratio across the spectrum of COVID-19 disease and show a strong age gradient in risk of death.

FUNDING

UK Medical Research Council.

Estimates of the severity of coronavirus disease 2019: a model-based analysis

BACKGROUND

In the face of rapidly changing data, a range of case fatality ratio estimates for coronavirus disease 2019 (COVID-19) have been produced that differ substantially in magnitude. We aimed to provide robust estimates, accounting for censoring and ascertainment biases.

METHODS

We collected individual-case data for patients who died from COVID-19 in Hubei, mainland China (reported by national and provincial health commissions to Feb 8, 2020), and for cases outside of mainland China (from government or ministry of health websites and media reports for 37 countries, as well as Hong Kong and Macau, until Feb 25, 2020). These individual-case data were used to estimate the time between onset of symptoms and outcome (death or discharge from hospital). We next obtained age-stratified estimates of the case fatality ratio by relating the aggregate distribution of cases to the observed cumulative deaths in China, assuming a constant attack rate by age and adjusting for demography and age-based and location-based under-ascertainment. We also estimated the case fatality ratio from individual line-list data on 1334 cases identified outside of mainland China. Using data on the prevalence of PCR-confirmed cases in international residents repatriated from China, we obtained age-stratified estimates of the infection fatality ratio. Furthermore, data on age-stratified severity in a subset of 3665 cases from China were used to estimate the proportion of infected individuals who are likely to require hospitalisation.

FINDINGS

Using data on 24 deaths that occurred in mainland China and 165 recoveries outside of China, we estimated the mean duration from onset of symptoms to death to be 17·8 days (95% credible interval [CrI] 16·9-19·2) and to hospital discharge to be 24·7 days (22·9-28·1). In all laboratory confirmed and clinically diagnosed cases from mainland China (n=70 117), we estimated a crude case fatality ratio (adjusted for censoring) of 3·67% (95% CrI 3·56-3·80). However, after further adjusting for demography and under-ascertainment, we obtained a best estimate of the case fatality ratio in China of 1·38% (1·23-1·53), with substantially higher ratios in older age groups (0·32% [0·27-0·38] in those aged <60 years vs 6·4% [5·7-7·2] in those aged ≥60 years), up to 13·4% (11·2-15·9) in those aged 80 years or older. Estimates of case fatality ratio from international cases stratified by age were consistent with those from China (parametric estimate 1·4% [0·4-3·5] in those aged <60 years [n=360] and 4·5% [1·8-11·1] in those aged ≥60 years [n=151]). Our estimated overall infection fatality ratio for China was 0·66% (0·39-1·33), with an increasing profile with age. Similarly, estimates of the proportion of infected individuals likely to be hospitalised increased with age up to a maximum of 18·4% (11·0-37·6) in those aged 80 years or older.

INTERPRETATION

These early estimates give an indication of the fatality ratio across the spectrum of COVID-19 disease and show a strong age gradient in risk of death.

FUNDING

UK Medical Research Council.

Estimating the number of undetected COVID-19 cases among travellers from mainland China

Background: As of August 2021, every region of the world has been affected by the COVID-19 pandemic, with more than 196,000,000 cases worldwide. Methods: We analysed COVID-19 cases among travellers from mainland China to different regions and countries, comparing the region- and country-specific rates of detected and confirmed cases per flight volume to estimate the relative sensitivity of surveillance in different regions and countries. Results: Although travel restrictions from Wuhan City and other cities across China may have reduced the absolute number of travellers to and from China, we estimated that up to 70% (95% CI: 54% - 80%) of imported cases could remain undetected relative to the sensitivity of surveillance in Singapore. The percentage of undetected imported cases rises to 75% (95% CI 66% - 82%) when comparing to the surveillance sensitivity in multiple countries. Conclusions: Our analysis shows that a large number of COVID-19 cases remain undetected across the world.  These undetected cases potentially resulted in multiple chains of human-to-human transmission outside mainland China.

Estimating the number of undetected COVID-19 cases among travellers from mainland China

Background: As of August 2021, every region of the world has been affected by the COVID-19 pandemic, with more than 196,000,000 cases worldwide. Methods: We analysed COVID-19 cases among travellers from mainland China to different regions and countries, comparing the region- and country-specific rates of detected and confirmed cases per flight volume to estimate the relative sensitivity of surveillance in different regions and countries. Results: Although travel restrictions from Wuhan City and other cities across China may have reduced the absolute number of travellers to and from China, we estimated that up to 70% (95% CI: 54% - 80%) of imported cases could remain undetected relative to the sensitivity of surveillance in Singapore. The percentage of undetected imported cases rises to 75% (95% CI 66% - 82%) when comparing to the surveillance sensitivity in multiple countries. Conclusions: Our analysis shows that a large number of COVID-19 cases remain undetected across the world.  These undetected cases potentially resulted in multiple chains of human-to-human transmission outside mainland China.

Estimating the number of undetected COVID-19 cases among travellers from mainland China

Background: Since the start of the COVID-19 epidemic in late 2019, there have been more than 152 affected regions and countries with over 110,000 confirmed cases outside mainland China. Methods: We analysed COVID-19 cases among travellers from mainland China to different regions and countries, comparing the region- and country-specific rates of detected and confirmed cases per flight volume to estimate the relative sensitivity of surveillance in different regions and countries. Results: Although travel restrictions from Wuhan City and other cities across China may have reduced the absolute number of travellers to and from China, we estimated that more than two thirds (70%, 95% CI: 54% - 80%, compared to Singapore; 75%, 95% CI: 66% - 82%, compared to multiple countries) of cases exported from mainland China have remained undetected. Conclusions: These undetected cases potentially resulted in multiple chains of human-to-human transmission outside mainland China.

HIV-1 Transmission Clustering and Phylodynamics Highlight the Important Role of Young Men Who Have Sex with Men

More persons living with HIV reside in the Southern United States than in any other region, yet little is known about HIV molecular epidemiology in the South. We used cluster and phylodynamic analyses to evaluate HIV transmission patterns in middle Tennessee. We performed cross-sectional analyses of HIV-1 pol sequences and clinical data collected from 2001 to 2015 among persons attending the Vanderbilt Comprehensive Care Clinic. Transmission clusters were identified using maximum likelihood phylogenetics and patristic distance differences. Demographic, risk behavior, and clinical factors were assessed evaluating “active” clusters (clusters including sequences sampled 2011–2015) and associations estimated with logistic regression. Transmission risk ratios for men who have sex with men (MSM) were estimated with phylodynamic models. Among 2915 persons (96% subtype-B sequences), 963 (33%) were members of 292 clusters (distance ≤1.5%, size range 2–39). Most clusters (62%, n = 690 persons) were active, either being newly identified (n = 80) or showing expansion on existing clusters (n = 101). Correlates of active clustering among persons with sequences collected during 2011–2015 included MSM risk and ≤30 years of age. Active clusters were significantly more concentrated in MSM and younger persons than historical clusters. Young MSM (YMSM) (≤26.4 years) had high estimated transmission risk [risk ratio = 4.04 (2.85–5.65) relative to older MSM] and were much more likely to transmit to YMSM. In this Tennessee cohort, transmission clusters over time were more concentrated by MSM and younger age, with high transmission risk among and between YMSM, highlighting the importance of interventions among this group. Detecting active clusters could help direct interventions to disrupt ongoing transmission chains.

Estimating the contribution of key populations towards the spread of HIV in Dakar, Senegal

INTRODUCTION

Key populations including female sex workers (FSW) and men who have sex with men (MSM) bear a disproportionate burden of HIV. However, the role of focusing prevention efforts on these groups for reducing a country's HIV epidemic is debated. We estimate the extent to which HIV transmission among FSW and MSM contributes to overall HIV transmission in Dakar, Senegal, using a dynamic assessment of the population attributable fraction (PAF).

METHODS

A dynamic transmission model of HIV among FSW, their clients, MSM and the lower-risk adult population was parameterized and calibrated within a Bayesian framework using setting-specific demographic, behavioural, HIV epidemiological and antiretroviral treatment (ART) coverage data for 1985 to 2015. We used the model to estimate the 10-year PAF of commercial sex between FSW and their clients, and sex between men, to overall HIV transmission (defined as the percentage of new infections prevented when these modes of transmission are removed). In addition, we estimated the prevention benefits associated with historical increases in condom use and ART uptake, and impact of further increases in prevention and treatment.

RESULTS

The model projections suggest that unprotected sex between men contributed to 42% (2.5 to 97.5th percentile range 24 to 59%) of transmissions between 1995 and 2005, increasing to 64% (37 to 79%) from 2015 to 2025. The 10-year PAF of commercial sex is smaller, diminishing from 21% (7 to 39%) in 1995 to 14% (5 to 35%) in 2015. Without ART, 49% (32 to 71%) more HIV infections would have occurred since 2000, when ART was initiated, whereas without condom use since 1985, 67% (27 to 179%) more HIV infections would have occurred, and the overall HIV prevalence would have been 60% (29 to 211%) greater than what it is now. Further large decreases in HIV incidence (68%) can be achieved by scaling up ART in MSM to 74% coverage and reducing their susceptibility to HIV by two-thirds through any prevention modality.

CONCLUSIONS

Unprotected sex between men may be an important contributor to HIV transmission in Dakar, due to suboptimal coverage of evidence-informed interventions. Although existing interventions have effectively reduced HIV transmission among adults, it is crucial that further strategies address the unmet need among MSM.

Share2Quit: Online Social Network Peer Marketing of Tobacco Cessation Systems

Introduction

Although technology-assisted tobacco interventions (TATIs) are effective, they are underused due to recruitment challenges. We tested whether we could successfully recruit smokers to a TATI using peer marketing through a social network (Facebook).

Methods

We recruited smokers on Facebook using online advertisements. These recruited smokers (seeds) and subsequent waves of smokers (peer recruits) were provided the Share2Quit peer recruitment Facebook app and other tools. Smokers were incentivized for up to seven successful peer recruitments and had 30 days to recruit from date of registration. Successful peer recruitment was defined as a peer recruited smoker completing the registration on the TATI following a referral. Our primary questions were (1) whether smokers would recruit other smokers and (2) whether peer recruitment would extend the reach of the intervention to harder-to-reach groups, including those not ready to quit and minority smokers.

Results

Overall, 759 smokers were recruited (seeds: 190; peer recruits: 569). Fifteen percent (n = 117) of smokers successfully recruited their peers (seeds: 24.7%; peer recruits: 7.7%) leading to four recruitment waves. Compared to seeds, peer recruits were less likely to be ready to quit (peer recruits 74.2% vs. seeds 95.1%), more likely to be male (67.1% vs. 32.9%), and more likely to be African American (23.8% vs. 10.8%) (p < .01 for all comparisons).

Conclusions

Peer marketing quadrupled our engaged smokers and enriched the sample with not-ready-to-quit and African American smokers. Peer recruitment is promising, and our study uncovered several important challenges for future research.

Implications

This study demonstrates the successful recruitment of smokers to a TATI using a Facebook-based peer marketing strategy. Smokers on Facebook were willing and able to recruit other smokers to a TATI, yielding a large and diverse population of smokers.

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.

Design and methods of a social network isolation study for reducing respiratory infection transmission: The eX-FLU cluster randomized trial

BACKGROUND

Social networks are increasingly recognized as important points of intervention, yet relatively few intervention studies of respiratory infection transmission have utilized a network design. Here we describe the design, methods, and social network structure of a randomized intervention for isolating respiratory infection cases in a university setting over a 10-week period.

METHODOLOGY/PRINCIPAL FINDINGS

590 students in six residence halls enrolled in the eX-FLU study during a chain-referral recruitment process from September 2012-January 2013. Of these, 262 joined as "seed" participants, who nominated their social contacts to join the study, of which 328 "nominees" enrolled. Participants were cluster-randomized by 117 residence halls. Participants were asked to respond to weekly surveys on health behaviors, social interactions, and influenza-like illness (ILI) symptoms. Participants were randomized to either a 3-Day dorm room isolation intervention or a control group (no isolation) upon illness onset. ILI cases reported on their isolation behavior during illness and provided throat and nasal swab specimens at onset, day-three, and day-six of illness. A subsample of individuals (N=103) participated in a sub-study using a novel smartphone application, iEpi, which collected sensor and contextually-dependent survey data on social interactions. Within the social network, participants were significantly positively assortative by intervention group, enrollment type, residence hall, iEpi participation, age, gender, race, and alcohol use (all P<0.002).

CONCLUSIONS/SIGNIFICANCE

We identified a feasible study design for testing the impact of isolation from social networks in a university setting. These data provide an unparalleled opportunity to address questions about isolation and infection transmission, as well as insights into social networks and behaviors among college-aged students. Several important lessons were learned over the course of this project, including feasible isolation durations, the need for extensive organizational efforts, as well as the need for specialized programmers and server space for managing survey and smartphone data.

Dynamic Variation in Sexual Contact Rates in a Cohort of HIV-Negative Gay Men

Human immunodeficiency virus (HIV) transmission models that include variability in sexual behavior over time have shown increased incidence, prevalence, and acute-state transmission rates for a given population risk profile. This raises the question of whether dynamic variation in individual sexual behavior is a real phenomenon that can be observed and measured. To study this dynamic variation, we developed a model incorporating heterogeneity in both between-person and within-person sexual contact patterns. Using novel methodology that we call iterated filtering for longitudinal data, we fitted this model by maximum likelihood to longitudinal survey data from the Centers for Disease Control and Prevention's Collaborative HIV Seroincidence Study (1992-1995). We found evidence for individual heterogeneity in sexual behavior over time. We simulated an epidemic process and found that inclusion of empirically measured levels of dynamic variation in individual-level sexual behavior brought the theoretical predictions of HIV incidence into closer alignment with reality given the measured per-act probabilities of transmission. The methods developed here provide a framework for quantifying variation in sexual behaviors that helps in understanding the HIV epidemic among gay men.

Phylodynamic analysis of ebola virus in the 2014 sierra leone epidemic

Background: The Ebola virus (EBOV) epidemic in Western Africa is the largest in recorded history and control efforts have so far failed to stem the rapid growth in the number of infections. Mathematical models serve a key role in estimating epidemic growth rates and the reproduction number (R0) from surveillance data and, recently, molecular sequence data. Phylodynamic analysis of existing EBOV time-stamped sequence data may provide independent estimates of the unobserved number of infections, reveal recent epidemiological history, and provide insight into selective pressures acting upon viral genes. Methods: We fit a series mathematical models of infectious disease dynamics to phylogenies estimated from 78 whole EBOV genomes collected from distinct patients in May and June of 2014 in Sierra Leone, and perform evolutionary analysis on these genomes combined with closely related EBOV genomes from previous outbreaks. Two analyses are conducted with values of the latent period that have been used in recent modelling efforts. We also examined the EBOV sequences for evidence of possible episodic adaptive molecular evolution during the 2014 outbreak. Results: We find evidence for adaptive evolution affecting L and GP protein coding regions of the EBOV genome, which is unlikely to bias molecular clock and phylodynamic analyses. We estimate R0=2.40 (95% HPD:1.54-3.87 ) if the mean latent period is 5.3 days, and R0=3.81, (95% HPD:2.47-6.3) if the mean latent period is 12.7 days. The estimated coefficient of variation (CV) of the number of transmissions per infected host is very high, and a large proportion of infections yield no transmissions. Conclusions: Estimates of R0 are sensitive to the unknown latent infectious period which can not be reliably estimated from genetic data alone. EBOV phylogenies show significant evidence for superspreading and extreme variance in the number of transmissions per infected individual during the early epidemic in Sierra Leone.

Innovative recruitment using online networks: lessons learned from an online study of alcohol and other drug use utilizing a web-based, respondent-driven sampling (webRDS) strategy

OBJECTIVE

We used a web version of Respondent-Driven Sampling (webRDS) to recruit a sample of young adults (ages 18-24) and examined whether this strategy would result in alcohol and other drug (AOD) prevalence estimates comparable to national estimates (National Survey on Drug Use and Health [NSDUH]).

METHOD

We recruited 22 initial participants (seeds) via Facebook to complete a web survey examining AOD risk correlates. Sequential, incentivized recruitment continued until our desired sample size was achieved. After correcting for webRDS clustering effects, we contrasted our AOD prevalence estimates (past 30 days) to NSDUH estimates by comparing the 95% confidence intervals of prevalence estimates.

RESULTS

We found comparable AOD prevalence estimates between our sample and NSDUH for the past 30 days for alcohol, marijuana, cocaine, Ecstasy (3,4-methylenedioxymethamphetamine, or MDMA), and hallucinogens. Cigarette use was lower than NSDUH estimates.

CONCLUSIONS

WebRDS may be a suitable strategy to recruit young adults online. We discuss the unique strengths and challenges that may be encountered by public health researchers using webRDS methods.

Web-based peer-driven chain referrals for smoking cessation

BACKGROUND

We are testing web-based respondent-driven sampling (RDS) chain referrals to recruit smokers to the Decide2Quit.org (D2Q) web-assisted tobacco intervention.

METHODS

Using an online survey of smokers, we assessed the potential of recruiting 1200 smokers in 9 months using RDS chain referrals. RDS is a complex sample design, and many factors can influence its success. We conducted simulations to determine the design of optimal RDS chains.

RESULTS

Smokers (n=48) were mostly female (72%) and between ages 30-60 (82%). Estimation of smokers in their network: 1-5 (40%), 6-10 (24%), and 10-20 (22%), with mean number of intimate family (2.2, SD=2.1) and close friend smokers (3.7, SD=3.8). Most smokers (82%) were willing to refer to D2Q and thought their friends (mean=5.0, SD=4.4, range=0-20) would be open to referral. Simulations suggested that with a quota of 3 and 10 seeds, 99.9% of the sample would be achieved in 107 days if the acceptance probability was 0.5. Acceptance probability of 25% would necessitate an increased quota.

CONCLUSIONS

Our study suggests that it is possible to recruit smokers using RDS.

Disease transmission in territorial populations: the small-world network of Serengeti lions

Territoriality in animal populations creates spatial structure that is thought to naturally buffer disease invasion. Often, however, territorial populations also include highly mobile, non-residential individuals that potentially serve as disease superspreaders. Using long-term data from the Serengeti Lion Project, we characterize the contact network structure of a territorial wildlife population and address the epidemiological impact of nomadic individuals. As expected, pride contacts are dominated by interactions with neighbouring prides and interspersed by encounters with nomads as they wander throughout the ecosystem. Yet the pride-pride network also includes occasional long-range contacts between prides, making it surprisingly small world and vulnerable to epidemics, even without nomads. While nomads increase both the local and global connectivity of the network, their epidemiological impact is marginal, particularly for diseases with short infectious periods like canine distemper virus. Thus, territoriality in Serengeti lions may be less protective and non-residents less important for disease transmission than previously considered.

Epidemic thresholds in dynamic contact networks

The reproductive ratio, R0, is a fundamental quantity in epidemiology, which determines the initial increase in an infectious disease in a susceptible host population. In most epidemic models, there is a specific value of R0, the epidemic threshold, above which epidemics are possible, but below which epidemics cannot occur. As the complexity of an epidemic model increases, so too does the difficulty of calculating epidemic thresholds. Here we derive the reproductive ratio and epidemic thresholds for susceptible-infected-recovered (SIR) epidemics in a simple class of dynamic random networks. As in most epidemiological models, R0 depends on two basic epidemic parameters, the transmission and recovery rates. We find that R0 also depends on social parameters, namely the degree distribution that describes heterogeneity in the numbers of concurrent contacts and the mixing parameter that gives the rate at which contacts are initiated and terminated. We show that social mixing fundamentally changes the epidemiological landscape and, consequently, that static network approximations of dynamic networks can be inadequate.