Experimental evolution of Staphylococcus aureus in macrophages: dissection of a conditional adaptive trait promoting intracellular survival.

Staphylococcus aureus is a major pathogen associated with important diseases in humans and animals. Macrophages are a key component of the innate immune response to S. aureus infection and play a major role in disease outcomes. To investigate the adaptive evolution of S. aureus in response to macrophages, we developed an experimental infection assay. S. aureus strains representing major human epidemic clones were passaged many times in a macrophage cell line, accumulating mutations in an array of genomic loci. Phenotypic analysis revealed the emergence of a lineage exhibiting increased survival in macrophages and human blood, and resistance to vancomycin. The evolved lineage exhibited a previously undescribed small colony variant (SCV) phenotype characterized by hyper-pigmentation, which resulted from a missense mutation in rsbW. Notably, the novel SCV was a conditional adaptive trait that was unstable in nutrient-replete conditions in vitro, rapidly converting from hyper-pigmented SCV to a non-pigmented large colony variant via spontaneous sigB deletion events. Importantly, we identified similar deletions in the genome sequences of a limited number of clinical S. aureus isolates from public databases, indicating that related events may occur during clinical infection. Experimental infection of zebrafish did not reveal a difference in virulence between parent and novel SCV but demonstrated an in vivo fitness cost for the compensatory sigB deletion events. Taken together, we report an experimental evolutionary approach for investigating bacterial innate immune cell interactions, revealing a conditional adaptation that promotes S. aureus survival in macrophages and resistance to vancomycin. IMPORTANCE: Staphylococcus aureus is an important human bacterial pathogen. The host response to S. aureus involves the production of innate immune cells such as macrophages which are important for fighting infection. Here we report a new model of experimental evolution for studying how S. aureus can evade killing by macrophages. We identified a novel adaptive phenotype that promotes survival in macrophages and blood and resistance to antibiotics. The phenotype is lost rapidly upon growth in nutrient-rich conditions via disruption of the alternative sigma factor sigB, revealing a conditional niche-specific fitness advantage. Genomic analysis of clinical isolates suggests similar adaptations may occur during human infections. Our model may be used broadly to identify adaptations of S. aureus to the innate immune response.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock.

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

A graph-based approach for the visualisation and analysis of bacterial pangenomes.

BACKGROUND: The advent of low cost, high throughput DNA sequencing has led to the availability of thousands of complete genome sequences for a wide variety of bacterial species. Examining and interpreting genetic variation on this scale represents a significant challenge to existing methods of data analysis and visualisation. RESULTS: Starting with the output of standard pangenome analysis tools, we describe the generation and analysis of interactive, 3D network graphs to explore the structure of bacterial populations, the distribution of genes across a population, and the syntenic order in which those genes occur, in the new open-source network analysis platform, Graphia. Both the analysis and the visualisation are scalable to datasets of thousands of genome sequences. CONCLUSIONS: We anticipate that the approaches presented here will be of great utility to the microbial research community, allowing faster, more intuitive, and flexible interaction with pangenome datasets, thereby enhancing interpretation of these complex data.

SARS-CoV-2 Omicron is an immune escape variant with an altered cell entry pathway.

Vaccines based on the spike protein of SARS-CoV-2 are a cornerstone of the public health response to COVID-19. The emergence of hypermutated, increasingly transmissible variants of concern (VOCs) threaten this strategy. Omicron (B.1.1.529), the fifth VOC to be described, harbours multiple amino acid mutations in spike, half of which lie within the receptor-binding domain. Here we demonstrate substantial evasion of neutralization by Omicron BA.1 and BA.2 variants in vitro using sera from individuals vaccinated with ChAdOx1, BNT162b2 and mRNA-1273. These data were mirrored by a substantial reduction in real-world vaccine effectiveness that was partially restored by booster vaccination. The Omicron variants BA.1 and BA.2 did not induce cell syncytia in vitro and favoured a TMPRSS2-independent endosomal entry pathway, these phenotypes mapping to distinct regions of the spike protein. Impaired cell fusion was determined by the receptor-binding domain, while endosomal entry mapped to the S2 domain. Such marked changes in antigenicity and replicative biology may underlie the rapid global spread and altered pathogenicity of the Omicron variant.

Evolutionary and Functional Analysis of Coagulase Positivity among the Staphylococci.

The bacterial genus Staphylococcus comprises a large group of pathogenic and nonpathogenic species associated with an array of host species. Staphylococci are differentiated into coagulase-positive or coagulase-negative groups based on the capacity to promote clotting of plasma, a phenotype historically associated with the ability to cause disease. However, the genetic basis of this important diagnostic and pathogenic trait across the genus has not been examined to date. Here, we selected 54 representative staphylococcal species and subspecies to examine coagulation of plasma derived from six representative host species. In total, 13 staphylococcal species mediated coagulation of plasma from at least one host species including one previously identified as coagulase negative (Staphylococcus condimenti). Comparative genomic analysis revealed that coagulase activity correlated with the presence of a gene (vwb) encoding the von Willebrand binding protein (vWbp) whereas only the Staphylococcus aureus complex contained a gene encoding staphylocoagulase (Coa), the classical mediator of coagulation. Importantly, S. aureus retained vwb-dependent coagulase activity in an S. aureus strain deleted for coa whereas deletion of vwb in Staphylococcus pseudintermedius resulted in loss of coagulase activity. Whole-genome-based phylogenetic reconstruction of the Staphylococcus genus revealed that the vwb gene has been acquired on at least four different occasions during the evolution of the Staphylococcus genus followed by allelic diversification via mutation and recombination. Allelic variants of vWbp from selected coagulase-positive staphylococci mediated coagulation in a host-dependent manner indicative of host-adaptive evolution. Taken together, we have determined the genetic and evolutionary basis of staphylococcal coagulation, revealing vWbp to be its archetypal determinant. IMPORTANCE The ability of some species of staphylococci to promote coagulation of plasma is a key pathogenic and diagnostic trait. Here, we provide a comprehensive analysis of the coagulase positivity of the staphylococci and its evolutionary genetic basis. We demonstrate that the von Willebrand binding protein rather than staphylocoagulase is the archetypal coagulation factor of the staphylococci and that the vwb gene has been acquired several times independently during the evolution of the staphylococci. Subsequently, vwb has undergone adaptive diversification to facilitate host-specific functionality. Our findings provide important insights into the evolution of pathogenicity among the staphylococci and the genetic basis for a defining diagnostic phenotype.

Phylogenetic Networks and Parameters Inferred from HIV Nucleotide Sequences of High-Risk and General Population Groups in Uganda: Implications for Epidemic Control.

Phylogenetic inference is useful in characterising HIV transmission networks and assessing where prevention is likely to have the greatest impact. However, estimating parameters that influence the network structure is still scarce, but important in evaluating determinants of HIV spread. We analyzed 2017 HIV pol sequences (728 Lake Victoria fisherfolk communities (FFCs), 592 female sex workers (FSWs) and 697 general population (GP)) to identify transmission networks on Maximum Likelihood (ML) phylogenetic trees and refined them using time-resolved phylogenies. Network generative models were fitted to the observed degree distributions and network parameters, and corrected Akaike Information Criteria and Bayesian Information Criteria values were estimated. 347 (17.2%) HIV sequences were linked on ML trees (maximum genetic distance ≤4.5%, ≥95% bootstrap support) and, of these, 303 (86.7%) that consisted of pure A1 (n = 168) and D (n = 135) subtypes were analyzed in BEAST v1.8.4. The majority of networks (at least 40%) were found at a time depth of ≤5 years. The waring and yule models fitted best networks of FFCs and FSWs respectively while the negative binomial model fitted best networks in the GP. The network structure in the HIV-hyperendemic FFCs is likely to be scale-free and shaped by preferential attachment, in contrast to the GP. The findings support the targeting of interventions for FFCs in a timely manner for effective epidemic control. Interventions ought to be tailored according to the dynamics of the HIV epidemic in the target population and understanding the network structure is critical in ensuring the success of HIV prevention programs.

Radical genome remodelling accompanied the emergence of a novel host-restricted bacterial pathogen.

The emergence of new pathogens is a major threat to public and veterinary health. Changes in bacterial habitat such as a switch in host or disease tropism are typically accompanied by genetic diversification. Staphylococcus aureus is a multi-host bacterial species associated with human and livestock infections. A microaerophilic subspecies, Staphylococcus aureus subsp. anaerobius, is responsible for Morel's disease, a lymphadenitis restricted to sheep and goats. However, the evolutionary history of S. aureus subsp. anaerobius and its relatedness to S. aureus are unknown. Population genomic analyses of clinical S. aureus subsp. anaerobius isolates revealed a highly conserved clone that descended from a S. aureus progenitor about 1000 years ago before differentiating into distinct lineages that contain African and European isolates. S. aureus subsp. anaerobius has undergone limited clonal expansion, with a restricted population size, and an evolutionary rate 10-fold slower than S. aureus. The transition to its current restricted ecological niche involved acquisition of a pathogenicity island encoding a ruminant host-specific effector of abscess formation, large chromosomal re-arrangements, and the accumulation of at least 205 pseudogenes, resulting in a highly fastidious metabolism. Importantly, expansion of ~87 insertion sequences (IS) located largely in intergenic regions provided distinct mechanisms for the control of expression of flanking genes, including a novel mechanism associated with IS-mediated anti-anti-sense decoupling of ancestral gene repression. Our findings reveal the remarkable evolutionary trajectory of a host-restricted bacterial pathogen that resulted from extensive remodelling of the S. aureus genome through an array of diverse mechanisms in parallel.

The Molecular Epidemiology and Transmission Dynamics of HIV Type 1 in a General Population Cohort in Uganda.

The General Population Cohort (GPC) in south-western Uganda has a low HIV-1 incidence rate (<1%). However, new infections continue to emerge. In this research, 3796 HIV-1 pol sequences (GPC: n = 1418, non-GPC sites: n = 1223, Central Uganda: n = 1010 and Eastern Uganda: n = 145) generated between 2003-2015 were analysed using phylogenetic methods with demographic data to understand HIV-1 transmission in this cohort and inform the epidemic response. HIV-1 subtype A1 was the most prevalent strain in the GPC area (GPC and non-GPC sites) (39.8%), central (45.9%) and eastern (52.4%) Uganda. However, in the GPC alone, subtype D was the predominant subtype (39.1%). Of the 524 transmission clusters identified by Cluster Picker, all large clusters (≥5 individuals, n = 8) involved individuals from the GPC. In a multivariate analysis, clustering was strongly associated with being female (adjusted Odds Ratio, aOR = 1.28; 95% CI, 1.06-1.54), being >25 years (aOR = 1.52; 95% CI, 1.16-2.0) and being a resident in the GPC (aOR = 6.90; 95% CI, 5.22-9.21). Phylogeographic analysis showed significant viral dissemination (Bayes Factor test, BF > 3) from the GPC without significant viral introductions (BF < 3) into the GPC. The findings suggest localized HIV-1 transmission in the GPC. Intensifying geographically focused combination interventions in the GPC would contribute towards controlling HIV-1 infections.

Population Analysis of Staphylococcus aureus Reveals a Cryptic, Highly Prevalent Superantigen SElW That Contributes to the Pathogenesis of Bacteremia.

Staphylococcal superantigens (SAgs) are a family of secreted toxins that stimulate T cell activation and are associated with an array of diseases in humans and livestock. Most SAgs produced by Staphylococcus aureus are encoded by mobile genetic elements, such as pathogenicity islands, bacteriophages, and plasmids, in a strain-dependent manner. Here, we carried out a population genomic analysis of >800 staphylococcal isolates representing the breadth of S. aureus diversity to investigate the distribution of all 26 identified SAg genes. Up to 14 SAg genes were identified per isolate with the most common gene selw (encoding a putative SAg, SElW) identified in 97% of isolates. Most isolates (62.5%) have a full-length open reading frame of selw with an alternative TTG start codon that may have precluded functional characterization of SElW to date. Here, we demonstrate that S. aureus uses the TTG start codon to translate a potent SAg SElW that induces Vß-specific T cell proliferation, a defining feature of classical SAgs. SElW is the only SAg predicted to be expressed by isolates of the CC398 lineage, an important human and livestock epidemic clone. Deletion of selw in a representative CC398 clinical isolate, S. aureus NM001, resulted in complete loss of T cell mitogenicity in vitro, and in vivo expression of SElW by S. aureus increased the bacterial load in the liver during bloodstream infection of SAg-sensitive HLA-DR4 transgenic mice. Overall, we report the characterization of a novel, highly prevalent, and potent SAg that contributes to the pathogenesis of S. aureus infection.IMPORTANCEStaphylococcus aureus is an important human and animal pathogen associated with an array of diseases, including life-threatening necrotizing pneumonia and infective endocarditis. The success of S. aureus as a pathogen has been linked in part to its ability to manipulate the host immune response through the secretion of toxins and immune evasion molecules. The staphylococcal superantigens (SAgs) have been studied for decades, but their role in S. aureus pathogenesis is not well understood, and an appreciation for how SAgs manipulate the host immune response to promote infection may be crucial for the development of novel intervention strategies. Here, we characterized a widely prevalent, previously cryptic, staphylococcal SAg, SElW, that contributes to the severity of S. aureus infections caused by an important epidemic clone of S. aureus CC398. Our findings add to the understanding of staphylococcal SAg diversity and function and provide new insights into the capacity of S. aureus to cause disease.

Acapsular Staphylococcus aureus with a non-functional agr regains capsule expression after passage through the bloodstream in a bacteremia mouse model.

Selection pressures exerted on Staphylococcus aureus by host factors during infection may lead to the emergence of regulatory phenotypes better adapted to the infection site. Traits convenient for persistence may be fixed by mutation thus turning these mutants into microevolution endpoints. The feasibility that stable, non-encapsulated S. aureus mutants can regain expression of key virulence factors for survival in the bloodstream was investigated. S. aureus agr mutant HU-14 (IS256 insertion in agrC) from a patient with chronic osteomyelitis was passed through the bloodstream using a bacteriemia mouse model and derivative P3.1 was obtained. Although IS256 remained inserted in agrC, P3.1 regained production of capsular polysaccharide type 5 (CP5) and staphyloxanthin. Furthermore, P3.1 expressed higher levels of asp23/SigB when compared with parental strain HU-14. Strain P3.1 displayed decreased osteoclastogenesis capacity, thus indicating decreased adaptability to bone compared with strain HU-14 and exhibited a trend to be more virulent than parental strain HU-14. Strain P3.1 exhibited the loss of one IS256 copy, which was originally located in the HU-14 noncoding region between dnaG (DNA primase) and rpoD (sigA). This loss may be associated with the observed phenotype change but the mechanism remains unknown. In conclusion, S. aureus organisms that escape the infected bone may recover the expression of key virulence factors through a rapid microevolution pathway involving SigB regulation of key virulence factors.

Pervasive and non-random recombination in near full-length HIV genomes from Uganda.

Recombination is an important feature of HIV evolution, occurring both within and between the major branches of diversity (subtypes). The Ugandan epidemic is primarily composed of two subtypes, A1 and D, that have been co-circulating for 50 years, frequently recombining in dually infected patients. Here, we investigate the frequency of recombinants in this population and the location of breakpoints along the genome. As part of the PANGEA-HIV consortium, 1,472 consensus genome sequences over 5 kb have been obtained from 1,857 samples collected by the MRC/UVRI & LSHTM Research unit in Uganda, 465 (31.6 per cent) of which were near full-length sequences (>8 kb). Using the subtyping tool SCUEAL, we find that of the near full-length dataset, 233 (50.1 per cent) genomes contained only one subtype, 30.8 per cent A1 (n = 143), 17.6 per cent D (n = 82), and 1.7 per cent C (n = 8), while 49.9 per cent (n = 232) contained more than one subtype (including A1/D (n = 164), A1/C (n = 13), C/D (n = 9); A1/C/D (n = 13), and 33 complex types). K-means clustering of the recombinant A1/D genomes revealed a section of envelope (C2gp120-TMgp41) is often inherited intact, whilst a generalized linear model was used to demonstrate significantly fewer breakpoints in the gag-pol and envelope C2-TM regions compared with accessory gene regions. Despite similar recombination patterns in many recombinants, no clearly supported circulating recombinant form (CRF) was found, there was limited evidence of the transmission of breakpoints, and the vast majority (153/164; 93 per cent) of the A1/D recombinants appear to be unique recombinant forms. Thus, recombination is pervasive with clear biases in breakpoint location, but CRFs are not a significant feature, characteristic of a complex, and diverse epidemic.

Metagenomic sequencing of clinical samples reveals a single widespread clone of Lawsonia intracellularis responsible for porcine proliferative enteropathy.

Lawsonia intracellularis is a Gram-negative obligate intracellular bacterium that is the aetiological agent of proliferative enteropathy (PE), a common intestinal disease of major economic importance in pigs and other animal species. To date, progress in understanding the biology of L. intracellularis for improved disease control has been hampered by the inability to culture the organism in vitro. In particular, our understanding of the genomic diversity and population structure of clinical L. intercellularis is very limited. Here, we utilized a metagenomic shotgun approach to directly sequence and assemble 21 L. intracellularis genomes from faecal and ileum samples of infected pigs and horses across three continents. Phylogenetic analysis revealed a genetically monomorphic clonal lineage responsible for infections in pigs, with distinct subtypes associated with infections in horses. The genome was highly conserved, with 94 % of genes shared by all isolates and a very small accessory genome made up of only 84 genes across all sequenced strains. In part, the accessory genome was represented by regions with a high density of SNPs, indicative of recombination events importing novel gene alleles. In summary, our analysis provides the first view of the population structure for L. intracellularis, revealing a single major lineage associated with disease of pigs. The limited diversity and broad geographical distribution suggest the recent emergence and clonal expansion of an important livestock pathogen.

Phylogeography of HIV-1 suggests that Ugandan fishing communities are a sink for, not a source of, virus from general populations.

Although fishing communities (FCs) in Uganda are disproportionately affected by HIV-1 relative to the general population (GP), the transmission dynamics are not completely understood. We earlier found most HIV-1 transmissions to occur within FCs of Lake Victoria. Here, we test the hypothesis that HIV-1 transmission in FCs is isolated from networks in the GP. We used phylogeography to reconstruct the geospatial viral migration patterns in 8 FCs and 2 GP cohorts and a Bayesian phylogenetic inference in BEAST v1.8.4 to analyse the temporal dynamics of HIV-1 transmission. Subtype A1 (pol region) was most prevalent in the FCs (115, 45.1%) and GP (177, 50.4%). More recent HIV transmission pairs from FCs were found at a genetic distance (GD) <1.5% than in the GP (Fisher's exact test, p = 0.001). The mean time depth for pairs was shorter in FCs (5 months) than in the GP (4 years). Phylogeographic analysis showed strong support for viral migration from the GP to FCs without evidence of substantial viral dissemination to the GP. This suggests that FCs are a sink for, not a source of, virus strains from the GP. Targeted interventions in FCs should be extended to include the neighbouring GP for effective epidemic control.

A high HIV-1 strain variability in London, UK, revealed by full-genome analysis: Results from the ICONIC project.

BACKGROUND & METHODS: The ICONIC project has developed an automated high-throughput pipeline to generate HIV nearly full-length genomes (NFLG, i.e. from gag to nef) from next-generation sequencing (NGS) data. The pipeline was applied to 420 HIV samples collected at University College London Hospitals NHS Trust and Barts Health NHS Trust (London) and sequenced using an Illumina MiSeq at the Wellcome Trust Sanger Institute (Cambridge). Consensus genomes were generated and subtyped using COMET, and unique recombinants were studied with jpHMM and SimPlot. Maximum-likelihood phylogenetic trees were constructed using RAxML to identify transmission networks using the Cluster Picker. RESULTS: The pipeline generated sequences of at least 1Kb of length (median = 7.46Kb, IQR = 4.01Kb) for 375 out of the 420 samples (89%), with 174 (46.4%) being NFLG. A total of 365 sequences (169 of them NFLG) corresponded to unique subjects and were included in the down-stream analyses. The most frequent HIV subtypes were B (n = 149, 40.8%) and C (n = 77, 21.1%) and the circulating recombinant form CRF02_AG (n = 32, 8.8%). We found 14 different CRFs (n = 66, 18.1%) and multiple URFs (n = 32, 8.8%) that involved recombination between 12 different subtypes/CRFs. The most frequent URFs were B/CRF01_AE (4 cases) and A1/D, B/C, and B/CRF02_AG (3 cases each). Most URFs (19/26, 73%) lacked breakpoints in the PR+RT pol region, rendering them undetectable if only that was sequenced. Twelve (37.5%) of the URFs could have emerged within the UK, whereas the rest were probably imported from sub-Saharan Africa, South East Asia and South America. For 2 URFs we found highly similar pol sequences circulating in the UK. We detected 31 phylogenetic clusters using the full dataset: 25 pairs (mostly subtypes B and C), 4 triplets and 2 quadruplets. Some of these were not consistent across different genes due to inter- and intra-subtype recombination. Clusters involved 70 sequences, 19.2% of the dataset. CONCLUSIONS: The initial analysis of genome sequences detected substantial hidden variability in the London HIV epidemic. Analysing full genome sequences, as opposed to only PR+RT, identified previously undetected recombinants. It provided a more reliable description of CRFs (that would be otherwise misclassified) and transmission clusters.

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain: A phylogeographic approach.

Phylogenetic studies are a valuable tool to understand viral transmission patterns and the role of immigration in HIV-1 spread. We analyzed the spatio-temporal relationship of different HIV-1 non-B subtype variants over time using phylogenetic analysis techniques. We collected 693 pol (PR+RT) sequences that were sampled from 2005 to 2012 from naïve patients in different hospitals in southern Spain. We used REGA v3.0 to classify them into subtypes and recombinant forms, which were confirmed by phylogenetic analysis through maximum likelihood (ML) using RAxML. For the main HIV-1 non-B variants, publicly available, genetically similar sequences were sought using HIV-BLAST. The presence of HIV-1 lineages circulating in our study population was established using ML and Bayesian inference (BEAST v1.7.5) and transmission networks were identified. We detected 165 (23.4%) patients infected with HIV-1 non-B variants: 104 (63%) with recombinant viruses in pol: CRF02_AG (71, 43%), CRF14_BG (8, 4.8%), CRF06_cpx (5, 3%) and nine other recombinant forms (11, 6.7%) and unique recombinants (9, 5.5%). The rest (61, 37%) were infected with non-recombinant subtypes: A1 (30, 18.2%), C (7, [4.2%]), D (3, [1.8%]), F1 (9, 5.5%) and G (12, 7.3%). Most patients infected with HIV-1 non-B variants were men (63%, p < 0.001) aged over 35 (73.5%, p < 0.001), heterosexuals (92.2%, p < 0.001), from Africa (59.5%, p < 0.001) and living in the El Ejido area (62.4%, p<0.001). We found lineages of epidemiological relevance (mainly within Subtype A1), imported primarily through female sex workers from East Europe. We detected 11 transmission clusters of HIV-1 non-B Subtypes, which included patients born in Spain in half of them. We present the phylogenetic profiles of the HIV-1 non-B variants detected in southern Spain, and explore their putative geographical origins. Our data reveals a high HIV-1 genetic diversity likely due to the import of viral lineages that circulate in other countries. The highly immigrated El Ejido area acts as a gateway through which different subtypes are introduced into other regions, hence the importance of setting up epidemiological control measures to prevent future outbreaks.

Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic.

HIV molecular epidemiology studies analyse viral pol gene sequences due to their availability, but whole genome sequencing allows to use other genes. We aimed to determine what gene(s) provide(s) the best approximation to the real phylogeny by analysing a simulated epidemic (created as part of the PANGEA_HIV project) with a known transmission tree. We sub-sampled a simulated dataset of 4662 sequences into different combinations of genes (gag-pol-env, gag-pol, gag, pol, env and partial pol) and sampling depths (100%, 60%, 20% and 5%), generating 100 replicates for each case. We built maximum-likelihood trees for each combination using RAxML (GTR + Γ), and compared their topologies to the corresponding true tree's using CompareTree. The accuracy of the trees was significantly proportional to the length of the sequences used, with the gag-pol-env datasets showing the best performance and gag and partial pol sequences showing the worst. The lowest sampling depths (20% and 5%) greatly reduced the accuracy of tree reconstruction and showed high variability among replicates, especially when using the shortest gene datasets. In conclusion, using longer sequences derived from nearly whole genomes will improve the reliability of phylogenetic reconstruction. With low sample coverage, results can be highly variable, particularly when based on short sequences.

Phylodynamic and Phylogeographic Profiles of Subtype B HIV-1 Epidemics in South Spain.

BACKGROUND: Since 1982, HIV-1 epidemics have evolved to different scenarios in terms of transmission routes, subtype distribution and characteristics of transmission clusters. We investigated the evolutionary history of HIV-1 subtype B in south Spain. PATIENTS & METHODS: We studied all newly diagnosed HIV-1 subtype B patients in East Andalusia during the 2005-2012 period. For the analysis, we used the reverse transcriptase and protease sequences from baseline resistance, and the Trugene® HIV Genotyping kit (Siemens, Barcelona, Spain). Subtyping was done with REGA v3.0. The maximum likelihood trees constructed with RAxML were used to study HIV-1 clustering. Phylogeographic and phylodynamic profiles were studied by Bayesian inference methods with BEAST v1.7.5 and SPREAD v1.0.6. RESULTS: Of the 493 patients infected with HIV-1 subtype B, 234 grouped into 55 clusters, most of which were small (44 clusters ≤ 5 patients, 31 with 2 patients, 13 with 3). The rest (133/234) were grouped into 11 clusters with ≥ 5 patients, and most (82%, 109/133) were men who have sex with men (MSM) grouped into 8 clusters. The association with clusters was more frequent in Spanish (p = 0.02) men (p< 0.001), MSM (p<0.001) younger than 35 years (p = 0.001) and with a CD4+ T-cell count above 350 cells/ul (p<0.001). We estimated the date of HIV-1 subtype B regional epidemic diversification around 1970 (95% CI: 1965-1987), with an evolutionary rate of 2.4 (95%CI: 1.7-3.1) x 10-3 substitutions/site/year. Most clusters originated in the 1990s in MSMs. We observed exponential subtype B HIV-1 growth in 1980-1990 and 2005-2008. The most significant migration routes for subtype B went from inland cities to seaside locations. CONCLUSIONS: We provide the first data on the phylodynamic and phylogeographic profiles of HIV-1 subtype B in south Spain. Our findings of transmission clustering among MSMs should alert healthcare managers to enhance preventive measures in this risk group in order to prevent future outbreaks.

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.

Analysis of the history and spread of HIV-1 in Uganda using phylodynamics.

HIV prevalence has decreased in Uganda since the 1990s, but remains substantial within high-risk groups. Here, we reconstruct the history and spread of HIV subtypes A1 and D in Uganda and explore the transmission dynamics in high-risk populations. We analysed HIV pol sequences from female sex workers in Kampala (n = 42), Lake Victoria fisher-folk (n = 46) and a rural clinical cohort (n = 74), together with publicly available sequences from adjacent regions in Uganda (n = 412) and newly generated sequences from samples taken in Kampala in 1986 (n = 12). Of the sequences from the three Ugandan populations, 60 (37.1 %) were classified as subtype D, 54 (33.3 %) as subtype A1, 31 (19.1 %) as A1/D recombinants, six (3.7 %) as subtype C, one (0.6 %) as subtype G and 10 (6.2 %) as other recombinants. Among the A1/D recombinants we identified a new candidate circulating recombinant form. Phylodynamic and phylogeographic analyses using BEAST indicated that the Ugandan epidemics originated in 1960 (1950-1968) for subtype A1 and 1973 (1970-1977) for D, in rural south-western Uganda with subsequent spread to Kampala. They also showed extensive interconnection with adjacent countries. The sequence analysis shows both epidemics grew exponentially during the 1970s-1980s and decreased from 1992, which agrees with HIV prevalence reports in Uganda. Inclusion of sequences from the 1980s indicated the origin of both epidemics was more recent than expected and substantially narrowed the confidence intervals in comparison to previous estimates. We identified three transmission clusters and ten pairs, none of them including patients from different populations, suggesting active transmission within a structured transmission network.