Increasing HIV type 1 polymorphic diversity but no resistance to antiretroviral drugs in untreated patients from Central African Republic: a 2005 study

Phylogenetic evidence of extensive spatial mixing of diverse HIV-1 group M lineages within Cameroon but not between its neighbours

From the perspective of developing relevant interventions for treating HIV and controlling its spread, it is particularly important to comprehensively understand the underlying diversity of the virus, especially in countries where the virus has been present and evolving since the cross-species transmission event that triggered the global pandemic. Here, we generate and phylogenetically analyse sequences derived from the gag-protease (2010 bp; n = 115), partial integrase (345 bp; n = 36), and nef (719 bp; n = 321) genes of HIV-1 group M (HIV-1M) isolates sampled between 2000 and 2022 from two cosmopolitan cities and 40 remote villages of Cameroon. While 52.4% of all sequenced viruses belonged to circulating recombinant form (CRF) 02_AG (CRF02_AG), the remainder were highly diverse, collectively representing seven subtypes and sub-subtypes, eight CRFs, and 36 highly divergent lineages that fall outside the established HIV-1M classification. Additionally, in 77 samples for which at least two genes were typed, 31% of the studied viruses apparently had fragments from viruses belonging to different clades. Furthermore, we found that the distribution of HIV-1M populations is similar between different regions of Cameroon. In contrast, HIV-1M demographics in Cameroon differ significantly from those in its neighbouring countries in the Congo Basin (CB). In phylogenetic trees, viral sequences cluster according to the countries where they were sampled, suggesting that while there are minimal geographical or social barriers to viral dissemination throughout Cameroon, there is strongly impeded dispersal of HIV-1M lineages between Cameroon and other locations of the CB. This suggests that the apparent stability of highly diverse Cameroonian HIV-1M populations may be attributable to the extensive mixing of human populations within the country and the concomitant trans-national movements of major lineages with very similar degrees of fitness; coupled with the relatively infrequent inter-national transmission of these lineages from neighbouring countries in the CB.

New Genomes from the Congo Basin Expand History of CRF01_AE Origin and Dissemination

Although the first HIV circulating recombinant form (CRF01_AE) is the predominant strain in many Asian countries, it is uncommonly found in the Congo Basin from where it first originated. To fill the gap in the evolutionary history of this important strain, we sequenced near complete genomes from HIV samples with subgenomic CRF01_AE regions collected in Cameroon and the Democratic Republic of the Congo from 2001 to 2006. HIV genomes were generated from N = 13 plasma specimens by next-generation sequencing of metagenomic libraries prepared with spiked primers targeting HIV, followed by Sanger gap-filling. Genome sequences were aligned to reference strains, including Asian and African CRF01_AE sequences, and evaluated by phylogenetic and recombinant analysis to identify four CRF01_AE strains from Cameroon. We also identified two CRF02, one CRF27, and six unique recombinant form genomes (01|A1|G, 01|02|F|U, F|G|01, A1|D|01, F|G|01, and A1|G|01). Phylogenetic analysis, including the four new African CRF01_AE genomes, placed these samples as a bridge between basal Central African Republic CRF01_AE strains and all Asian, European, and American CRF01_AE strains. Molecular dating confirmed previous estimates indicating that the most recent common CRF01_AE ancestor emerged in the early 1970s (1968–1970) and spread beyond Africa around 1980 to Asia. The new sequences and analysis presented in this study expand the molecular history of the CRF01_AE clade, and are illustrated in an interactive Next Strain phylogenetic tree, map, and timeline at (https://nextstrain.org/community/EduanWilkinson/hiv-1_crf01).

Unravelling the complicated evolutionary and dissemination history of HIV-1M subtype A lineages

Subtype A is one of the rare HIV-1 group M (HIV-1M) lineages that is both widely distributed throughout the world and persists at high frequencies in the Congo Basin (CB), the site where HIV-1M likely originated. This, together with its high degree of diversity suggests that subtype A is amongst the fittest HIV-1M lineages. Here we use a comprehensive set of published near full-length subtype A sequences and A-derived genome fragments from both circulating and unique recombinant forms (CRFs/URFs) to obtain some insights into how frequently these lineages have independently seeded HIV-1M sub-epidemics in different parts of the world. We do this by inferring when and where the major subtype A lineages and subtype A-derived CRFs originated. Following its origin in the CB during the 1940s, we track the diversification and recombination history of subtype A sequences before and during its dissemination throughout much of the world between the 1950s and 1970s. Collectively, the timings and numbers of detectable subtype A recombination and dissemination events, the present broad global distribution of the sub-epidemics that were seeded by these events, and the high prevalence of subtype A sequences within the regions where these sub-epidemics occurred, suggest that ancestral subtype A viruses (and particularly sub-subtype A1 ancestral viruses) may have been genetically predisposed to become major components of the present epidemic.

Phylodynamic applications in 21st century global infectious disease research

BACKGROUND

Phylodynamics, the study of the interaction between epidemiological and pathogen evolutionary processes within and among populations, was originally defined in the context of rapidly evolving viruses and used to characterize transmission dynamics. The concept of phylodynamics has evolved since the early 21st century, extending its reach to slower-evolving pathogens, including bacteria and fungi, and to the identification of influential factors in disease spread and pathogen population dynamics.

RESULTS

The phylodynamic approach has now become a fundamental building block for the development of comparative phylogenetic tools capable of incorporating epidemiological surveillance data with molecular sequences into a single statistical framework. These innovative tools have greatly enhanced scientific investigations of the temporal and geographical origins, evolutionary history, and ecological risk factors associated with the growth and spread of viruses such as human immunodeficiency virus (HIV), Zika, and dengue and bacteria such as Methicillin-resistant Staphylococcus aureus.

CONCLUSIONS

Capitalizing on an extensive review of the literature, we discuss the evolution of the field of infectious disease epidemiology and recent accomplishments, highlighting the advancements in phylodynamics, as well as the challenges and limitations currently facing researchers studying emerging pathogen epidemics across the globe.

Time-scale of minor HIV-1 complex circulating recombinant forms from Central and West Africa

Background

Several HIV-1 circulating recombinant forms with a complex mosaic structure (CRFs_cpx) circulate in central and western African regions. Here we reconstruct the evolutionary history of some of these complex CRFs (09_cpx, 11_cpx, 13_cpx and 45_cpx) and further investigate the dissemination dynamic of the CRF11_cpx clade by using a Bayesian coalescent-based method.

Results

The analysis of two HIV-1 datasets comprising 181 pol (36 CRF09_cpx, 116 CRF11_cpx, 20 CRF13_cpx and 9 CRF45_cpx) and 125 env (12 CRF09_cpx, 67 CRF11_cpx, 17 CRF13_cpx and 29 CRF45_cpx) sequences pointed to quite consistent onset dates for CRF09_cpx (~1966: 1958–1979), CRF11_cpx (~1957: 1950–1966) and CRF13_cpx (~1965: 1958–1973) clades; while some divergence was found for the estimated date of origin of CRF45_cpx clade [pol = 1970 (1964–1976); env = 1960 (1952–1969)]. Phylogeographic reconstructions indicate that the HIV-1 CRF11_cpx clade most probably emerged in Cameroon and from there it was first disseminated to the Central Africa Republic and Chad in the early 1970s and to other central and western African countries from the early 1980s onwards. Demographic reconstructions suggest that the CRF11_cpx epidemic grew between 1960 and 1990 with a median exponential growth rate of 0.27 year⁻¹, and stabilized after.

Conclusions

These results reveal that HIV-1 CRFs_cpx clades have been circulating in Central Africa for a period comparable to other much more prevalent HIV-1 group M lineages. Cameroon was probably the epicenter of dissemination of the CRF11_cpx clade that seems to have experienced a long epidemic growth phase before stabilization. The epidemic growth of the CRF11_cpx clade was roughly comparable to other HIV-1 group M lineages circulating in Central Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0824-8) contains supplementary material, which is available to authorized users.

Predominance of CRF06_cpx and Transmitted HIV Resistance in Algeria: Update 2013-2014

Since 2008, no data on HIV diversity or the transmission rate of HIV resistance mutations in naive patients have been presented for Algeria, a country of MENA region. Between 2013 and 2014, we studied 152 samples including 89 naive patients. The current study describes the change in HIV diversity in Algeria with the predominance of CRF06_cpx and the huge increase of transmitted HIV resistance, which now reaches 15%.

Spatiotemporal dynamics of the HIV-1 subtype G epidemic in West and Central Africa

The human immunodeficiency virus type 1 (HIV-1) subtype G is the second most prevalent HIV-1 clade in West Africa, accounting for nearly 30% of infections in the region. There is no information about the spatiotemporal dynamics of dissemination of this HIV-1 clade in Africa. To this end, we analyzed a total of 305 HIV-1 subtype G pol sequences isolated from 11 different countries from West and Central Africa over a period of 20 years (1992 to 2011). Evolutionary, phylogeographic and demographic parameters were jointly estimated from sequence data using a Bayesian coalescent-based method. Our analyses indicate that subtype G most probably emerged in Central Africa in 1968 (1956–1976). From Central Africa, the virus was disseminated to West and West Central Africa at multiple times from the middle 1970s onwards. Two subtype G strains probably introduced into Nigeria and Togo between the middle and the late 1970s were disseminated locally and to neighboring countries, leading to the origin of two major western African clades (G_WA-I and G_WA-II). Subtype G clades circulating in western and central African regions displayed an initial phase of exponential growth followed by a decline in growth rate since the early/middle 1990s; but the mean epidemic growth rate of G_WA-I (0.75 year⁻¹) and G_WA-II (0.95 year⁻¹) clades was about two times higher than that estimated for central African lineages (0.47 year⁻¹). Notably, the overall evolutionary and demographic history of G_WA-I and G_WA-II clades was very similar to that estimated for the CRF06_cpx clade circulating in the same region. These results support the notion that the spatiotemporal dissemination dynamics of major HIV-1 clades circulating in western Africa have probably been shaped by the same ecological factors.

Towards the embodiment of biosocial resistance? How to account for the unexpected effects of antiretroviral scale-up in the Central African Republic

At the fringes of the unprecedented medication scale-up in the treatment of HIV, many African countries have experienced dramatic antiretroviral drug stock-outs. Usually considered the result of irrational decisions on behalf of local politicians, programme managers and even patients (who are stigmatised as immoral), these problems seem not to be so exceptional. However, ethnographic attention to the social consequences of the presence and absence of antiretroviral drugs in the Central African Republic (CAR) suggests that these stock-outs entail far more than logistical failures. In 2010 and 2011 in the CAR, major antiretroviral treatment (ARV) stock-outs resulted in the renewal of 'therapeutic' social ties and also significant social resistance and defiance. While this paper explores reasons for the shortage, its focus is on subsequent popular reactions to it, particularly among people who are HIV-positive and dependent on ARVs. The exceptional and ambiguous consequences of these drug stock-outs raise new concerns relevant to the politics of global public health.

Phylodynamics of the HIV-1 epidemic in Cuba

Previous studies have shown that the HIV-1 epidemic in Cuba displayed a complex molecular epidemiologic profile with circulation of several subtypes and circulating recombinant forms (CRF); but the evolutionary and population history of those viral variants remains unknown. HIV-1 pol sequences of the most prevalent Cuban lineages (subtypes B, C and G, CRF18_cpx, CRF19_cpx, and CRFs20/23/24_BG) isolated between 1999 and 2011 were analyzed. Maximum-likelihood analyses revealed multiple introductions of subtype B (n≥66), subtype C (n≥10), subtype G (n≥8) and CRF18_cpx (n≥2) viruses in Cuba. The bulk of HIV-1 infections in this country, however, was caused by dissemination of a few founder strains probably introduced from North America/Europe (clades B_CU-I and B_CU-II), east Africa (clade C_CU-I) and central Africa (clades G_CU, CRF18_CU and CRF19_CU), or locally generated (clades CRFs20/23/24_BG). Bayesian-coalescent analyses show that the major HIV-1 founder strains were introduced into Cuba during 1985–1995; whereas the CRFs_BG strains emerged in the second half of the 1990s. Most HIV-1 Cuban clades appear to have experienced an initial period of fast exponential spread during the 1990s and early 2000s, followed by a more recent decline in growth rate. The median initial growth rate of HIV-1 Cuban clades ranged from 0.4 year⁻¹ to 1.6 year⁻¹. Thus, the HIV-1 epidemic in Cuba has been a result of the successful introduction of a few viral strains that began to circulate at a rather late time of the AIDS pandemic, but then were rapidly disseminated through local transmission networks.

Virological response and resistance profiles after 24 months of first-line antiretroviral treatment in adults living in Bangui, Central African Republic

The rate of virological failure was assessed in 386 adult patients attending the Centre National Hospitalier Universitaire of Bangui, the capital city of the Central African Republic (CAR), receiving their first-line antiretroviral (ARV) drug regimen for 24 months, according to the World Health Organization (WHO) recommendations. In addition, genotypic resistance testing was carried out in 45 of 145 randomly selected patients whose plasma HIV-1 RNA load was detectable. Overall, 28.5% of ARV-treated patients were in virological failure (e.g., HIV-1 RNA >3.7 log(10) copies/ml). Twenty-four percent of patients in virological failure showed wild-type viruses, likely indicating poor adherence. Even after excluding the M184V mutation, all 76% of patients in virological failure displayed viruses harboring at least one major drug resistance mutation to nucleoside reverse transcriptase inhibitors (NRTI), non-NRTI, or protease inhibitors. Whereas the second-line regimen proposed by the 2010 WHO recommendations, including zidovudine, tenofovir, lopinavir, and atazanavir, could be effective in more than 90% of patients in virological failure with resistant viruses, the remaining patients showed genotypic profiles highly predictive of resistance to the usual WHO second-line regimen, including complex genotypic profiles diagnosed only by genotypic resistance tests in some patients. In conclusion, our observations highlight the high frequency of therapeutic failure in ARV-treated adults in this study, as well as the urgent and absolute need for improving viral load assessment in the CAR to prevent and/or, from now on, to monitor therapeutic failure.

Phylodynamics of the HIV-1 CRF02_AG clade in Cameroon

Evolutionary analyses have revealed an origin of pandemic HIV-1 group M in the Congo River basin in the first part of the XX century, but the patterns of historical viral spread in or around its epicentre remain largely unexplored. Here, we combine epidemiologic and molecular sequence data to investigate the spatiotemporal patterns of the CRF02_AG clade. By explicitly integrating prevalence counts and genetic population size estimates we date the epidemic emergence of CRF02_AG at 1973.1 (1972.1, 1975.3, 95% CI). To infer the phylogeographic signature of this clade at a regional scale, we analyze pol and env time-stamped sequence data from 10 countries using a Bayesian phylogeographic approach based on an asymmetric discretized diffusion model. Our data confirms a spatial origin of CRF02_AG in the Democratic Republic of Congo (DRC) and suggests that viral dissemination to Cameroon occurred at an early stage of the evolutionary history of CRF02_AG. We find considerable support for epidemiological linkage between neighbour countries. Compilation of ethnographic data suggested that well-supported viral migration did not reflect sustained human migratory flows. Finally, using sequence data from 15 locations in Cameroon, we use relaxed random walk models to explore the spatiotemporal dynamics of CRF02_AG at a finer geographical detail. Phylogeographic dispersal in continuous space reveals that at least two distinct CRF02_AG lineages are circulating in overlapping regions that are evolving at different evolutionary and diffusion rates. In conclusion, by combining molecular and epidemiological data, our results provide a time scale for CRF02_AG, early 70s, place its spatial root in the DRC within the putative root of group-M diversity and propose a scenario of chance-exportation events for the spatiotemporal patterns of a successful HIV-1 lineage both at a regional and country-scale.

First data on HIV-1 resistance mutations to antiretroviral drugs in Central African Republic

In a background of high genomic HIV-1 variability with a predominance of CRF11_cpx and CRF22_01A1, we have studied the emergence of resistance mutations in isolates from Central African patients at failure of d4T-AZT/3TC/NVP-EFV plus two at failure of a PI-including regimen; the resistance mutations observed are those which are expected on HIV-1 subtype B.

Transmitted drug resistance in nonsubtype B HIV-1 infection

HIV-1 nonsubtype B variants account for the majority of HIV infections worldwide. Drug resistance in individuals who have never undergone antiretroviral therapy can lead to early failure and limited treatment options and, therefore, is an important concern. Evaluation of reported transmitted drug resistance (TDR) is challenging owing to varying definitions and study designs, and is further complicated by HIV-1 subtype diversity. In this article, we discuss the importance of various mutation lists for TDR definition, summarize TDR in nonsubtype B HIV-1 and highlight TDR reporting and interpreting challenges in the context of HIV-1 diversity. When examined carefully, TDR in HIV-1 non-B protease and reverse transcriptase is still relatively low in most regions. Whether it will increase with time and therapy access, as observed in subtype-B-predominant regions, remains to be determined.

HIV type 1 group M subtype G in Cameroon: five genome sequences

Near full-length viral genome sequences were obtained for five putative subtype G candidates identified in HIV-infected Cameroonian blood donors, based on partial genome sequences for the gag, pol, and env regions. Phylogenetic analysis of the genome sequences shows that all five strains are pure subtype G with no indication of intersubtype recombination. The Cameroon subtype G sequences did not form a geographically based subcluster and were intermixed within the subtype G branch with isolates from several different countries. HIV-1 group M subtype G accounts for only 4.5% of HIV infections in Cameroon. However, genome segments of subtype G are present in 67% of all infections and 80% of infections due to intersubtype recombinant strains in Cameroon. The addition of five subtype G genome sequences to the HIV database may contribute to a better understanding of the origins and classification of HIV-1 subtypes and CRFs.

High HIV Type 1 prevalence and wide genetic diversity with dominance of recombinant strains but low level of antiretroviral drug-resistance mutations in untreated patients in northeast Gabon, Central Africa

The northeast of Gabon, central Africa is characterized by high population density and a high rate of immigration from the surrounding countries. To determine the prevalence, circulating subtypes, and antiretroviral resistance mutations of HIV-1, 810 blood samples were collected from the general population of the two main cities (Oyem and Makokou) of this region. Of these, 61 (7.5%) were found to be positive for HIV-1. Analysis of the env (gp120), pol, and gag (p24) sequences as well as phylogenetic analyses showed at least eight different viral lineages. The most prevalent strains were CRF02 recombinants, followed by subtypes A, D, and C. The remaining strains were found to be F, J, G, and also, for the first time in Gabon, the recombinant form CRF11cpx. Analysis of antiretroviral drug-resistance mutations in protease and reverse transcriptase from this untreated population showed a low level of specific mutations. These mutations were associated with subtype polymorphism rather than with resistance to antiretroviral drugs. The wide diversity and the emergence of recombinant strains are in accordance with the rapid spread of new HIV strains in the population and, thus, the dynamic evolution of the epidemic.

Identification of new CRF43_02G and CRF25_cpx in Saudi Arabia based on full genome sequence analysis of six HIV type 1 isolates

Recently, we reported a high level of HIV-1 strain diversity in patients at the King Faisal Specialist Hospital and Research Center, Jeddah, Saudi Arabia. Based on phylogenetic analysis of gag p24, pol integrase, and env gp41 sequences, subtypes A, B, C, D, and G, and CRF02_AG, as well as unique recombinant forms were identified. Subtype G accounted for 25% of the infections in the Saudi population and this high prevalence was unexpected. Although subtype G is found in west central Africa, pure subtype G strains are uncommon. To further characterize the subtype G infections in Saudi Arabia, six strains that appeared to be pure subtype G were selected for full genome sequencing. Near full-length genomes were obtained using RT-PCR amplification to generate overlapping fragments from viral RNA extracted from plasma. The six strains are not subtype G throughout their entire genome. Four isolates have a recombinant structure composed of CRF02_AG and subtype G and share three identical breakpoints. This recombinant form defines a new CRF designated CRF43_02G. The remaining two isolates are CRF25_cpx, a circulating recombinant form identified in Cameroon composed of subtypes A and G and unclassified segments. Reanalysis of the previously reported Saudi HIV-1 partial genome sequences revealed additional isolates classified as CRF43_02G and CRF25_cpx and one isolate was reclassified to CRF22_01A. Identification of CRF43_02G in Saudi Arabia could indicate a transmission network within the country. Alternatively, the new CRF could have been introduced from an external source where this CRF is not yet recognized.

Human immunodeficiency virus type 1 seroprevalence and antiretroviral drug resistance-associated mutations in miners in Gabon, central Africa

Miners in sub-Saharan African are known to have an extremely high prevalence of HIV-1 infection. We therefore evaluated the prevalence of HIV-1 infection among manganese miners in Gabon, central Africa and examined the diversity of HIV-1 strains by characterizing the polymorphism of the pol gene in order to observe drug resistance-associated mutations. In 857 samples tested, the HIV-1 prevalence was 2.9%. By pol sequence analysis, we showed that all the HIV-1 strains belonged to group M, with a majority of CRF02_AG (57%) followed by subtype A (9%) and CRF01_AE or subtype B (4%). The remaining HIV-1 strains demonstrated discordant genomic results and exhibited a mosaic pol genome (30%). Most of the mutations detected in pol coding regions corresponded to the subtype polymorphism, with no specific antiretroviral drug resistance. To avoid the rapid emergence of resistant viruses in this part of central Africa, continuous surveillance of the circulation of drug-resistant viruses must be maintained to guide treatment strategies.

Rapid spread and genetic diversification of HIV type 1 subtype C in a rural area of southern Mozambique

In this study, we analyzed the human immunodeficiency type 1 (HIV-1) viruses circulating between 1999 and 2004 in antiretroviral-naive women from a rural area of southern Mozambique. Nucleotide sequencing of the HIV-1 long terminal repeat (LTR) U3, envelope (env) C2V3C3, and protease (pr) genomic regions was performed from women sera samples collected in 1999 and 2004. Phylogenetic analysis revealed that all amplified sequences belonged to subtype C. Although env sequences were predominantly CCR5-tropic (R5), CXCR4-tropic (X4) variants were also identified (13%). Both 1999 and 2004 sequences were widely dispersed across multiple clusters and were related to different reference sequences from neighboring countries. Sequences from 2004 showed significantly more nucleotide genetic diversity than sequences from 1999. Importantly, genetic diversification was also observed at the pr and env amino acid level, suggesting that positive selection forces were implicated in the viral diversification. These results indicate the rapid spread and diversification of subtype C virus in Mozambique where HIV-1 prevalence in the Manhiça antenatal clinic reached 23% in 2004.