The molecular population genetics of HIV-1 group O

CGRWDL: alignment-free phylogeny reconstruction method for viruses based on chaos game representation weighted by dynamical language model

Traditional alignment-based methods meet serious challenges in genome sequence comparison and phylogeny reconstruction due to their high computational complexity. Here, we propose a new alignment-free method to analyze the phylogenetic relationships (classification) among species. In our method, the dynamical language (DL) model and the chaos game representation (CGR) method are used to characterize the frequency information and the context information of k-mers in a sequence, respectively. Then for each DNA sequence or protein sequence in a dataset, our method converts the sequence into a feature vector that represents the sequence information based on CGR weighted by the DL model to infer phylogenetic relationships. We name our method CGRWDL. Its performance was tested on both DNA and protein sequences of 8 datasets of viruses to construct the phylogenetic trees. We compared the Robinson-Foulds (RF) distance between the phylogenetic tree constructed by CGRWDL and the reference tree by other advanced methods for each dataset. The results show that the phylogenetic trees constructed by CGRWDL can accurately classify the viruses, and the RF scores between the trees and the reference trees are smaller than that with other methods.

Replicative fitness and pathogenicity of primate lentiviruses in lymphoid tissue, primary human and chimpanzee cells: relation to possible jumps to humans

BACKGROUND

Simian immunodeficiency viruses (SIV) have been jumping between non-human primates in West/Central Africa for thousands of years and yet, the HIV-1 epidemic only originated from a primate lentivirus over 100 years ago.

METHODS

This study examined the replicative fitness, transmission, restriction, and cytopathogenicity of 22 primate lentiviruses in primary human lymphoid tissue and both primary human and chimpanzee peripheral blood mononuclear cells.

FINDINGS

Pairwise competitions revealed that SIV from chimpanzees (cpz) had the highest replicative fitness in human or chimpanzee peripheral blood mononuclear cells, even higher fitness than HIV-1 group M strains responsible for worldwide epidemic. The SIV strains belonging to the "HIV-2 lineage" (including SIVsmm, SIVmac, SIVagm) had the lowest replicative fitness. SIVcpz strains were less inhibited by human restriction factors than the "HIV-2 lineage" strains. SIVcpz efficiently replicated in human tonsillar tissue but did not deplete CD4+ T-cells, consistent with the slow or nonpathogenic disease observed in most chimpanzees. In contrast, HIV-1 isolates and SIV of the HIV-2 lineage were pathogenic to the human tonsillar tissue, almost independent of the level of virus replication.

INTERPRETATION

Of all primate lentiviruses, SIV from chimpanzees appears most capable of infecting and replicating in humans, establishing HIV-1. SIV from other Old World monkeys, e.g. the progenitor of HIV-2, replicate slowly in humans due in part to restriction factors. Nonetheless, many of these SIV strains were more pathogenic than SIVcpz. Either SIVcpz evolved into a more pathogenic virus while in humans or a rare SIVcpz, possibly extinct in chimpanzees, was pathogenic immediately following the jump into human.

FUNDING

Support for this study to E.J.A. was provided by the NIH/NIAID R01 AI49170 and CIHR project grant 385787. Infrastructure support was provided by the NIH CFAR AI36219 and Canadian CFI/Ontario ORF 36287. Efforts of J.A.B. and N.J.H. was provided by NIH AI099473 and for D.H.C., by VA and NIH AI AI080313.

Disease-causing human viruses: novelty and legacy

About 270 viruses are known to infect humans. Some of these viruses have been known for centuries, whereas others have recently emerged. During their evolutionary history, humans have moved out of Africa to populate the world. In historical times, human migrations resulted in the displacement of large numbers of people. All these events determined the movement and dispersal of human-infecting viruses. Technological advances have resulted in the characterization of the genetic variability of human viruses, both in extant and in archaeological samples. Field studies investigated the diversity of viruses hosted by other animals. In turn, these advances provided insight into the evolutionary history of human viruses back in time and defined the key events through which they originated and spread.

Contrasting patterns in the early stage of SARS-CoV-2 evolution between humans and minks

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor binding domain and receptor binding motif. An excess of high frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies the virus may not be pre-adapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.

Site-Specific Evolutionary Rate Shifts in HIV-1 and SIV

Site-specific evolutionary rate shifts are defined as protein sites, where the rate of substitution has changed dramatically across the phylogeny. With respect to a given clade, sites may either undergo a rate acceleration or a rate deceleration, reflecting a site that was conserved and became variable, or vice-versa, respectively. Sites displaying such a dramatic evolutionary change may point to a loss or gain of function at the protein site, reflecting adaptation, or they may indicate epistatic interactions among sites. Here, we analyzed full genomes of HIV and SIV-1 and identified 271 rate-shifting sites along the HIV-1/SIV phylogeny. The majority of rate shifts occurred at long branches, often corresponding to cross-species transmission branches. We noted that in most proteins, the number of rate accelerations and decelerations was equal, and we suggest that this reflects epistatic interactions among sites. However, several accessory proteins were enriched for either accelerations or decelerations, and we suggest that this may be a signature of adaptation to new hosts. Interestingly, the non-pandemic HIV-1 group O clade exhibited a substantially higher number of rate-shift events than the pandemic group M clade. We propose that this may be a reflection of the height of the species barrier between gorillas and humans versus chimpanzees and humans. Our results provide a genome-wide view of the constraints operating on proteins of HIV-1 and SIV.

Phylogenetic Analysis of HIV-1 Genomes Based on the Position-Weighted K-mers Method

HIV-1 viruses, which are predominant in the family of HIV viruses, have strong pathogenicity and infectivity. They can evolve into many different variants in a very short time. In this study, we propose a new and effective alignment-free method for the phylogenetic analysis of HIV-1 viruses using complete genome sequences. Our method combines the position distribution information and the counts of the k-mers together. We also propose a metric to determine the optimal k value. We name our method the Position-Weighted k-mers (PWkmer) method. Validation and comparison with the Robinson-Foulds distance method and the modified bootstrap method on a benchmark dataset show that our method is reliable for the phylogenetic analysis of HIV-1 viruses. PWkmer can resolve within-group variations for different known subtypes of Group M of HIV-1 viruses. This method is simple and computationally fast for whole genome phylogenetic analysis.

Structural perspectives on HIV-1 Vif and APOBEC3 restriction factor interactions

Human immunodeficiency virus (HIV) is a retroviral pathogen that targets human immune cells such as CD4+ T cells, macrophages, and dendritic cells. The human apolipoprotein B mRNA- editing catalytic polypeptide 3 (APOBEC3 or A3) cytidine deaminases are a key class of intrinsic restriction factors that inhibit replication of HIV. When HIV-1 enters the cell, the immune system responds by inducing the activation of the A3 family proteins, which convert cytosines to uracils in single-stranded DNA replication intermediates, neutralizing the virus. HIV counteracts this intrinsic immune response by encoding a protein termed viral infectivity factor (Vif). Vif targets A3 to an E3 ubiquitin ligase complex for poly-ubiquitination and proteasomal degradation. Vif is unique in that it can recognize and counteract multiple A3 restriction factor substrates. Structural biology studies have provided significant insights into the overall architectures and functions of Vif and A3 proteins; however, a structure of the Vif-A3 complex has remained elusive. In this review, we summarize and reanalyze experimental data from recent structural, biochemical, and functional studies to provide key perspectives on the residues involved in Vif-A3 protein-protein interactions.

Evolutionary dynamics of Ceratonova species (Cnidaria: Myxozoa) reveal different host adaptation strategies

The myxozoan parasite Ceratonova shasta is an important pathogen that infects multiple species of Pacific salmonids. Ongoing genetic surveillance has revealed stable host-parasite relationships throughout the parasite's endemic range. We applied Bayesian phylogenetics to test specific hypotheses about the evolution of these host-parasite relationships within the well-studied Klamath River watershed in Oregon and California, USA. The results provide statistical support that different genotypes of C. shasta are distinct lineages of one species, which is related to two other Ceratonova species in the same ecosystems; Ceratonova X in speckled dace and C. gasterostea in threespine stickleback. Furthermore, we found strong support for the hypothesis that C. shasta type 0 in native steelhead trout and type I in Chinook salmon each evolved with a specialist host adaptation strategy, while C. shasta type II in coho salmon resulted from a generalist host adaptation strategy. Inferred date and host species of the most recent common ancestor of extant Klamath basin types indicate that it occurred between 14,000 and 21,000 years ago, and most likely infected a native steelhead or rainbow trout host.

The epidemic emergence of HIV: what novel enabling factors were involved?

Humans acquired retroviruses from simians, mainly through bushmeat handling. All epidemically successful HIV groups started to spread in early 20th century, contrasting with the antiquity of T-cell lymphotropic viruses, implying that novel enabling factors were involved in HIV emergence. Here we review the Parenteral Serial Transmission and the Enhanced Heterosexual Transmission hypotheses for the adaptation and early spread of HIV. Epidemic start roughly coincides in time with peak genital ulcer disease in cities, suggesting a major role for sexual transmission. Only ill-adapted and rare HIV groups emerged after approximately 1950, when injections and transfusions attained their maximal levels, suggesting that if parenteral serial transmission was necessary for HIV adaptation, it had to be complemented by sexual transmission for HIV to reach epidemic potential.

[Formula: see text]

Male Circumcision and the Epidemic Emergence of HIV-2 in West Africa

BACKGROUND

Epidemic HIV-2 (groups A and B) emerged in humans circa 1930-40. Its closest ancestors are SIVsmm infecting sooty mangabeys from southwestern Côte d'Ivoire. The earliest large-scale serological surveys of HIV-2 in West Africa (1985-91) show a patchy spread. Côte d'Ivoire and Guinea-Bissau had the highest prevalence rates by then, and phylogeographical analysis suggests they were the earliest epicenters. Wars and parenteral transmission have been hypothesized to have promoted HIV-2 spread. Male circumcision (MC) is known to correlate negatively with HIV-1 prevalence in Africa, but studies examining this issue for HIV-2 are lacking.

METHODS

We reviewed published HIV-2 serosurveys for 30 cities of all West African countries and obtained credible estimates of real prevalence through Bayesian estimation. We estimated past MC rates of 218 West African ethnic groups, based on ethnographic literature and fieldwork. We collected demographic tables specifying the ethnic partition in cities. Uncertainty was incorporated by defining plausible ranges of parameters (e.g. timing of introduction, proportion circumcised). We generated 1,000 sets of past MC rates per city using Latin Hypercube Sampling with different parameter combinations, and explored the correlation between HIV-2 prevalence and estimated MC rate (both logit-transformed) in the 1,000 replicates.

RESULTS AND CONCLUSIONS

Our survey reveals that, in the early 20th century, MC was far less common and geographically more variable than nowadays. HIV-2 prevalence in 1985-91 and MC rates in 1950 were negatively correlated (Spearman rho = -0.546, IQR: -0.553--0.546, p≤0.0021). Guinea-Bissau and Côte d'Ivoire cities had markedly lower MC rates. In addition, MC was uncommon in rural southwestern Côte d'Ivoire in 1930.The differential HIV-2 spread in West Africa correlates with different historical MC rates. We suggest HIV-2 only formed early substantial foci in cities with substantial uncircumcised populations. Lack of MC in rural areas exposed to bushmeat may have had a role in successful HIV-2 emergence.

Genomic Analysis of the Emergence, Evolution, and Spread of Human Respiratory RNA Viruses

The emergence and reemergence of rapidly evolving RNA viruses-particularly those responsible for respiratory diseases, such as influenza viruses and coronaviruses-pose a significant threat to global health, including the potential of major pandemics. Importantly, recent advances in high-throughput genome sequencing enable researchers to reveal the genomic diversity of these viral pathogens at much lower cost and with much greater precision than they could before. In particular, the genome sequence data generated allow inferences to be made on the molecular basis of viral emergence, evolution, and spread in human populations in real time. In this review, we introduce recent computational methods that analyze viral genomic data, particularly in combination with metadata such as sampling time, geographic location, and virulence. We then outline the insights these analyses have provided into the fundamental patterns and processes of evolution and emergence in human respiratory RNA viruses, as well as the major challenges in such genomic analyses.

Comparative Microbial Genomics and Forensics

Forensic science concerns the application of scientific techniques to questions of a legal nature and may also be used to address questions of historical importance. Forensic techniques are often used in legal cases that involve crimes against persons or property, and they increasingly may involve cases of bioterrorism, crimes against nature, medical negligence, or tracing the origin of food- and crop-borne disease. Given the rapid advance of genome sequencing and comparative genomics techniques, we ask how these might be used to address cases of a forensic nature, focusing on the use of microbial genome sequence analysis. Such analyses rely on the increasingly large numbers of microbial genomes present in public databases, the ability of individual investigators to rapidly sequence whole microbial genomes, and an increasing depth of understanding of their evolution and function. Suggestions are made as to how comparative microbial genomics might be applied forensically and may represent possibilities for the future development of forensic techniques. A particular emphasis is on the nascent field of genomic epidemiology, which utilizes rapid whole-genome sequencing to identify the source and spread of infectious outbreaks. Also discussed is the application of comparative microbial genomics to the study of historical epidemics and deaths and how the approaches developed may also be applicable to more recent and actionable cases.

Reconstructing the HIV-1 CRF02_AG and CRF06_cpx epidemics in Burkina Faso and West Africa using early samples

BACKGROUND

HIV-1 circulating recombinant forms (CRFs) represent viral recombinant lineages that play a significant role in the global epidemic. Two of them dominate the epidemic in Burkina Faso: CRF06_cpx (first described in this country) and CRF02_AG. We reconstructed the phylodynamics of both recombinant viruses in Burkina Faso and throughout West Africa.

METHODS

We analysed CRF06_cpx and CRF02_AG sequences (protease/gp41) from early samples collected in Burkina Faso in 1986 together with other GenBank sequences (1984-2013) in 4 datasets: African CRF06_cpx (210/60); down-sampled CRF06_cpx (146/45); Burkina Faso CRF02_AG (130/39) and West/Central African CRF02_AG (691/298). For each dataset, we analysed both protease and gp41 jointly using the BEAST multilocus analysis and conducted phylogeographic analysis to reconstruct the early migration routes between countries.

RESULTS

The time to the most recent common ancestor (tMRCA) of CRF06_cpx was 1979 (1973-1983) for protease and 1981 (1978-1983) for gp41. The gp41 analysis inferred the origin of CRF06_cpx (or at least its parental subtype G lineage) in the Democratic Republic of Congo but migrated to Burkina Faso soon after (1982). Both genes showed that CRF06_cpx radiated to the rest of West Africa predominantly after around 1990. These results were robust to the oversampling of Burkina Faso sequences as they were confirmed in the down-sampled dataset. The tMRCA of the Burkina Faso CRF02_AG lineage was 1979 (1977-1983) for protease and 1980 (1978-1981) for gp41. However, we reconstructed its presence in West Africa much earlier (mid-1960s), with an initial origin in Cameroon and/or Nigeria, and its phylogeographic analysis revealed much interconnection within the region with a lack of country-specific phylogenetic patterns, which prevents tracking its exact migration routes.

CONCLUSIONS

Burkina Faso presents a relatively young HIV epidemic, with the diversification of the current in-country CRF02_AG and CRF06_cpx lineages taking place around 1980. This country represents the main source of CRF06_cpx in West Africa. The CRF02_AG epidemic started at least a decade earlier and showed much interchange between West African countries (especially involving coastal countries) suggesting great population mobility and an extensive viral spread in the region.

First evidence of transmission of an HIV-1 M/O intergroup recombinant virus

OBJECTIVE

Despite the genetic divergence between HIV-1 groups M and O, HIV-1  M/O intergroup recombinants were reported. Actually, there is no data on the transmissibility of such recombinant forms. During a surveillance of HIV genetic diversity in Cameroon, we investigated the possible direct transmission of an HIV-1  M/O recombinant virus in an HIV-infected couple.

METHODS

Consecutive samples obtained from the couple were analysed for detection of dual HIV-1 groups M and O infections, and HIV-1  M/O recombinant forms. Analyses were performed using a serological and molecular algorithm based on HIV serotyping and group-specific PCRs targeting the polymerase and envelope genes. Pattern characterization of the strains found in both patients was based on complete genome sequencing. Phylogenetic and similarity profile analyses were performed to investigate the genetic relationship between viruses from both spouses and the previously described recombinant forms.

RESULTS

The sero-molecular algorithm data showed a group O serotype confirmed by molecular analysis in the envelope regions, whereas molecular tests identified HIV-1 group M in the polymerase. Phylogenetic analyses and similarity profiles of the full-length genome sequences showed that both spouses were infected with a unique recombinant virus having two recombination breakpoints in the vpr gene and LTR region. No phylogenetic link was found with the previous M/O recombinants.

CONCLUSION

We provide, for the first time, molecular evidence of direct transmission of an HIV-1  M/O recombinant, highlighting the potential spread of these divergent viruses. The importance of HIV-1 recombination on genetic evolution and public health when implying divergent strains as group O has to be carefully considered.

The Two-Phase Emergence of Non Pandemic HIV-1 Group O in Cameroon

Unlike the pandemic form of HIV-1 (group M), group O viruses are endemic in west central Africa, especially in Cameroon. However, little is known about group O's genetic evolution, and why this highly divergent lineage has not become pandemic. Using a unique and large set of group O sequences from samples collected from 1987 to 2012, we find that this lineage has evolved in successive slow and fast phases of diversification, with a most recent common ancestor estimated to have existed around 1930 (1914-1944). The most rapid periods of diversification occurred in the 1950s and in the 1980s, and could be linked to favourable epidemiological contexts in Cameroon. Group O genetic diversity reflects this two-phase evolution, with two distinct populations potentially having different viral properties. The currently predominant viral population emerged in the 1980s, from an ancient population which had first developed in the 1950s, and is characterized by higher growth and evolutionary rates, and the natural presence of the Y181C resistance mutation, thought to confer a phenotypic advantage. Our findings show that although this evolutionary pattern is specific to HIV-1 group O, it paralleled the early spread of HIV-1 group M in the Democratic Republic of Congo. Both viral lineages are likely to have benefited from similar epidemiological contexts. The relative role of virological and social factors in the distinct epidemic histories of HIV-1 group O and M needs to be reassessed.

Short Communication: HIV-1 Subtype B in the Dominican Republic: Evolution and Molecular Epidemiology

The Caribbean region has the world second highest incidence rate of acquired immunodeficiency syndrome. The island of Hispaniola is composed of two sovereign nations: the Dominican Republic and Haiti. Together, they account for more than 85% of HIV/AIDS cases in the Caribbean; and the Dominican Republic alone has approximately 46,000 (33,000-59,000) HIV-1-infected adults and children. Despite this, the magnitude of the genetic variability and evolution of the HIV-1 virus in the Dominican Republic is unclear. In the current study, we analyzed 195 reverse transcriptase (RT) sequences obtained from the Los Alamos HIV database. The data were used to assess the course of the viral epidemic over time in the Dominican Republic, using a coalescent approach. Based on the data, we estimated that the timing of the most recent common ancestor (tMRCA) of local HIV-1 subtype B emerged in 1963, approximately. In addition, the Bayesian analysis provided new information that suggests that the epidemic in the Dominican Republic experienced a significant decrease in relative genetic diversity in the past 2 decades. The results suggest that adherence to antiretroviral therapy, adequate prevention campaigns, and better access to health care may be altering the virus's evolution in the Dominican Republic.

Eight challenges in phylodynamic inference

The field of phylodynamics, which attempts to enhance our understanding of infectious disease dynamics using pathogen phylogenies, has made great strides in the past decade. Basic epidemiological and evolutionary models are now well characterized with inferential frameworks in place. However, significant challenges remain in extending phylodynamic inference to more complex systems. These challenges include accounting for evolutionary complexities such as changing mutation rates, selection, reassortment, and recombination, as well as epidemiological complexities such as stochastic population dynamics, host population structure, and different patterns at the within-host and between-host scales. An additional challenge exists in making efficient inferences from an ever increasing corpus of sequence data.

Molecular epidemiology, phylogeny and evolution of the filarial nematode Wuchereria bancrofti

Wuchereria bancrofti (Wb) is the most widely distributed of the three nematodes known to cause lymphatic filariasis (LF), the other two being Brugia malayi and Brugia timori. Current tools available to monitor LF are limited to diagnostic tests targeting DNA repeats, filarial antigens, and anti-filarial antibodies. While these tools are useful for detection and surveillance, elimination programs have yet to take full advantage of molecular typing for inferring infection history, strain fingerprinting, and evolution. To date, molecular typing approaches have included whole mitochondrial genomes, genotyping, targeted sequencing, and random amplified polymorphic DNA (RAPDs). These studies have revealed much about Wb biology. For example, in one study in Papua New Guinea researchers identified 5 major strains that were widespread and many minor strains some of which exhibit geographic stratification. Genome data, while rare, has been utilized to reconstruct evolutionary relationships among taxa of the Onchocercidae (the clade of filarial nematodes) and identify gene synteny. Their phylogeny reveals that speciation from the common ancestor of both B. malayi and Wb occurred around 5-6 millions years ago with shared ancestry to other filarial nematodes as recent as 15 million years ago. These discoveries hold promise for gene discovery and identifying drug targets in species that are more amenable to in vivo experiments. Continued technological developments in whole genome sequencing and data analysis will likely replace many other forms of molecular typing, multiplying the amount of data available on population structure, genetic diversity, and phylogenetics. Once widely available, the addition of population genetic data from genomic studies should hasten the elimination of LF parasites like Wb. Infectious disease control programs have benefited greatly from population genetics data and recently from population genomics data. However, while there is currently a surplus of data for diseases like malaria and HIV, there is a scarcity of this data for filarial nematodes. With the falling cost of genome sequencing, research on filarial nematodes could benefit from the addition of population genetics statistics and phylogenetics especially in dealing with elimination programs. A comprehensive review focusing on population genetics of filarial nematode does not yet exist. Here our goal is to provide a current overview of the molecular epidemiology of W. bancrofti (Wb) the primary causative agent of LF. We begin by reviewing studies utilizing molecular typing techniques with specific focus on genomic and population datasets. Next, we used whole mitochondrial genome data to construct a phylogeny and examine the evolutionary history of the Onchocercidae. Then, we provide a perspective to aid in understanding how population genetic techniques translate to modern epidemiology. Finally, we introduce the concept of genomic epidemiology and provide some examples that will aid in future studies of Wb.

Nef proteins of epidemic HIV-1 group O strains antagonize human tetherin

Most simian immunodeficiency viruses use their Nef protein to antagonize the host restriction factor tetherin. A deletion in human tetherin confers Nef resistance, representing a hurdle to successful zoonotic transmission. HIV-1 group M evolved to utilize the viral protein U (Vpu) to counteract tetherin. Although HIV-1 group O has spread epidemically in humans, it has not evolved a Vpu-based tetherin antagonism. Here we show that HIV-1 group O Nef targets a region adjacent to this deletion to inhibit transport of human tetherin to the cell surface, enhances virion release, and increases viral resistance to inhibition by interferon-α. The Nef protein of the inferred common ancestor of group O viruses is also active against human tetherin. Thus, Nef-mediated antagonism of human tetherin evolved prior to the spread of HIV-1 group O and likely facilitated secondary virus transmission. Our results may explain the epidemic spread of HIV-1 group O.

HIV epidemiology. The early spread and epidemic ignition of HIV-1 in human populations

Thirty years after the discovery of HIV-1, the early transmission, dissemination, and establishment of the virus in human populations remain unclear. Using statistical approaches applied to HIV-1 sequence data from central Africa, we show that from the 1920s Kinshasa (in what is now the Democratic Republic of Congo) was the focus of early transmission and the source of pre-1960 pandemic viruses elsewhere. Location and dating estimates were validated using the earliest HIV-1 archival sample, also from Kinshasa. The epidemic histories of HIV-1 group M and nonpandemic group O were similar until ~1960, after which group M underwent an epidemiological transition and outpaced regional population growth. Our results reconstruct the early dynamics of HIV-1 and emphasize the role of social changes and transport networks in the establishment of this virus in human populations.

HIV epidemiology. The early spread and epidemic ignition of HIV-1 in human populations

Thirty years after the discovery of HIV-1, the early transmission, dissemination, and establishment of the virus in human populations remain unclear. Using statistical approaches applied to HIV-1 sequence data from central Africa, we show that from the 1920s Kinshasa (in what is now the Democratic Republic of Congo) was the focus of early transmission and the source of pre-1960 pandemic viruses elsewhere. Location and dating estimates were validated using the earliest HIV-1 archival sample, also from Kinshasa. The epidemic histories of HIV-1 group M and nonpandemic group O were similar until ~1960, after which group M underwent an epidemiological transition and outpaced regional population growth. Our results reconstruct the early dynamics of HIV-1 and emphasize the role of social changes and transport networks in the establishment of this virus in human populations.

Phylogenetic studies of transmission dynamics in generalized HIV epidemics: an essential tool where the burden is greatest?

Efficient and effective HIV prevention measures for generalized epidemics in sub-Saharan Africa have not yet been validated at the population level. Design and impact evaluation of such measures requires fine-scale understanding of local HIV transmission dynamics. The novel tools of HIV phylogenetics and molecular epidemiology may elucidate these transmission dynamics. Such methods have been incorporated into studies of concentrated HIV epidemics to identify proximate and determinant traits associated with ongoing transmission. However, applying similar phylogenetic analyses to generalized epidemics, including the design and evaluation of prevention trials, presents additional challenges. Here we review the scope of these methods and present examples of their use in concentrated epidemics in the context of prevention. Next, we describe the current uses for phylogenetics in generalized epidemics and discuss their promise for elucidating transmission patterns and informing prevention trials. Finally, we review logistic and technical challenges inherent to large-scale molecular epidemiological studies of generalized epidemics and suggest potential solutions.

Non-M variants of human immunodeficiency virus type 1

The AIDS pandemic that started in the early 1980s is due to human immunodeficiency virus type 1 (HIV-1) group M (HIV-M), but apart from this major group, many divergent variants have been described (HIV-1 groups N, O, and P and HIV-2). The four HIV-1 groups arose from independent cross-species transmission of the simian immunodeficiency viruses (SIVs) SIVcpz, infecting chimpanzees, and SIVgor, infecting gorillas. This, together with human adaptation, accounts for their genomic, phylogenetic, and virological specificities. Nevertheless, the natural course of non-M HIV infection seems similar to that of HIV-M. The virological monitoring of infected patients is now possible with commercial kits, but their therapeutic management remains complex. All non-M variants were principally described for patients linked to Cameroon, where HIV-O accounts for 1% of all HIV infections; only 15 cases of HIV-N infection and 2 HIV-P infections have been reported. Despite improvements in our knowledge, many fascinating questions remain concerning the origin, genetic evolution, and slow spread of these variants. Other variants may already exist or may arise in the future, calling for close surveillance. This review provides a comprehensive, up-to-date summary of the current knowledge on these pathogens, including the historical background of their discovery; the latest advances in the comprehension of their origin and spread; and clinical, therapeutic, and laboratory aspects that may be useful for the management and the treatment of patients infected with these divergent viruses.

Within-host and between-host evolutionary rates across the HIV-1 genome

Background

HIV evolves rapidly at the epidemiological level but also at the within-host level. The virus’ within-host evolutionary rates have been argued to be much higher than its between-host evolutionary rates. However, this conclusion relies on analyses of a short portion of the virus envelope gene. Here, we study in detail these evolutionary rates across the HIV genome.

Results

We build phylogenies using a relaxed molecular clock assumption to estimate evolutionary rates in different regions of the HIV genome. We find that these rates vary strongly across the genome, with higher rates in the envelope gene (env). Within-host evolutionary rates are consistently higher than between-host rates throughout the HIV genome. This difference is significantly more pronounced in env. Finally, we find weak differences between overlapping and non-overlapping regions.

Conclusions

We provide a genome-wide overview of the differences in the HIV rates of molecular evolution at the within- and between-host levels. Contrary to hepatitis C virus, where differences are only located in the envelope gene, within-host evolutionary rates are higher than between-host evolutionary rates across the whole HIV genome. This supports the hypothesis that HIV strains that are less adapted to the host have an advantage during transmission. The most likely mechanism for this is storage and then preferential transmission of viruses in latent T-cells. These results shed a new light on the role of the transmission bottleneck in the evolutionary dynamics of HIV.

Origins of HIV and the AIDS pandemic

Acquired immunodeficiency syndrome (AIDS) of humans is caused by two lentiviruses, human immunodeficiency viruses types 1 and 2 (HIV-1 and HIV-2). Here, we describe the origins and evolution of these viruses, and the circumstances that led to the AIDS pandemic. Both HIVs are the result of multiple cross-species transmissions of simian immunodeficiency viruses (SIVs) naturally infecting African primates. Most of these transfers resulted in viruses that spread in humans to only a limited extent. However, one transmission event, involving SIVcpz from chimpanzees in southeastern Cameroon, gave rise to HIV-1 group M-the principal cause of the AIDS pandemic. We discuss how host restriction factors have shaped the emergence of new SIV zoonoses by imposing adaptive hurdles to cross-species transmission and/or secondary spread. We also show that AIDS has likely afflicted chimpanzees long before the emergence of HIV. Tracing the genetic changes that occurred as SIVs crossed from monkeys to apes and from apes to humans provides a new framework to examine the requirements of successful host switches and to gauge future zoonotic risk.

Viral phylodynamics

Viral phylodynamics is defined as the study of how epidemiological, immunological, and evolutionary processes act and potentially interact to shape viralphylogenies. Since the coining of the term in 2004, research on viral phylodynamics has focused on transmission dynamics in an effort to shed light on how these dynamics impact viral genetic variation. Transmission dynamics can be considered at the level of cells within an infected host, individual hosts within a population, or entire populations of hosts. Many viruses, especially RNA viruses, rapidly accumulate genetic variation because of short generation times and high mutation rates. Patterns of viral genetic variation are therefore heavily influenced by how quickly transmission occurs and by which entities transmit to one another. Patterns of viral genetic variation will also be affected by selection acting on viral phenotypes. Although viruses can differ with respect to many phenotypes, phylodynamic studies have to date tended to focus on a limited number of viral phenotypes. These include virulence phenotypes, phenotypes associated with viral transmissibility, cell or tissue tropism phenotypes, and antigenic phenotypes that can facilitate escape from host immunity. Due to the impact that transmission dynamics and selection can have on viral genetic variation, viral phylogenies can therefore be used to investigate important epidemiological, immunological, and evolutionary processes, such as epidemic spread[2], spatio-temporal dynamics including metapopulation dynamics[3], zoonotic transmission, tissue tropism[4], and antigenic drift[5]. The quantitative investigation of these processes through the consideration of viral phylogenies is the central aim of viral phylodynamics.

Estimating the rate of intersubtype recombination in early HIV-1 group M strains

West Central Africa has been implicated as the epicenter of the HIV-1 epidemic, and almost all group M subtypes can be found there. Previous analysis of early HIV-1 group M sequences from Kinshasa in the Democratic Republic of Congo, formerly Zaire, revealed that isolates from a number of individuals fall in different positions in phylogenetic trees constructed from sequences from opposite ends of the genome as a result of recombination between viruses of different subtypes. Here, we use discrete ancestral trait mapping to develop a procedure for quantifying HIV-1 group M intersubtype recombination across phylogenies, using individuals' gag (p17) and env (gp41) subtypes. The method was applied to previously described HIV-1 group M sequences from samples obtained in Kinshasa early in the global radiation of HIV. Nine different p17 and gp41 intersubtype recombinant combinations were present in the data set. The mean number of excess ancestral subtype transitions (NEST) required to map individuals' p17 subtypes onto the gp14 phylogeny samples, compared to the number required to map them onto the p17 phylogenies, and vice versa, indicated that excess subtype transitions occurred at a rate of approximately 7 × 10(-3) to 8 × 10(-3) per lineage per year as a result of intersubtype recombination. Our results imply that intersubtype recombination may have occurred in approximately 20% of lineages evolving over a period of 30 years and confirm intersubtype recombination as a substantial force in generating HIV-1 group M diversity.

Intracontinental spread of human invasive Salmonella Typhimurium pathovariants in sub-Saharan Africa

A highly invasive form of non-typhoidal Salmonella (iNTS) disease has recently been documented in many countries in sub-Saharan Africa. The most common Salmonella enterica serovar causing this disease is Typhimurium (Salmonella Typhimurium). We applied whole-genome sequence-based phylogenetic methods to define the population structure of sub-Saharan African invasive Salmonella Typhimurium isolates and compared these to global Salmonella Typhimurium populations. Notably, the vast majority of sub-Saharan invasive Salmonella Typhimurium isolates fell within two closely related, highly clustered phylogenetic lineages that we estimate emerged independently ∼52 and ∼35 years ago in close temporal association with the current HIV pandemic. Clonal replacement of isolates from lineage I by those from lineage II was potentially influenced by the use of chloramphenicol for the treatment of iNTS disease. Our analysis suggests that iNTS disease is in part an epidemic in sub-Saharan Africa caused by highly related Salmonella Typhimurium lineages that may have occupied new niches associated with a compromised human population and antibiotic treatment.

Enhanced heterosexual transmission hypothesis for the origin of pandemic HIV-1

HIV-1 M originated from SIVcpz endemic in chimpanzees from southeast Cameroon or neighboring areas, and it started to spread in the early 20th century. Here we examine the factors that may have contributed to simian-to-human transmission, local transmission between humans, and export to a city. The region had intense ape hunting, social disruption, commercial sex work, STDs, and traffic to/from Kinshasa in the period 1899-1923. Injection treatments increased sharply around 1930; however, their frequency among local patients was far lower than among modern groups experiencing parenteral HIV-1 outbreaks. Recent molecular datings of HIV-1 M fit better the period of maximal resource exploitation and trade links than the period of high injection intensity. We conclude that although local parenteral outbreaks might have occurred, these are unlikely to have caused massive transmission. World War I led to additional, and hitherto unrecognized, risks of HIV-1 emergence. We propose an Enhanced Heterosexual Transmission Hypothesis for the origin of HIV-1 M, featuring at the time and place of its origin a coincidence of favorable co-factors (ape hunting, social disruption, STDs, and mobility) for both cross-species transmission and heterosexual spread. Our hypothesis does not exclude a role for parenteral transmission in the initial viral adaptation.

Interplay between HIV-1 and Host Genetic Variation: A Snapshot into Its Impact on AIDS and Therapy Response

As of February 2012, 50 circulating recombinant forms (CRFs) have been reported for HIV-1 while one CRF for HIV-2. Also according to HIV sequence compendium 2011, the HIV sequence database is replete with 414,398 sequences. The fact that there are CRFs, which are an amalgamation of sequences derived from six or more subtypes (CRF27_cpx (cpx refers to complex) is a mosaic with sequences from 6 different subtypes besides an unclassified fragment), serves as a testimony to the continual divergent evolution of the virus with its approximate 1% per year rate of evolution, and this phenomena per se poses tremendous challenge for vaccine development against HIV/AIDS, a devastating disease that has killed 1.8 million patients in 2010. Here, we explore the interaction between HIV-1 and host genetic variation in the context of HIV/AIDS and antiretroviral therapy response.

Simian retroviruses in African apes

It is now well established that simian immunodeficiency viruses (SIVs) from chimpanzees (SIVcpz) and gorillas (SIVgor) from west Central Africa are at the origin of HIV-1/AIDS. Apes are also infected with other retroviruses, notably simian T-cell lymphotropic viruses (STLVs) and simian foamy viruses (SFVs), that can be transmitted to humans. We discuss the actual knowledge on SIV, STLV and SFV infections in chimpanzees, gorillas, and bonobos. We especially elaborate on how the recent development of non-invasive methods has allowed us to identify the reservoirs of the HIV-1 ancestors in chimpanzees and gorillas, and increased our knowledge of the natural history of SIV infections in chimpanzees. Multiple cross-species events with retroviruses from apes to humans have occurred, but only one transmission of SIVcpz from chimpanzees in south-eastern Cameroon spread worldwide, and is responsible for the actual HIV pandemic. Frequent SFV transmissions have been recently reported, but no human-to-human transmission has been documented yet. Because humans are still in contact with apes, identification of pathogens in wild ape populations can signal which pathogens may be cause risk for humans, and allow the development of serological and molecular assays with which to detect transmissions to humans. Finally, non-invasive sampling also allows the study of the impact of retroviruses and other pathogens on the health and survival of endangered species such as chimpanzees, gorillas, and bonobos.

Genomic analysis of hepatitis B virus reveals antigen state and genotype as sources of evolutionary rate variation

Hepatitis B virus (HBV) genomes are small, semi-double-stranded DNA circular genomes that contain alternating overlapping reading frames and replicate through an RNA intermediary phase. This complex biology has presented a challenge to estimating an evolutionary rate for HBV, leading to difficulties resolving the evolutionary and epidemiological history of the virus. Here, we re-examine rates of HBV evolution using a novel data set of 112 within-host, transmission history (pedigree) and among-host genomes isolated over 20 years from the indigenous peoples of the South Pacific, combined with 313 previously published HBV genomes. We employ Bayesian phylogenetic approaches to examine several potential causes and consequences of evolutionary rate variation in HBV. Our results reveal rate variation both between genotypes and across the genome, as well as strikingly slower rates when genomes are sampled in the Hepatitis B e antigen positive state, compared to the e antigen negative state. This Hepatitis B e antigen rate variation was found to be largely attributable to changes during the course of infection in the preCore and Core genes and their regulatory elements.

Genetic diversity and molecular epidemiology of HIV transmission

The high genetic diversity of HIV is one of its most significant features, as it has consequences in global distribution, vaccine design, therapy success, disease progression, transmissibility and viral load testing. Studying HIV diversity helps to understand its origins, migration patterns, current distribution and transmission events. New advances in sequencing technologies based on the parallel acquisition of data are now used to characterize within-host and population processes in depth. Additionally, we have seen similar advances in statistical methods designed to model the past history of lineages (the phylodynamic framework) to ultimately gain better insights into the evolutionary history of HIV. We can, for example, estimate population size changes, lineage dispersion over geographic areas and epidemiological parameters solely from sequence data. In this article, we review some of the evolutionary approaches used to study transmission patterns and processes in HIV and the insights gained from such studies.

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.

The origin and diversity of the HIV-1 pandemic

This review examines the enormous progress that has been made in the past decade in understanding the origin of HIV, HIV genetic variability, and the impact of global HIV diversity on the pandemic. Multiple zoonotic transmissions of simian immunodeficiency virus (SIV) have resulted in different HIV lineages in humans. In addition, the high mutation and recombination rates during viral replication result in a great genetic variability of HIV within individuals, as well as within populations, upon which evolutionary selection pressures act. The global HIV pandemic is examined in the context of HIV evolution, and the global diversity of HIV subtypes and recombinants is discussed in detail. Finally, the impact of HIV diversity on pathogenesis, transmission, diagnosis, treatment, the immune response, and vaccine development is reviewed.

Viral evolution in deep time: lentiviruses and mammals

Lentiviruses are a distinctive genus of retroviruses that cause chronic, persistent infections in mammals, including humans. The emergence of pandemic HIV type-1 (HIV-1) infection during the late 20th century shaped a view of lentiviruses as 'modern' viruses. However, recent research has revealed an entirely different perspective, elucidating aspects of an evolutionary relationship with mammals that extends across many millions of years. Such deep evolutionary history is likely to be typical of many host-virus systems, fundamentally underpinning their interactions in the present day. For this reason, establishing the deep history of virus and host interaction is key to developing a fully informed approach to tackling viral diseases. Here, I use the example of lentiviruses to illustrate how paleovirological, geographic and genetic calibrations allow observations of virus and host interaction across a wide range of temporal and spatial scales to be integrated into a coherent ecological and evolutionary framework.

HIV-1 Group P is unable to antagonize human tetherin by Vpu, Env or Nef

Background

A new subgroup of HIV-1, designated Group P, was recently detected in two unrelated patients of Cameroonian origin. HIV-1 Group P phylogenetically clusters with SIVgor suggesting that it is the result of a cross-species transmission from gorillas. Until today, HIV-1 Group P has only been detected in two patients, and its degree of adaptation to the human host is largely unknown. Previous data have shown that pandemic HIV-1 Group M, but not non-pandemic Group O or rare Group N viruses, efficiently antagonize the human orthologue of the restriction factor tetherin (BST-2, HM1.24, CD317) suggesting that primate lentiviruses may have to gain anti-tetherin activity for efficient spread in the human population. Thus far, three SIV/HIV gene products (vpu, nef and env) are known to have the potential to counteract primate tetherin proteins, often in a species-specific manner. Here, we examined how long Group P may have been circulating in humans and determined its capability to antagonize human tetherin as an indicator of adaptation to humans.

Results

Our data suggest that HIV-1 Group P entered the human population between 1845 and 1989. Vpu, Env and Nef proteins from both Group P viruses failed to counteract human or gorilla tetherin to promote efficient release of HIV-1 virions, although both Group P Nef proteins moderately downmodulated gorilla tetherin from the cell surface. Notably, Vpu, Env and Nef alleles from the two HIV-1 P strains were all able to reduce CD4 cell surface expression.

Conclusions

Our analyses of the two reported HIV-1 Group P viruses suggest that zoonosis occurred in the last 170 years and further support that pandemic HIV-1 Group M strains are better adapted to humans than non-pandemic or rare Group O, N and P viruses. The inability to antagonize human tetherin may potentially explain the limited spread of HIV-1 Group P in the human population.

Analysis of the human immunodeficiency virus type 1 M group Vpu domains involved in antagonizing tetherin

Zoonosis of chimpanzee simian immunodeficiency virus cpz to humans has given rise to both pandemic (M) and non-pandemic (O, N and P) groups of human immunodeficiency virus type-1 (HIV). These lentiviruses encode accessory proteins, including Vpu, which has been shown to reduce CD4 levels on the cell surface, as well as increase virion release from the cell by antagonizing tetherin (CD317, BST2). Here, we confirm that O group Vpus (Ca9 and BCF06) are unable to counteract tetherin or downregulate the protein from the cell surface, although they are still able to reduce cell-surface CD4 levels. We hypothesize that this inability to antagonize tetherin may have contributed to O group viruses failing to achieve pandemic levels of human-to-human transmission. Characterization of chimeric O/M group Vpus and Vpu mutants demonstrate that the Vpu-tetherin interaction is complex, involving several domains. We identify specific residues within the transmembrane proximal region that, along with the transmembrane domain, are crucial for tetherin counteraction and enhanced virion release. We have also shown that the critical domains are responsible for the localization of M group Vpu to the trans-Golgi network, where it relocalizes tetherin to counteract its function. This work sheds light on the acquisition of anti-tetherin activity and the molecular details of pandemic HIV infection in humans.

Phylogenetic and epidemic modeling of rapidly evolving infectious diseases

Epidemic modeling of infectious diseases has a long history in both theoretical and empirical research. However the recent explosion of genetic data has revealed the rapid rate of evolution that many populations of infectious agents undergo and has underscored the need to consider both evolutionary and ecological processes on the same time scale. Mathematical epidemiology has applied dynamical models to study infectious epidemics, but these models have tended not to exploit--or take into account--evolutionary changes and their effect on the ecological processes and population dynamics of the infectious agent. On the other hand, statistical phylogenetics has increasingly been applied to the study of infectious agents. This approach is based on phylogenetics, molecular clocks, genealogy-based population genetics and phylogeography. Bayesian Markov chain Monte Carlo and related computational tools have been the primary source of advances in these statistical phylogenetic approaches. Recently the first tentative steps have been taken to reconcile these two theoretical approaches. We survey the Bayesian phylogenetic approach to epidemic modeling of infection diseases and describe the contrasts it provides to mathematical epidemiology as well as emphasize the significance of the future unification of these two fields.

Reviewing the history of HIV-1: spread of subtype B in the Americas

The dispersal of HIV-1 subtype B (HIV-1B) is a reflection of the movement of human populations in response to social, political, and geographical issues. The initial dissemination of HIV-1B outside Africa seems to have included the passive involvement of human populations from the Caribbean in spreading the virus to the United States. However, the exact pathways taken during the establishment of the pandemic in the Americas remain unclear. Here, we propose a geographical scenario for the dissemination of HIV-1B in the Americas, based on phylogenetic and genetic statistical analyses of 313 available sequences of the pol gene from 27 countries. Maximum likelihood and Bayesian inference methods were used to explore the phylogenetic relationships between HIV-1B sequences, and molecular variance estimates were analyzed to infer the genetic structure of the viral population. We found that the initial dissemination and subsequent spread of subtype B in the Americas occurred via a single introduction event in the Caribbean around 1964 (1950–1967). Phylogenetic trees present evidence of several primary outbreaks in countries in South America, directly seeded by the Caribbean epidemic. Cuba is an exception insofar as its epidemic seems to have been introduced from South America. One clade comprising isolates from different countries emerged in the most-derived branches, reflecting the intense circulation of the virus throughout the American continents. Statistical analysis supports the genetic compartmentalization of the virus among the Americas, with a close relationship between the South American and Caribbean epidemics. These findings reflect the complex establishment of the HIV-1B pandemic and contribute to our understanding between the migration process of human populations and virus diffusion.

Molecular epidemiology of HIV type 1 CRF02_AG in Cameroon and African patients living in Italy

HIV-1 CRF02_AG accounts for >50% of infected individuals in Cameroon. CRF02_AG prevalence has been increasing both in Africa and Europe, particularly in Italy because of migrations from the sub-Saharan region. This study investigated the molecular epidemiology of CRF02_AG in Cameroon by employing Bayesian phylodynamics and analyzed the relationship between HIV-1 CRF02_AG isolates circulating in Italy and those prevalent in Africa to understand the link between the two epidemics. Among 291 Cameroonian reverse transcriptase sequences analyzed, about 70% clustered within three distinct clades, two of which shared a most recent common ancestor, all related to sequences from Western Africa. The major Cameroonian clades emerged during the mid-1970s and slowly spread during the next 30 years. Little or no geographic structure was detected within these clades. One of the major driving forces of the epidemic was likely the high accessibility between locations in Southern Cameroon contributing to the mobility of the population. The remaining Cameroonian sequences and the new strains isolated from Italian patients were interspersed mainly within West and Central African sequences in the tree, indicating a continuous exchange of CRF02_AG viral strains between Cameroon and other African countries, as well as multiple independent introductions in the Italian population. The evaluation of the spread of CRF02_AG may provide significant insight about the future dynamics of the Italian and European epidemic.

Dates of HIV infection can be estimated for seroprevalent patients by coalescent analysis of serial next-generation sequencing data

OBJECTIVE

To reconstruct dates of HIV infection by the coalescent analysis of longitudinal next-generation sequencing (NGS) data.

DESIGN

The coalescent predicts the time that has elapsed since the most recent common ancestor (MRCA) of a population. Because HIV tends to undergo severe bottlenecks upon transmission, the MRCA may be a good predictor of the time of infection. NGS provides an efficient means for performing large-scale clonal sequencing of HIV populations within patients, and the ideal raw material for coalescent analysis.

METHODS

Baseline and follow-up plasma samples were obtained from 19 individuals enrolled into the Montréal Primary HIV Infection cohort. Dates of infection were initially estimated at baseline from nongenetic data (clinical and serological markers and patient questionnaires). HIV RNA was extracted and seven regions of the genome were amplified, subjected to parallel-tagged 454 pyrosequencing, and analyzed using the software package BEAST.

RESULTS

Mean estimates of the time to the MRCA per patient were significantly correlated with nongenetic estimates (Spearman's ρ = 0.65, P = 4.4 × 10(-3)). The median absolute difference between coalescent and nongenetic date estimates was smallest (median 29.4 days) for highly variable regions of the HIV genome such as env V3, and greater (median 114.9 days) for more conserved regions such as pol.

CONCLUSION

This application of NGS represents an important advancement, not only because accurate estimates of dates of infection can be derived retrospectively from archived specimens, but also because each analysis is patient-specific and, therefore, robust to variation in rates of HIV evolution.

A possible geographic origin of endemic hepatitis C virus 6a in Hong Kong: evidences for the association with Vietnamese immigration

Background

Hepatitis C virus (HCV) 6a accounts for 23.6% of all HCV infections of the general population and 58.5% of intravenous drug users in Hong Kong. However, the geographical origin of this highly predominant HCV subgenotype is largely unknown. This study explores a hypothesis for one possible transmission route of HCV 6a to Hong Kong.

Methods

NS5A sequences derived from 26 HCV 6a samples were chosen from a five year period (1999–2004) from epidemiologically unrelated patients from Hong Kong. Partial-NS5A sequences (513-bp from nt 6728 to 7240) were adopted for Bayesian coalescent analysis to reconstruct the evolutionary history of HCV infections in Hong Kong using the BEAST v1.3 program. A rooted phylogenetic tree was drawn for these sequences by alignment with reference Vietnamese sequences. Demographic data were accessed from “The Statistic Yearbooks of Hong Kong”.

Results

Bayesian coalescent analysis showed that the rapid increase in 6a infections, which had increased more than 90-fold in Hong Kong from 1986 to 1994 correlated to two peaks of Vietnamese immigration to Hong Kong from 1978 to 1997. The second peak, which occurred from 1987 through 1997, overlapped with the rapid increase of HCV 6a occurrence in Hong Kong. Phylogenetic analyses have further revealed that HCV 6a strains from Vietnam may be ancestral to Hong Kong counterparts.

Conclusions

The high predominance of HCV 6a infections in Hong Kong was possibly associated with Vietnamese immigration during 1987–1997.

Inference for nonlinear epidemiological models using genealogies and time series

Phylodynamics - the field aiming to quantitatively integrate the ecological and evolutionary dynamics of rapidly evolving populations like those of RNA viruses - increasingly relies upon coalescent approaches to infer past population dynamics from reconstructed genealogies. As sequence data have become more abundant, these approaches are beginning to be used on populations undergoing rapid and rather complex dynamics. In such cases, the simple demographic models that current phylodynamic methods employ can be limiting. First, these models are not ideal for yielding biological insight into the processes that drive the dynamics of the populations of interest. Second, these models differ in form from mechanistic and often stochastic population dynamic models that are currently widely used when fitting models to time series data. As such, their use does not allow for both genealogical data and time series data to be considered in tandem when conducting inference. Here, we present a flexible statistical framework for phylodynamic inference that goes beyond these current limitations. The framework we present employs a recently developed method known as particle MCMC to fit stochastic, nonlinear mechanistic models for complex population dynamics to gene genealogies and time series data in a Bayesian framework. We demonstrate our approach using a nonlinear Susceptible-Infected-Recovered (SIR) model for the transmission dynamics of an infectious disease and show through simulations that it provides accurate estimates of past disease dynamics and key epidemiological parameters from genealogies with or without accompanying time series data.

Human immunodeficiency virus (HIV) immunopathogenesis and vaccine development: a review

The development of a safe, effective and globally affordable HIV vaccine offers the best hope for the future control of the HIV-1 pandemic. Since 1987, scores of candidate HIV-1 vaccines have been developed which elicited varying degrees of protective responses in nonhuman primate models, including DNA vaccines, subunit vaccines, live vectored recombinant vaccines and various prime-boost combinations. Four of these candidate vaccines have been tested for efficacy in human volunteers, but, to the exception of the recent RV144 Phase III trial in Thailand, which elicited a modest but statistically significant level of protection against infection, none has shown efficacy in preventing HIV-1 infection or in controlling virus replication and delaying progression of disease in humans. Protection against infection was observed in the RV144 trial, but intensive research is needed to try to understand the protective immune mechanisms at stake. Building-up on the results of the RV144 trial and deciphering what possibly are the immune correlates of protection are the top research priorities of the moment, which will certainly accelerate the development of an highly effective vaccine that could be used in conjunction with other HIV prevention and treatment strategies. This article reviews the state of the art of HIV vaccine development and discusses the formidable scientific challenges met in this endeavor, in the context of a better understanding of the immunopathogenesis of the disease.

Analysis of the origin and evolutionary history of HIV-1 CRF28_BF and CRF29_BF reveals a decreasing prevalence in the AIDS epidemic of Brazil

BACKGROUND

HIV-1 subtype B and subtype F are prevalent in the AIDS epidemic of Brazil. Recombinations between these subtypes have generated at least four BF circulating recombinant forms (CRFs). CRF28_BF and CRF29_BF are among the first two BF recombinants being identified in Brazil and they contributed significantly to the epidemic. However, the evolution and demographic histories of the CRFs are unclear.

METHODOLOGY/PRINCIPAL FINDINGS

A collection of gag and pol sequences sampled within Brazil was screened for CRF28_BF-like and CRF29_BF-like recombination patterns. A Bayesian coalescent framework was employed to delineate the phylogenetic, divergence time and population dynamics of the virus having CRF28_BF-like and CRF29_BF-like genotype. These recombinants were phylogenetically related to each other and formed a well-supported monophyletic clade dated to 1988-1989. The effective number of infections by these recombinants grew exponentially over a five-year period after their emergence, but then decreased toward the present following a logistic model of population growth. The demographic pattern of both recombinants closely resembles those previously reported for CRF31_BC.

CONCLUSIONS

We revealed that HIV-1 recombinants of the CRF28_BF/CRF29_BF clade are still circulating in the Brazilian population. These recombinants did not exhibit a strong founder effect and showed a decreasing prevalence in the AIDS epidemic of Brazil. Our data suggested that multiple URFs may also play a role in shaping the epidemic of recombinant BF HIV-1 in the region.