Development and application of a high-throughput HIV type 1 genotyping assay to identify CRF02_AG in West/West Central Africa

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

HIV-1 genetic diversity and demographic characteristics in Bulgaria

HIV-1 strain diversity in Bulgaria is extensive and includes contributions from nearly all major subtypes and the Circulating Recombinant Forms (CRF): 01_AE, 02_AG, and 05_DF. Prior to this study, HIV-1 sequence information from Bulgaria has been based solely on the pro-RT gene, which represent less than 15% of the viral genome. To further characterize HIV-1 in Bulgaria, assess participant risk behaviors, and strengthen knowledge of circulating strains in the region, the study “Genetic Subtypes of HIV-1 in Bulgaria (RV240)” was conducted. This study employed the real time-PCR based Multi-region Hybridization Assay (MHA) B/non-B and HIV-1 sequencing to survey 215 of the approximately 1100 known HIV-1 infected Bulgarian adults (2008–2009) and determine if they were infected with subtype B HIV-1. The results indicated a subtype B prevalence of 40% and demonstrate the application of the MHA B/non-B in an area containing broad HIV-1 strain diversity. Within the assessed risk behaviors, the proportion of subtype B infection was greatest in men who have sex with men and lowest among those with drug use risk factors. During this study, 15 near full-length genomes and 22 envelope sequences were isolated from study participants. Phylogenetic analysis shows the presence of subtypes A1, B, C, F1, and G, CRF01_AE, CRF02_AG, CRF05_DF, and one unique recombinant form (URF). These sequences also show the presence of two strain groups containing participants with similar risk factors. Previous studies in African and Asian cohorts have shown that co-circulation of multiple subtypes can lead to viral recombination within super-infected individuals and the emergence of new URFs. The low prevalence of URFs in the presence of high subtype diversity in this study, may be the result of successful infection prevention and control programs. Continued epidemiological monitoring and support of infection prevention programs will help maintain control of the HIV-1 epidemic in Bulgaria.

Rev-RRE Functional Activity Differs Substantially Among Primary HIV-1 Isolates

The HIV-1 replication cycle requires the nucleocytoplasmic export of intron-containing viral RNAs, a process that is ordinarily restricted. HIV overcomes this by means of the viral Rev protein, which binds to an RNA secondary structure called the Rev response element (RRE) present in all unspliced or incompletely spliced viral RNA transcripts. The resulting mRNP complex is exported through interaction with cellular factors. The Rev-RRE binding interaction is increasingly understood to display remarkable structural plasticity, but little is known about how Rev-RRE sequence differences affect functional activity. To study this issue, we utilized a lentiviral vector assay in which vector titer is dependent on the activity of selected Rev-RRE pairs. We found that Rev-RRE functional activity varies significantly (up to 24-fold) between naturally occurring viral isolates. The activity differences of the Rev-RRE cognate pairs track closely with Rev, but not with RRE activity. This variation in Rev activity is not correlated with differences in Rev steady state protein levels. These data suggest that Rev sequence differences drive substantial variation in Rev-RRE functional activity between patients. Such variation may play a role in viral adaptation to different immune milieus within and between patients and may be significant in the establishment of latency. The identification of differences in Rev-RRE functional activity in naturally occurring isolates may also permit more efficient production of lentiviral vectors.

Significant contribution of subtype G to HIV-1 genetic complexity in Nigeria identified by a newly developed subtyping assay specific for subtype G and CRF02_AG

While abundant sequence information is available from human immunodeficiency virus type 1 (HIV-1) subtypes A, B, C and CRF01_AE for HIV-1 vaccine design, sequences from West Africa are less represented. We sought to augment our understanding of HIV-1 variants circulating in 6 Nigerian cities as a step to subsequent HIV-1 vaccine development.The G/CRF02_AG multi-region hybridization assay (MHA) was developed to differentiate subtype G, CRF02_AG and their recombinants from other subtypes based on 7 HIV-1 segments. Plasma from 224 HIV-1 infected volunteers enrolled in a cohort examining HIV-1 prevalence, risk factor, and subtype from Makurdi (30), Abuja (18), Enugu (11), Kaduna (12), Tafa (95), and Ojo/Lagos (58) was analyzed using MHA. HIV-1 genomes from 42 samples were sequenced to validate the MHA and fully explore the recombinant structure of G and CRF02_AG variants.The sensitivity and specificity of MHA varied between 73-100% and 90-100%, respectively. The subtype distribution as identified by MHA among 224 samples revealed 38% CRF02_AG, 28% G, and 26% G/CRF02_AG recombinants while 8% remained nontypeable strains. In envelope (env) gp120, 38.84% of the samples reacted to a G probe while 31.25% reacted to a CRF02 (subtype A) probe. Full genome characterization of 42 sequences revealed the complexity of Nigerian HIV-1 variants.CRF02_AG, subtype G, and their recombinants were the major circulating HIV-1 variants in 6 Nigerian cities. High proportions of samples reacted to a G probe in env gp120 confirms that subtype G infections are abundant and should be considered in strategies for global HIV-1 vaccine development.

An effective tool for identifying HIV-1 subtypes B, C, CRF01_AE, their recombinant forms, and dual infections in Southeast Asia by the multi-region subtype specific PCR (MSSP) assay

The RV144 Thai vaccine trial has been the only vaccine study to show efficacy in preventing HIV infection. Ongoing molecular surveillance of HIV-1 in Southeast Asia is vital for vaccine development and evaluation. In this study a novel tool, the multi-region subtype specific PCR (MSSP) assay, that was able to identify subtypes B, C, CRF01_AE for Thailand, other Southeast Asian countries, India and China is described. The MSSP assay is based on a nested PCR strategy and amplifies eight short regions distributed along the HIV-1 genome using subtype-specific primers. A panel of 41 clinical DNA samples obtained primarily from opiate users in northern Thailand was used to test the assay performance. The MSSP assay provided 73-100% sensitivity and 100% specificity for the three subtypes in each genome region. The assay was then field-tested on 337 sera from HIV infected northern Thai drug users collected between 1999 and 2002. Subtype distribution was CRF01_AE 77.4% (n=261), subtype B 3.3% (n=11), CRF01_AE/B recombinant 12.2% (n=41), CRF01_AE/C recombinant 0.6% (n=2), and non-typeable 6.5% (n=22). The MSSP assay is a simple, cost-effective, and accurate genotyping tool for laboratory settings with limited resources and is sensitive enough to capture the recombinant genomes and dual infections.

Molecular detection of HIV-1 subtype B, CRF01_AE, CRF33_01B, and newly emerging recombinant lineages in Malaysia

A molecular genotyping assay for human immunodeficiency virus type 1 (HIV-1) circulating in Southeast Asia is difficult to design because of the high level of genetic diversity. We developed a multiplex real-time polymerase chain reaction (PCR) assay to detect subtype B, CRF01_AE, CRF33_01B, and three newly described circulating recombinant forms, (CRFs) (CRF53_01B, CRF54_01B, and CRF58_01B). A total of 785 reference genomes were used for subtype-specific primers and TaqMan probes design targeting the gag, pol, and env genes. The performance of this assay was compared and evaluated with direct sequencing and phylogenetic analysis. A total of 180 HIV-infected subjects from Kuala Lumpur, Malaysia were screened and 171 samples were successfully genotyped, in agreement with the phylogenetic data. The HIV-1 genotype distribution was as follows: subtype B (16.7%); CRF01_AE (52.8%); CRF33_01B (24.4%); CRF53_01B (1.1%); CRF54_01B (0.6%); and CRF01_AE/B unique recombinant forms (4.4%). The overall accuracy of the genotyping assay was over 95.0%, in which the sensitivities for subtype B, CRF01_AE, and CRF33_01B detection were 100%, 100%, and 97.7%, respectively. The specificity of genotyping was 100%, inter-subtype specificities were > 95% and the limit of detection of 10(3) copies/mL for plasma. The newly developed real-time PCR assay offers a rapid and cost-effective alternative for large-scale molecular epidemiological surveillance for HIV-1.

Novel multiregion hybridization assay for the identification of the most prevalent genetic forms of the human immunodeficiency virus type 1 circulating in Portugal

The most efficient method for HIV-1 genetic characterization involves full-genome sequencing, but the associated costs, technical features, and low throughput preclude it from being routinely used for the analysis of large numbers of viral strains. Multiregion hybridization assays (MHA) represent an alternative for a consistent genetic analysis of large numbers of viral strains. Classically, MHA rely on the amplification by real-time PCR of several regions scattered along the HIV-1 genome, and on their characterization with clade-specific TaqMan probes (also known as hydrolysis probes). In this context, the aim of our study was the development of a technical variant of an MHA (vMHA(B/G/02)) for genotyping the most prevalent genetic forms of HIV-1 circulating in Portugal. Different sets of primers were designed for universal and clade-specific amplifications of several sections of the viral genome: gag, pol(Pr), pol(RT), vpu, env(gp120), and env(gp41). vMHA(B/G/02) was implemented using a real-time PCR-based approach, with detection dependent on the use of SYBR Green I. As an alternative, a technically less demanding strategy based on conventional PCR and agarose gel analysis of the reaction products was also developed. This method performed with overall good sensitivity and specificity (>91%) when a convenience sample of 45 plasma-derived HIV-1 strains was analyzed. Apart from the detection of subtype B, G, CRF02_AG, and CRF14_BG viruses, several unique B/G recombinant were also detected. Curiously, recombinant viruses including CRF02_AG sequences were not detected in the group of samples analyzed.

A novel multiregion hybridization assay reveals high frequency of dual inter-subtype infections among HIV-positive individuals in Cameroon, West Central Africa

In West and West Central Africa, multiple subtypes, circulating recombinant forms (CRF), and high proportions of unique recombinant forms (URF) are documented. The predominance of recombinants strongly suggests that dual infections occur frequently. In the present study, we adapted the multi-region hybridization assay (MHA), previously developed to identify dual infections in geographic regions where few HIV-1 variants circulate, to identify HIV-1 variants and dual infections. We designed clade-specific probes in three genomic regions (gag p17, vpu, nef) to detect eight different variants that are common in this part of Africa (A, B/D, C, F, G, CRF02_AG, CRF06_cpx, CRF22_01A1). The assay was validated with 163 samples representing the corresponding HIV-1 variants. Depending on the genomic regions, the global sensitivity of the assay ranged from 86% to 94%, and the global specificity was between 85% and 96%. The assay was then applied on 156 antiretroviral treatment-naive patients from Cameroon. The MHA assay identified 79%, 85% and 90% of the strains in nef, gag and vpu regions, respectively. The subtype/CRF distribution and the proportion of inter-region recombinants obtained by the new MHA assay were in accordance with known subtype/CRF distribution in Cameroon. Moreover, the MHA assay identified 35 (22.4%) patients as dually infected, from which 20 were reactive in more than one region and/or with concordant multigenomic recombination pattern. Despite the high genetic diversity, we successfully developed an hybridization assay allowing identification of eight common HIV-1 variants circulating in West and West Central Africa. We documented high rates of dual infection in a low-risk population group, illustrating that the global evolution of HIV diversity is driven by dual infections. This assay could become a useful screening tool for the global surveillance and monitoring of inter-subtype/CRF dual infections in West and West Central Africa.

CRF22_01A1 is involved in the emergence of new HIV-1 recombinants in Cameroon

Cameroon is a West African country where high genetic diversity of HIV-1 has been reported. The predominant CRF02_AG is involved in the emergence of more complex intersubtype recombinants. In this study, we sequenced the full-length genome of a novel unique recombinant form of HIV-1, 02CAMLT04 isolated in blood donors in urban Cameroon. Phylogenetic tree and bootscan analysis showed that 02CAMLT04 was complex and seemed to be a secondary recombinant derived from CRF02_AG and CRF22_01A1. The genomic composition of 02CAMLT04 strain showed that it is composed of 3 segments; 24% of the genome is classified as CRF02_AG, spanning most of the envelope gene. The remaining 76% of the genome is classified as CRF22_01A1. In addition, the sequence analysis of 13 full-length sequences from HIV-1-positive specimens received from Cameroon between 2002 and 2010 indicated that 5 specimens are pure CRF22_01A1 viruses, and 6 others have homology with CRF22_01A1 sequences in either gag, pol, or env region, whereas 6% of strains contain portions of CRF22_01A1. Further study demonstrated that CRF22_01A1 is a primary prevalence strain co-circulating in Cameroon and is involved in complex intersubtype recombination events with subtypes (D or F), subsubtypes (A1 or F2), and CRFs (CRF01_AE or CRF02_AG). Our studies show that novel recombinants between CRF22_01A1 and other clades and recombinant forms may be emerging in Cameroon that could contribute to the future global diversity of HIV-1 in this region and worldwide.

Identification and genetic characterization of a novel CRF22_01A1 recombinant form of HIV type 1 in Cameroon

Cameroon is a country in West Central Africa in which all four groups of HIV-1 (M, N, O, and P), some circulating recombinant forms (CRFs) and unique recombinant forms (URFs) are prevalent. The CRF22 was initially identified through a novel URF strain, 01CM53122, and later defined from two additional sequences; however, the genomic properties of CRF22 have never been demonstrated in detail. In this study, we describe the characterization of five CRF22_01A1 strains, 02CMLT72, 01CM1867LE, 01CM001BBY, 02CM3097MN, and 02CM1917LE, identified in Cameroon without apparent epidemiological links. A typical CRF22_01A1 strain contains five fragments that can be assigned to the CRF01_AE and subsubtype A1 radiations. Forty-eight percent of the genome is classified as CRF01_AE, spanning the entire region of the gag gene, part of the pol gene, and accessory genes as well as the beginning and the end of the env gene and nef gene. Fifty-two percent of the genome is subsubtype A1 including regions mostly in the pol, vif, and env genes. The five CRF22_01A1 viruses formed a deep branch outside the groups of CRF01_AE and displayed similar mosaic structure but were moderately different from the original strain of CRF22_01A1, 01CM53122. Further analysis of the 01CM53122 genome showed that this virus represents a diverse set of mosaic genomes from CRF22_01A1, including a 446-nt segment of 01CM53122 in the env region, but unlike other CRF22 strains, clustered with CRF01_AE rather than the A1 sequence, suggesting that the 01CM53122 strain is a recombinant of CRF22_01A1 and CRF01_AE.

HIV-1 recombinants with multiple parental strains in low-prevalence, remote regions of Cameroon: evolutionary relics?

Background

The HIV pandemic disseminated globally from Central West Africa, beginning in the second half of the twentieth century. To elucidate the virologic origins of the pandemic, a cross-sectional study was conducted of the genetic diversity of HIV-1 strains in villagers in 14 remote locations in Cameroon and in hospitalized and STI patients. DNA extracted from PBMC was PCR amplified from HIV(+) subjects. Partial pol amplicons (N = 164) and nearly full virus genomes (N = 78) were sequenced. Among the 3956 rural villagers studied, the prevalence of HIV infection was 4.9%; among the hospitalized and clinic patients, it was 8.6%.

Results

Virus genotypes fell into two distinctive groups. A majority of the genotyped strains (109/164) were the circulating recombinant form (CRF) known to be endemic in West Africa and Central West Africa, CRF02_AG. The second most common genetic form (9/164) was the recently described CRF22_01A1, and the rest were a collection of 4 different subtypes (A2, D, F2, G) and 6 different CRFs (-01, -11, -13, -18, -25, -37). Remarkably, 10.4% of HIV-1 genomes detected (17/164) were heretofore undescribed unique recombinant forms (URF) present in only a single person. Nearly full genome sequencing was completed for 78 of the viruses of interest. HIV genetic diversity was commonplace in rural villages: 12 villages each had at least one newly detected URF, and 9 villages had two or more.

Conclusions

These results show that while CRF02_AG dominated the HIV strains in the rural villages, the remainder of the viruses had tremendous genetic diversity. Between the trans-species transmission of SIV_cpz and the dispersal of pandemic HIV-1, there was a time when we hypothesize that nascent HIV-1 was spreading, but only to a limited extent, recombining with other local HIV-1, creating a large variety of recombinants. When one of those recombinants began to spread widely (i.e. became epidemic), it was recognized as a subtype. We hypothesize that the viruses in these remote Cameroon villages may represent that pre-epidemic stage of viral evolution.

Thermodynamically modulated partially double-stranded linear DNA probe design for homogeneous real-time PCR

Real-time PCR assays have recently been developed for diagnostic and research purposes. Signal generation in real-time PCR is achieved with probe designs that usually depend on exonuclease activity of DNA polymerase (e.g. TaqMan probe) or oligonucleotide hybridization (e.g. molecular beacon). Probe design often needs to be specifically tailored either to tolerate or to differentiate between sequence variations. The conventional probe technologies offer limited flexibility to meet these diverse requirements. Here, we introduce a novel partially double-stranded linear DNA probe design. It consists of a hybridization probe 5′-labeled with a fluorophore and a shorter quencher oligo of complementary sequence 3′-labeled with a quencher. Fluorescent signal is generated when the hybridization probe preferentially binds to amplified targets during PCR. This novel class of probe can be thermodynamically modulated by adjusting (i) the length of hybridization probe, (ii) the length of quencher oligo, (iii) the molar ratio between the two strands and (iv) signal detection temperature. As a result, pre-amplification signal, signal gain and the extent of mismatch discrimination can be reliably controlled and optimized. The applicability of this design strategy was demonstrated in the Abbott RealTime HIV-1 assay.

The complexity of circulating HIV type 1 strains in Oyo state, Nigeria

Multiple HIV-1 subtypes and circulating recombinant forms (CRFs) are known to circulate in West Africa. We undertook a survey of HIVs in Oyo state, in southwestern Nigeria. We analyzed 71 samples from Ibadan, the capital city, and 33 samples from Saki, 100 miles west of Ibadan. We sequenced part of the gag gene and the envelope C2V3 region from 102 and 89 samples, respectively. In the 87 samples for which both genes were sequenced, subtype G and CRF02_AG were found in equal proportions (32.2% each). Other samples included CRF06_cpx (8.0%), subtype A (2.3%), C (1.1%), unclassified (1.1%), or discordant sequences suggesting the presence of a large number of recombinants involving CRF02_AG and/or subtype G (20.7%) or other subtypes (2.3%). The subtype/CRF designation was concordant in two gene fragments in the majority of samples evaluated. However, we observed differences in subtype distribution between the two locations with a predominance of subtype G in Ibadan and CRF02 in Saki. This is the first in-depth analysis of HIV variability at a state level in Nigeria. Our analysis revealed a significant level of viral heterogeneity and a geographical difference in subtype distribution, and demonstrated that CRF02_AG does not account for the majority of circulating strains.

Identification of a novel circulating recombinant form (CRF) 36_cpx in Cameroon that combines two CRFs (01_AE and 02_AG) with ancestral lineages of subtypes A and G

An array of CRFs have been identified in Cameroon, the most notable being CRF02_AG. HIV-1 in the East Province of Cameroon is particularly diverse: in a recent study, we found a high proportion of unique recombinant forms (URFs). Herein we describe the analysis of the full-length sequences of two of these URFs, which, after preliminary analysis of gag, pol, and env fragments, appeared to be a novel CRF. This novel strain, CRF36_cpx, contains fragments that can be assigned to the CRF01_AE, CRF02_AG, and subtype A and G radiations. Forty percent of the genome can be classified as CRF02_AG, including regions in gag, pol, env, and the accessory genes. Twenty-seven percent is CRF01_AE, comprising the majority of gag, the beginning of env, and the end of env into the 3' LTR. Twenty percent of the genome can be assigned to subtype A, with segments in pol and env. The remaining 13% of the sequence is classifiable as subtype G, in pol and vpu. The subtype A and G lineages formed by the CRF36_cpx sequences are unique and appear ancestral in nature. CRF36_cpx is both the first to combine more than one CRF and the first to include fragments of CRF02_AG. The ancestral sequences present in CRF36_cpx represent a link to extinct strains, and, potentially, insight into the evolution of HIV-1.

Identification of conserved lentiviral sequences as landmarks of genomic flexibility

Considering that recombinations produce quasispecies in lentivirus spreading, we identified and localized highly conserved sequences that may play an important role in viral ontology. Comparison of entire genomes, including 237 human, simian and non-primate mammal lentiviruses and 103 negative control viruses, led to identify 28 Conserved Lentiviral Sequences (CLSs). They were located mainly in the structural genes forming hot spots particularly in the gag and pol genes and to a lesser extent in LTRs and regulatory genes. The CLS pattern was the same throughout the different HIV-1 subtypes, except for some HIV-1-O strains. Only CLS 3 and 4 were detected in both negative control HTLV-1 oncornaviruses and D-particle-forming simian viruses, which are not immunodeficiency inducers and display a genetic stability. CLSs divided the virus genomes into domains allowing us to distinguish sequence families leading to the notion of ‘species self’ besides that of ‘lentiviral self’. Most of acutely localized CLSs in HIV-1s (82%) corresponded to wide recombination segments being currently reported. To cite this article: M.L.J. Moncany et al., C. R. Biologies 329 (2006).

Molecular epidemiology of HIV Type 1 in preparation for a Phase III prime-boost vaccine trial in Thailand and a new approach to HIV Type 1 genotyping

To characterize HIV-1 genotypes in candidate populations for a prime-boost phase III vaccine trial in Thailand, specimens from prevalent and incident HIV-1 infections from a family planning clinic population in Rayong Province and a community cohort in Chon Buri Province, collected from 1998 to 2001, were genotyped. A new multiregion hybridization assay, MHAbce, capable of distinguishing HIV-1 CRF01_AE, subtype B, and subtype C and their recombinants, was developed and applied to prevalent infections. Most incident and selected prevalent infections were studied by complete genome sequencing. By MHAbce, 168 of 194 prevalent infections were genotyped. Of these, 90.5% were CRF01_AE, 2.4% were subtype B, and 7.2% showed discordant or dual probe reactivity, indicative of recombination or dual infection, respectively. Among 23 incident infections, 20 were sequenced. Eighteen CRF01_AE, one subtype B, and one CRF01/B recombinant strains were seen. Two CRF01/B and one CRF01/C recombinant were identified among selected prevalent infections. These results indicate that incident and prevalent HIV-1 infections in Rayong and Chon Buri during 1998-2001 were 90% CRF01_AE, 3% subtype B, and 7% either recombinant or dual. This study frames the genetic diversity of HIV-1 in these cohorts in their preparatory phase for the ongoing ALVACHIV (vCP1521) prime, AIDSVAX B/E boost, phase III vaccine trial and will provide a benchmark for interpretation and analysis.

Circulating and unique recombinant forms of HIV type 1 containing subsubtype A2

HIV-1 strains containing subsubtype A2 are relatively rare in the pandemic but have been repeatedly identified in Kenya, where candidate vaccines based in part on subtype A, but not A2 strains, may be evaluated. Among the most recent is CRF16_A2D, a circulating recombinant form (CRF) whose prototypes are complete or partial HIV-1 sequences from Kenya, Korea, and Argentina. Using samples from blood bank discards in Kenya and complete genome sequencing, this report further documents CRF16_A2D and related recombinants and identifies a second CRF, CRF21_A2D. The two A2-containing CRFs, and two recombinants related to CRF16_A2D, share common structural elements but appear to have been independently derived. Concerted selection may have influenced the emergence and spread of certain A2-containing strains in Kenya. The second complete subtype C sequence from Kenya is also reported here. Monitoring of A2-containing recombinants and subtype C strains, both relatively rare in Kenya, may be informative in the course of cohort development and evaluation of candidate vaccines.

A jumping profile Hidden Markov Model and applications to recombination sites in HIV and HCV genomes

BACKGROUND

Jumping alignments have recently been proposed as a strategy to search a given multiple sequence alignment A against a database. Instead of comparing a database sequence S to the multiple alignment or profile as a whole, S is compared and aligned to individual sequences from A. Within this alignment, S can jump between different sequences from A, so different parts of S can be aligned to different sequences from the input multiple alignment. This approach is particularly useful for dealing with recombination events.

RESULTS

We developed a jumping profile Hidden Markov Model (jpHMM), a probabilistic generalization of the jumping-alignment approach. Given a partition of the aligned input sequence family into known sequence subtypes, our model can jump between states corresponding to these different subtypes, depending on which subtype is locally most similar to a database sequence. Jumps between different subtypes are indicative of intersubtype recombinations. We applied our method to a large set of genome sequences from human immunodeficiency virus (HIV) and hepatitis C virus (HCV) as well as to simulated recombined genome sequences.

CONCLUSION

Our results demonstrate that jumps in our jumping profile HMM often correspond to recombination breakpoints; our approach can therefore be used to detect recombinations in genomic sequences. The recombination breakpoints identified by jpHMM were found to be significantly more accurate than breakpoints defined by traditional methods based on comparing single representative sequences.

In-depth analysis of a heterosexually acquired human immunodeficiency virus type 1 superinfection: evolution, temporal fluctuation, and intercompartment dynamics from the seronegative window period through 30 months postinfection

Human immunodeficiency virus type 1 (HIV-1) superinfection refers to the acquisition of another strain by an already infected individual. Here we report a comprehensive genetic analysis of an HIV-1 superinfection acquired heterosexually. The infected individual was in a high-risk cohort in Tanzania, was exposed to multiple subtypes, and was systematically evaluated every 3 months with a fluorescent multi-region genotyping assay. The subject was identified in the window period and was first infected with a complex ACD recombinant strain, became superinfected 6 to 9 months later with an AC recombinant, and was monitored for >2.5 years. The plasma viral load exceeded 400,000 copies/ml during the first 9 months of infection but resolved to the set point of 67,000 copies/ml by 3 months after superinfection; the CD4 cell count was 377 cells/mul at 30 months. Viral diversity was evaluated with techniques designed to fully sample the quasi-species, permitting direct observation of the evolution, temporal fluctuation, and intercompartment dynamics of the initial and superinfecting strains and recombinants derived from them. Within 3 months of superinfection, seven different molecular forms were detected in gag and six were detected in env. The proportions of forms fluctuated widely over time in plasma and peripheral blood mononuclear cells, illustrating how challenging the detection of dually infected individuals can be. Strain-specific nested PCR confirmed that the superinfecting strain was not present until the 9 month follow-up. This study further defines the parameters and dynamics of superinfection and will foster appropriate studies and approaches to gain a more complete understanding of risk factors for superinfection and its impact on clinical progression, epidemiology, and vaccine design.

Biologic and genetic characterization of a panel of 60 human immunodeficiency virus type 1 isolates, representing clades A, B, C, D, CRF01_AE, and CRF02_AG, for the development and assessment of candidate vaccines

A critical priority for human immunodeficiency virus type 1 (HIV-1) vaccine development is standardization of reagents and assays for evaluation of immune responses elicited by candidate vaccines. To provide a panel of viral reagents from multiple vaccine trial sites, 60 international HIV-1 isolates were expanded in peripheral blood mononuclear cells and characterized both genetically and biologically. Ten isolates each from clades A, B, C, and D and 10 isolates each from CRF01_AE and CRF02_AG were prepared from individuals whose HIV-1 infection was evaluated by complete genome sequencing. The main criterion for selection was that the candidate isolate was pure clade or pure circulating recombinant. After expansion in culture, the complete envelope (gp160) of each isolate was verified by sequencing. The 50% tissue culture infectious dose and p24 antigen concentration for each viral stock were determined; no correlation between these two biologic parameters was found. Syncytium formation in MT-2 cells and CCR5 or CXCR4 coreceptor usage were determined for all isolates. Isolates were also screened for neutralization by soluble CD4, a cocktail of monoclonal antibodies, and a pool of HIV-1-positive patient sera. The panel consists of 49 nonsyncytium-inducing isolates that use CCR5 as a major coreceptor and 11 syncytium-inducing isolates that use only CXCR4 or both coreceptors. Neutralization profiles suggest that the panel contains both neutralization-sensitive and -resistant isolates. This collection of HIV-1 isolates represents the six major globally prevalent strains, is exceptionally large and well characterized, and provides an important resource for standardization of immunogenicity assessment in HIV-1 vaccine trials.