Abstract
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.
Introduction
HIV continues to be a significant cause of morbidity and mortality in the United States (US) (). Public health officials and providers are interested in inferring transmission links between individuals with HIV to inform and improve treatment and prevention approaches (). In the absence of reliable patient contact histories, phylogenetic analysis of HIV sequence data can and has been used to infer transmission clusters (), under the assumption that two individuals sharing a most recent common ancestor in a phylogeny are more likely to share or lead to an epidemiological link in the real, unobservable transmission network. The application of molecular epidemiology and cluster inference techniques in public health interventions to disrupt transmission was delineated as one of the four key pillars for ending the HIV epidemic in the US ().
While historically phylogenetic informativeness of the HIV pol genomic region was suggested and contested (; ), its use is now widespread in cluster inference, often due to the availability of sequences from guideline-recommended routine drug resistance testing, typically performed by commercial Sanger sequencing (). In a recent review of HIV cluster inference, 98 out of 105 (93%) analyzed the pol region ().
The increasing availability of NGS technology has led to longer (across more genes) and deeper (multiple reads that correspond to multiple within-host genomes) sequencing of HIV, and data sets more routinely cover nearly the whole genome at great depth (). Recent evidence suggests improvements in both phylogenetic analysis and cluster inference from longer near-whole-genome HIV sequences obtained with NGS. For example, found that the accuracy of phylogenetic reconstruction and cluster inference on simulated sequences improved with longer genomic regions (with the best accuracy from a gag-pol-env concatenation). similarly studied effects on cluster inference of using longer genomic regions from near-whole-genome publicly available Sanger sequences and found that the proportion of sequences in clusters increased with longer sequences. Even before the availability of NGS, using longer regions of the HIV genome was shown to improve phylogenetic reconstruction. In one of the earliest studies of HIV sequence data with a known HIV transmission network, found that combining data from the gag and env regions improved the accuracy of phylogenetic reconstruction.
While potential advantages of longer NGS sequences in inferring clusters have been examined (; ), advantages of deeper sequencing are less investigated, and whether it can improve HIV molecular clustering inference is unknown. This is due to limitations in established practices of inferring HIV phylogenies across hosts. Researchers often rely on a single consensus sequence for each host that discard all but the majority variant at each site, since most molecular epidemiology approaches require a single fully-resolved sequence per individual in the phylogeny. Accordingly, researchers studying HIV transmission networks discard available information on within-host variation, known to impact phylogenetic inference (; ).
The consensus approach, which to date has been employed with Sanger sequencing data in multiple studies of HIV molecular epidemiology (), carries an underlying statistical assumption of low relative entropy (). For HIV, this is equivalent to the strong assumption that a consensus sequence adequately captures all relevant information about HIV diversity within an individual and that variation within hosts has no information about relationships across hosts. While many researchers understand that this assumption is likely wrong and intra-host variation is relevant for phylogenetic analysis of HIV [for a recent review, see ], in practice researchers have faced limitations in data collection that prevent measuring intra-host variation or in available analysis methods that preserve intra-host variation during alignment and phylogeny. With the advent of long read sequencing technologies for full HIV genomes, obtaining fully resolved sequences that represent the within-host viral population will be possible, but methods to incorporate intra-host variation for transmission cluster analysis will still need to be developed.
Two previous studies have accounted for within-host variation in deeply-sequenced NGS data with coalescent evolutionary models (; ), but such models still assume a consensus sequence as the observed data. Two other studies introduced methods to use deeply-sequenced HIV data without assuming a consensus, for a different but related epidemiological goal of estimating transmission directionality and identifying multiple infections (; ). Methods also exist that combine haplotype estimation from deeply-sequenced NGS data and phylogenetics () as a way to incorporate within-host diversity, but available haplotyping methods have a high computational cost and results are often not sufficiently accurate for cluster analysis (). Additionally, all aforementioned methods that do not rely on a consensus incorporate within-host diversity by including multiple sequences or tips per sample, which presents difficulties with summarizing or collapsing the resulting phylogenetic tree in order to identify transmission clusters and measure cluster certainty.
In this study, we develop a new method we call profile sampling that incorporates within-host HIV genome variation into phylogenetic analyses used to identify transmission clusters. We examine if, and to what extent, incorporation of within-host variation available from deeply-sequenced Illumina-based NGS data provides improved phylogenetic inference and clustering relative to traditional consensus-sequence-based approaches. We focus our analyses on all newly HIV-diagnosed individuals during six months from the largest HIV center in Rhode Island, US.
Materials and Methods
Data Collection and Sequencing
HIV-1 pol Sanger sequences (HXB2 positions 2253-3554), available through clinical care, were collected from the 37 adults (18 years) newly-diagnosed with HIV-1 during the first six months of 2013 and treated at The Miriam Hospital Immunology Center in Providence, Rhode Island, US. Patients at this Center represent ∼80% of the state’s HIV epidemic.
In addition, blood was obtained from consenting participants and processed to isolate RNA from plasma (n = 27), and proviral DNA from whole blood (n = 4) or peripheral blood mononuclear cells (PBMC; n = 6). Using Sanger sequencing and Illumina-based NGS, near-whole-genome viral sequences were obtained from one compartment; plasma for participants with detectable viral load and proviral DNA for participants with undetectable viral load or unsuccessful plasma genotyping. The study was approved by the Institutional Review Board at Lifespan, which is the parent health network of The Miriam Hospital.
Total nucleic acids were extracted and an in-house genotyping assay was used to generate the near-whole genome sequence (wgs), based on previously published methods (; ). For each sample, two cDNA templates were generated by SuperscriptIII First Strand Synthesis System (Thermofisher, Carlsbad, CA, United States), followed by eight separate nested PCR reactions; these eight amplicons span the near-whole HIV genome. Final amplicon products were sequenced by Sanger using the 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA, United States) and by NGS using Nextera XT DNA Library Prep chemistry (Illumina, San Diego, CA, United States) to generate multiplexed libraries for Illumina’s MiSeq platform with 250 base paired-end reads. Sanger consensus sequences were generated manually using Sequencher version 5.2.4 (Gene Codes, Ann Arbor, MI, United States) to confirm degenerate nucleotides. NGS data were processed and demultiplexed using BaseSpace cloud application (Illumina, San Diego, CA, United States). NGS consensus sequences were called at a 20% threshold.
Profile Sampling
We introduce a new approach for incorporating within-host viral variation into phylogenetic analysis, called profile sampling (Figure 1). Profile sampling builds upon existing methods of phylogenetic and cluster inference by also sampling from within-host viral diversity. We start by aligning each individual’s HIV NGS reads using the hivmmer pipeline (), which we developed and now extended to support near-whole-genome HIV data and perform codon-aware alignment within each gene (hivmmer version 0.2.1). A key feature of this pipeline is its use of profile hidden Markov models (HMMs) to model and align collections of HIV sequences. Profile HMMs have been abundantly used for biological sequence analyses and are particularly well-suited to modeling variation in populations of sequences (). Briefly, hivmmer performs quality control and error correction in overlapping regions of read pairs using PEAR version 0.9.11 (), translates them into possible reading frames, aligns them in amino acid space to profile HMMs of all HIV-1 group M reference sequences () using the profile HMM alignment tool HMMER version 3.1b2 (), and produces a codon frequency table across the near-whole HIV genome. We refer to this resulting codon frequency table as the individual’s HIV profile.
FIGURE 1
Subsequently, we sample 500 fully-resolved sequences from each of the 37 individuals’ HIV profile according to the frequency of observed codons in the profile, for a total of 37 × 500 = 18,500 profile-sampled sequences. These sequences do not correspond to real strains present in the biological sample, but do capture the empirical distribution of within-host variation at the individual codon level. We note that the sequences do not capture linkage across codons, which is important for the detection and elimination of recombinant HIV strains as part of quality control, but is unessential for phylogenetic analyses. We then collate the 18,500 sequences into 500 profile-sampled data sets, by sampling without replacement so that each data set has 37 sequences (one for each individual) and can be used in a phylogenetic analysis with existing methods.
We also use the 18,500 sequences to estimate the within-host diversity for the 37 individuals as the average percent difference in nucleotides across all pairwise comparisons of each individual’s 500 profile-sampled nucleotide sequences. These pairwise differences are calculated using the Hamming distance () [also called p-distance (; )].
Phylogenetic Inference
For profile sampling, we perform phylogenetic inference of wgs (HXB2 positions 790-9417) on each of the 500 profile-sampled data sets by estimating a multiple sequence alignment with OMM_MACSE version 10.02 () and a maximum-likelihood phylogeny with the GTRCAT model and 100 rapid bootstrap replicates using RAxML version 8.2.12 (), with HIV-1 group O (GenBank accession L20587.1) as the outgroup. We perform this same phylogenetic inference on three clinically-relevant sub-genomic regions: protease and reverse transcriptase at the beginning of the pol gene (“prrt”, HXB2 positions 2253-3554), int gene (HXB2 positions 4230-5096), and env gene (HXB2 positions 6225-8790). The prrt and int regions are routinely sequenced in clinical care to detect drug resistance and inform clinical anti-retroviral therapy choices. The env region is sequenced to genotypically infer viral tropism and co-receptor usage. Cluster inference is performed on all phylogenies using Cluster Picker () with a threshold of 99% bootstrap support.
In addition to profile sampling, we infer phylogenies with similar tools, regions and parameters for the NGS consensus sequences at the 20% threshold and the Sanger sequences. We perform cluster inference on all consensus phylogenies using a similar method as for profile sampling. We do not impose a genetic distance threshold because empirically-justified thresholds that are comparable across the near-whole-genome and the int and env regions have not to our knowledge been established. This approach of using only bootstrap criteria for cluster detection is consistent with methods commonly used in the broader HIV cluster analysis literature ().
We investigate the impact of within-host diversity on phylogenetic topology and evolutionary distance estimates in the sub-genomic and near-whole-genome regions. To examine variation in topology, we first calculate pairwise geodesic distance (; ) among the 500 phylogenies from the profile samples, as well as phylogenies from the NGS consensus and Sanger sequences. Then we perform multi-dimensional scaling on the resulting distance matrix to visualize topological space in two dimensions. Next, to examine variation in estimated evolutionary distance, we sum the branch lengths within each phylogeny across all branches and across only tip branches and visualize the distribution of these branch length sums. These analyses establish to what extent phylogenies from consensus sequences (which are point estimates) summarize the underlying variation in two important aspects of the phylogeny: estimated topology and estimated evolutionary distance.
Finally, we examine the clusters that are detected in phylogenies of NGS consensus sequences versus Sanger sequences, and across the four genomic regions and their profile sampling support, using the frequency that a cluster appears across the 500 profile-sampled phylogenies. We refer to this value as the profile-sampled support and note that it can be conceived as analogical to the conventional bootstrap support for evaluating robustness of an individual phylogeny’s topology but extends that feature to evaluating robustness of cluster detection using within-host sequence variation.
All analysis source code is available from https://github.com/kantorlab/hiv-profile-sampling.
Results
Profile Sampling Estimates of Within-Host Diversity
Figure 2 shows the estimated within-host percent diversity in each examined genomic region across individuals, ordered by env, which we expected a priori to be the most variable region. The largest estimated diversity is in env for individual MC28 (3.9%), and env has the overall largest estimated diversity range (0.2–3.9%, mean 1.5%). The other regions have ranges of 0.2–1.9% (mean 0.9%) for prrt, 0.1–2.0% (mean 0.9%) for int, and 0.2–2.6% (mean 1.2%) for wgs. Such within-host estimations are not feasible with conventional consensus Sanger or NGS approaches, although methods such as phyloscanner and HAPHPIPE that utilize deeper NGS sequencing to build phylogenies with multiple tips per sample, are able to also quantify within-host diversity (; ).
FIGURE 2
Phylogenetic Estimates Are Sensitive to Within-Host Diversity
Figure 3 demonstrates multi-dimensional scaling on the profile-sampled phylogenies and the phylogenies from Sanger and NGS consensus sequences within each genomic region. The profile sampling approach reveals for each genomic region a multi-modal topological space in which phylogenies inferred from both Sanger and NGS consensus sequences are outliers; a result that is confirmed by multi-dimensional scaling across all regions (Figure 4). A key difference between the consensus and profile-sampled sequences is that consensus sequences contain ambiguous nucleotides at sites with ≥2 nucleotides by Sanger Sequencing base calling or with ≥20% frequency for NGS. In contrast, profile-sampled sequences by construction have no ambiguous sites, and ambiguity is instead incorporated into analyses through frequency of the ambiguous nucleotides across the 500 samples.
FIGURE 3
FIGURE 4
Figure 5 shows the distribution of branch length sums across compared phylogenies. Overall, estimates are larger in env and wgs, and smaller when restricting to only tip branches. In some cases, consensus phylogenies provide an adequate summary of the distribution (as in the phylogeny of the NGS consensus sequence for tip branches for wgs). In other cases, consensus phylogenies have estimates that are outliers in the distribution (as in the phylogenies from NGS and Sanger consensus sequences for all branches in wgs and env).
FIGURE 5
Taken together, the heterogeneity between the phylogenetic results from profile sampling and consensus-inferred point estimates demonstrate that within-host virus sequence diversity impacts the inference of virus phylogeny across individuals, and that the consensus approach to handling ambiguity and collapsing within-host sequence variation can obscure both the magnitude and effect of these impacts.
Profile-Sampled Cluster Support Differs by Sequencing Depth and Genomic Region
Combining the results of all examined methods (Sanger, NGS consensus, NGS profile sampling) and genomic regions (prrt, int, env, wgs) there were overall 12 identified clusters among the 37 participants. Seven clusters had two members, four had three members, and one had five members. Figure 6 demonstrates comparison of cluster detection by examined methods and genomic regions. Some clusters had consistently high support (>75%) across all regions (e.g., MC25/MC26/MC52 and MC14/MC59). Other clusters had higher support in certain regions (e.g., MC17/MC20/MC21 in env and wgs). Eight clusters across different regions were detected by profile sampling but not by consensus methods, while all clusters detected by consensus methods were detected by profile sampling. One larger cluster, MC37/MC41/MC47/MC53/MC56, was detected only by profile sampling with the wgs dataset.
FIGURE 6
By providing previously-unavailable cluster support that considers within-host “deep” viral variation, profile sampling in the wgs dataset allowed detection of the largest (all 12) overall number of clusters. The clusters detected in wgs also had the highest overall profile-sampled support, as compared to the other genomic regions. The median profile-sampled support was 99.8% for wgs, 77.6% for env, 41.4% for int, and 51.3% for prrt. The profile-sampled support for wgs was significantly larger than for int (p-value = 0.005, Dunn’s test of multiple comparisons using paired rank sums with Holm-Bonferroni correction) and prrt (p-value = 0.010), but not significantly larger than for env (p-value = 0.092).
The phylogenies of NGS consensus sequences detected only six clusters in prrt, seven in int, seven in env, and seven in wgs (Figure 7). The phylogenies of Sanger sequences detected only four clusters in prrt, six in int, seven in env, and seven in wgs (Figure 8). Only one cluster (MC25/MC26/MC52) was consistently detected across phylogenies from NGS and Sanger consensus sequences, and across all regions.
FIGURE 7
FIGURE 8
The median profile-sampled support for clusters detected by Sanger consensus sequences (green and yellow cells, Figure 6) was 60.2%, not different than for those detected by NGS consensus sequences (72.5%; orange and yellow cells in Figure 6; p-value = 0.415, Wilcoxon signed-rank test). Totaling the clusters detected across the four regions, phylogenies of Sanger consensus sequences detected fewer clusters (27) than phylogenies of NGS consensus sequences (31) or profile sampling (43); and detected fewer clusters in each region except env. Cluster support values were higher for clusters detected by phylogenies of both NGS and Sanger sequences (yellow cells, Figure 6; median cluster support 98.5%) than those detected only by one or the other (orange or green cells, Figure 6; median cluster support 68.1%).
Discussion
Current phylogenetic approaches to inference of HIV transmission clusters utilize consensus sequences to summarize within-host sequence variation. We introduce a different summarization strategy, profile sampling, that preserves the within-host sequence variation provided by the deeper sequencing that is now widely available. In a dataset of all newly HIV-diagnosed individuals over six months at the largest HIV center in Rhode Island, United States, deeper sequencing provided by NGS and incorporated by the newly-introduced profile sampling captured within-host diversity, revealing clusters detected by profile sampling but not by consensus approaches, including one larger cluster found only with profile sampling of the wgs. This suggests that routinely used consensus sequence approaches discard potentially relevant information present in NGS data, and that considering this additional information in phylogenetic analysis may improve the robustness of HIV cluster detection. Profile sampling can thus provide a new quantitative measure of cluster confidence with potential applications to public health activities. Such activities could be better justified in scenarios where clusters triggering them have high cluster support from deep-sequenced data, though this was not addressed here and needs to be demonstrated.
Profile sampling complements well-established bootstrapping methods, and in some senses is orthogonal to them. Phylogenetic inference depends on a sequence alignment, where each row corresponds to a single host and each column corresponds to a given genomic position. Bootstrapping resamples columns of the matrix with replacement, giving an indication of how consistent signal is across genomic positions. Profile sampling, on the other hand, resamples each site in the alignment given the sequence diversity observed within each host. This gives an indication of how consistent signal is across HIV genomes within each host. Variation in within-host diversity could be due to a variety of biological factors, like viral mutation, effective viral population size, and time since infection, as well as technical factors like sequencing depth and sequencing errors. Not accommodating this variation could lead to overconfidence in results or missed clinically relevant phylogenetic signals.
Although the standard practice of collapsing within-host diversity into a single consensus sequence simplifies downstream analyses, the results presented here demonstrate that this practice discards potentially relevant biological results and may mislead phylogenetic analysis and resulting epidemiological consequences. For example, public health activities triggered by phylogenetic inference of HIV molecular clustering to inform and improve prevention and treatment interventions can be affected (). In our data, clusters vary in their profile-sampled support, and consensus approaches can fail to detect clusters supported by deep-sequenced data, as in the case of the largest cluster, which was detected by profile sampling, not by consensus approaches. As data acquisition increasingly shifts to NGS approaches (), it is important to compare results of larger datasets from new methods to the more common conventional Sanger pol consensus sequences.
Much of the enthusiasm about shifting from Sanger sequencing to NGS has been due to reducing costs and the ability to more easily collect data on the entire HIV genome rather than few genes. Our results suggest that much benefit from NGS may also come from its greater sequencing depth, capturing viral sequence variation within individuals. This benefit can only be realized, though, if this variation is propagated to phylogenetic analyses, such as by profile sampling introduced here, rather than being collapsed to a consensus sequence, as is conventionally done. We suggest creating a profile that captures that variation, performing multiple phylogenetic analyses on sequences sampled from the profile, and then summarizing the phylogenetic analyses. This summary method can also take advantage of the output from long-read sequencing technologies which are able to provide fully resolved sequences from the viral population and which are starting to replace short-read NGS sequencing in a number of HIV labs. Future tools could incorporate the variation directly into the phylogenetic inference process itself ().
In our comparison of cluster inference across genomic regions, we found that fewer clusters were detected overall in prrt and int compared to env and wgs. Prior studies of clustering from Sanger consensus sequences present mixed results on prevalence of clustering across genomic regions. Some studies found concordant clustering across gag-env () and gag-pol-env (; ), while others found fewer clusters in pol than in env (), or in gag-env than in pol (). The additional information available in deep-sequenced NGS data, along with cluster support measures provided by profile sampling, could help resolve differences, as suggested here. In our data, not only were more clusters detected in the near-whole length genome, but those clusters also had higher cluster support as measured by deep-sequenced NGS data. While prior studies demonstrated better accuracy of cluster inference on simulated NGS sequences when using wgs () and that the proportion of Sanger sequences in clusters increased with longer sequence regions (), we have demonstrated here that deeply-sequenced, near-whole length NGS data can be used with profile sampling to detect clusters undetectable by consensus approaches. The impact of this approach for public health remains to be determined.
One limitation of our study is the small number of participants and the short timeframe they were enrolled in. However, participants represent a dense temporal sampling, and comprise all newly HIV-diagnosed individuals in a six-month period at the largest HIV center in Rhode Island, in which 80% of the state’s people with HIV are cared for. The overall size of the HIV epidemic in Rhode Island was estimated as 2,396 individuals in 2016 (), but NGS data for this population are not currently available beyond those presented here. In future work, we will apply profile sampling to larger NGS data sets, to assess cluster inference concordance between Sanger and NGS data, and its impact on public health actions to halt HIV transmission. Additionally, we do not know the real number of clusters or the true transmission chains, a limitation with all studies on HIV transmission networks. Our construction of HIV profiles from NGS data is also limited by the accuracy of the NGS assays themselves. The codon frequencies in the profiles may be biased measures of the true within-host diversity because of biases in PCR amplification, as well as a variety of technical factors related to next-generation sequencing and analysis [see for a detailed discussion]. Sequencing protocols such as Primer ID (; ) have been introduced to reduce and correct for these biases, and should be considered in the future.
In conclusion, the true HIV transmission network is unknown, but phylogenetic analysis and cluster inference are promising tools for aiding clinicians and public health officials in better understanding and in disrupting HIV transmission (). Most current phylogenetic approaches do not fully utilize the information on within-host diversity available in deep-sequenced, near-whole-genome NGS data. As NGS data sets are increasingly available and become more representative of HIV epidemics, we suggest that the additional information they measure has the potential to improve the robustness of HIV molecular cluster inference, the impact of which needs to be further investigated.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this article are not readily available because of patient privacy concerns. Requests to access the datasets should be directed to RK at rkantor@brown.edu.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Review Board at Lifespan. The patients/participants provided their written informed consent to participate in this study.
Author contributions
AG, MH, PC, CD, and RK contributed to the conception of the study. LL, MD’A, PC, and RK contributed to data collection and organization. AG and CL created the method. AG and MH contributed to implementation, analysis and visualization. AG, MH, CD, and RK contributed to writing the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was supported by the National Institutes of Health (R01AI136058, K24AI134359, P30AI042853 and P20GM109035) and a Brown University DEANS Award.
Acknowledgments
We thank Dr. Mia Coetzer for assistance with next-generation sequencing and feedback on an earlier draft of the manuscript. This research was conducted using computational resources and services at the Center for Computation and Visualization, Brown University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AllamO.SamaraniS.AhmadA. (2011). Hammering out HIV-1 incidence with Hamming distance.AIDS252047–2048. 10.1097/QAD.0b013e32834bac66
2
BendallM. L.GibsonK. M.SteinerM. C.RentiaU.Pérez-LosadaM.CrandallK. A. (2021). HAPHPIPE: haplotype reconstruction and phylodynamics for deep sequencing of intrahost viral populations.Mol. Biol. Evol.381677–1690. 10.1093/molbev/msaa315
3
BilleraL. J.HolmesS. P.VogtmannK. (2001). Geometry of the space of phylogenetic trees.Adv. Appl. Math.27733–767.
4
Di GiallonardoF.TöpferA.ReyM.PrabhakaranS.DuportY.LeemannC.et al (2014). Full-length haplotype reconstruction to infer the structure of heterogeneous virus populations.Nucleic Acids Res.42:e115. 10.1093/nar/gku537
5
EddyS. R. (2004). What is a hidden Markov model?Nat. Biotechnol.221315–1316.
6
EddyS. R. (2011). Accelerated profile HMM searches.PLoS Comput. Biol.7:e1002195. 10.1371/journal.pcbi.1002195
7
EnglishS.KatzourakisA.BonsallD.FlanaganP.DudaA.FidlerS.et al (2011). Phylogenetic analysis consistent with a clinical history of sexual transmission of HIV-1 from a single donor reveals transmission of highly distinct variants.Retrovirology8:54. 10.1186/1742-4690-8-54
8
FauciA. S.LaneH. C. (2020). Four decades of HIV/AIDS – much accomplished, much to do.N. Engl. J. Med.3831–4. 10.1056/NEJMp1916753
9
FauciA. S.RedfieldR. R.SigounasG.WeahkeeM. D.GiroirB. P. (2019). Ending the HIV epidemic: a plan for the United States.JAMA321:844. 10.1001/jama.2019.1343
10
GiardinaF.Romero-SeversonE. O.AlbertJ.BrittonT.LeitnerT. (2017). Inference of transmission network structure from HIV phylogenetic trees.PLoS Comput. Biol.13:e1005316. 10.1371/journal.pcbi.1005316
11
GuangA.ZapataF.HowisonM.LawrenceC. E.DunnC. W. (2016). An integrated perspective on phylogenetic workflows.Trends Ecol. Evol.31116–126. 10.1016/j.tree.2015.12.007
12
HanZ.LeungT. W.ZhaoJ.WangM.FanL.LiK.et al (2009). A HIV-1 heterosexual transmission chain in Guangzhou, China: a molecular epidemiological study.Virol. J.6:148. 10.1186/1743-422X-6-148
13
HassanA. S.PybusO. G.SandersE. J.AlbertJ.EsbjörnssonJ. (2017). Defining HIV-1 transmission clusters based on sequence data.AIDS311211–1222. 10.1097/QAD.0000000000001470
14
HogbenM.CollinsD.HootsB.O’ConnorK. (2016). Partner services in sexually transmitted disease prevention programs: a review.Sex. Transm. Dis.43S53–S62. 10.1097/OLQ.0000000000000328
15
HowisonM.CoetzerM.KantorR. (2019). Measurement error and variant-calling in deep Illumina sequencing of HIV.Bioinformatics352029–2035. 10.1093/bioinformatics/bty919
16
HuéS.ClewleyJ. P.CaneP. A.PillayD. (2004). HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy.AIDS18719–728. 10.1097/00002030-200403260-00002
17
JabaraC. B.JonesC. D.RoachJ.AndersonJ. A.SwanstromR. (2011). Accurate sampling and deep sequencing of the HIV-1 protease gene using a Primer ID.Proc. Natl. Acad. Sci. U.S.A.10820166–20171. 10.1073/pnas.1110064108
18
KantorR. (2021). Next generation sequencing for HIV-1 drug resistance testing-a special issue walkthrough.Viruses13:340. 10.3390/v13020340
19
KapaataA.LyagobaF.SsemwangaD.MagamboB.NanyonjoM.LevinJ.et al (2013). HIV-1 subtype distribution trends and evidence of transmission clusters among incident cases in a rural clinical cohort in Southwest Uganda, 2004–2010.AIDS Res. Hum. Retroviruses29520–527. 10.1089/AID.2012.0170
20
KayeM.ChiboD.BirchC. (2008). Phylogenetic investigation of transmission pathways of drug-resistant HIV-1 utilizing pol sequences derived from resistance genotyping.J. Acquir. Immune Defic. Syndr.499–16. 10.1097/QAI.0b013e318180c8af
21
LeitnerT. (2019). Phylogenetics in HIV transmission: taking within-host diversity into account.Curr. Opin. HIV AIDS14181–187. 10.1097/COH.0000000000000536
22
LeitnerT.EscanillaD.FranzenC.UhlenM.AlbertJ. (1996). Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.Proc. Natl. Acad. Sci. U.S.A.9310864–10869. 10.1073/pnas.93.20.10864
23
Los Alamos National Lab (2020). HIV Database [Internet].Available online at: http://www.hiv.lanl.gov/(accessed February 9, 2020).
24
MaldarelliF.KearneyM.PalmerS.StephensR.MicanJ.PolisM. A.et al (2013). HIV populations are large and accumulate high genetic diversity in a nonlinear fashion.J. Virol.8710313–10323. 10.1128/JVI.01225-12
25
NadaiY.EyzaguirreL. M.ConstantineN. T.SillA. M.CleghornF.BlattnerW. A.et al (2008). Protocol for nearly full-length sequencing of HIV-1 RNA from plasma.PLoS One3:e1420. 10.1371/journal.pone.0001420
26
NdiayeH. D.TchiakpeE.VidalN.NdiayeO.DiopA. K.PeetersM.et al (2013). HIV type 1 subtype c remains the predominant subtype in men having sex with men in Senegal.AIDS Res. Hum. Retroviruses291265–1272. 10.1089/aid.2013.0140
27
NovitskyV.MoyoS.LeiQ.DeGruttolaV.EssexM. (2015). Importance of viral sequence length and number of variable and informative sites in analysis of HIV clustering.AIDS Res. Hum. Retroviruses31531–542. 10.1089/AID.2014.0211
28
OwenM.ProvanJ. S. (2011). A fast algorithm for computing geodesic distances in tree space.IEEE/ACM Trans. Comput. Biol. Bioinformatics82–13. 10.1109/TCBB.2010.3
29
Panel on Antiretroviral Guidelines for Adults and Adolescents (2020). Guidelines for the Use of Antiretroviral Agents in Adults and Adolescents with HIV.Available online at: https://clinicalinfo.hiv.gov/sites/default/files/guidelines/documents/AdultandAdolescentGL.pdf(accessed February 9, 2020).
30
PetersP. J.PontonesP.HooverK. W.PatelM. R.GalangR. R.ShieldsJ.et al (2016). HIV infection linked to injection use of Oxymorphone in Indiana, 2014–2015.N. Engl. J. Med.375229–239. 10.1056/NEJMoa1515195
31
Ragonnet-CroninM.HodcroftE.HuéS.FearnhillE.DelpechV.BrownA. J.et al (2013). Automated analysis of phylogenetic clusters.BMC Bioinformatics14:317. 10.1186/1471-2105-14-317
32
RanwezV.DouzeryE. J. P.CambonC.ChantretN.DelsucF. (2018). MACSE v2: toolkit for the alignment of coding sequences accounting for frameshifts and stop codons.Mol. Biol. Evol.352582–2584. 10.1093/molbev/msy159
33
Rhode Island Department of Health (2019). HIV Progress [Internet].Available online at: https://health.ri.gov/data/hiv/(accessed December 12, 2019).
34
Romero-SeversonE.SkarH.BullaI.AlbertJ.LeitnerT. (2014). Timing and order of transmission events is not directly reflected in a pathogen phylogeny.Mol. Biol. Evol.312472–2482. 10.1093/molbev/msu179
35
SkumsP.ZelikovskyA.SinghR.GusslerW.DimitrovaZ.KnyazevS.et al (2018). QUENTIN: reconstruction of disease transmissions from viral quasispecies genomic data.Bioinformatics34163–170. 10.1093/bioinformatics/btx402
36
StamatakisA. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.Bioinformatics301312–1313. 10.1093/bioinformatics/btu033
37
StürmerM.PreiserW.GuteP.NisiusG.DoerrH. W. (2004). Phylogenetic analysis of HIV-1 transmission: pol gene sequences are insufficient to clarify true relationships between patient isolates.AIDS182109–2113. 10.1097/00002030-200411050-00002
38
VoelkerdingK. V.DamesS. A.DurtschiJ. D. (2009). Next-generation sequencing: from basic research to diagnostics.Clin. Chem.55641–658. 10.1373/clinchem.2008.112789
39
WymantC.HallM.RatmannO.BonsallD.GolubchikT.de CesareM.et al (2018). PHYLOSCANNER: inferring transmission from within- and between-host pathogen genetic diversity.Mol. Biol. Evol.25719–733. 10.1093/molbev/msx304
40
YebraG.HodcroftE. B.Ragonnet-CroninM. L.PillayD.BrownA. J. L.et alPANGEA_HIV Consortium (2016). Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic.Sci. Rep.6:39489. 10.1038/srep39489
41
ZhangJ.KobertK.FlouriT.StamatakisA. (2014). PEAR: a fast and accurate Illumina paired-end read merger.Bioinformatics30614–620. 10.1093/bioinformatics/btt593
42
ZhouS.SizemoreS.MoeserM.ZimmermanS.SamoffE.MobleyV.et al (2020). Near real-time identification of recent HIV transmissions, transmitted drug resistance mutations, and transmission networks by MPID-NGS in North Carolina.J. Infect. Dis.223876–884. 10.1093/infdis/jiaa417
Summary
Keywords
HIV, cluster inference, profile sampling, phylogenetics, next generation sequencing (NGS), near-whole-genome, consensus sequence, transmission disruption
Citation
Guang A, Howison M, Ledingham L, D’Antuono M, Chan PA, Lawrence C, Dunn CW and Kantor R (2022) Incorporating Within-Host Diversity in Phylogenetic Analyses for Detecting Clusters of New HIV Diagnoses. Front. Microbiol. 12:803190. doi: 10.3389/fmicb.2021.803190
Received
27 October 2021
Accepted
22 December 2021
Published
17 February 2022
Volume
12 - 2021
Edited by
Michael M. Thomson, Instituto de Salud Carlos III (ISCIII), Spain
Reviewed by
Marcos Perez Losada, George Washington University, United States; Brittany Rife Magalis, University of Florida, United States
Updates
Copyright
© 2022 Guang, Howison, Ledingham, D’Antuono, Chan, Lawrence, Dunn and Kantor.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: August Guang, august_guang@brown.edu
†These authors have contributed equally to this work and share first authorship
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.