jpHMM at GOBICS: a web server to detect genomic recombinations in HIV-1

Identification of a novel first-generation HIV-1 circulating recombinant form (CRF152_DG) among people living with HIV in Karachi, Pakistan

The objective of this study was to characterize a novel circulating recombinant form of human immunodeficiency virus type 1 (HIV-1) among people living with HIV in Karachi, Pakistan. We conducted near-full-length genome (NFLG) sequencing on eight samples exhibiting D/G recombination signals in the pol gene region. We successfully obtained NFLG sequences (790-9,614; with reference to the HXB2 genome) from four of the eight samples and then conducted phylogenetic and recombination analyses on them. The four NFLG sequences from our study and one DG unique recombinant form previously identified in the United Kingdom (GenBank accession: MF109700) formed a distinct monophyletic cluster with an Shimodaira-Hasegawa approximate likelihood ratio test node support value of 100%. Bootscan analyses of the five NFLG sequences of DG recombinants showed that all five NFLGs shared the same unique mosaic pattern of recombination breakpoints between D and G clades, with two D fragments in the pol and vif regions inserted into a G backbone. Subregion phylogenetic analyses confirmed these sequences to be a novel circulating recombinant form (CRF) composed of subtypes D and G. The DG recombinant sequences were eventually designated as CRF152_DG by the Los Alamos HIV Sequence Database staff.

IMPORTANCE

In Pakistan, the genetic diversity of human immunodeficiency virus type 1 (HIV-1) is becoming increasingly complex, compared to the early years of the epidemic that started after the detection of the first cases of HIV-1 in 1987 in Karachi. Based on the available molecular studies, two dominant HIV-1 clades, sub-subtype A1 and CRF02_AG, have been found to co-circulate with other clades, namely B, C, D, G, CRF01_AE, CRF35_A1D, and CRF56_cpx, in various urban areas of Pakistan. Several novel recombinant forms have also been detected. This first report of CRF152_DG highlights the complex nature of the HIV epidemic in Pakistan and emphasizes the importance of continual molecular surveillance (ideally based on whole-genome sequences) of HIV.

Molecular Epidemiology of HIV-1 in Eastern Europe and Russia

The HIV epidemic in Eastern Europe and Russia is large and not well-controlled. To describe the more recent molecular epidemiology of HIV-1, transmitted drug resistance, and the relationship between the epidemics in this region, we sequenced the protease and reverse transcriptase genes of HIV-1 from 812 people living with HIV from Ukraine (n = 191), Georgia (n = 201), and Russia (n = 420) before the initiation of antiretroviral therapy. In 190 Ukrainian patients, the integrase gene sequence was also determined. The most reported route of transmission was heterosexual contact, followed by intravenous drug use, and men having sex with men (MSM). Several pre-existing drug resistance mutations were found against non-nucleoside reverse transcriptase inhibitors (RTIs) (n = 103), protease inhibitors (n = 11), and nucleoside analogue RTIs (n = 12), mostly polymorphic mutations or revertants. In the integrase gene, four strains with accessory integrase strand transfer inhibitor mutations were identified. Sub-subtype A6 caused most of the infections (713/812; 87.8%) in all three countries, including in MSM. In contrast to earlier studies, no clear clusters related to the route of transmission were identified, indicating that, within the region, the exchange of viruses among the different risk groups may occur more often than earlier reported.

First complete-genome documentation of HIV-1 intersubtype superinfection with transmissions of diverse recombinants over time to five recipients

Human immunodeficiency virus type 1 (HIV-1) recombinants in the world are believed to be generated through recombination between distinct HIV-1 strains among coinfection or superinfection cases. However, direct evidence to support transmission of HIV-1 recombinants from a coinfected/superinfected donor to putative recipient is lacking. Here, we report on the origin and evolutionary relationship between a set of recombinants from a CRF01_AE/CRF07_BC superinfected putative donor and diverse CRF01_AE/CRF07_BC recombinants from five putative recipients. Interviews on sociodemographic characteristics and sexual behaviors for these six HIV-1-infected men who have sex with men showed that they had similar ways of partner seeking: online dating sites and social circles. Phylogenetic and recombination analyses demonstrated that the near-full-length genome sequences from six patients formed a monophyletic cluster different from known HIV-1 genotypes in maximum likelihood phylogenetic trees, were all composed of CRF01_AE and CRF07_BC fragments with two common breakpoints on env, and shared 4-7 breakpoints with each other. Moreover, 3' half-genomes of recombinant strains from five recipients had identical/similar recombinant structures with strains at longitudinal samples from the superinfected donor. Recombinants from the donor were paraphyletic, whereas five recipients were monophyletic or polyphyletic in the maximum clade credibility tree. Bayesian analyses confirmed that the estimated time to the most recent common ancestor (tMRCA) of CRF01_AE and CRF07_BC strains of the donor was 2009.2 and 2010.7, respectively, and all were earlier than the emergence of recombinants from five recipients. Our results demonstrated that the closely related unique recombinant forms of HIV-1 might be the descendent of a series of recombinants generated gradually in a superinfected patient. This finding highlights the importance of early initiation of antiretroviral therapy as well as tracing and testing of partners in patients with multiple HIV-1 infection.

Temporal evolution and global spread of hepatitis B virus genotype G

Hepatitis B virus (HBV) infection is considered a major health problem in the world. HBV is classified into genotypes A to J disseminated worldwide. Genotypes A, D and F are the most frequent in the Western World, B and C are predominant in the East, and E, F, H and J are infrequent and restricted to specific regions. HBV-G is a rare genotype, but it has been detected in different continents. This study aimed to report the temporal evolution and global spread of HBV-G comparing whole-genome sequences of this genotype from different regions in the world. Bayesian coalescent analysis was performed to estimate the time to the most recent common ancestor (tMRCA) and the population dynamics in the last decades. The results demonstrated that tMRCA of all HBV-Gs dated back to 1855 (95% highest posterior density interval [HPD 95%]: 1778 - 1931). This genotype has a possible origin in North America and it was disseminated to other continents (South and Central America, Europe, Asia and Africa) more than one century later (around the 1970s). The viral population demonstrated constant spreading from 1855 to the 1980s, followed by an increase in the 1990s and reached a plateau after the 2000s. Wide spreading at the beginning of the 1990s was probably associated with the dissemination by highly sexual active groups and injecting drug users. In conclusion, the present study demonstrated that HBV-G was originated in the 19th century with main events of spread at the end of the 20th century.

Molecular Epidemiology of HIV-1 Infected Migrants Followed up in Portugal: Trends between 2001-2017

Migration is associated with HIV-1 vulnerability. Objectives: To identify long-term trends in HIV-1 molecular epidemiology and antiretroviral drug resistance (ARV) among migrants followed up in Portugal Methods: 5177 patients were included between 2001 and 2017. Rega, Scuel, Comet, and jPHMM algorithms were used for subtyping. Transmitted drug resistance (TDR) and Acquired drug resistance (ADR) were defined as the presence of surveillance drug resistance mutations (SDRMs) and as mutations of the IAS-USA 2015 algorithm, respectively. Statistical analyses were performed. Results: HIV-1 subtypes infecting migrants were consistent with the ones prevailing in their countries of origin. Over time, overall TDR significantly increased and specifically for Non-nucleoside reverse transcriptase inhibitor (NNRTIs) and Nucleoside reverse transcriptase inhibitor (NRTIs). TDR was higher in patients from Mozambique. Country of origin Mozambique and subtype B were independently associated with TDR. Overall, ADR significantly decreased over time and specifically for NRTIs and Protease Inhibitors (PIs). Age, subtype B, and viral load were independently associated with ADR. Conclusions: HIV-1 molecular epidemiology in migrants suggests high levels of connectivity with their country of origin. The increasing levels of TDR in migrants could indicate an increase also in their countries of origin, where more efficient surveillance should occur.

Homogenous HIV-1 subtype B from the Brazilian Amazon with infrequent diverse BF1 recombinants, subtypes F1 and C among blood donors

In the last decade a growing HIV/AIDS epidemic with increased incidence and AIDS-related mortality has been reported in Northern Brazil from which molecular data are scarce. Also, apparently healthy, adult blood donors, recently diagnosed with HIV-1 represent important sentinel populations for molecular studies. This cross-sectional study describes HIV-1 subtypes in blood donors from three reference public blood centers located in three States in Northern Brazil. HIV-1 pol sequencing (protease/PR, reverse transcriptase/RT) was performed on plasma samples of HIV-1 positive donors from HEMOAM, Manaus, Amazonas (n = 198), HEMERON, Porto Velho, Rondônia (n = 20) and HEMORAIMA, Boa Vista, Roraima (n = 9) collected from 2011-2017. HIV-1 subtypes were identified by REGA, phylogenetic inference; recombinant viruses were characterized by SIMPLOT. Young, single, males predominated, around half was first-time donors. Syphilis co-infection was detected in 17% (39 out of 227), 8% (18 out of 227) was anti-HBc positive. Subtype B represented ≥ 90% in Amazonas, Rondônia and Roraima, subtype C (3.1%) was found in Amazonas and Rondônia; subtype F1 (0.9%) and BF1 recombinants (5.3%) were only detected in Amazonas. Subtype B sequences from Amazonas (n = 179), Rondônia (n = 18) and Roraima (n = 9) were combined with viral strains representative of the BPANDEMIC (n = 300) and BCARIBBEAN/BCAR (n = 200) lineages. The BPANDEMIC lineage predominated (78%) although BCAR lineages were frequent in Roraima (56%) and Amazonas (22%). Subtype C and subtype F1 sequences identified here clustered within Brazilian CBR and F1BR lineages, respectively. Twelve BF1 mosaics showed 11 different recombination profiles: six were singleton unique-recombinant-forms/URFs, one displays a CRF28/29_BF-like recombinant pattern and the remaining four BF1 isolates branched with other Brazilian BF1 viruses previously described and may represent putative new CRF_BF1 from Northern Brazil. Our study shows a highly homogeneous molecular pattern with prevalent subtype B, followed by BF1, and sporadic subtype C and F1 in blood donors from the Northern region. Surveillance studies are important to monitor HIV-1 diversity which can reveal patterns of viral dissemination, especially in a highly endemic, remote and geographically isolated region as Northern Brazil.

Detection of human immunodeficiency virus Type 1 phylogenetic clusters with multidrug resistance mutations among 2011 to 2017 blood donors from the highly endemic Northern Brazilian Amazon

BACKGROUND

This study describes transmitted drug resistance (TDR) in blood donors diagnosed with human immunodeficiency virus Type 1 (HIV-1) infection from 2011 to 2017 in three reference public blood centers from the Northern Brazilian Amazon.

STUDY DESIGN AND METHODS

This was a cross-sectional study on HIV-positive blood donors from HEMOAM, Manaus, Amazonas, AM (n = 198); HEMERON, Porto Velho, Rondônia, RO (n = 20); and HEMORAIMA, Boa Vista, Roraima, RR (n = 9). HIV-1 pol sequences (protease, reverse transcriptase) were analyzed for drug resistance mutations (DRMs) using the Calibrated Population Resistance tool (Stanford). TDR/DRM clusters were investigated by phylogenetic analysis after removing positions associated with drug resistance of Subtype B sequences from untreated and treated subjects from Northern Brazil.

RESULTS

Transmitted drug resistance/DRM in blood donors was 11% (25 of 227), all of them from HEMOAM. Most blood donors with TDR/DRM had multiple and similar DRMs. Nonnucleoside reverse transcriptase inhibitor (NNRTI) mutations predominated (10.1%), followed by nucleoside reverse transcriptase inhibitor (NRTI) mutations (5.3%) and protease inhibitor mutations (0.4%). Dual-class NNRTI/NRTI mutations represented 4.8%. Three highly supported Subtype B monophyletic clades mostly composed by individuals from Amazonas with TDR/DRM mutations were identified. The largest transmission cluster contained 10 sequences, eight from HEMOAM and two sequences described previously (one from a treated subject from Amazonas and the other one from Roraima). This cluster was characterized by NRTI (D67N, T69D, T215S/F/L, K219Q) and NNRTI (K101H, K103 N, G190A) mutations. The other two transmission clades comprised only three and two sequences from HEMOAM sharing the E138A NNRTI mutation.

CONCLUSIONS

The identification of transmission clusters of multidrug-resistant viruses in blood donors from Amazonas highlight the need of continued monitoring of TDR/DRM and the importance of pretreatment genotyping in the highly endemic Amazonas state.

Antiretroviral resistance, genotypic characterization and origin of Human Immunodeficiency Virus among the infected wives of Intravenous drug users in Manipur

Increasing incidence of drug resistance is ascertained to be the main obstacles in limiting the virus among the human immunodeficiency virus (HIV) infected individuals. This study investigates the drug resistance mutations (DRMs), genetic variants and origin of transmitted drug resistance of HIV-1 among the HIV-1 infected wives of intravenous drug users (IDUs) in Manipur. 44 HIV pol gene sequences were generated from 56 blood samples by viral gene amplification and sequencing. Sequences were then analysed for drug resistance, genetic variants and origin. The result revealed that among the treatment naive cases, 35.7% had Transmitted Drug Resistance Mutations (TDRMs) while among treatment experienced cases, 50% had Acquired Drug Resistant Mutations (ADRMs). These TDRMs and ADRMs conferred resistance to nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs) and/or protease inhibitors (PIs). Majority of the isolated HIV-1 sequences (77.3%) were subtype C while 9.1% was discordant subtype, 6.8% was subtype B, 4.5% was CRF_01AE and 2.3% was URF_BC. TDRM strains were found to be introduced from Myanmar, Vietnam and mainland India. This study also reveals the appearance of CRF_01AE for the first time in Manipur. The finding of this study indicates high prevalence of drug resistant mutations and complex molecular epidemiology in Manipur.

Genetic characterization of a novel unique recombinant form of HIV-1 originating from subtype A1 and C in Guangdong Province, China

The dominant human immunodeficiency virus type 1 (HIV-1) subtypes are CRF01_AE, CRF07_BC and CRF08_BC in Guangdong Province, China. In this study, we report a unique recombinant form (URF) of HIV-1 that was identified in an HIV/hepatitis B virus (HBV)/hepatitis C virus (HCV) triple-infected patient who was an intravenous drug user (IDU) in Heyuan City, Guangdong Province. The near full-length genome was amplified, and the PCR products were sequenced by Sanger's method. The Recombination Identification Program (RIP 3.0) and jpHMM online tools showed that four subtype C fragments were inserted into the A1 backbone genome in the gag, pol, vpr and nef gene regions. In the phylogenetic tree analysis, the subtype A1 and C fragments clustered with HIV-1 A1 and C reference sequences, respectively. No similar breakpoints between our strain and the other strains in the Los Alamos HIV database were observed. The results of evolutionary analysis using BEAST software showed that the subtype A1 fragment originated from Guangzhou City, China; however, the subtype C fragment originated from East Africa. This is the first report of HIV-1 URF A1C in Guangdong Province, China. The identification of this URF suggested that further dynamic surveillance of new recombinant forms is needed.

Molecular Epidemiology of HIV-1 Virus in Puerto Rico: Novel Cases of HIV-1 Subtype C, D, and CRF-24BG

HIV-1 subtype B virus is the most prevalent subtype in Puerto Rico (PR), accounting for about 90% of infection in the island. Recently, other subtypes and circulating recombinant forms (CRFs), including F(12_BF), A (01_BF), and CRF-39 BF-like, have been identified. The purpose of this study is to assess the distribution of drug resistance mutations and subtypes in PR. A total of 846 nucleotide sequences from the period comprising 2013 through 2017 were obtained from our "HIV Genotyping" test file. Phylogenetic and molecular epidemiology analyses were performed to evaluate the evolutionary dynamics and prevalence of drug resistance mutations. According to our results, we detected a decrease in the prevalence of protease inhibitor, nucleoside reverse transcriptase inhibitor (NRTI), and non-NRTI (NNRTI) resistance mutations over time. In addition, we also detected recombinant forms and, for the first time, identified subtypes C, D, and CRF-24BG in PR. Recent studies suggest that non-subtypes B are associated with a high risk of treatment failure and disease progression. The constant monitoring of viral evolution and drug resistance mutation dynamics is important to establish appropriate efforts for controlling viral expansion.

Genetic characterization and antiretroviral resistance mutations among treatment-naive HIV-infected individuals in Jiaxing, China

The aim of this study was to characterize HIV-1 genotypes and antiretroviral resistance mutations among treatment-naive HIV-infected individuals in Jiaxing, China. The HIV-1 partial polymerase (pol) genes in 93 of the 99 plasma samples were successfully amplified and analyzed. Phylogenetic analysis revealed the existence of five HIV-1 genotypes, of which the most prevalent genotype was CRF01_AE (38.7%), followed by CRF07_BC (34.4%), CRF08_BC (16.1%), subtype B/B' (5.4%), and CRF55_01B (2.1%). Besides, three types of unique recombination forms (URFs) were also observed, including C/F2/A1, CRF01_AE/B, and CRF08_BC/CRF07_BC. Among 93 amplicons, 46.2% had drug resistance-associated mutations, including 23.7% for protease inhibitors (PIs) mutations, 1.1% for nucleoside reverse transcriptase inhibitors (NRTIs) mutations, and 20.4% for non-nucleoside reverse transcriptase inhibitors (NNRTIs) mutations. Six (6.5%) out of 93 treatment-naive subjects were identified to be resistant to one or more NNRTIs, while resistance to NRTIs or PIs was not observed. Our study showed the genetic diversity of HIV-1 strains circulating in Jiaxing and a relative high proportion of antiretroviral resistance mutations among treatment-naive patients, indicating a serious challenge for HIV prevention and treatment program.

Decoding noises in HIV computational genotyping

Lack of a consistent and reliable genotyping system can critically impede HIV genomic research on pathogenesis, fitness, virulence, drug resistance, and genomic-based healthcare and treatment. At present, mis-genotyping, i.e., background noises in molecular genotyping, and its impact on epidemic surveillance is unknown. For the first time, we present a comprehensive assessment of HIV genotyping quality. HIV sequence data were retrieved from worldwide published records, and subjected to a systematic genotyping assessment pipeline. Results showed that mis-genotyped cases occurred at 4.6% globally, with some regional and high-risk population heterogeneities. Results also revealed a consistent mis-genotyping pattern in gp120 in all studied populations except the group of men who have sex with men. Our study also suggests novel virus diversities in the mis-genotyped cases. Finally, this study reemphasizes the importance of implementing a standardized genotyping pipeline to avoid genotyping disparity and to advance our understanding of virus evolution in various epidemiological settings.

Characteristics of HIV-infected U.S. Army soldiers linked in molecular transmission clusters, 2001-2012

Objective

Recent surveillance data suggests the United States (U.S.) Army HIV epidemic is concentrated among men who have sex with men. To identify potential targets for HIV prevention strategies, the relationship between demographic and clinical factors and membership within transmission clusters based on baseline pol sequences of HIV-infected Soldiers from 2001 through 2012 were analyzed.

Methods

We conducted a retrospective analysis of baseline partial pol sequences, demographic and clinical characteristics available for all Soldiers in active service and newly-diagnosed with HIV-1 infection from January 1, 2001 through December 31, 2012. HIV-1 subtype designations and transmission clusters were identified from phylogenetic analysis of sequences. Univariate and multivariate logistic regression models were used to evaluate and adjust for the association between characteristics and cluster membership.

Results

Among 518 of 995 HIV-infected Soldiers with available partial pol sequences, 29% were members of a transmission cluster. Assignment to a southern U.S. region at diagnosis and year of diagnosis were independently associated with cluster membership after adjustment for other significant characteristics (p<0.10) of age, race, year of diagnosis, region of duty assignment, sexually transmitted infections, last negative HIV test, antiretroviral therapy, and transmitted drug resistance. Subtyping of the pol fragment indicated HIV-1 subtype B infection predominated (94%) among HIV-infected Soldiers.

Conclusion

These findings identify areas to explore as HIV prevention targets in the U.S. Army. An increased frequency of current force testing may be justified, especially among Soldiers assigned to duty in installations with high local HIV prevalence such as southern U.S. states.

Enhanced surveillance of HIV-1 drug resistance in recently infected MSM in the UK

OBJECTIVES

To determine the prevalence of inferred low-frequency HIV-1 transmitted drug resistance (TDR) in MSM in the UK and its predicted effect on first-line therapy.

METHODS

The HIV-1 pol gene was amplified from 442 newly diagnosed MSM identified as likely recently infected by serological avidity testing in 2011-13. The PCR products were sequenced by next-generation sequencing with a mutation frequency threshold of >2% and TDR mutations defined according to the 2009 WHO surveillance drug resistance mutations list.

RESULTS

The majority (75.6%) were infected with subtype B and 6.6% with rare complex or unique recombinant forms. At a mutation frequency threshold of >20%, 7.2% (95% CI 5.0%-10.1%) of the sequences had TDR and this doubled to 15.8% (95% CI 12.6%-19.6%) at >2% mutation frequency (P < 0.0001). The majority (26/42, 62%) of low-frequency variants were against PIs. The most common mutations detected at >20% and 2%-20% mutation frequency differed for each drug class, these respectively being: L90M (n = 7) and M46IL (n = 10) for PIs; T215rev (n = 9) and D67GN (n = 4) for NRTIs; and K103N (n = 5) and G190E (n = 2) for NNRTIs. Combined TDR was more frequent in subtype B than non-B (OR = 0.38; 95% CI = 0.17-0.88; P = 0.024) and had minimal predicted effect on recommended first-line therapies.

CONCLUSIONS

The data suggest differences in the types of low-frequency compared with majority TDR variants that require a better understanding of the origins and clinical significance of low-frequency variants. This will better inform diagnostic and treatment strategies.

Rare HIV-1 Subtype J Genomes and a New H/U/CRF02_AG Recombinant Genome Suggests an Ancient Origin of HIV-1 in Angola

Angola has an extremely diverse HIV-1 epidemic fueled in part by the frequent interchange of people with the Democratic Republic of Congo (DRC) and Republic of Congo (RC). Characterization of HIV-1 strains circulating in Angola should help to better understand the origin of HIV-1 subtypes and recombinant forms and their transmission dynamics. In this study we characterize the first near full-length HIV-1 genomic sequences from HIV-1 infected individuals from Angola. Samples were obtained in 1993 from three HIV-1 infected patients living in Cabinda, Angola. Near full-length genomic sequences were obtained from virus isolates. Maximum likelihood phylogenetic tree inference and analyses of potential recombination patterns were performed to evaluate the sequence classifications and origins. Phylogenetic and recombination analyses revealed that one virus was a pure subtype J, another mostly subtype J with a small uncertain region, and the final virus was classified as a H/U/CRF02_AG recombinant. Consistent with their epidemiological data, the subtype J sequences were more closely related to each other than to other J sequences previously published. Based on the env gene, taxa from Angola occur throughout the global subtype J phylogeny. HIV-1 subtypes J and H are present in Angola at low levels since at least 1993. Low transmission efficiency and/or high recombination potential may explain their limited epidemic success in Angola and worldwide. The high diversity of rare subtypes in Angola suggests that Angola was part of the early establishment of the HIV-1 pandemic.

Los Alamos Hepatitis C Virus Sequence and Human Immunology Databases: An Expanding Resource for Antiviral Research

The hepatitis C virus (HCV) resource at Los Alamos (hcv.lanl.gov) provides access to multiple databases: one containing annotated sequences and the other a repository of immunogenic epitopes. They are derived from databases originally developed for HIV research (hiv.lanl.gov). HCV and HIV are RNA viruses with relatively compact genomes (around 10 kb) that are extraordinarily variable, both within and between hosts. This diversity requires methods to track and exclude variants from an individual infection or from epidemiologically related infections, and tools to analyse the variation. The HCV immunology database contains a curated inventory of immunogenic epitopes and information about their interaction with the host immune system, with associated retrieval and analysis tools. This interactive resource provides flexible retrieval tools for sequences, epitopes, clinical information, and metadata, as well as utilities for scientific data analysis, to investigators with internet access and a web browser. This paper describes the types of data and the services that these databases offer, the tools they provide, and their configuration and use. Examples of applications to clonal analysis for drug-resistance mutations are shown.

Implementation of an HIV-1 Triple-Target NAT Assay in the Routine Screening at Three German Red Cross Blood Centres

BACKGROUND

Blood product safety was significantly improved by the introduction of NAT testing in the late 1990s, resulting in a strong decrease of transfusion-transmitted infections (TTIs). Due to the occurrence of HIV-1 NAT test failures as a consequence of mismatch mutations in the amplicon regions of mono-target NAT assays, the Paul Ehrlich Institute mandated the implementation of multi-target NAT assays for HIV-1 in 2014. Commercial suppliers mostly developed dual-target NAT assays, with only one implementing a triple-target NAT assay.

METHODS

The HIV-1 triple-target NAT assay v3 (GFE Blut) was tested on mutated specimens and synthetic DNA bearing mutations that resulted in sample underquantification or false-negative test results. In addition, data from 2 years routine testing at three German Red Cross Blood centres were analysed.

RESULTS

The HIV-1 triple-target PCR could compensate for all mutations tested and could compensate the loss of one amplicon without a significant loss of sensitivity. Data from 2 years routine testing showed a solid performance.

CONCLUSION

The HIV-1 triple-target v3 assay (GFE Blut) can compensate mutations in target sequences better than a dual-target assay and is applicable to high-throughput screening, thus increasing blood product safety.

Geographic origin and evolutionary history of China's two predominant HIV-1 circulating recombinant forms, CRF07_BC and CRF08_BC

To determine the origin and evolutionary history of two predominant and closely-related circulating recombinant forms (CRF07_BC and CRF08_BC), recombinant structures and phylogenies of 7 unique recombinant forms comprised of subtypes of B’ (Thai B linage) and C (designated URFs_BC) from archival specimens of injection drug users (IDUs) collected in 1996 to 1998 from western Yunnan and 4 circulating recombinant forms with B’/C recombinants recently identified (designated nCRFs_BC) in China were compared with those of CRF07_BC and CRF08_BC. The results showed that 5 of 7 URFs_BC and all the nCRFs_BC shared recombination breakpoints with CRF07_BC and/or CRF08_BC. Yunnan URFs_BC consistently occupied the basal branch positions compared with CRF07_BC, CRF08_BC, and nCRFs_BC in phylogenetic trees. The estimated most recent common ancestors (tMRCA) for Yunnan URFs_BC were from ~1987, approximately half a decade earlier than those for CRF07_BC (~1994) and CRF08_BC (~1992). Discrete phylogeographic and spatial diffusion analysis revealed that both CRF07_BC and CRF08 BC came from western Yunnan in the early 1990s. Our results provide compelling evidence for western Yunnan as the geographic origin of CRF07_BC and CRF08_BC, which emerged from a swarm of URFs_BC by a series of recombination events in western Yunnan in the early 1990s.

HIV-1 Thai B strain has spread out of former plasma donors into general population through sexual contact in Henan, China

Henan, China is characterized by the outbreak of HIV epidemic of Thai B strain in former plasma donors in 1990s. After the forbidden of paid blood donation, whether Thai B strain will spread out of former plasma donors into sexual transmitted population is unknown. To answer the question, phylogenetic analysis was used to explore relationships of HIV strains circulating in those two populations in the study. HIV-1 sero-positive drug-naïve patients infected through sexual contact were enrolled into the study. Full length gag and pol genes were amplified with nested RT-PCR followed by sequencing and phylogenetic analysis. The genotypes of anti-HIV drug resistance were also analyzed with available pol genes. HIV subtypes were determined in 249 individuals from 288 participants. Subtype B was dominant (202/249, 81.1%), followed by CRF01_AE (25/249, 10.0%), CRF07_BC (14/249, 5.6%), C (4/249, 1.6%), URF (3/249, 1.2%), and CRF08_BC (1/249, 0.4%). Most of subtype B strains belong to Thailand B lineage. All of Thai B strains identified in sexual transmitted population intermixed with those from former blood donors in phylogenetic tree, suggesting close phylogenetic relationship between strains epidemic in those two populations. TDR was identified in 9.9% individuals. Thai B strain has spread out of former blood donors in Henan province. The finding will contribute to understanding the distribution and evolution of HIV-1 in Henan province and also provide clue to behavior change intervention.

The HAND Database: a gateway to understanding the role of HIV in HIV-associated neurocognitive disorders

Background

Despite an augmented research effort and scale-up of highly active antiretroviral therapy, a high prevalence of HIV-1-associated neurocognitive disorders (HAND) persists in the HIV-infected population. Nearly 50 % of all HIV-1-infected individuals suffer from a neurocognitive disorder due to neural and synaptodendritic damage. Challenges in HAND research, including limited availability of brain tissue from HIV patients, variation in HAND study protocols, and virus genotyping inconsistency and errors, however, have resulted in studies with insufficient power to delineate molecular mechanisms underlying HAND pathogenesis. There exists, therefore, a great need for a reliable and centralized resource specific to HAND research, particularly for epidemiological study and surveillance in resource-limited countries where severe forms of HAND persist.

Description

To address the aforementioned imperative need, here we present the HAND Database, a resource containing well-curated and up-to-date HAND virus information and associated clinical and epidemiological data. This database provides information on 5,783 non-redundant HIV-1 sequences from global HAND research published to date, representing a total of 163 unique individuals that have been assessed for HAND. A user-friendly interface allows for flexible searching, filtering, browsing, and downloading of data. The most comprehensive database of its kind, the HAND Database not only bolsters current HAND research by increasing sampling power and reducing study biases caused by protocol variation and genotyping inconsistency, it allows for comparison between HAND studies across different dimensions. Development of the HAND Database has also revealed significant knowledge gaps in HIV-driven neuropathology. These gaps include inadequate sequencing of viral genes beyond env, lack of HAND viral data from HIV epidemiologically important regions including Asian and Sub-Saharan African countries, and biased sampling toward the male gender, all factors that impede efforts toward providing an improved quality of life to HIV-infected individuals, and toward elimination of viruses in the brain.

Conclusion

Our aim with the HAND database is to provide researchers in both the HIV and neuroscience fields a comprehensive and rigorous data source toward better understanding virus compartmentalization and to help in design of improved strategies against HAND viruses. We also expect this resource, which will be updated on a regular basis, to be useful as a reliable reference for further HAND epidemiology studies. The HAND Database is freely available and accessible online at http://www.handdatabase.org.

Genetic Characterization of a Unique Recombinant Originating from CRF55_01B, CRF01_AE, and CRF07_BC in Shenzhen, China

Multiple subtypes were found to be epidemic in the Shenzhen men who have sex with men (MSM) population, which always predicts the emergence of a unique recombinant. In 2012, CRF55_01B was first reported, which later was proven to have originated in MSM in Shenzhen city. In this study, we reported a unique recombinant form (URF) of HIV-1 identified in a man who has had sex with men in Shenzhen city. The strain showed a genomic schematic map similar to CRF55_01B with subtype C segments inserted in the gag and pol genes. The full-length genome was amplified in two halves with 1-kb overlap regions. The PCR products were cloned and sequenced. A recombination detection program showed that two subtype C fragments and two subtype B fragments were inserted into the CRF01_AE backbone genome in the gag and pol regions. In the phylogenetic tree, the subtype C fragments clustered with CRF07_BC variants and the other segments grouped with CRF55_01B strains except for one segment that clustered with CRF01_AE. Similar breakpoints between our strain and CRF65_cpx were also observed. The data suggested that the URF strain might be the recombinant form of CRF55_01B, CRF01_AE, and CRF07_BC. This is the first report of a third generation of recombination of HIV-1 that originated from CRF55_01B in China. The identification of the URF indicated the severity of the HIV epidemic in Shenzhen MSM and the urgent need for epidemiological surveillance of the new recombination.

Sequencing and phylogenetic analysis of near full-length HIV-1 subtypes A, B, G and unique recombinant AC and AD viral strains identified in South Africa

By the end of 2012, more than 6.1 million people were infected with HIV-1 in South Africa. Subtype C was responsible for the majority of these infections and more than 300 near full-length genomes (NFLGs) have been published. Currently very few non-subtype C isolates have been identified and characterized within the country, particularly full genome non-C isolates. Seven patients from the Tygerberg Virology (TV) cohort were previously identified as possible non-C subtypes and were selected for further analyses. RNA was isolated from five individuals (TV047, TV096, TV101, TV218, and TV546) and DNA from TV016 and TV1057. The NFLGs of these samples were amplified in overlapping fragments and sequenced. Online subtyping tools REGA version 3 and jpHMM were used to screen for subtypes and recombinants. Maximum likelihood (ML) phylogenetic analysis (phyML) was used to infer subtypes and SimPlot was used to confirm possible intersubtype recombinants. We identified three subtype B (TV016, TV047, and TV1057) isolates, one subtype A1 (TV096), one subtype G (TV546), one unique AD (TV101), and one unique AC (TV218) recombinant form. This is the first NFLG of subtype G that has been described in South Africa. The subtype B sequences described also increased the NFLG subtype B sequences in Africa from three to six. There is a need for more NFLG sequences, as partial HIV-1 sequences may underrepresent viral recombinant forms. It is also necessary to continue monitoring the evolution and spread of HIV-1 in South Africa, because understanding viral diversity may play an important role in HIV-1 prevention strategies.

Molecular epidemiology of HIV-1 in Jilin province, northeastern China: emergence of a new CRF07_BC transmission cluster and intersubtype recombinants

Objective

To investigate the HIV-1 molecular epidemiology among newly diagnosed HIV-1 infected persons living in the Jilin province of northeastern China.

Methods

Plasma samples from 189 newly diagnosed HIV-1 infected patients were collected between June 2010 and August 2011 from all nine cities of Jilin province. HIV-1 nucleotide sequences of gag P17–P24 and env C2–C4 gene regions were amplified using a multiplex RT-PCR method and sequenced. Phylogenetic and recombination analyses were used to determine the HIV-1 genotypes.

Results

Based on all sequences generated, the subtype/CFR distribution was as follows: CRF01_AE (58.1%), CRF07_BC (13.2%), subtype B’ (13.2%), recombinant viruses (8.1%), subtype B (3.7%), CRF02_AG (2.9%), subtype C (0.7%). In addition to finding CRF01_AE strains from previously reported transmission clusters 1, 4 and 5, a new transmission cluster was described within the CRF07_BC radiation. Among 11 different recombinants identified, 10 contained portions of gene regions from the CRF01_AE lineage. CRF02_AG was found to form a transmission cluster of 4 in local Jilin residents.

Conclusions

Our study presents a molecular epidemiologic investigation describing the complex structure of HIV-1 strains co-circulating in Jilin province. The results highlight the critical importance of continuous monitoring of HIV-infections, along with detailed socio-demographic data, in order to design appropriate prevention measures to limit the spread of new HIV infections.

A model-based information sharing protocol for profile Hidden Markov Models used for HIV-1 recombination detection

BACKGROUND

In many applications, a family of nucleotide or protein sequences classified into several subfamilies has to be modeled. Profile Hidden Markov Models (pHMMs) are widely used for this task, modeling each subfamily separately by one pHMM. However, a major drawback of this approach is the difficulty of dealing with subfamilies composed of very few sequences. One of the most crucial bioinformatical tasks affected by the problem of small-size subfamilies is the subtyping of human immunodeficiency virus type 1 (HIV-1) sequences, i.e., HIV-1 subtypes for which only a small number of sequences is known.

RESULTS

To deal with small samples for particular subfamilies of HIV-1, we introduce a novel model-based information sharing protocol. It estimates the emission probabilities of the pHMM modeling a particular subfamily not only based on the nucleotide frequencies of the respective subfamily but also incorporating the nucleotide frequencies of all available subfamilies. To this end, the underlying probabilistic model mimics the pattern of commonality and variation between the subtypes with regards to the biological characteristics of HI viruses. In order to implement the proposed protocol, we make use of an existing HMM architecture and its associated inference engine.

CONCLUSIONS

We apply the modified algorithm to classify HIV-1 sequence data in the form of partial HIV-1 sequences and semi-artificial recombinants. Thereby, we demonstrate that the performance of pHMMs can be significantly improved by the proposed technique. Moreover, we show that our algorithm performs significantly better than Simplot and Bootscanning.

Near full-length genome sequence of a novel HIV type 1 second-generation recombinant form (CRF01_AE/CRF07_BC) identified among men who have sex with men in Jilin, China

We report here a novel HIV-1 second-generation recombinant form (CRF01_AE/CRF07_BC) composed of CRF01_AE and CRF07_BC, identified among men who have sex with men (MSM) in Jilin, with four breakpoints observed in the pol, vif, and vpr genes. The CRF01_AE regions of the recombinant were clustered with the CRF01_AE lineage, which is mainly circulating among MSM in northern China, with the support of 100% bootstrap value, indicating that the parental origin of the CRF01_AE regions was from MSM, in which recombination events may be more likely to occur. To the best of our knowledge, this is the first detection of a novel HIV-1 second-generation recombinant form (CRF01AE/CRF07_BC) in Jilin, which indicates active transmission networks of HIV-1 infection among MSM in the region. Therefore, it is necessary to continue monitoring the molecular epidemiology of HIV-1 among MSM in Jilin to obtain a better understanding of the transmission and potential public health impact of HIV-1 among MSM in the region.

Bayesian inference of the evolution of HBV/E

Despite its wide spread and high prevalence in sub-Saharan Africa, hepatitis B virus genotype E (HBV/E) has a surprisingly low genetic diversity, indicating an only recent emergence of this genotype in the general African population. Here, we performed extensive phylogeographic analyses, including Bayesian MCMC modeling. Our results indicate a mutation rate of 1.9×10⁻⁴ substitutions per site and year (s/s/y) and confirm a recent emergence of HBV/E, most likely within the last 130 years, and only after the transatlantic slave-trade had come to an end. Our analyses suggest that HBV/E originated from the area of Nigeria, before rapidly spreading throughout sub-Saharan Africa. Interestingly, viral strains found in Haiti seem to be the result of multiple introductions only in the second half of the 20th century, corroborating an absence of a significant number of HBV/E strains in West Africa several centuries ago. Our results confirm that the hyperendemicity of HBV(E) in today's Africa is a recent phenomenon and likely the result of dramatic changes in the routes of viral transmission in a relatively recent past.

Molecular characterization of unique intersubtype HIV type 1 A1/C recombinant strain circulating in Pune, India

An increasing number of circulating recombinant forms (CRFs) and unique recombinant forms (URFs) all over the world has necessitated being vigilant about new recombinants. Since the first report of a recombinant virus with an A1/C mosaic in 1998 more and more B/C and A/C recombinant viruses are being reported from India. Here we report the identification and characterization of a unique HIV-1 A1/C recombinant circulating in Western India. Analysis of the full-length genome using RIP, SimPlot, and jpHMM@Gobics has confirmed its mosaic structure with insertion of subtype A1 in the backbone of subtype C at three positions: gag-pol (1973±15-2617±47), pol-vif (4879±37-5582±32), and gp41 (8437±106-8811±8); however, RIP and SimPlot showed one more small insertion in integrase (4343-4519). Phylogenetic analysis confirmed that the recombinant virus has an insertion of clade A1 in the backbone of subtype C, which has come from Indian subtype C.

Recombination form and epidemiology of HIV-1 unique recombinant strains identified in Yunnan, China

Several studies identified HIV-1 recombination in some distinct areas in Yunnan, China. However, no comprehensive studies had been fulfilled in the whole province up to now. To illustrate the epidemiology and recombination form of Unique Recombinant Forms (URFs) circulating in Yunnan, 788 HIV-1 positive individuals residing in 15 prefectures of Yunnan were randomly enrolled into the study. Full-length gag and pol genes were amplified and sequenced. Maximum likelihood tree was constructed for phylogenetic analysis. Recombinant breakpoints and genomic schematics were identified with online software jpHMM. 63 (10.2%) unique recombinant strains were identified from 617 strains with subtypes. The URFs distributed significantly differently among prefectures (Pearson chi-square test, P<0.05). IDUs contained more URFs than sexual transmitted population (Pearson chi-square test, P<0.05). Two main recombinant forms were identified by considering the presence of CRF01_AE segments in full length gag-pol genes, which were B′/C and B′/C/CRF01-AE recombinants. Three clusters were identified in the ML tree which contained more than three sequences and supported by high bootstrap values. One CRF was identified. Many of URFs contained identical breakpoints. The results will contribute to our understanding on HIV recombination and provide clues to the identification of potential CRFs in China.

Characterization of primary isolates of HIV type 1 CRF28_BF, CRF29_BF, and unique BF recombinants circulating in São Paulo, Brazil

We report for the first time the genetic and biological characterization of 10 HIV-1 primary isolates representing CRF28_BF and CRF29_BF together with additional unique BF recombinant forms (URFs) obtained by PBMC cocultivation. Recombination is an important factor promoting the increase in the genetic diversity of HIV-1. Notably, more than 20% of HIV-1 sequences worldwide were recombinants. Several recombinant viruses were reported in Brazil, and six circulating recombinant forms (CRFs) have been identified (CRF28_BF, CRF29_BF, CRF31_BC, CRF39_BF, CRF40_BF, and CRF46_BF). CRF28_BF and CRF29_BF were found to infect almost 30% of the patients in São Paulo State. The near full-length genomes of these 10 primary isolates were amplified by nested PCR in three overlapping segments, purified, and sequenced. Three samples were related to CRF28_BF, three to CRF29_BF, and four were unique recombinant forms (URFs), as determined by their breakpoint profile determined with the jpHMM program. Additionally, the coreceptor usage of these isolates was investigated in vitro using GHOST assays, which revealed three dual-tropic (X4/R5) viruses, four lymphotropic (X4) viruses, and three macrophage-tropic (R5) viruses with different V3-loop motifs, which challenges the notion that GWGR-carrying viruses are macrophage-tropic only. In sum, we report a much-anticipated well-characterized panel of viruses representing CRF28_BF, CRF29_BF, and URFs from São Paulo State, Brazil.

A clustering phenomenon among HCV-1a strains among patients coinfected with HIV from Buenos Aires, Argentina

The human immunodeficiency virus (HIV) and hepatitis C virus (HCV) share the same transmission routes which lead to high coinfection rates. Among HIV-infected individuals such rates reached 21% in Argentina, being HCV-1a the most predominant subtype. In this work, 25 HCV subtype 1a (HCV-1a) strains from Argentinean patients coinfected with HIV were studied based on E2 and NS5A sequences. Phylogenetic analyses indicated that 12 strains were highly related to each other, constituting a highly supported (posterior probability = 0.95) monophyletic group that we called "M." The remaining HCV strains (group dispersed or "D") were interspersed along the phylogenetic trees. When comparing both groups of HCV-1a, 10 amino acid differences were located in functional domains of E2 and NS5A proteins that appeared to affect eventually the peptides binding to MHC-I molecules thus favoring immune escape and contributing to the divergence of HCV genotypes. Bayesian coalescent analyses for HCV-1a cluster M isolates indicated that the time to the most recent common ancestor (tMRCA) overlaps with the age estimated recently for the HIV-BF epidemic in Argentina. Furthermore, the genomic characterization based on pol gene analysis from HIV viremic patients showed that most HIV isolates from patients coinfected with HCV-1a cluster M were BF recombinants with identical recombination patterns. In conclusion, these results suggest the presence of an HCV-1a monophyletic cluster with a potential HIV co-transmission by phylogenetic analyses.

jpHMM: recombination analysis in viruses with circular genomes such as the hepatitis B virus

jpHMM is a very accurate and widely used tool for recombination detection in genomic sequences of HIV-1. Here, we present an extension of jpHMM to analyze recombinations in viruses with circular genomes such as the hepatitis B virus (HBV). Sequence analysis of circular genomes is usually performed on linearized sequences using linear models. Since linear models are unable to model dependencies between nucleotides at the 5'- and 3'-end of a sequence, this can result in inaccurate predictions of recombination breakpoints and thus in incorrect classification of viruses with circular genomes. The proposed circular jpHMM takes into account the circularity of the genome and is not biased against recombination breakpoints close to the 5'- or 3'-end of the linearized version of the circular genome. It can be applied automatically to any query sequence without assuming a specific origin for the sequence coordinates. We apply the method to genomic sequences of HBV and visualize its output in a circular form. jpHMM is available online at http://jphmm.gobics.de for download and as a web server for HIV-1 and HBV sequences.

BF integrase genes of HIV-1 circulating in São Paulo, Brazil, with a recurrent recombination region

Although some studies have shown diversity in HIV integrase (IN) genes, none has focused particularly on the gene evolving in epidemics in the context of recombination. The IN gene in 157 HIV-1 integrase inhibitor-naïve patients from the São Paulo State, Brazil, were sequenced tallying 128 of subtype B (23 of which were found in non-B genomes), 17 of subtype F (8 of which were found in recombinant genomes), 11 integrases were BF recombinants, and 1 from subtype C. Crucially, we found that 4 BF recombinant viruses shared a recurrent recombination breakpoint region between positions 4900 and 4924 (relative to the HXB2) that includes 2 gRNA loops, where the RT may stutter. Since these recombinants had independent phylogenetic origin, we argue that these results suggest a possible recombination hotspot not observed so far in BF CRF in particular, or in any other HIV-1 CRF in general. Additionally, 40% of the drug-naïve and 45% of the drug-treated patients had at least 1 raltegravir (RAL) or elvitegravir (EVG) resistance-associated amino acid change, but no major resistance mutations were found, in line with other studies. Importantly, V151I was the most common minor resistance mutation among B, F, and BF IN genes. Most codon sites of the IN genes had higher rates of synonymous substitutions (dS) indicative of a strong negative selection. Nevertheless, several codon sites mainly in the subtype B were found under positive selection. Consequently, we observed a higher genetic diversity in the B portions of the mosaics, possibly due to the more recent introduction of subtype F on top of an ongoing subtype B epidemics and a fast spread of subtype F alleles among the B population.

Detection of viral sequence fragments of HIV-1 subfamilies yet unknown

Background

Methods of determining whether or not any particular HIV-1 sequence stems - completely or in part - from some unknown HIV-1 subtype are important for the design of vaccines and molecular detection systems, as well as for epidemiological monitoring. Nevertheless, a single algorithm only, the Branching Index (BI), has been developed for this task so far. Moving along the genome of a query sequence in a sliding window, the BI computes a ratio quantifying how closely the query sequence clusters with a subtype clade. In its current version, however, the BI does not provide predicted boundaries of unknown fragments.

Results

We have developed Unknown Subtype Finder (USF), an algorithm based on a probabilistic model, which automatically determines which parts of an input sequence originate from a subtype yet unknown. The underlying model is based on a simple profile hidden Markov model (pHMM) for each known subtype and an additional pHMM for an unknown subtype. The emission probabilities of the latter are estimated using the emission frequencies of the known subtypes by means of a (position-wise) probabilistic model for the emergence of new subtypes. We have applied USF to SIV and HIV-1 sequences formerly classified as having emerged from an unknown subtype. Moreover, we have evaluated its performance on artificial HIV-1 recombinants and non-recombinant HIV-1 sequences. The results have been compared with the corresponding results of the BI.

Conclusions

Our results demonstrate that USF is suitable for detecting segments in HIV-1 sequences stemming from yet unknown subtypes. Comparing USF with the BI shows that our algorithm performs as good as the BI or better.

Molecular characterization of HIV type 1 among HIV-infected respondents in a cohort being prepared for HIV Phase III vaccine clinical trials, Western Kenya

Kenya is one of the sub-Saharan African countries affected by HIV-1 infection and AIDS. We investigated HIV-1 genetic diversity in 130 individuals from Busia, Bungoma, and Kakamega in western Kenya as part of an HIV-1 vaccine feasibility study in preparation for Phase III efficacy clinical trials. After RNA extraction the partial gag (484 bp) and env (1297 bp) regions were amplified and directly sequenced. Phylogenetic analysis was done using MEGA version 4 and recombinants were identified using the jpHMM tool and phylogenetic analysis. HIV-1 sequences were amplified from 122 of the 130 samples, 118 (90.8%) from the gag region and 78 (60 %) from the env region and 74 samples (56.9%) from both the gag and env regions. Of these sequenced on both regions, 51.4% were subtype A, 9.4% subtype D, 1.4% subtype C, 4.1% subtype G, and 33.7% were discordant and thus possible recombinants, including A1/C, A1/D, A1/A2, and A2/C. The jpHMM tool indicated a further two samples with CD and BD breakpoints within the env gene and one within the gag gene (A1C). An additional sample had an A1D breakpoint in the gag gene, but the envelope was not amplified. HIV-1 subtype diversity in western Kenya should be considered in vaccines designed for clinical trials in this region and this genetic diversity should be continuously monitored.

A classification approach for genotyping viral sequences based on multidimensional scaling and linear discriminant analysis

BACKGROUND

Accurate classification into genotypes is critical in understanding evolution of divergent viruses. Here we report a new approach, MuLDAS, which classifies a query sequence based on the statistical genotype models learned from the known sequences. Thus, MuLDAS utilizes full spectra of well characterized sequences as references, typically of an order of hundreds, in order to estimate the significance of each genotype assignment.

RESULTS

MuLDAS starts by aligning the query sequence to the reference multiple sequence alignment and calculating the subsequent distance matrix among the sequences. They are then mapped to a principal coordinate space by multidimensional scaling, and the coordinates of the reference sequences are used as features in developing linear discriminant models that partition the space by genotype. The genotype of the query is then given as the maximum a posteriori estimate. MuLDAS tests the model confidence by leave-one-out cross-validation and also provides some heuristics for the detection of 'outlier' sequences that fall far outside or in-between genotype clusters. We have tested our method by classifying HIV-1 and HCV nucleotide sequences downloaded from NCBI GenBank, achieving the overall concordance rates of 99.3% and 96.6%, respectively, with the benchmark test dataset retrieved from the respective databases of Los Alamos National Laboratory.

CONCLUSIONS

The highly accurate genotype assignment coupled with several measures for evaluating the results makes MuLDAS useful in analyzing the sequences of rapidly evolving viruses such as HIV-1 and HCV. A web-based genotype prediction server is available at http://www.muldas.org/MuLDAS/.

HIV classification using the coalescent theory

MOTIVATION

Existing coalescent models and phylogenetic tools based on them are not designed for studying the genealogy of sequences like those of HIV, since in HIV recombinants with multiple cross-over points between the parental strains frequently arise. Hence, ambiguous cases in the classification of HIV sequences into subtypes and circulating recombinant forms (CRFs) have been treated with ad hoc methods in lack of tools based on a comprehensive coalescent model accounting for complex recombination patterns.

RESULTS

We developed the program ARGUS that scores classifications of sequences into subtypes and recombinant forms. It reconstructs ancestral recombination graphs (ARGs) that reflect the genealogy of the input sequences given a classification hypothesis. An ARG with maximal probability is approximated using a Markov chain Monte Carlo approach. ARGUS was able to distinguish the correct classification with a low error rate from plausible alternative classifications in simulation studies with realistic parameters. We applied our algorithm to decide between two recently debated alternatives in the classification of CRF02 of HIV-1 and find that CRF02 is indeed a recombinant of Subtypes A and G.

AVAILABILITY

ARGUS is implemented in C++ and the source code is available at http://gobics.de/software.

The role of recombination in the emergence of a complex and dynamic HIV epidemic

Background

Inter-subtype recombinants dominate the HIV epidemics in three geographical regions. To better understand the role of HIV recombinants in shaping the current HIV epidemic, we here present the results of a large-scale subtyping analysis of 9435 HIV-1 sequences that involve subtypes A, B, C, G, F and the epidemiologically important recombinants derived from three continents.

Results

The circulating recombinant form CRF02_AG, common in West Central Africa, appears to result from recombination events that occurred early in the divergence between subtypes A and G, followed by additional recent recombination events that contribute to the breakpoint pattern defining the current recombinant lineage. This finding also corrects a recent claim that G is a recombinant and a descendant of CRF02, which was suggested to be a pure subtype. The BC and BF recombinants in China and South America, respectively, are derived from recent recombination between contemporary parental lineages. Shared breakpoints in South America BF recombinants indicate that the HIV-1 epidemics in Argentina and Brazil are not independent. Therefore, the contemporary HIV-1 epidemic has recombinant lineages of both ancient and more recent origins.

Conclusions

Taken together, we show that these recombinant lineages, which are highly prevalent in the current HIV epidemic, are a mixture of ancient and recent recombination. The HIV pandemic is moving towards having increasing complexity and higher prevalence of recombinant forms, sometimes existing as "families" of related forms. We find that the classification of some CRF designations need to be revised as a consequence of (1) an estimated > 5% error in the original subtype assignments deposited in the Los Alamos sequence database; (2) an increasing number of CRFs are defined while they do not readily fit into groupings for molecular epidemiology and vaccine design; and (3) a dynamic HIV epidemic context.

Use of massive parallel pyrosequencing for near full-length characterization of a unique HIV Type 1 BF recombinant associated with a fatal primary infection

Near full length genome characterization of a BF recombinant from a patient who died from multiorgan failure during HIV-1 seroconversion is reported. Massive parallel pyrosequencing was used with the shotgun approach. Intrahost genetic variability along the whole genome was calculated and coreceptor usage of viral quasispecies was predicted. A consensus sequence was established to perform subtype assignment, phylogenetic analysis, and recombination tests. The sequence clustered with two recently described BF unique recombinant forms from Brazil, consistent with the recombination pattern, yielding breakpoints located at the same positions, with the exception of the second env breakpoint. The actual prevalence of recombinant forms is probably underestimated if partial genomic regions are considered. Here the first full length BF recombinant from Italy is described, together with an evaluation of quasispecies heterogeneity. Our data provide evidence that next generation sequencing may provide a major contribution to HIV-1 molecular epidemiology and to the comprehension of intrapatient heterogeneity.

jpHMM: improving the reliability of recombination prediction in HIV-1

Previously, we developed jumping profile hidden Markov model (jpHMM), a new method to detect recombinations in HIV-1 genomes. The jpHMM predicts recombination breakpoints in a query sequence and assigns to each position of the sequence one of the major HIV-1 subtypes. Since incorrect subtype assignment or recombination prediction may lead to wrong conclusions in epidemiological or vaccine research, information about the reliability of the predicted parental subtypes and breakpoint positions is valuable. For this reason, we extended the output of jpHMM to include such information in terms of ‘uncertainty’ regions in the recombination prediction and an interval estimate of the breakpoint. Both types of information are computed based on the posterior probabilities of the subtypes at each query sequence position. Our results show that this extension strongly improves the reliability of the jpHMM recombination prediction. The jpHMM is available online at http://jphmm.gobics.de/.

Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning

Background

Recombination has a profound impact on the evolution of viruses, but characterizing recombination patterns in molecular sequences remains a challenging endeavor. Despite its importance in molecular evolutionary studies, identifying the sequences that exhibit such patterns has received comparatively less attention in the recombination detection framework. Here, we extend a quartet-mapping based recombination detection method to enable identification of recombinant sequences without prior specifications of either query and reference sequences. Through simulations we evaluate different recombinant identification statistics and significance tests. We compare the quartet approach with triplet-based methods that employ additional heuristic tests to identify parental and recombinant sequences.

Results

Analysis of phylogenetic simulations reveal that identifying the descendents of relatively old recombination events is a challenging task for all methods available, and that quartet scanning performs relatively well compared to the triplet based methods. The use of quartet scanning is further demonstrated by analyzing both well-established and putative HIV-1 recombinant strains. In agreement with recent findings, we provide evidence that the presumed circulating recombinant CRF02_AG is a 'pure' lineage, whereas the presumed parental lineage subtype G has a recombinant origin. We also demonstrate HIV-1 intrasubtype recombination, confirm the hybrid origin of SIV in chimpanzees and further disentangle the recombinant history of SIV lineages in a primate immunodeficiency virus data set.

Conclusion

Quartet scanning makes a valuable addition to triplet-based methods for identifying recombinant sequences without prior specifications of either query and reference sequences. The new method is available in the VisRD v.3.0 package .

Discrepancies in assignment of subtype/recombinant forms by genotyping programs for HIV type 1 drug resistance testing may falsely predict superinfection

With the growing diversity of the HIV pandemic, routine genotyping is an important tool for monitoring viral subtype as well as drug resistance. In this regard, numerous subtyping tools and drug resistance algorithms are available online. However, there are discrepancies in the use of these online tools in the designation of HIV-1 subtypes or recombinant forms that may have an impact on drug susceptibility profiles. Indeed, inconsistencies in some of these tools may lead to a false designation of dual infection and/or superinfection. In this case study, we evaluated the sequence diversity of an infection that was referred to us as a potential case of superinfection as a result of variations in designation of subtype. We evaluated sequences using five different online tools and finally determined by phylogenetic analysis that the sequence was a unique A1/C intersubtype recombinant at baseline and not a case of superinfection.

Near full-length genomic analysis of a novel subtype A1/C recombinant HIV type 1 isolate from Canada

The rapid introduction of non-B HIV-1 subtypes into Quebec, mostly from persons from regions where HIV prevalence is high and where different variants circulate, means that recombination must now be an important consideration in the epidemiologic surveillance of HIV infection. No circulating recombinant form (CRF), currently involving exclusively subtypes A1 and C, exists in the Los Alamos HIV database. This study presents a near full-length genomic analysis of a novel HIV-1 recombinant involving subtypes A1 and C. Bootscanning revealed that the recombinant structure involves three breakpoints that separate the genome into four regions, alternating between subtypes A1 and C. The intersubtype recombinant breakpoint in the pol gene was at midpoint between the protease and reverse transcriptase open reading frames. This is the first report of a recombinant involving subtypes A1 and C in Canada, the epidemiologic significance of which is not yet known.

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.