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

Sources of Human Immunodeficiency Virus Infections Among Men Who Have Sex With Men With a Migration Background: A Viral Phylogenetic Case Study in Amsterdam, The Netherlands

BACKGROUND

Men and women with a migration background comprise an increasing proportion of incident human immunodeficiency virus (HIV) cases across Western Europe.

METHODS

To characterize sources of transmission in local transmission chains, we used partial HIV consensus sequences with linked demographic and clinical data from the opt-out AIDS Therapy Evaluation in the Netherlands (ATHENA) cohort of people with HIV in the Netherlands and identified phylogenetically and epidemiologically possible HIV transmission pairs in Amsterdam. We interpreted these in the context of estimated infection dates, and quantified population-level sources of transmission to foreign-born and Dutch-born Amsterdam men who have sex with men (MSM) within Amsterdam transmission chains.

RESULTS

We estimate that Dutch-born MSM were the predominant sources of infections among all Amsterdam MSM who acquired their infection locally in 2010-2021, and among almost all foreign-born Amsterdam MSM subpopulations. Stratifying by 2-year intervals indicated time trends in transmission dynamics, with a majority of infections originating from foreign-born MSM since 2016, although uncertainty ranges remained wide.

CONCLUSIONS

Native-born MSM have predominantly driven HIV transmissions in Amsterdam in 2010-2021. However, in the context of rapidly declining incidence in Amsterdam, the contribution from foreign-born MSM living in Amsterdam is increasing, with some evidence that most local transmissions have been from foreign-born Amsterdam MSM since 2016.

Population sequencing for diversity and transmission analyses

Genomic diversity in a pathogen population is the foundation for evolution and adaptations in virulence, drug resistance, pathogenesis, and immune evasion. Characterizing, analyzing, and understanding population-level diversity is also essential for epidemiological and forensic tracking of sources and revealing detailed pathways of transmission and spread. For bacteria, culturing, isolating, and sequencing the large number of individual colonies required to adequately sample diversity can be prohibitively time-consuming and expensive. While sequencing directly from a mixed population will show variants among reads, they cannot be linked to reveal allele combinations associated with particular traits or phylogenetic inheritance patterns. Here, we describe the theory and method of how population sequencing directly from a mixed sample can be used in conjunction with sequencing a very small number of colonies to describe the phylogenetic diversity of a population without haplotype reconstruction. To demonstrate the utility of population sequencing in capturing phylogenetic diversity, we compared isogenic clones to population sequences of Burkholderia pseudomallei from the sputum of a single patient. We also analyzed population sequences of Staphylococcus aureus derived from different people and different body sites. Sequencing results confirm our ability to capture and characterize phylogenetic diversity in our samples. Our analyses of B. pseudomallei populations led to the surprising discovery that the pathogen population is highly structured in sputum, suggesting that for some pathogens, sputum sampling may preserve structuring in the lungs and thus present a non-invasive alternative to understanding colonization, movement, and pathogen/host interactions. Our analyses of S. aureus samples show how comparing phylogenetic diversity across populations can reveal directionality of transmission between hosts and across body sites, demonstrating the power and utility for characterizing the spread of disease and identification of reservoirs at the finest levels. We anticipate that population sequencing and analysis can be broadly applied to accelerate research in a broad range of fields reliant on a foundational understanding of population diversity.

Assessing transmission attribution risk from simulated sequencing data in HIV molecular epidemiology

BACKGROUND

HIV molecular epidemiology (ME) is the analysis of sequence data together with individual-level clinical, demographic, and behavioral data to understand HIV epidemiology. The use of ME has raised concerns regarding identification of the putative source in direct transmission events. This could result in harm ranging from stigma to criminal prosecution in some jurisdictions. Here we assessed the risks of ME using simulated HIV genetic sequencing data.

METHODS

We simulated social networks of men-who-have-sex-with-men, calibrating the simulations to data from San Diego. We used these networks to simulate consensus and next-generation sequence (NGS) data to evaluate the risks of identifying direct transmissions using different HIV sequence lengths, and population sampling depths. To identify the source of transmissions, we calculated infector probability and used phyloscanner software for the analysis of consensus and NGS data, respectively.

RESULTS

Consensus sequence analyses showed that the risk of correctly inferring the source (direct transmission) within identified transmission pairs was very small and independent of sampling depth. Alternatively, NGS analyses showed that identification of the source of a transmission was very accurate, but only for 6.5% of inferred pairs. False positive transmissions were also observed, where one or more unobserved intermediaries were present when compared to the true network.

CONCLUSION

Source attribution using consensus sequences rarely infers direct transmission pairs with high confidence but is still useful for population studies. In contrast, source attribution using NGS data was much more accurate in identifying direct transmission pairs, but for only a small percentage of transmission pairs analyzed.

HIV transmission dynamics and population-wide drug resistance in rural South Africa

Despite expanded antiretroviral therapy (ART) in South Africa, HIV-1 transmission persists. Integrase strand transfer inhibitors (INSTI) and long-acting injectables offer potential for superior viral suppression, but pre-existing drug resistance could threaten their effectiveness. In a community-based study in rural KwaZulu-Natal, prior to widespread INSTI usage, we enroled 18,025 individuals to characterise HIV-1 drug resistance and transmission networks to inform public health strategies. HIV testing and reflex viral load quantification were performed, with deep sequencing (20% variant threshold) used to detect resistance mutations. Phylogenetic and geospatial analyses characterised transmission clusters. One-third of participants were HIV-positive, with 21.7% having detectable viral loads; 62.1% of those with detectable viral loads were ART-naïve. Resistance to older reverse transcriptase (RT)-targeting drugs was found, but INSTI resistance remained low (<1%). Non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, particularly to rilpivirine (RPV) even in ART-naïve individuals, was concerning. Twenty percent of sequenced individuals belonged to transmission clusters, with geographic analysis highlighting higher clustering in peripheral and rural areas. Our findings suggest promise for INSTI-based strategies in this setting but underscore the need for RPV resistance screening before implementing long-acting cabotegravir (CAB) + RPV. The significant clustering emphasises the importance of geographically targeted interventions to effectively curb HIV-1 transmission.

Consideration of within-patient diversity highlights transmission pathways and antimicrobial resistance gene variability in vancomycin-resistant Enterococcus faecium

BACKGROUND

WGS is increasingly being applied to healthcare-associated vancomycin-resistant Enterococcus faecium (VREfm) outbreaks. Within-patient diversity could complicate transmission resolution if single colonies are sequenced from identified cases.

OBJECTIVES

Determine the impact of within-patient diversity on transmission resolution of VREfm.

MATERIALS AND METHODS

Fourteen colonies were collected from VREfm positive rectal screens, single colonies were collected from clinical samples and Illumina WGS was performed. Two isolates were selected for Oxford Nanopore sequencing and hybrid genome assembly to generate lineage-specific reference genomes. Mapping to closely related references was used to identify genetic variations and closely related genomes. A transmission network was inferred for the entire genome set using Phyloscanner.

RESULTS AND DISCUSSION

In total, 229 isolates from 11 patients were sequenced. Carriage of two or three sequence types was detected in 27% of patients. Presence of antimicrobial resistance genes and plasmids was variable within genomes from the same patient and sequence type. We identified two dominant sequence types (ST80 and ST1424), with two putative transmission clusters of two patients within ST80, and a single cluster of six patients within ST1424. We found transmission resolution was impaired using fewer than 14 colonies.

CONCLUSIONS

Patients can carry multiple sequence types of VREfm, and even within related lineages the presence of mobile genetic elements and antimicrobial resistance genes can vary. VREfm within-patient diversity could be considered in future to aid accurate resolution of transmission networks.

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

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

Inference of Infectious Disease Transmission through a Relaxed Bottleneck Using Multiple Genomes Per Host

In recent times, pathogen genome sequencing has become increasingly used to investigate infectious disease outbreaks. When genomic data is sampled densely enough amongst infected individuals, it can help resolve who infected whom. However, transmission analysis cannot rely solely on a phylogeny of the genomes but must account for the within-host evolution of the pathogen, which blurs the relationship between phylogenetic and transmission trees. When only a single genome is sampled for each host, the uncertainty about who infected whom can be quite high. Consequently, transmission analysis based on multiple genomes of the same pathogen per host has a clear potential for delivering more precise results, even though it is more laborious to achieve. Here, we present a new methodology that can use any number of genomes sampled from a set of individuals to reconstruct their transmission network. Furthermore, we remove the need for the assumption of a complete transmission bottleneck. We use simulated data to show that our method becomes more accurate as more genomes per host are provided, and that it can infer key infectious disease parameters such as the size of the transmission bottleneck, within-host growth rate, basic reproduction number, and sampling fraction. We demonstrate the usefulness of our method in applications to real datasets from an outbreak of Pseudomonas aeruginosa amongst cystic fibrosis patients and a nosocomial outbreak of Klebsiella pneumoniae.

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

BACKGROUND

In the last decade, universally available antiretroviral therapy (ART) has led to greatly improved health and survival of people living with HIV in sub-Saharan Africa, but new infections continue to appear. The design of effective prevention strategies requires the demographic characterisation of individuals acting as sources of infection, which is the aim of this study.

METHODS

Between 2014 and 2018, the HPTN 071 PopART study was conducted to quantify the public health benefits of ART. Viral samples from 7124 study participants in Zambia were deep-sequenced as part of HPTN 071-02 PopART Phylogenetics, an ancillary study. We used these sequences to identify likely transmission pairs. After demographic weighting of the recipients in these pairs to match the overall HIV-positive population, we analysed the demographic characteristics of the sources to better understand transmission in the general population.

FINDINGS

We identified a total of 300 likely transmission pairs. 178 (59·4%) were male to female, with 130 (95% CI 110-150; 43·3%) from males aged 25-40 years. Overall, men transmitted 2·09-fold (2·06-2·29) more infections per capita than women, a ratio peaking at 5·87 (2·78-15·8) in the 35-39 years source age group. 40 (26-57; 13·2%) transmissions linked individuals from different communities in the trial. Of 288 sources with recorded information on drug resistance mutations, 52 (38-69; 18·1%) carried viruses resistant to first-line ART.

INTERPRETATION

HIV-1 transmission in the HPTN 071 study communities comes from a wide range of age and sex groups, and there is no outsized contribution to new infections from importation or drug resistance mutations. Men aged 25-39 years, underserved by current treatment and prevention services, should be prioritised for HIV testing and ART.

FUNDING

National Institute of Allergy and Infectious Diseases, US President's Emergency Plan for AIDS Relief, International Initiative for Impact Evaluation, Bill & Melinda Gates Foundation, National Institute on Drug Abuse, and National Institute of Mental Health.

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

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

Signatures of transmission in within-host M. tuberculosis variation

Background

Because M. tuberculosis evolves slowly, transmission clusters often contain multiple individuals with identical consensus genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared M. tuberculosis variation could help overcome this problem. Previous studies have reported M. tuberculosis diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear.

Methods

To evaluate the transmission information present in within-host M. tuberculosis variation, we re-analyzed publicly available sequence data from three household transmission studies, using household membership as a proxy for transmission linkage between donor-recipient pairs.

Findings

We found moderate levels of minority variation present in M. tuberculosis sequence data from cultured isolates that varied significantly across studies (mean: 6, 7, and 170 minority variants above a 1% minor allele frequency threshold, outside of PE/PPE genes). Isolates from household members shared more minority variants than did isolates from unlinked individuals in the three studies (mean 98 shared minority variants vs. 10; 0.8 vs. 0.2, and 0.7 vs. 0.2, respectively). Shared within-host variation was significantly associated with household membership (OR: 1.51 [1.30,1.71], for one standard deviation increase in shared minority variants). Models that included shared within-host variation improved the accuracy of predicting household membership in all three studies as compared to models without within-host variation (AUC: 0.95 versus 0.92, 0.99 versus 0.95, and 0.93 versus 0.91).

Interpretation

Within-host M. tuberculosis variation persists through culture and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlights the need to optimize approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference.

Funding

NIAID: 5K01AI173385.

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

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

Phylogenetic identification of influenza virus candidates for seasonal vaccines

The seasonal influenza (flu) vaccine is designed to protect against those influenza viruses predicted to circulate during the upcoming flu season, but identifying which viruses are likely to circulate is challenging. We use features from phylogenetic trees reconstructed from hemagglutinin (HA) and neuraminidase (NA) sequences, together with a support vector machine, to predict future circulation. We obtain accuracies of 0.75 to 0.89 (AUC 0.83 to 0.91) over 2016-2020. We explore ways to select potential candidates for a seasonal vaccine and find that the machine learning model has a moderate ability to select strains that are close to future populations. However, consensus sequences among the most recent 3 years also do well at this task. We identify similar candidate strains to those proposed by the World Health Organization, suggesting that this approach can help inform vaccine strain selection.

Inferring bacterial transmission dynamics using deep sequencing genomic surveillance data

Identifying and interrupting transmission chains is important for controlling infectious diseases. One way to identify transmission pairs - two hosts in which infection was transmitted from one to the other - is using the variation of the pathogen within each single host (within-host variation). However, the role of such variation in transmission is understudied due to a lack of experimental and clinical datasets that capture pathogen diversity in both donor and recipient hosts. In this work, we assess the utility of deep-sequenced genomic surveillance (where genomic regions are sequenced hundreds to thousands of times) using a mouse transmission model involving controlled spread of the pathogenic bacterium Citrobacter rodentium from infected to naïve female animals. We observe that within-host single nucleotide variants (iSNVs) are maintained over multiple transmission steps and present a model for inferring the likelihood that a given pair of sequenced samples are linked by transmission. In this work we show that, beyond the presence and absence of within-host variants, differences arising in the relative abundance of iSNVs (allelic frequency) can infer transmission pairs more precisely. Our approach further highlights the critical role bottlenecks play in reserving the within-host diversity during transmission.

Phylogenetic inference of pneumococcal transmission from cross-sectional data, a pilot study

Background: Inference on pneumococcal transmission has mostly relied on longitudinal studies which are costly and resource intensive. Therefore, we conducted a pilot study to test the ability to infer who infected whom from cross-sectional pneumococcal sequences using phylogenetic inference. Methods: Five suspected transmission pairs, for which there was epidemiological evidence of who infected whom, were selected from a household study. For each pair, Streptococcus pneumoniae full genomes were sequenced from nasopharyngeal swabs collected on the same day. The within-host genetic diversity of the pneumococcal population was used to infer the transmission direction and then cross-validated with the direction suggested by the epidemiological records. Results: The pneumococcal genomes clustered into the five households from which the samples were taken. The proportion of concordantly inferred transmission direction generally increased with increasing minimum genome fragment size and single nucleotide polymorphisms. We observed a larger proportion of unique polymorphic sites in the source bacterial population compared to that of the recipient in four of the five pairs, as expected in the case of a transmission bottleneck. The only pair that did not exhibit this effect was also the pair that had consistent discordant transmission direction compared to the epidemiological records suggesting potential misdirection as a result of false-negative sampling. Conclusions: This pilot provided support for further studies to test if the direction of pneumococcal transmission can be reliably inferred from cross-sectional samples if sequenced with sufficient depth and fragment length.

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low-frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

Limited emergence of resistance to integrase strand transfer inhibitors (INSTIs) in ART-experienced participants failing dolutegravir-based antiretroviral therapy: a cross-sectional analysis of a Northeast Nigerian cohort

BACKGROUND

Due to the high prevalence of resistance to NNRTI-based ART since 2018, consolidated recommendations from the WHO have indicated dolutegravir as the preferred drug of choice for HIV treatment globally. There is a paucity of resistance outcome data from HIV-1 non-B subtypes circulating across West Africa.

AIMS

We characterized the mutational profiles of persons living with HIV from a cross-sectional cohort in North-East Nigeria failing a dolutegravir-based ART regimen.

METHODS

WGS of plasma samples collected from 61 HIV-1-infected participants following virological failure of dolutegravir-based ART were sequenced using the Illumina platform. Sequencing was successfully completed for samples from 55 participants. Following quality control, 33 full genomes were analysed from participants with a median age of 40 years and median time on ART of 9 years. HIV-1 subtyping was performed using SNAPPy.

RESULTS

Most participants had mutational profiles reflective of exposure to previous first- and second-line ART regimens comprised NRTIs and NNRTIs. More than half of participants had one or more drug resistance-associated mutations (DRMs) affecting susceptibility to NRTIs (17/33; 52%) and NNRTIs (24/33; 73%). Almost a quarter of participants (8/33; 24.4%) had one or more DRMs affecting tenofovir susceptibility. Only one participant, infected with HIV-1 subtype G, had evidence of DRMs affecting dolutegravir susceptibility-this was characterized by the T66A, G118R, E138K and R263K mutations.

CONCLUSIONS

This study found a low prevalence of resistance to dolutegravir; the data are therefore supportive of the continual rollout of dolutegravir as the primary first-line regimen for ART-naive participants and the preferred switch to second-line ART across the region. However, population-level, longer-term data collection on dolutegravir outcomes are required to further guide implementation and policy action across the region.

Validation of Dried Blood Spots for Capturing Hepatitis C Virus Diversity for Genomic Surveillance

Dried blood spots (DBS) have emerged as a promising alternative to traditional venous blood for HCV testing. However, their capacity to accurately reflect the genetic diversity of HCV remains poorly understood. We employed deep sequencing and advanced phylogenetic analyses on paired plasma and DBS samples to evaluate the suitability of DBS for genomic surveillance. Results demonstrated that DBS captured equivalent viral diversity compared to plasma with no phylogenetic discordance observed. The ability of DBS to accurately reflect the profile of viral genetic diversity suggests it may be a promising avenue for future surveillance efforts to curb HCV outbreaks.

Public health use of HIV phylogenetic data in sub-Saharan Africa: ethical issues

Phylogenetic analyses of HIV are an increasingly accurate method of clarifying population-level patterns of transmission and linking individuals or groups with transmission events. Viral genetic data may be used by public health agencies to guide policy interventions focused on clusters of transmission or segments of the population in which transmission is concentrated. Analyses of HIV phylogenetics in high-income countries have often found that clusters of transmission play a significant role in HIV epidemics. In sub-Saharan Africa, HIV phylogenetic analyses to date suggest that clusters of transmission play a relatively minor role in local epidemics. Such analyses could nevertheless be used to guide priority setting and HIV public health programme design in Africa for sub-populations in which transmission events are more concentrated. Phylogenetic analysis raises ethical issues, in part due to the range of potential benefits and potential harms (ie, risks). Potential benefits include (1) improving knowledge of transmission patterns, (2) informing the design of focused public health interventions for subpopulations in which transmission is concentrated, (3) identifying and responding to clusters of transmission, (4) reducing stigma (in some cases) and (5) informing estimates of the (cost-)effectiveness of HIV treatment programmes. Potential harms include (1) privacy infringements, (2) increasing stigma (in some cases), (3) reducing trust in public health programmes, and (4) increased prosecution of legal cases where HIV transmission, homosexuality or sex work is criminalised. This paper provides analysis of relevant issues with a focus on sub-Saharan Africa in order to inform consultations regarding ethical best practice for HIV phylogenetics.

Quantifying the value of viral genomics when inferring who infected whom in the 2014-16 Ebola virus outbreak in Guinea

Transmission trees can be established through detailed contact histories, statistical or phylogenetic inference, or a combination of methods. Each approach has its limitations, and the extent to which they succeed in revealing a 'true' transmission history remains unclear. In this study, we compared the transmission trees obtained through contact tracing investigations and various inference methods to identify the contribution and value of each approach. We studied eighty-six sequenced cases reported in Guinea between March and November 2015. Contact tracing investigations classified these cases into eight independent transmission chains. We inferred the transmission history from the genetic sequences of the cases (phylogenetic approach), their onset date (epidemiological approach), and a combination of both (combined approach). The inferred transmission trees were then compared to those from the contact tracing investigations. Inference methods using individual data sources (i.e. the phylogenetic analysis and the epidemiological approach) were insufficiently informative to accurately reconstruct the transmission trees and the direction of transmission. The combined approach was able to identify a reduced pool of infectors for each case and highlight likely connections among chains classified as independent by the contact tracing investigations. Overall, the transmissions identified by the contact tracing investigations agreed with the evolutionary history of the viral genomes, even though some cases appeared to be misclassified. Therefore, collecting genetic sequences during outbreak is key to supplement the information contained in contact tracing investigations. Although none of the methods we used could identify one unique infector per case, the combined approach highlighted the added value of mixing epidemiological and genetic information to reconstruct who infected whom.

Dissemination of carbapenemase-producing Enterobacterales in Ireland from 2012 to 2017: a retrospective genomic surveillance study

The spread of carbapenemase-producing Enterobacterales (CPE) is of major public health concern. The transmission dynamics of CPE in hospitals, particularly at the national level, are not well understood. Here, we describe a retrospective nationwide genomic surveillance study of CPE in Ireland between 2012 and 2017. We sequenced 746 national surveillance CPE samples obtained between 2012 and 2017. After clustering the sequences, we used thresholds based on pairwise SNPs, and reported within-host diversity along with epidemiological data to infer recent putative transmissions. All clusters in circulating clones, derived from high-resolution phylogenies, of a species (Klebsiella pneumoniae, Escherichia coli, Klebsiella oxytoca, Enterobacter cloacae, Enterobacter hormaechei and Citrobacter freundii) were individually examined for evidence of transmission. Antimicrobial resistance trends over time were also assessed. We identified 352 putative transmission events in six species including widespread and frequent transmissions in three species. We detected putative outbreaks in 4/6 species with three hospitals experiencing prolonged outbreaks. The bla OXA-48 gene was the main cause of carbapenem resistance in Ireland in almost all species. An expansion in the number of sequence types carrying bla OXA-48 was an additional cause of the increasing prevalence of carbapenemase-producing K. pneumoniae and E. coli.

Low Viral Diversity Limits the Effectiveness of Sequence-Based Transmission Inference for SARS-CoV-2

Tracking the spread of infection amongst individuals within and between communities has been a major challenge during viral outbreaks. With the unprecedented scale of viral sequence data collection during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, the possibility of using phylogenetics to reconstruct past transmission events has been explored as a more rigorous alternative to traditional contact tracing; however, the reliability of sequence-based inference of transmission networks has yet to be directly evaluated. E. E. Bendall, G. Paz-Bailey, G. A. Santiago, C. A. Porucznik, et al. (mSphere 7:e00400-22, 2022, https://doi.org/10.1128/mSphere.00400-22) evaluate the potential of this technique by applying best practices sequence comparison methods to three geographically distinct cohorts that include known transmission pairs and demonstrate that linked pairs are often indistinguishable from unrelated samples. This study clearly establishes how low viral diversity limits the utility of genomic methods of epidemiological inference for SARS-CoV-2.

Challenges in Harnessing Shared Within-Host Severe Acute Respiratory Syndrome Coronavirus 2 Variation for Transmission Inference

Background

The limited variation observed among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) consensus sequences makes it difficult to reconstruct transmission linkages in outbreak settings. Previous studies have recovered variation within individual SARS-CoV-2 infections but have not yet measured the informativeness of within-host variation for transmission inference.

Methods

We performed tiled amplicon sequencing on 307 SARS-CoV-2 samples, including 130 samples from 32 individuals in 14 households and 47 longitudinally sampled individuals, from 4 prospective studies with household membership data, a proxy for transmission linkage.

Results

Consensus sequences from households had limited diversity (mean pairwise distance, 3.06 single-nucleotide polymorphisms [SNPs]; range, 0-40). Most (83.1%, 255 of 307) samples harbored at least 1 intrahost single-nucleotide variant ([iSNV] median, 117; interquartile range [IQR], 17-208), above a minor allele frequency threshold of 0.2%. Pairs in the same household shared significantly more iSNVs (mean, 1.20 iSNVs; 95% confidence interval [CI], 1.02-1.39) than did pairs in different households infected with the same viral clade (mean, 0.31 iSNVs; 95% CI, .28-.34), a signal that decreases with increasingly stringent minor allele frequency thresholds. The number of shared iSNVs was significantly associated with an increased odds of household membership (adjusted odds ratio, 1.35; 95% CI, 1.23-1.49). However, the poor concordance of iSNVs detected across sequencing replicates (24.8% and 35.0% above a 0.2% and 1% threshold) confirms technical concerns that current sequencing and bioinformatic workflows do not consistently recover low-frequency within-host variants.

Conclusions

Shared within-host variation may augment the information in consensus sequences for predicting transmission linkages. Improving sensitivity and specificity of within-host variant identification will improve the informativeness of within-host variation.

Community Insights in Phylogenetic HIV Research: The CIPHR Project Protocol

Inferring HIV transmission networks from HIV sequences is gaining popularity in the field of HIV molecular epidemiology. However, HIV sequences are often analyzed at distance from those affected by HIV epidemics, namely without the involvement of communities most affected by HIV. These remote analyses often mean that knowledge is generated in absence of lived experiences and socio-economic realities that could inform the ethical application of network-derived information in 'real world' programmes. Procedures to engage communities are noticeably absent from the HIV molecular epidemiology literature. Here we present our team's protocol for engaging community activists living in Nairobi, Kenya in a knowledge exchange process - The CIPHR Project (Community Insights in Phylogenetic HIV Research). Drawing upon a community-based participatory approach, our team will (1) explore the possibilities and limitations of HIV molecular epidemiology for key population programmes, (2) pilot a community-based HIV molecular study, and (3) co-develop policy guidelines on conducting ethically safe HIV molecular epidemiology. Critical dialogue with activist communities will offer insight into the potential uses and abuses of using such information to sharpen HIV prevention programmes. The outcome of this process holds importance to the development of policy frameworks that will guide the next generation of the global response.

Evaluating the potential of whole-genome sequencing for tracing transmission routes in experimental infections and natural outbreaks of bovine respiratory syncytial virus

Bovine respiratory syncytial virus (BRSV) is a major cause of respiratory disease in cattle. Genomic sequencing can resolve phylogenetic relationships between virus populations, which can be used to infer transmission routes and potentially inform the design of biosecurity measures. Sequencing of short (<2000 nt) segments of the 15 000-nt BRSV genome has revealed geographic and temporal clustering of BRSV populations, but insufficient variation to distinguish viruses collected from herds infected close together in space and time. This study investigated the potential for whole-genome sequencing to reveal sufficient genomic variation for inferring transmission routes between herds. Next-generation sequencing (NGS) data were generated from experimental infections and from natural outbreaks in Jämtland and Uppsala counties in Sweden. Sufficient depth of coverage for analysis of consensus and sub-consensus sequence diversity was obtained from 47 to 20 samples respectively. Few (range: 0-6 polymorphisms across the six experiments) consensus-level polymorphisms were observed along experimental transmissions. A much higher level of diversity (146 polymorphic sites) was found among the consensus sequences from the outbreak samples. The majority (144/146) of polymorphisms were between rather than within counties, suggesting that consensus whole-genome sequences show insufficient spatial resolution for inferring direct transmission routes, but might allow identification of outbreak sources at the regional scale. By contrast, within-sample diversity was generally higher in the experimental than the outbreak samples. Analyses to infer known (experimental) and suspected (outbreak) transmission links from within-sample diversity data were uninformative. In conclusion, analysis of the whole-genome sequence of BRSV from experimental samples discriminated between circulating isolates from distant areas, but insufficient diversity was observed between closely related isolates to aid local transmission route inference.

Using genetic data to identify transmission risk factors: Statistical assessment and application to tuberculosis transmission

Identifying host factors that influence infectious disease transmission is an important step toward developing interventions to reduce disease incidence. Recent advances in methods for reconstructing infectious disease transmission events using pathogen genomic and epidemiological data open the door for investigation of host factors that affect onward transmission. While most transmission reconstruction methods are designed to work with densely sampled outbreaks, these methods are making their way into surveillance studies, where the fraction of sampled cases with sequenced pathogens could be relatively low. Surveillance studies that use transmission event reconstruction then use the reconstructed events as response variables (i.e., infection source status of each sampled case) and use host characteristics as predictors (e.g., presence of HIV infection) in regression models. We use simulations to study estimation of the effect of a host factor on probability of being an infection source via this multi-step inferential procedure. Using TransPhylo-a widely-used method for Bayesian estimation of infectious disease transmission events-and logistic regression, we find that low sensitivity of identifying infection sources leads to dilution of the signal, biasing logistic regression coefficients toward zero. We show that increasing the proportion of sampled cases improves sensitivity and some, but not all properties of the logistic regression inference. Application of these approaches to real world data from a population-based TB study in Botswana fails to detect an association between HIV infection and probability of being a TB infection source. We conclude that application of a pipeline, where one first uses TransPhylo and sparsely sampled surveillance data to infer transmission events and then estimates effects of host characteristics on probabilities of these events, should be accompanied by a realistic simulation study to better understand biases stemming from imprecise transmission event inference.

SOPHIE: Viral outbreak investigation and transmission history reconstruction in a joint phylogenetic and network theory framework

Genomic epidemiology is now widely used for viral outbreak investigations. Still, this methodology faces many challenges. First, few methods account for intra-host viral diversity. Second, maximum parsimony principle continues to be employed for phylogenetic inference of transmission histories, even though maximum likelihood or Bayesian models are usually more consistent. Third, many methods utilize case-specific data, such as sampling times or infection exposure intervals. This impedes study of persistent infections in vulnerable groups, where such information has a limited use. Finally, most methods implicitly assume that transmission events are independent, although common source outbreaks violate this assumption. We propose a maximum likelihood framework, SOPHIE, based on the integration of phylogenetic and random graph models. It infers transmission networks from viral phylogenies and expected properties of inter-host social networks modeled as random graphs with given expected degree distributions. SOPHIE is scalable, accounts for intra-host diversity, and accurately infers transmissions without case-specific epidemiological data.

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

Characterizing viral within-host diversity in fast and non-equilibrium demo-genetic dynamics

High-throughput sequencing has opened the route for a deep assessment of within-host genetic diversity that can be used, e.g., to characterize microbial communities and to infer transmission links in infectious disease outbreaks. The performance of such characterizations and inferences cannot be analytically assessed in general and are often grounded on computer-intensive evaluations. Then, being able to simulate within-host genetic diversity across time under various demo-genetic assumptions is paramount to assess the performance of the approaches of interest. In this context, we built an original model that can be simulated to investigate the temporal evolution of genotypes and their frequencies under various demo-genetic assumptions. The model describes the growth and the mutation of genotypes at the nucleotide resolution conditional on an overall within-host viral kinetics, and can be tuned to generate fast non-equilibrium demo-genetic dynamics. We ran simulations of this model and computed classic diversity indices to characterize the temporal variation of within-host genetic diversity (from high-throughput amplicon sequences) of virus populations under three demographic kinetic models of viral infection. Our results highlight how demographic (viral load) and genetic (mutation, selection, or drift) factors drive variations in within-host diversity during the course of an infection. In particular, we observed a non-monotonic relationship between pathogen population size and genetic diversity, and a reduction of the impact of mutation on diversity when a non-specific host immune response is activated. The large variation in the diversity patterns generated in our simulations suggests that the underlying model provides a flexible basis to produce very diverse demo-genetic scenarios and test, for instance, methods for the inference of transmission links during outbreaks.

A Systematic Molecular Epidemiology Screen Reveals Numerous Human Immunodeficiency Virus (HIV) Type 1 Superinfections in the Swiss HIV Cohort Study

BACKGROUND

Studying human immunodeficiency virus type 1 (HIV-1) superinfection is important to understand virus transmission, disease progression, and vaccine design. But detection remains challenging, with low sampling frequencies and insufficient longitudinal samples.

METHODS

Using the Swiss HIV Cohort Study (SHCS), we developed a molecular epidemiology screening for superinfections. A phylogeny built from 22 243 HIV-1 partial polymerase sequences was used to identify potential superinfections among 4575 SHCS participants with longitudinal sequences. A subset of potential superinfections was tested by near-full-length viral genome sequencing (NFVGS) of biobanked plasma samples.

RESULTS

Based on phylogenetic and distance criteria, 325 potential HIV-1 superinfections were identified and categorized by their likelihood of being detected as superinfections due to sample misidentification. NFVGS was performed for 128 potential superinfections; of these, 52 were confirmed by NFVGS, 15 were not confirmed, and for 61 sampling did not allow confirming or rejecting superinfection because the sequenced samples did not include the relevant time points causing the superinfection signal in the original screen. Thus, NFVGS could support 52 of 67 adequately sampled potential superinfections.

CONCLUSIONS

This cohort-based molecular approach identified, to our knowledge, the largest population of confirmed superinfections, showing that, while rare with a prevalence of 1%-7%, superinfections are not negligible events.

Using phylogenetics to infer HIV-1 transmission direction between known transmission pairs

Identifying the transmission direction between individuals provides unparalleled power to understand infectious disease epidemiology. With epidemiological and clinical information typically unavailable to infer transmission direction, phylogenetic analysis of pathogen sequence data offers an alternative approach. While the success of this phylogenetic analysis varies, the reasons remain unknown. We analyze sequence data from over 100 transmission pairs for which both the transmission direction of HIV is known and detailed additional information is available. We find that easily quantifiable phylogenetic and sampling characteristics discriminate whether a phylogenetically inferred transmission direction is correct. Our analysis highlights that while phylogenetic approaches to infer transmission direction are unsuitable for individual-level analysis, such as forensic investigations, confidence in source attribution can be incorporated in population-level analyses.

Inferring the transmission direction between linked individuals living with HIV provides unparalleled power to understand the epidemiology that determines transmission. Phylogenetic ancestral-state reconstruction approaches infer the transmission direction by identifying the individual in whom the most recent common ancestor of the virus populations originated. While these methods vary in accuracy, it is unclear why. To evaluate the performance of phylogenetic ancestral-state reconstruction to determine the transmission direction of HIV-1 infection, we inferred the transmission direction for 112 transmission pairs where transmission direction and detailed additional information were available. We then fit a statistical model to evaluate the extent to which epidemiological, sampling, genetic, and phylogenetic factors influenced the outcome of the inference. Finally, we repeated the analysis under real-life conditions with only routinely available data. We found that whether ancestral-state reconstruction correctly infers the transmission direction depends principally on the phylogeny's topology. For example, under real-life conditions, the probability of identifying the correct transmission direction increases from 32%—when a monophyletic–monophyletic or paraphyletic–polyphyletic tree topology is observed and when the tip closest to the root does not agree with the state at the root—to 93% when a paraphyletic–monophyletic topology is observed and when the tip closest to the root agrees with the root state. Our results suggest that documenting larger differences in relative intrahost diversity increases our confidence in the transmission direction inference of linked pairs for population-level studies of HIV. These findings provide a practical starting point to determine our confidence in transmission direction inference from ancestral-state reconstruction.

Apathogenic proxies for transmission dynamics of a fatal virus

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

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5-50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

Estimating the potential to prevent locally acquired HIV infections in a UNAIDS Fast-Track City, Amsterdam

Background:

More than 300 cities including the city of Amsterdam in the Netherlands have joined the UNAIDS Fast-Track Cities initiative, committing to accelerate their HIV response and end the AIDS epidemic in cities by 2030. To support this commitment, we aimed to estimate the number and proportion of Amsterdam HIV infections that originated within the city, from Amsterdam residents. We also aimed to estimate the proportion of recent HIV infections during the 5-year period 2014–2018 in Amsterdam that remained undiagnosed.

Methods:

We located diagnosed HIV infections in Amsterdam using postcode data (PC4) at time of registration in the ATHENA observational HIV cohort, and used HIV sequence data to reconstruct phylogeographically distinct, partially observed Amsterdam transmission chains. Individual-level infection times were estimated from biomarker data, and used to date the phylogenetically observed transmission chains as well as to estimate undiagnosed proportions among recent infections. A Bayesian Negative Binomial branching process model was used to estimate the number, size, and growth of the unobserved Amsterdam transmission chains from the partially observed phylogenetic data.

Results:

Between 1 January 2014 and 1 May 2019, there were 846 HIV diagnoses in Amsterdam residents, of whom 516 (61%) were estimated to have been infected in 2014–2018. The rate of new Amsterdam diagnoses since 2014 (104 per 100,000) remained higher than the national rates excluding Amsterdam (24 per 100,000), and in this sense Amsterdam remained a HIV hotspot in the Netherlands. An estimated 14% [12–16%] of infections in Amsterdan MSM in 2014–2018 remained undiagnosed by 1 May 2019, and 41% [35–48%] in Amsterdam heterosexuals, with variation by region of birth. An estimated 67% [60–74%] of Amsterdam MSM infections in 2014–2018 had an Amsterdam resident as source, and 56% [41–70%] in Amsterdam heterosexuals, with heterogeneity by region of birth. Of the locally acquired infections, an estimated 43% [37–49%] were in foreign-born MSM, 41% [35–47%] in Dutch-born MSM, 10% [6–18%] in foreign-born heterosexuals, and 5% [2–9%] in Dutch-born heterosexuals. We estimate the majority of Amsterdam MSM infections in 2014–2018 originated in transmission chains that pre-existed by 2014.

Conclusions:

This combined phylogenetic, epidemiologic, and modelling analysis in the UNAIDS Fast-Track City Amsterdam indicates that there remains considerable potential to prevent HIV infections among Amsterdam residents through city-level interventions. The burden of locally acquired infection remains concentrated in MSM, and both Dutch-born and foreign-born MSM would likely benefit most from intensified city-level interventions.

Funding:

This study received funding as part of the H-TEAM initiative from Aidsfonds (project number P29701). The H-TEAM initiative is being supported by Aidsfonds (grant number: 2013169, P29701, P60803), Stichting Amsterdam Dinner Foundation, Bristol-Myers Squibb International Corp. (study number: AI424-541), Gilead Sciences Europe Ltd (grant number: PA-HIV-PREP-16-0024), Gilead Sciences (protocol numbers: CO-NL-276-4222, CO-US-276-1712, CO-NL-985-6195), and M.A.C AIDS Fund.

The bounded coalescent model: Conditioning a genealogy on a minimum root date

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online.

Within-host diversity improves phylogenetic and transmission reconstruction of SARS-CoV-2 outbreaks

Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches to reconstruct outbreaks exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is stable among repeated serial samples from the same host, is transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

Longitudinal deep sequencing informs vector selection and future deployment strategies for transmissible vaccines

Vaccination is a powerful tool in combating infectious diseases of humans and companion animals. In most wildlife, including reservoirs of emerging human diseases, achieving sufficient vaccine coverage to mitigate disease burdens remains logistically unattainable. Virally vectored "transmissible" vaccines that deliberately spread among hosts are a potentially transformative, but still theoretical, solution to the challenge of immunising inaccessible wildlife. Progress towards real-world application is frustrated by the absence of frameworks to guide vector selection and vaccine deployment prior to major in vitro and in vivo investments in vaccine engineering and testing. Here, we performed deep sequencing on field-collected samples of Desmodus rotundus betaherpesvirus (DrBHV), a candidate vector for a transmissible vaccine targeting vampire bat-transmitted rabies. We discovered 11 strains of DrBHV that varied in prevalence and geographic distribution across Peru. The phylogeographic structure of DrBHV strains was predictable from both host genetics and landscape topology, informing long-term DrBHV-vectored vaccine deployment strategies and identifying geographic areas for field trials where vaccine spread would be naturally contained. Multistrain infections were observed in 79% of infected bats. Resampling of marked individuals over 4 years showed within-host persistence kinetics characteristic of latency and reactivation, properties that might boost individual immunity and lead to sporadic vaccine transmission over the lifetime of the host. Further, strain acquisitions by already infected individuals implied that preexisting immunity and strain competition are unlikely to inhibit vaccine spread. Our results support the development of a transmissible vaccine targeting a major source of human and animal rabies in Latin America and show how genomics can enlighten vector selection and deployment strategies for transmissible vaccines.

Deep-sequence phylogenetics to quantify patterns of HIV transmission in the context of a universal testing and treatment trial - BCPP/ Ya Tsie trial

Background: Mathematical models predict that community-wide access to HIV testing-and-treatment can rapidly and substantially reduce new HIV infections. Yet several large universal test-and-treat HIV prevention trials in high-prevalence epidemics demonstrated variable reduction in population-level incidence.

Methods: To elucidate patterns of HIV spread in universal test-and-treat trials we quantified the contribution of geographic-location, gender, age and randomized-HIV-intervention to HIV transmissions in the 30-community Ya Tsie trial in Botswana. We sequenced HIV viral whole genomes from 5,114 trial participants among the 30 trial communities.

Results: Deep-sequence phylogenetic analysis revealed that most inferred HIV transmissions within the trial occurred within the same or between neighboring communities, and between similarly-aged partners. Transmissions into intervention communities from control communities were more common than the reverse post-baseline (30% [12.2 - 56.7] versus 3% [0.1 - 27.3]) than at baseline (7% [1.5 - 25.3] versus 5% [0.9 - 22.9]) compatible with a benefit from treatment-as-prevention.

Conclusion: Our findings suggest that population mobility patterns are fundamental to HIV transmission dynamics and to the impact of HIV control strategies.

Funding: This study was supported by the National Institute of General Medical Sciences (U54GM088558); the Fogarty International Center (FIC) of the U.S. National Institutes of Health (D43 TW009610); and the President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention (CDC) (Cooperative agreements U01 GH000447 and U2G GH001911).

Incorporating Within-Host Diversity in Phylogenetic Analyses for Detecting Clusters of New HIV Diagnoses

Background

Phylogenetic analyses of HIV sequences are used to detect clusters and inform public health interventions. Conventional approaches summarize within-host HIV diversity with a single consensus sequence per host of the pol gene, obtained from Sanger or next-generation sequencing (NGS). There is growing recognition that this approach discards potentially important information about within-host sequence variation, which can impact phylogenetic inference. However, whether alternative summary methods that incorporate intra-host variation impact phylogenetic inference of transmission network features is unknown.

Methods

We introduce profile sampling, a method to incorporate within-host NGS sequence diversity into phylogenetic HIV cluster inference. We compare this approach to Sanger- and NGS-derived pol and near-whole-genome consensus sequences and evaluate its potential benefits in identifying molecular clusters among all newly-HIV-diagnosed individuals over six months at the largest HIV center in Rhode Island.

Results

Profile sampling cluster inference demonstrated that within-host viral diversity impacts phylogenetic inference across individuals, and that consensus sequence approaches can obscure both magnitude and effect of these impacts. Clustering differed between Sanger- and NGS-derived consensus and profile sampling sequences, and across gene regions.

Discussion

Profile sampling can incorporate within-host HIV diversity captured by NGS into phylogenetic analyses. This additional information can improve robustness of cluster detection.

A highly virulent variant of HIV-1 circulating in the Netherlands

We discovered a highly virulent variant of subtype-B HIV-1 in the Netherlands. One hundred nine individuals with this variant had a 0.54 to 0.74 log₁₀ increase (i.e., a ~3.5-fold to 5.5-fold increase) in viral load compared with, and exhibited CD4 cell decline twice as fast as, 6604 individuals with other subtype-B strains. Without treatment, advanced HIV-CD4 cell counts below 350 cells per cubic millimeter, with long-term clinical consequences-is expected to be reached, on average, 9 months after diagnosis for individuals in their thirties with this variant. Age, sex, suspected mode of transmission, and place of birth for the aforementioned 109 individuals were typical for HIV-positive people in the Netherlands, which suggests that the increased virulence is attributable to the viral strain. Genetic sequence analysis suggests that this variant arose in the 1990s from de novo mutation, not recombination, with increased transmissibility and an unfamiliar molecular mechanism of virulence.

Parasites as conservation tools

Parasite success typically depends on a close relationship with one or more hosts; therefore, attributes of parasitic infection have the potential to provide indirect details of host natural history and are biologically relevant to animal conservation. Characterization of parasite infections has been useful in delineating host populations and has served as a proxy for assessment of environmental quality. In other cases, the utility of parasites is just being explored, for example, as indicators of host connectivity. Innovative studies of parasite biology can provide information to manage major conservation threats by using parasite assemblage, prevalence, or genetic data to provide insights into the host. Overexploitation, habitat loss and fragmentation, invasive species, and climate change are major threats to animal conservation, and all of these can be informed by parasites.

Many but small HIV-1 non-B transmission chains in the Netherlands

OBJECTIVE

The aim of this study was to investigate introductions and spread of different HIV-1 subtypes in the Netherlands.

DESIGN

We identified distinct HIV-1 transmission chains in the Netherlands within the global epidemic context through viral phylogenetic analysis of partial HIV-1 polymerase sequences from individuals enrolled in the ATHENA national HIV cohort of all persons in care since 1996, and publicly available international background sequences.

METHODS

Viral lineages circulating in the Netherlands were identified through maximum parsimony phylogeographic analysis. The proportion of HIV-1 infections acquired in-country among heterosexuals and MSM was estimated from phylogenetically observed, national transmission chains using a branching process model that accounts for incomplete sampling.

RESULTS

As of 1 January 2019, 2589 (24%) of 10 971 (41%) HIV-1 sequenced individuals in ATHENA had non-B subtypes (A1, C, D, F, G) or circulating recombinant forms (CRF01AE, CRF02AG, CRF06-cpx). The 1588 heterosexuals were in 1224, and 536 MSM in 270 phylogenetically observed transmission chains. After adjustments for incomplete sampling, most heterosexual (75%) and MSM (76%) transmission chains were estimated to include only the individual introducing the virus (size = 1). Onward transmission occurred mostly in chains size 2-5 amongst heterosexuals (62%) and in chains size at least 10 amongst MSM (64%). Considering some chains originated in-country from other risk-groups, 40% (95% confidence interval: 36-44) of non-B-infected heterosexuals and 62% (95% confidence interval: 49-73) of MSM-acquired infection in-country.

CONCLUSION

Although most HIV-1 non-B introductions showed no or very little onward transmission, a considerable proportion of non-B infections amongst both heterosexuals and MSM in the Netherlands have been acquired in-country.

OUP accepted manuscript

Escherichia coli has a rich history as biology's ‘rock star’, driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance—including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

TNet: Transmission Network Inference Using Within-Host Strain Diversity and its Application to Geographical Tracking of COVID-19 Spread

The inference of disease transmission networks is an important problem in epidemiology. One popular approach for building transmission networks is to reconstruct a phylogenetic tree using sequences from disease strains sampled from infected hosts and infer transmissions based on this tree. However, most existing phylogenetic approaches for transmission network inference are highly computationally intensive and cannot take within-host strain diversity into account. Here, we introduce a new phylogenetic approach for inferring transmission networks, TNet, that addresses these limitations. TNet uses multiple strain sequences from each sampled host to infer transmissions and is simpler and more accurate than existing approaches. Furthermore, TNet is highly scalable and able to distinguish between ambiguous and unambiguous transmission inferences. We evaluated TNet on a large collection of 560 simulated transmission networks of various sizes and diverse host, sequence, and transmission characteristics, as well as on 10 real transmission datasets with known transmission histories. Our results show that TNet outperforms two other recently developed methods, phyloscanner and SharpTNI, that also consider within-host strain diversity. We also applied TNet to a large collection of SARS-CoV-2 genomes sampled from infected individuals in many countries around the world, demonstrating how our inference framework can be adapted to accurately infer geographical transmission networks. TNet is freely available from https://compbio.engr.uconn.edu/software/TNet/.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

Drug resistance mutations in HIV: new bioinformatics approaches and challenges

Drug resistance mutations appear in HIV under treatment pressure. Resistant variants can be transmitted to treatment-naive individuals, which can lead to rapid virological failure and can limit treatment options. Consequently, quantifying the prevalence, emergence and transmission of drug resistance is critical to effectively treating patients and to shape health policies. We review recent bioinformatics developments and in particular describe: (1) the machine learning approaches intended to predict and explain the level of resistance of HIV variants from their sequence data; (2) the phylogenetic methods used to survey the emergence and dynamics of resistant HIV transmission clusters; (3) the impact of deep sequencing in studying within-host and between-host genetic diversity of HIV variants, notably regarding minority resistant variants.

Quantifying transmission dynamics of acute hepatitis C virus infections in a heterogeneous population using sequence data

Opioid substitution and syringes exchange programs have drastically reduced hepatitis C virus (HCV) spread in France but HCV sexual transmission in men having sex with men (MSM) has recently arisen as a significant public health concern. The fact that the virus is transmitting in a heterogeneous population, with different transmission routes, makes prevalence and incidence rates poorly informative. However, additional insights can be gained by analyzing virus phylogenies inferred from dated genetic sequence data. By combining a phylodynamics approach based on Approximate Bayesian Computation (ABC) and an original transmission model, we estimate key epidemiological parameters of an ongoing HCV epidemic among MSMs in Lyon (France). We show that this new epidemic is largely independent of the previously observed non-MSM HCV epidemics and that its doubling time is ten times lower (0.44 years versus 4.37 years). These results have practical implications for HCV control and illustrate the additional information provided by virus genomics in public health.

Identification of Genetically Related HCV Infections Among Self-Described Injecting Partnerships

BACKGROUND

The current opioid epidemic across the United States has fueled a surge in the rate of new hepatitis C virus (HCV) infections among young persons who inject drugs (PWIDs). Paramount to interrupting transmission is targeting these high-risk populations and understanding the underlying network structures facilitating transmission within these communities.

METHODS

Deep sequencing data were obtained for 52 participants from 32 injecting partnerships enrolled in the U-Find-Out (UFO) Partner Study, which is a prospective study of self-described injecting dyad partnerships from a large community-based study of HCV infection in young adult PWIDs from San Francisco. Phylogenetically linked transmission events were identified using traditional genetic-distance measures and viral deep sequence phylogenies reconstructed to determine the statistical support of inferences and the direction of transmission within partnerships.

RESULTS

Using deep sequencing data, we found that 12 of 32 partnerships were genetically similar and clustered. Three additional phylogenetic clusters were found describing novel putative transmission links outside of the injecting relationship. Transmission direction was inferred correctly for 5 partnerships with the incorrect transmission direction inferred in more than 50% of cases. Notably, we observed that phylogenetic linkage was most often associated with a lower number of network partners and involvement in a sexual relationship.

CONCLUSIONS

Deep sequencing of HCV among self-described injecting partnerships demonstrates that the majority of transmission events originate from outside of the injecting partnership. Furthermore, these findings caution that phylogenetic methods may be unable to routinely infer the direction of transmission among PWIDs especially when transmission events occur in rapid succession within high-risk networks.