RecDraw: a software package for the representation of HIV-1 recombinant structures

Streamlined Subpopulation, Subtype, and Recombination Analysis of HIV-1 Half-Genome Sequences Generated by High-Throughput Sequencing

The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

High-throughput sequencing (HTS) has been widely used to characterize HIV-1 genome sequences. There are no algorithms currently that can directly determine genotype and quasispecies population using short HTS reads generated from long genome sequences without additional software. To establish a robust subpopulation, subtype, and recombination analysis workflow, we amplified the HIV-1 3′-half genome from plasma samples of 65 HIV-1-infected individuals and sequenced the entire amplicon (∼4,500 bp) by HTS. With direct analysis of raw reads using HIVE-hexahedron, we showed that 48% of samples harbored 2 to 13 subpopulations. We identified various subtypes (17 A1s, 4 Bs, 27 Cs, 6 CRF02_AGs, and 11 unique recombinant forms) and defined recombinant breakpoints of 10 recombinants. These results were validated with viral genome sequences generated by single genome sequencing (SGS) or the analysis of consensus sequence of the HTS reads. The HIVE-hexahedron workflow is more sensitive and accurate than just evaluating the consensus sequence and also more cost-effective than SGS.

IMPORTANCE The highly recombinogenic nature of human immunodeficiency virus type 1 (HIV-1) leads to recombination and emergence of quasispecies. It is important to reliably identify subpopulations to understand the complexity of a viral population for drug resistance surveillance and vaccine development. High-throughput sequencing (HTS) provides improved resolution over Sanger sequencing for the analysis of heterogeneous viral subpopulations. However, current methods of analysis of HTS reads are unable to fully address accurate population reconstruction. Hence, there is a dire need for a more sensitive, accurate, user-friendly, and cost-effective method to analyze viral quasispecies. For this purpose, we have improved the HIVE-hexahedron algorithm that we previously developed with in silico short sequences to analyze raw HTS short reads. The significance of this study is that our standalone algorithm enables a streamlined analysis of quasispecies, subtype, and recombination patterns from long HIV-1 genome regions without the need of additional sequence analysis tools. Distinct viral populations and recombination patterns identified by HIVE-hexahedron are further validated by comparison with sequences obtained by single genome sequencing (SGS).

New Subtype B Containing HIV-1 Circulating Recombinant of sub-Saharan Africa Origin in Nigerian Men Who Have Sex With Men

BACKGROUND

HIV-1 circulating recombinant forms (CRF) containing subtype B are uncommon in sub-Saharan Africa. Prevalent infections observed during enrollment of a prospective study of men who have sex with men (MSM) from Lagos, Nigeria, revealed the presence of a family of subtype B and CRF02_AG recombinants. This report describes the HIV-1 genetic diversity within a high-risk, high-prevalence, and previously undersampled cohort of Nigerian MSM.

METHODS

Between 2013 and 2016, 672 MSM were enrolled at the Lagos site of the TRUST/RV368 study. Prevalent HIV-1 infections were initially characterized by pol sequencing and phylogenetic subtyping analysis. Samples demonstrating the presence of subtype B were further characterized by near full-length sequencing, phylogenetic, and Bayesian analyses.

RESULTS

Within this cohort, HIV-1 prevalence was 59%. The major subtype was CRF02_AG (57%), followed by CRF02/B recombinants (15%), subtype G (13%), and smaller amounts of A1, B, and other recombinants. Nine clusters of closely related pol sequences indicate ongoing transmission events within this cohort. Among the CRF02_AG/B, a new CRF was identified and termed CRF95_02B. Shared risk factors and Bayesian phylogenetic inference of the new CRF95_02B and the similarly structured CRF56_cpx indicate a Nigerian or West African origin of CRF56_cpx before its observation in France.

CONCLUSION

With high HIV-1 prevalence, new strains, and multiple transmission networks, this cohort of Nigerian MSM represents a previously hidden reservoir of HIV-1 strains, including the newly identified CRF95_02B and closely related CRF56_cpx. These strains will need to be considered during vaccine selection and development to optimize the design of a globally effective HIV-1 vaccine.

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.

Increased predominance of HIV-1 CRF01_AE and its recombinants in the Philippines

The growth rate of new HIV infections in the Philippines was the fastest of any countries in the Asia-Pacific region between 2010 and 2016. To date, HIV-1 subtyping results in the Philippines have been determined by characterizing only partial viral genome sequences. It is not known whether recombination occurs in the majority of unsequenced genome regions. Near-full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes from plasma samples collected between 2015 and 2017 from 23 newly diagnosed infected individuals in the Philippines. Phylogenetic analysis showed that the newly characterized sequences were CRF01_AE (14), subtype B (3), CRF01/B recombinants (5) and a CRF01/CRF07/B recombinant (1). All 14 CRF01_AE formed a tight cluster, suggesting that they were derived from a single introduction. The time to the most recent common ancestor (tMRCA) for CRF01_AE in the Philippines was 1995 (1992-1998), about 10-15 years later than that of CRF01_AE in China and Thailand. All five CRF01/B recombinants showed distinct recombination patterns, suggesting ongoing recombination between the two predominant circulating viruses. The identification of partial CRF07_BC sequences in one CRF01/CRF07/B recombinant, not reported previously in the Philippines, indicated that CRF07_BC may have been recently introduced into that country from China, where CRF07_BC is prevalent. Our results show that the major epidemic strains may have shifted to an increased predominance of CRF01_AE and its recombinants, and that other genotypes such as CRF07_BC may have been introduced into the Philippines.

Molecular epidemiology of a primarily MSM acute HIV-1 cohort in Bangkok, Thailand and connections within networks of transmission in Asia

INTRODUCTION

Thailand plays a substantial role in global HIV-1 transmission of CRF01_AE. Worldwide, men who have sex with men (MSM) are at elevated risk for HIV-1 infection. Hence, understanding HIV-1 diversity in a primarily Thai MSM cohort with acute infection, and its connections to the broader HIV-1 transmission network in Asia is crucial for research and development of HIV-1 vaccines, treatment and cure.

METHODS

Subtypes and diversity of infecting viruses from individuals sampled from 2009 to 2015 within the RV254/SEARCH 010 cohort were assessed by multiregion hybridization assay (MHAbce), multiregion subtype-specific PCR assay (MSSPbce) and full-length single-genome sequencing (SGS). Phylogenetic analysis was performed by maximum likelihood. Pairwise genetic distances of envelope gp160 sequences obtained from the cohort and from Asia (Los Alamos National Laboratory HIV Database) were calculated to identify potential transmission networks.

RESULTS

MHAbce/MSSPbce results identified 81.6% CRF01_AE infecting strains in RV254. CRF01_AE/B recombinants and subtype B were found at 7.3% and 2.8% respectively. Western subtype B strains outnumbered Thai B' strains. Phylogenetic analysis revealed one C, one CRF01_AE/CRF02_AG recombinant and one CRF01_AE/B/C recombinant. Asian network analysis identified one hundred and twenty-three clusters, including five clusters of RV254 participants. None of the RV254 sequences clustered with non-RV254 sequences. The largest international cluster involved 15 CRF01_AE strains from China and Vietnam. The remaining clusters were mostly intracountry connections, of which 31.7% included Thai nodes and 43.1% included Chinese nodes.

CONCLUSION

While the majority of strains in Thailand are CRF01_AE and subtype B, emergence of unique recombinant forms (URFs) are found in a moderate fraction of new HIV-1 infections. Approaches to vaccine design and immunotherapeutics will need to monitor and consider the expanding proportion of recombinants and the increasing genetic diversity in the region. Identified HIV-1 transmission networks indicate ongoing spread of HIV-1 among MSM. As HIV-1 epidemics continue to expand in other Asian countries, transmission network analyses can inform strategies for prevention, intervention, treatment and cure.

HIV-1 Genetic Diversity Among Incident Infections in Mbeya, Tanzania

In preparation for vaccine trials, HIV-1 genetic diversity was surveyed between 2002 and 2006 through the Cohort Development study in the form of a retrospective and prospective observational study in and around the town of Mbeya in Tanzania's Southwest Highlands. This study describes the molecular epidemiology of HIV-1 strains obtained from 97 out of 106 incident HIV-1 infections identified in three subpopulations of participants (one rural, two urban) from the Mbeya area. Near full-genome or half-genome sequencing showed a subtype distribution of 40% C, 17% A1, 1% D, and 42% inter-subtype recombinants. Compared to viral subtyping results previously obtained from the retrospective phase of this study, the overall proportion of incident viral strains did not change greatly during the study course, suggesting maturity of the epidemic. A comparison to a current Phase I-II vaccine being tested in Africa shows ∼17% amino acid sequence difference between the gp120 of the vaccine and subtype C incident strains. Phylogenetic and recombinant breakpoint analysis of the incident strains revealed the emergence of CRF41_CD and many unique recombinants, as well as the presence of six local transmission networks most of which were confined to the rural subpopulation. In the context of vaccine cohort selection, these results suggest distinct infection transmission dynamics within these three geographically close subpopulations. The diversity and genetic sequences of the HIV-1 strains obtained during this study will greatly contribute to the planning, immunogen selection, and analysis of vaccine-induced immune responses observed during HIV-1 vaccine trials in Tanzania and neighboring countries.

Fast Dissemination of New HIV-1 CRF02/A1 Recombinants in Pakistan

A number of HIV-1 subtypes are identified in Pakistan by characterization of partial viral gene sequences. Little is known whether new recombinants are generated and how they disseminate since whole genome sequences for these viruses have not been characterized. Near full-length genome (NFLG) sequences were obtained by amplifying two overlapping half genomes or next generation sequencing from 34 HIV-1-infected individuals in Pakistan. Phylogenetic tree analysis showed that the newly characterized sequences were 16 subtype As, one subtype C, and 17 A/G recombinants. Further analysis showed that all 16 subtype A1 sequences (47%), together with the vast majority of sequences from Pakistan from other studies, formed a tight subcluster (A1a) within the subtype A1 clade, suggesting that they were derived from a single introduction. More in-depth analysis of 17 A/G NFLG sequences showed that five shared similar recombination breakpoints as in CRF02 (15%) but were phylogenetically distinct from the prototype CRF02 by forming a tight subcluster (CRF02a) while 12 (38%) were new recombinants between CRF02a and A1a or a divergent A1b viruses. Unique recombination patterns among the majority of the newly characterized recombinants indicated ongoing recombination. Interestingly, recombination breakpoints in these CRF02/A1 recombinants were similar to those in prototype CRF02 viruses, indicating that recombination at these sites more likely generate variable recombinant viruses. The dominance and fast dissemination of new CRF02a/A1 recombinants over prototype CRF02 suggest that these recombinant have more adapted and may become major epidemic strains in Pakistan.

Molecular evolution of the HIV-1 Thai epidemic between the time of RV144 immunogen selection to the execution of the vaccine efficacy trial

The RV144 HIV-1 vaccine trial (Thailand, 2003 to 2009), using immunogens genetically matched to the regional epidemic, demonstrated the first evidence of efficacy for an HIV-1 vaccine. Here we studied the molecular evolution of the HIV-1 epidemic from the time of immunogen selection to the execution of the efficacy trial. We studied HIV-1 genetic diversity among 390 volunteers who were deferred from enrollment in RV144 due to preexisting HIV-1 infection using a multiregion hybridization assay, full-genome sequencing, and phylogenetic analyses. The subtype distribution was 91.7% CRF01_AE, 3.5% subtype B, 4.3% B/CRF01_AE recombinants, and 0.5% dual infections. CRF01_AE strains were 31% more diverse than the ones from the 1990s Thai epidemic. Sixty-nine percent of subtype B strains clustered with the cosmopolitan Western B strains. Ninety-three percent of B/CRF01_AE recombinants were unique; recombination breakpoint analysis showed that these strains were highly embedded within the larger network that integrates recombinants from East/Southeast Asia. Compared to Thai sequences from the early 1990s, the distance to the RV144 immunogens increased 52% to 68% for CRF01_AE Env immunogens and 12% to 29% for subtype B immunogens. Forty-three percent to 48% of CRF01_AE sequences differed from the sequence of the vaccine insert in Env variable region 2 positions 169 and 181, which were implicated in vaccine sieve effects in RV144. In conclusion, compared to the molecular picture at the early stages of vaccine development, our results show an overall increase in the genetic complexity of viruses in the Thai epidemic and in the distance to vaccine immunogens, which should be considered at the time of the analysis of the trial results.

Disinformation squared: was the HIV-from-Fort-Detrick myth a Stasi success?

BACKGROUND

When in May 1983 the acquired immunodeficiency syndrome (AIDS) was first securely attributed to a virus, eventually called the human immunodeficiency virus (HIV), many controversies arose. Among these was one centering on HIV's origin. A startling hypothesis, called here the "HIV-from-Fort-Detrick myth," asserted that HIV had been a product, accidental or intentional, of bioweaponry research. While its earliest identifiable contributors were in the West, this myth's most dynamic propagators were in the East. The Soviet security service, the KGB, took "active measures" to create and disseminate AIDS disinformation beginning no later than July 1983 and ending no earlier than October 1987. The East German security service, a complex bureaucracy popularly known as "the Stasi," was involved, too, but how early, how deeply, how uniformly, how ably, and how successfully has not been clear. Following German reunification, claims arose attributing to the Stasi the masterful execution of ingenious elements in a disinformation campaign they helped shape and soon came to dominate. We have tested these claims.

QUESTION

Was the HIV-from-Fort-Detrick myth a Stasi success?

METHODS

Primary sources were documents and photographs assembled by the Ministry of State Security (MfS) of the German Democratic Republic (GDR or East Germany), the Ministry of Interior of the People's Republic of Bulgaria, and the United States Department of State; the estate of myth principals Jakob and Lilli Segal; the "AIDS box" in the estate of East German literary figure Stefan Heym; participant-observer recollections, interviews, and correspondence; and expert interviews. We examined secondary sources in light of primary sources.

FINDINGS

The HIV-from-Fort-Detrick myth had debuted in print in India in 1983 and had been described in publications worldwide prior to 1986, the earliest year for which we found any Stasi document mentioning the myth in any context. Many of the myth's exponents were seemingly independent conspiracy theorists. Its single most creative exponent was Jakob Segal, an idiosyncratic Soviet biologist long resident in, and long retired in, the GDR. Segal applied to the myth a thin but tenacious layer of plausibility. We could not exclude a direct KGB influence on him but found no evidence demonstrating it. The Stasi did not direct his efforts and had difficulty tracking his activities. The Stasi were prone to interpretive error and self-aggrandizement. They credited themselves with successes they did not achieve, and, in one instance, failed to appreciate that a major presumptive success had actually been a fiasco. Senior Stasi officers came to see the myth's propagation as an embarrassment threatening broader interests, especially the GDR's interest in being accepted as a scientifically sophisticated state. In 1986, 1988, and 1989, officers of HV A/X, the Stasi's disinformation and "active measures" department, discussed the myth in meetings with the Bulgarian secret service. In the last of these meetings, HV A/X officers tried to interest their Bulgarian counterparts in taking up, or taking over, the myth's propagation. Further efforts, if any, were obscured by collapse of the East German and Bulgarian governments.

CONCLUSION

No, the HIV-from-Fort-Detrick myth was not a Stasi success. Impressions to the contrary can be attributed to reliance on presumptions, boasts, and inventions. Presumptions conceding to the Stasi an extraordinary operational efficiency and an irresistible competence - qualities we could not confirm in this case - made the boasts and inventions more convincing than their evidentiary basis, had it been known, would have allowed. The result was disinformation about disinformation, a product we call "disinformation squared."