The Gap Procedure: for the identification of phylogenetic clusters in HIV-1 sequence data

From Alpha to Zeta: Identifying Variants and Subtypes of SARS-CoV-2 Via Clustering

The availability of millions of SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus-2) sequences in public databases such as GISAID (Global Initiative on Sharing All Influenza Data) and EMBL-EBI (European Molecular Biology Laboratory-European Bioinformatics Institute) (the United Kingdom) allows a detailed study of the evolution, genomic diversity, and dynamics of a virus such as never before. Here, we identify novel variants and subtypes of SARS-CoV-2 by clustering sequences in adapting methods originally designed for haplotyping intrahost viral populations. We asses our results using clustering entropy-the first time it has been used in this context. Our clustering approach reaches lower entropies compared with other methods, and we are able to boost this even further through gap filling and Monte Carlo-based entropy minimization. Moreover, our method clearly identifies the well-known Alpha variant in the U.K. and GISAID data sets, and is also able to detect the much less represented (<1% of the sequences) Beta (South Africa), Epsilon (California), and Gamma and Zeta (Brazil) variants in the GISAID data set. Finally, we show that each variant identified has high selective fitness, based on the growth rate of its cluster over time. This demonstrates that our clustering approach is a viable alternative for detecting even rare subtypes in very large data sets.

Factors influencing estimates of HIV-1 infection timing using BEAST

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3–9 days) and gag (IQR: 5–9 days), whilst the genome placed it at a median of 10 days (IQR: 4–19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

Molecular dating using phylogenetics allows us to estimate the date of an infection from time-stamped within-host sequences alone. There are large datasets of HIV-1 sequences, but genome and gene analyses are not often performed in parallel and rarely with the possibility to compare results against a known narrow window of infection. We showed that all but the longest genes are near-clonal in acute infection, with little information for dating purposes. For infections with single founders, we estimated the eclipse phase—the time between HIV-1 exposure and the first positive diagnostic test—to last between one and two weeks using env, gag, pol and near full-length genomes. This approach could be used to narrow the date of suspected infection in ongoing clinical trials for the prevention of HIV-1 infection.

Factors influencing HIV-1 phylogenetic clustering

PURPOSE OF REVIEW

A major goal of public health in relation to HIV/AIDS is to prevent new transmissions in communities. Phylogenetic techniques have improved our understanding of the structure and dynamics of HIV transmissions. However, there is still no consensus about phylogenetic methodology, sampling coverage, gene target and/or minimum fragment size.

RECENT FINDINGS

Several studies use a combined methodology, which includes both a genetic or patristic distance cut-off and a branching support threshold to identify phylogenetic clusters. However, the choice about these thresholds remains an inherently subjective process, which affects the results of these studies. There is still a lack of consensus about the genomic region and the size of fragments that should be used, although there seems to be emerging a consensus that using longer segments, allied with the use of a realistic model of evolution and a codon alignment, increases the likelihood of inferring true transmission clusters. The pol gene is still the most used genomic region, but recent studies have suggested that whole genomes and/or sequences from nef and gp41 are also good targets for cluster reconstruction.

SUMMARY

The development and application of standard methodologies for phylogenetic clustering analysis will advance our understanding of factors associated with HIV transmission. This will lead to the design of more precise public health interventions.

Molecular network-based intervention brings us closer to ending the HIV pandemic

Precise identification of HIV transmission among populations is a key step in public health responses. However, the HIV transmission network is usually difficult to determine. HIV molecular networks can be determined by phylogenetic approach, genetic distance-based approach, and a combination of both approaches. These approaches are increasingly used to identify transmission networks among populations, reconstruct the history of HIV spread, monitor the dynamics of HIV transmission, guide targeted intervention on key subpopulations, and assess the effects of interventions. Simulation and retrospective studies have demonstrated that these molecular network-based interventions are more cost-effective than random or traditional interventions. However, we still need to address several challenges to improve the practice of molecular network-guided targeting interventions to finally end the HIV epidemic. The data remain limited or difficult to obtain, and more automatic real-time tools are required. In addition, molecular and social networks must be combined, and technical parameters and ethnic issues warrant further studies.

Molecular dating and viral load growth rates suggested that the eclipse phase lasted about a week in HIV-1 infected adults in East Africa and Thailand

Most HIV-1 infected individuals do not know their infection dates. Precise infection timing is crucial information for studies that document transmission networks or drug levels at infection. To improve infection timing, we used the prospective RV217 cohort where the window when plasma viremia becomes detectable is narrow: the last negative visit occurred a median of four days before the first detectable HIV-1 viremia with an RNA test, referred below as diagnosis. We sequenced 1,280 HIV-1 genomes from 39 participants at a median of 4, 32 and 170 days post-diagnosis. HIV-1 infections were dated by using sequence-based methods and a viral load regression method. Bayesian coalescent and viral load regression estimated that infections occurred a median of 6 days prior to diagnosis (IQR: 9–3 and 11–4 days prior, respectively). Poisson-Fitter, which analyzes the distribution of hamming distances among sequences, estimated a median of 7 days prior to diagnosis (IQR: 15–4 days) based on sequences sampled 4 days post-diagnosis, but it did not yield plausible results using sequences sampled at 32 days. Fourteen participants reported a high-risk exposure event at a median of 8 days prior to diagnosis (IQR: 12 to 6 days prior). These different methods concurred that HIV-1 infection occurred about a week before detectable viremia, corresponding to 20 days (IQR: 34–15 days) before peak viral load. Together, our methods comparison helps define a framework for future dating studies in early HIV-1 infection.

HIV-1 infected individuals rarely know when they became infected but knowing when an infection occurred provides critical information regarding HIV-1 pathogenesis and epidemiology. Using a unique cohort in which infection was known to have occurred in a narrow interval, we investigated methods to estimate the timing of infections. Several methods suggested that HIV-1 infection typically occurs a median of one week before the infection can be detected by HIV-1 RNA testing. Going forward, we provide a strategy that can be used to elucidate the origin of an acute/early infection.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

OutbreakFinder: a visualization tool for rapid detection of bacterial strain clusters based on optimized multidimensional scaling

With the evolution of next generation sequencing (NGS) technologies, whole-genome sequencing of bacterial isolates is increasingly employed to investigate epidemiology. Phylogenetic analysis is the common method for using NGS data, usually for comparing closeness between bacterial isolates to detect probable outbreaks. However, interpreting a phylogenetic tree is not easy without training in evolutionary biology. Therefore, developing an easy-to-use tool that can assist people who wish to use a phylogenetic tree to investigate epidemiological relatedness is crucial. In this paper, we present a tool called OutbreakFinder that can accept a distance matrix in csv format; alignment files from Lyve-SET, Parsnp, and ClustalOmega; and a tree file in Newick format as inputs to compute a cluster-labeled two-dimensional plot based on multidimensional-scaling dimension reduction coupled with affinity propagation clustering. OutbreakFinder can be downloaded for free at https://github.com/skypes/Newton-method-MDS.

Assessing the role of transmission chains in the spread of HIV-1 among men who have sex with men in Quebec, Canada

BACKGROUND

Phylogenetics has been used to investigate HIV transmission among men who have sex with men. This study compares several methodologies to elucidate the role of transmission chains in the dynamics of HIV spread in Quebec, Canada.

METHODS

The Quebec Human Immunodeficiency Virus (HIV) genotyping program database now includes viral sequences from close to 4,000 HIV-positive individuals classified as Men who have Sex with Men (MSMs), collected between 1996 and early 2016. Assessment of chain expansion may depend on the partitioning scheme used, and so, we produce estimates from several methods: the conventional Bayesian and maximum likelihood-bootstrap methods, in combination with a variety of schemes for applying a maximum distance criterion, and two other algorithms, DM-PhyClus, a Bayesian algorithm that produces a measure of uncertainty for proposed partitions, and the Gap Procedure, a fast non-phylogenetic approach. Sequences obtained from individuals in the Primary HIV Infection (PHI) stage serve to identify incident cases. We focus on the period ranging from January 1st 2012 to February 1st 2016.

RESULTS AND CONCLUSION

The analyses reveal considerable overlap between chain estimates obtained from conventional methods, thus leading to similar estimates of recent temporal expansion. The Gap Procedure and DM-PhyClus suggest however moderately different chains. Nevertheless, all estimates stress that longer older chains are responsible for a sizeable proportion of the sampled incident cases among MSMs. Curbing the HIV epidemic will require strategies aimed specifically at preventing such growth.

DM-PhyClus: a Bayesian phylogenetic algorithm for infectious disease transmission cluster inference

Background

Conventional phylogenetic clustering approaches rely on arbitrary cutpoints applied a posteriori to phylogenetic estimates. Although in practice, Bayesian and bootstrap-based clustering tend to lead to similar estimates, they often produce conflicting measures of confidence in clusters. The current study proposes a new Bayesian phylogenetic clustering algorithm, which we refer to as DM-PhyClus (Dirichlet-Multinomial Phylogenetic Clustering), that identifies sets of sequences resulting from quick transmission chains, thus yielding easily-interpretable clusters, without using any ad hoc distance or confidence requirement.

Results

Simulations reveal that DM-PhyClus can outperform conventional clustering methods, as well as the Gap procedure, a pure distance-based algorithm, in terms of mean cluster recovery. We apply DM-PhyClus to a sample of real HIV-1 sequences, producing a set of clusters whose inference is in line with the conclusions of a previous thorough analysis.

Conclusions

DM-PhyClus, by eliminating the need for cutpoints and producing sensible inference for cluster configurations, can facilitate transmission cluster detection. Future efforts to reduce incidence of infectious diseases, like HIV-1, will need reliable estimates of transmission clusters. It follows that algorithms like DM-PhyClus could serve to better inform public health strategies.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2347-3) contains supplementary material, which is available to authorized users.

Genotypic and Phylogenetic Insights on Prevention of the Spread of HIV-1 and Drug Resistance in "Real-World" Settings

HIV continues to spread among vulnerable heterosexual (HET), Men-having-Sex with Men (MSM) and intravenous drug user (IDU) populations, influenced by a complex array of biological, behavioral and societal factors. Phylogenetics analyses of large sequence datasets from national drug resistance testing programs reveal the evolutionary interrelationships of viral strains implicated in the dynamic spread of HIV in different regional settings. Viral phylogenetics can be combined with demographic and behavioral information to gain insights on epidemiological processes shaping transmission networks at the population-level. Drug resistance testing programs also reveal emergent mutational pathways leading to resistance to the 23 antiretroviral drugs used in HIV-1 management in low-, middle- and high-income settings. This article describes how genotypic and phylogenetic information from Quebec and elsewhere provide critical information on HIV transmission and resistance, Cumulative findings can be used to optimize public health strategies to tackle the challenges of HIV in “real-world” settings.

Biased phylodynamic inferences from analysing clusters of viral sequences

Phylogenetic methods are being increasingly used to help understand the transmission dynamics of measurably evolving viruses, including HIV. Clusters of highly similar sequences are often observed, which appear to follow a ‘power law’ behaviour, with a small number of very large clusters. These clusters may help to identify subpopulations in an epidemic, and inform where intervention strategies should be implemented. However, clustering of samples does not necessarily imply the presence of a subpopulation with high transmission rates, as groups of closely related viruses can also occur due to non-epidemiological effects such as over-sampling. It is important to ensure that observed phylogenetic clustering reflects true heterogeneity in the transmitting population, and is not being driven by non-epidemiological effects. We qualify the effect of using a falsely identified ‘transmission cluster’ of sequences to estimate phylodynamic parameters including the effective population size and exponential growth rate under several demographic scenarios. Our simulation studies show that taking the maximum size cluster to re-estimate parameters from trees simulated under a randomly mixing, constant population size coalescent process systematically underestimates the overall effective population size. In addition, the transmission cluster wrongly resembles an exponential or logistic growth model 99% of the time. We also illustrate the consequences of false clusters in exponentially growing coalescent and birth-death trees, where again, the growth rate is skewed upwards. This has clear implications for identifying clusters in large viral databases, where a false cluster could result in wasted intervention resources.

Comparative Evaluation of Subtyping Tools for Surveillance of Newly Emerging HIV-1 Strains

HIV-1 non-B subtypes/circulating recombinant forms (CRFs) are increasing worldwide. Since subtype identification can be clinically relevant, we assessed the added value in HIV-1 subtyping using updated molecular phylogeny (Mphy) and the performance of routinely used automated tools. Updated Mphy (2015 updated reference sequences), used as a gold standard, was performed to subtype 13,116 HIV-1 protease/reverse transcriptase sequences and then compared with previous Mphy (reference sequences until 2014) and with COMET, REGA, SCUEAL, and Stanford subtyping tools. Updated Mphy classified subtype B as the most prevalent (73.4%), followed by CRF02_AG (7.9%), C (4.6%), F1 (3.4%), A1 (2.2%), G (1.6%), CRF12_BF (1.2%), and other subtypes (5.7%). A 2.3% proportion of sequences were reassigned as different subtypes or CRFs because of misclassification by previous Mphy. Overall, the tool most concordant with updated Mphy was Stanford-v8.1 (95.4%), followed by COMET (93.8%), REGA-v3 (92.5%), Stanford-old (91.1%), and SCUEAL (85.9%). All the tools had a high sensitivity (≥98.0%) and specificity (≥95.7%) for subtype B. Regarding non-B subtypes, Stanford-v8.1 was the best tool for C, D, and F subtypes and for CRFs 01, 02, 06, 11, and 36 (sensitivity, ≥92.6%; specificity, ≥99.1%). A1 and G subtypes were better classified by COMET (92.3%) and REGA-v3 (98.6%), respectively. Our findings confirm Mphy as the gold standard for accurate HIV-1 subtyping, although Stanford-v8.1, occasionally combined with COMET or REGA-v3, represents an effective subtyping approach in clinical settings. Periodic updating of HIV-1 reference sequences is fundamental to improving subtype characterization in the context of an effective epidemiological surveillance of non-B strains.

Defining HIV-1 transmission clusters based on sequence data

Understanding HIV-1 transmission dynamics is relevant to both screening and intervention strategies of HIV-1 infection. Commonly, HIV-1 transmission chains are determined based on sequence similarity assessed either directly from a sequence alignment or by inferring a phylogenetic tree. This review is aimed at both nonexperts interested in understanding and interpreting studies of HIV-1 transmission, and experts interested in finding the most appropriate cluster definition for a specific dataset and research question. We start by introducing the concepts and methodologies of how HIV-1 transmission clusters usually have been defined. We then present the results of a systematic review of 105 HIV-1 molecular epidemiology studies summarizing the most common methods and definitions in the literature. Finally, we offer our perspectives on how HIV-1 transmission clusters can be defined and provide some guidance based on examples from real life datasets.

Large cluster outbreaks sustain the HIV epidemic among MSM in Quebec

OBJECTIVE

HIV-1 epidemics among MSM remain unchecked despite advances in treatment and prevention paradigms. This study combined viral phylogenetic and behavioural risk data to better understand underlying factors governing the temporal growth of the HIV epidemic among MSM in Quebec (2002-2015).

METHODS

Phylogenetic analysis of pol sequences was used to deduce HIV-1 transmission dynamics (cluster size, size distribution and growth rate) in first genotypes of treatment-naïve MSM (2002-2015, n = 3901). Low sequence diversity of first genotypes (0-0.44% mixed base calls) was used as an indication of early-stage infection. Behavioural risk data were obtained from the Montreal rapid testing site and primary HIV-1-infection cohorts.

RESULTS

Phylogenetic analyses uncovered high proportion of clustering of new MSM infections. Overall, 27, 45, 48, 53 and 57% of first genotypes within one (singleton, n = 1359), 2-4 (n = 692), 5-9 (n = 367), 10-19 (n = 405) and 20+ (n = 1277) cluster size groups were early infections (<0.44% diversity). Thirty viruses within large 20+ clusters disproportionately fuelled the epidemic, representing 13, 25 and 42% of infections, first genotyped in 2004-2007 (n = 1314), 2008-2011 (n = 1356) and 2012-2015 (n = 1033), respectively. Of note, 35, 21 and 14% of MSM belonging to 20+, 2-19 and one (singleton) cluster groups were under 30 years of age, respectively. Half of persons seen at the rapid testing site (2009-2011, n = 1781) were untested in the prior year. Poor testing propensity was associated with fewer reported partnerships.

CONCLUSION

Addressing the heterogeneity in transmission dynamics among HIV-1-infected MSM populations may help guide testing, treatment and prevention strategies.

Impacts and shortcomings of genetic clustering methods for infectious disease outbreaks

For infectious diseases, a genetic cluster is a group of closely related infections that is usually interpreted as representing a recent outbreak of transmission. Genetic clustering methods are becoming increasingly popular for molecular epidemiology, especially in the context of HIV where there is now considerable interest in applying these methods to prioritize groups for public health resources such as pre-exposure prophylaxis. To date, genetic clustering has generally been performed with ad hoc algorithms, only some of which have since been encoded and distributed as free software. These algorithms have seldom been validated on simulated data where clusters are known, and their interpretation and similarities are not transparent to users outside of the field. Here, I provide a brief overview on the development and inter-relationships of genetic clustering methods, and an evaluation of six methods on data simulated under an epidemic model in a risk-structured population. The simulation analysis demonstrates that the majority of clustering methods are systematically biased to detect variation in sampling rates among subpopulations, not variation in transmission rates. I discuss these results in the context of previous work and the implications for public health applications of genetic clustering.