Publications
307 results found
Bonsall D, Golubchik T, Kosloff B, et al., 2018, HIV genotyping and phylogenetics in the HPTN 071 (PopART) study: validation of a high-throughput sequencing assay for viral load quantification, genotyping, resistance testing and high-resolution transmission networking, Publisher: JOHN WILEY & SONS LTD, Pages: 58-59
Raghwani J, Redd AD, Longosz AF, et al., 2018, Evolution of HIV-1 within untreated individuals and at the population scale in Uganda, PLOS PATHOGENS, Vol: 14, ISSN: 1553-7366
- Author Web Link
- Cite
- Citations: 16
Le Vu SOK, Ratmann O, Delpech V, et al., 2018, Comparison of cluster-based and source-attribution methods for estimating transmission risk using large HIV sequence databases, Epidemics, Vol: 23, Pages: 1-10, ISSN: 1755-4365
Phylogenetic clustering of HIV sequences from a random sample of patients can reveal epidemiological transmission patterns, but interpretation is hampered by limited theoretical support and statistical properties of clustering analysis remain poorly understood. Alternatively, source attribution methods allow fitting of HIV transmission models and thereby quantify aspects of disease transmission.A simulation study was conducted to assess error rates of clustering methods for detecting transmission risk factors. We modeled HIV epidemics among men having sex with men and generated phylogenies comparable to those that can be obtained from HIV surveillance data in the UK. Clustering and source attribution approaches were applied to evaluate their ability to identify patient attributes as transmission risk factors.We find that commonly used methods show a misleading association between cluster size or odds of clustering and covariates that are correlated with time since infection, regardless of their influence on transmission. Clustering methods usually have higher error rates and lower sensitivity than source attribution method for identifying transmission risk factors. But neither methods provide robust estimates of transmission risk ratios. Source attribution method can alleviate drawbacks from phylogenetic clustering but formal population genetic modeling may be required to estimate quantitative transmission risk factors.
Blanquart F, Lehtinen S, Lipsitch M, et al., 2018, The evolution of antibiotic resistance in a structured host population, JOURNAL OF THE ROYAL SOCIETY INTERFACE, Vol: 15, ISSN: 1742-5689
- Author Web Link
- Cite
- Citations: 31
van der Kuyl AC, Vink M, Zorgdrager F, et al., 2018, The evolution of subtype B HIV-1 <i>tat</i> in the Netherlands during 1985-2012, VIRUS RESEARCH, Vol: 250, Pages: 51-64, ISSN: 0168-1702
- Author Web Link
- Cite
- Citations: 6
Volz E, Le Vu S, Ratmann O, et al., 2018, Molecular epidemiology of HIV-1 subtype B reveals heterogeneous transmission risk: Implications for intervention and control, Journal of Infectious Diseases, Vol: 217, Pages: 1522-1529, ISSN: 0022-1899
BackgroundThe impact of HIV pre-exposure prophylaxis (PrEP) depends on infections averted by protecting vulnerable individuals as well as infections averted by preventing transmission by those who would have been infected if not receiving PrEP. Analysis of HIV phylogenies reveals risk factors for transmission, which we examine as potential criteria for allocating PrEP.MethodsWe analyzed 6912 HIV-1 partial pol sequences from men who have sex with men (MSM) in the United Kingdom combined with global reference sequences and patient-level metadata. Population genetic models were developed that adjust for stage of infection, global migration of HIV lineages, and changing incidence of infection through time. Models were extended to simulate the effects of providing susceptible MSM with PrEP.ResultsWe found that young age <25 years confers higher risk of HIV transmission (relative risk = 2.52 [95% confidence interval, 2.32–2.73]) and that young MSM are more likely to transmit to one another than expected by chance. Simulated interventions indicate that 4-fold more infections can be averted over 5 years by focusing PrEP on young MSM.ConclusionsConcentrating PrEP doses on young individuals can avert more infections than random allocation.
Stirrup OT, Dunn DT, Tostevin A, et al., 2018, Risk factors and outcomes for the Q151M and T69 insertion HIV-1 resistance mutations in historic UK data, AIDS Research and Therapy, Vol: 15, ISSN: 1742-6405
Background: The prevalence of HIV-1 resistance to antiretroviral therapies (ART) has declined in high-income countries over recent years, but drug resistance remains a substantial concern in many low and middle-income countries. The Q151M and T69 insertion (T69i) resistance mutations in the viral reverse transcriptase gene can reduce susceptibility to all nucleoside/tide analogue reverse transcriptase inhibitors, motivating the present study to investigate the risk factors and outcomes associated with these mutations. Methods: We considered all data in the UK HIV Drug Resistance Database for blood samples obtained in the period 1997-2014. Where available, treatment history and patient outcomes were obtained through linkage to the UK Collaborative HIV Cohort study. A matched case-control approach was used to assess risk factors associated with the appearance of each of the mutations in ART-experienced patients, and survival analysis was used to investigate factors associated with viral suppression. A further analysis using matched controls was performed to investigate the impact of each mutation on survival. Results: A total of 180 patients with Q151M mutation and 85 with T69i mutation were identified, almost entirely from before 2006. Occurrence of both the Q151M and T69i mutations was strongly associated with cumulative period of virological failure while on ART, and for Q151M there was a particular positive association with use of stavudine and negative association with use of boosted-protease inhibitors. Subsequent viral suppression was negatively associated with viral load at sequencing for both mutations, and for Q151M we found a negative association with didanosine use but a positive association with boosted-protease inhibitor use. The results obtained in these analyses were also consistent with potentially large associations with other drugs. Analyses were inconclusive regarding associations between the mutations and mortality, but mortality was high for pati
Croucher NJ, Apagyi KJ, Fraser C, 2018, Transformation asymmetry and the evolution of the bacterial accessory genome, Molecular Biology and Evolution, Vol: 35, Pages: 575-581, ISSN: 1537-1719
Bacterial transformation can insert or delete genomic islands (GIs), depending on the donor and recipient genotypes, if an homologous recombination spans the GI’s integration site and includes sufficiently long flanking homologous arms. Combining mathematical models of recombination with experiments using pneumococci found GI insertion rates declined geometrically with the GI’s size. The decrease in acquisition frequency with length (1.08x10−3 bp−1) was higher than a previous estimate of the analogous rate at which core genome recombinations terminated. Although most efficient for shorter GIs, transformation-mediated deletion frequencies did not vary consistently with GI length, with removal of 10 kb GIs approximately 50% as efficient as acquisition of base substitutions. Fragments of two kilobases, typical of transformation event sizes, could drive all these deletions independent of island length. The strong asymmetry of transformation, and its capacity to efficiently remove GIs, suggests non-mobile accessory loci will decline in frequency without preservation by selection.
Wymant C, Blanquart F, Golubchik T, et al., 2018, Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver, Virus Evolution, Vol: 4, ISSN: 2057-1577
Studying the evolution of viruses and their molecular epidemiology relies on accurate viral sequence data, so that small differences between similar viruses can be meaningfully interpreted. Despite its higher throughput and more detailed minority variant data, next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of large between- and within-host diversity, including frequent indels, may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to pre-process reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We used shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read whole-genome data produced with the Illumina platform, for sixty-five existing publicly available samples and fifty new samples. We show the systematic superiority of mapping to shiver's constructed reference compared with mapping the same reads to the closest of 3,249 real references: median values of 13 bases called differently and more accurately, 0 bases called differently and less accurately, and 205 bases of missing sequence recovered. We also successfully applied sh
Cornelissen M, Gall A, van der Kuyl A, et al., 2018, Workup of Human Blood Samples for Deep Sequencing of HIV-1 Genomes, VIRAL METAGENOMICS: METHODS AND PROTOCOLS, Vol: 1746, Pages: 55-61, ISSN: 1064-3745
Grandjean L, Gilman RH, Iwamoto T, et al., 2017, Convergent evolution and topologically disruptive polymorphisms among multidrug-resistant tuberculosis in Peru., PLoS ONE, Vol: 12, Pages: e0189838-e0189838, ISSN: 1932-6203
BACKGROUND: Multidrug-resistant tuberculosis poses a major threat to the success of tuberculosis control programs worldwide. Understanding how drug-resistant tuberculosis evolves can inform the development of new therapeutic and preventive strategies. METHODS: Here, we use novel genome-wide analysis techniques to identify polymorphisms that are associated with drug resistance, adaptive evolution and the structure of the phylogenetic tree. A total of 471 samples from different patients collected between 2009 and 2013 in the Lima suburbs of Callao and Lima South were sequenced on the Illumina MiSeq platform with 150bp paired-end reads. After alignment to the reference H37Rv genome, variants were called using standardized methodology. Genome-wide analysis was undertaken using custom written scripts implemented in R software. RESULTS: High quality homoplastic single nucleotide polymorphisms were observed in genes known to confer drug resistance as well as genes in the Mycobacterium tuberculosis ESX secreted protein pathway, pks12, and close to toxin/anti-toxin pairs. Correlation of homoplastic variant sites identified that many were significantly correlated, suggestive of epistasis. Variation in genes coding for ESX secreted proteins also significantly disrupted phylogenetic structure. Mutations in ESX genes in key antigenic epitope positions were also found to disrupt tree topology. CONCLUSION: Variation in these genes have a biologically plausible effect on immunogenicity and virulence. This makes functional characterization warranted to determine the effects of these polymorphisms on bacterial fitness and transmission.
Lehtinen S, Blanquart F, Lipsitch M, et al., 2017, On the evolutionary ecology of multidrug resistance in bacteria
<jats:title>Abstract</jats:title><jats:p>Resistance against different antibiotics appears on the same bacterial strains more often than expected by chance, leading to high frequencies of multidrug resistance. There are multiple explanations for this observation, but these tend to be specific to subsets of antibiotics and/or bacterial species, whereas the trend is pervasive. Here, we consider the question in terms of strain ecology: explaining why resistance to different antibiotics is often seen on the same strain requires an understanding of the competition between strains with different resistance profiles. This work builds on models originally proposed to explain another aspect of strain competition: the stable coexistence of antibiotic sensitivity and resistance observed in a number of bacterial species. We first demonstrate a partial structural similarity in these models of coexistence. We then generalise this unified underlying model to multidrug resistance and show that models with this structure predict high levels of association between resistance to different drugs and high multidrug resistance frequencies. We test predictions from this model in six bacterial datasets and find them to be qualitatively consistent with observed trends. The higher than expected frequencies of multidrug resistance are often interpreted as evidence that these strains are out-competing strains with lower resistance multiplicity. Our work provides an alternative explanation that is compatible with long-term stability in resistance frequencies.</jats:p><jats:sec><jats:title>Author summary</jats:title><jats:p>Antibiotic resistance is a serious public health concern, yet the ecology and evolution of drug resistance are not fully understood. This impacts our ability to design effective interventions to combat resistance. From a public health point of view, multidrug resistance is particularly problematic because resistance to different
Wymant C, Hall M, Ratmann O, et al., 2017, PHYLOSCANNER: inferring transmission from within- and between-host pathogen genetic diversity, Molecular Biology and Evolution, Vol: 35, Pages: 719-733, ISSN: 1537-1719
A central feature of pathogen genomics is that different infectious particles (virions, bacterial cells, etc.) within an infected individual may be genetically distinct, with patterns of relatedness amongst infectious particles being the result of both within-host evolution and transmission from one host to the next. Here we present a new software tool, phyloscanner, which analyses pathogen diversity from multiple infected hosts. phyloscanner provides unprecedented resolution into the transmission process, allowing inference of the direction of transmission from sequence data alone. Multiply infected individuals are also identified, as they harbour subpopulations of infectious particles that are not connected by within-host evolution, except where recombinant types emerge. Low-level contamination is flagged and removed. We illustrate phyloscanner on both viral and bacterial pathogens, namely HIV-1 sequenced on Illumina and Roche 454 platforms, HCV sequenced with the Oxford Nanopore MinION platform, and Streptococcus pneumoniae with sequences from multiple colonies per individual. phyloscanner is available from https://github.com/BDI-pathogens/phyloscanner.
Corander J, Fraser C, Gutmann MU, et al., 2017, Frequency-dependent selection in vaccine-associated pneumococcal population dynamics, Nature Ecology and Evolution, Vol: 1, Pages: 1950-1960, ISSN: 2397-334X
Many bacterial species are composed of multiple lineages distinguished by extensive variation in gene content. These often cocirculate in the same habitat, but the evolutionary and ecological processes that shape these complex populations are poorly understood. Addressing these questions is particularly important for Streptococcus pneumoniae, a nasopharyngeal commensal and respiratory pathogen, because the changes in population structure associated with the recent introduction of partial-coverage vaccines have substantially reduced pneumococcal disease. Here we show that pneumococcal lineages from multiple populations each have a distinct combination of intermediate-frequency genes. Functional analysis suggested that these loci may be subject to negative frequency-dependent selection (NFDS) through interactions with other bacteria, hosts or mobile elements. Correspondingly, these genes had similar frequencies in four populations with dissimilar lineage compositions. These frequencies were maintained following substantial alterations in lineage prevalences once vaccination programmes began. Fitting a multilocus NFDS model of post-vaccine population dynamics to three genomic datasets using Approximate Bayesian Computation generated reproducible estimates of the influence of NFDS on pneumococcal evolution, the strength of which varied between loci. Simulations replicated the stable frequency of lineages unperturbed by vaccination, patterns of serotype switching and clonal replacement. This framework highlights how bacterial ecology affects the impact of clinical interventions.Accessory loci are shown to have similar frequencies in diverse Streptococcus pneumoniae populations, suggesting negative frequency-dependent selection drives post-vaccination population restructuring.
Thomas RA, Burger R, Harper A, et al., 2017, Differences in health-related quality of life between HIV-positive and HIV-negative people in Zambia and South Africa: a cross-sectional baseline survey of the HPTN 071 (PopART) trial, The Lancet Global Health, Vol: 5, Pages: e1133-e1141, ISSN: 2214-109X
BackgroundThe life expectancy of HIV-positive individuals receiving antiretroviral therapy (ART) is approaching that of HIV-negative people. However, little is known about how these populations compare in terms of health-related quality of life (HRQoL). We aimed to compare HRQoL between HIV-positive and HIV-negative people in Zambia and South Africa.MethodsAs part of the HPTN 071 (PopART) study, data from adults aged 18–44 years were gathered between Nov 28, 2013, and March 31, 2015, in large cross-sectional surveys of random samples of the general population in 21 communities in Zambia and South Africa. HRQoL data were collected with a standardised generic measure of health across five domains. We used β-distributed multivariable models to analyse differences in HRQoL scores between HIV-negative and HIV-positive individuals who were unaware of their status; aware, but not in HIV care; in HIV care, but who had not initiated ART; on ART for less than 5 years; and on ART for 5 years or more. We included controls for sociodemographic variables, herpes simplex virus type-2 status, and recreational drug use.FindingsWe obtained data for 19 750 respondents in Zambia and 18 941 respondents in South Africa. Laboratory-confirmed HIV status was available for 19 330 respondents in Zambia and 18 004 respondents in South Africa; 4128 (21%) of these 19 330 respondents in Zambia and 4012 (22%) of 18 004 respondents in South Africa had laboratory-confirmed HIV. We obtained complete HRQoL information for 19 637 respondents in Zambia and 18 429 respondents in South Africa. HRQoL scores did not differ significantly between individuals who had initiated ART more than 5 years previously and HIV-negative individuals, neither in Zambia (change in mean score −0·002, 95% CI −0·01 to 0·001; p=0·219) nor in South Africa (0·000, −0·002 to 0·003; p=0·939). However, scores did differ between HIV-positive individu
Cobey S, Baskerville EB, Colijn C, et al., 2017, Host population structure and treatment frequency maintain balancing selection on drug resistance, JOURNAL OF THE ROYAL SOCIETY INTERFACE, Vol: 14, ISSN: 1742-5689
- Author Web Link
- Cite
- Citations: 23
Cornelissen M, Gall A, Vink M, et al., 2017, From clinical sample to complete genome: Comparing methods for the extraction of HIV-1 RNA for high-throughput deep sequencing, Virus Research, Vol: 239, Pages: 10-16, ISSN: 0168-1702
The BEEHIVE (Bridging the Evolution and Epidemiology of HIV in Europe) project aims to analyse nearly-complete viral genomes from >3000 HIV-1 infected Europeans using high-throughput deep sequencing techniques to investigate the virus genetic contribution to virulence. Following the development of a computational pipeline, including a new de novo assembler for RNA virus genomes, to generate larger contiguous sequences (contigs) from the abundance of short sequence reads that characterise the data, another area that determines genome sequencing success is the quality and quantity of the input RNA. A pilot experiment with 125 patient plasma samples was performed to investigate the optimal method for isolation of HIV-1 viral RNA for long amplicon genome sequencing. Manual isolation with the QIAamp Viral RNA Mini Kit (Qiagen) was superior over robotically extracted RNA using either the QIAcube robotic system, the mSample Preparation Systems RNA kit with automated extraction by the m2000sp system (Abbott Molecular), or the MagNA Pure 96 System in combination with the MagNA Pure 96 Instrument (Roche Diagnostics). We scored amplification of a set of four HIV-1 amplicons of ∼1.9, 3.6, 3.0 and 3.5 kb, and subsequent recovery of near-complete viral genomes.Subsequently, 616 BEEHIVE patient samples were analysed to determine factors that influence successful amplification of the genome in four overlapping amplicons using the QIAamp Viral RNA Kit for viral RNA isolation. Both low plasma viral load and high sample age (stored before 1999) negatively influenced the amplification of viral amplicons >3 kb. A plasma viral load of >100,000 copies/ml resulted in successful amplification of all four amplicons for 86% of the samples, this value dropped to only 46% for samples with viral loads of <20,000 copies/ml.
Blanquart F, Wymant C, Cornelissen M, et al., 2017, Correction: Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe., PLoS Biology, Vol: 15, ISSN: 1544-9173
[This corrects the article DOI: 10.1371/journal.pbio.2001855.].
Li LM, Grassly NC, Fraser C, 2017, Quantifying Transmission Heterogeneity Using Both Pathogen Phylogenies and Incidence Time Series., Molecular Biology and Evolution, Vol: 34, Pages: 2982-2995, ISSN: 1537-1719
Heterogeneity in individual-level transmissibility can be quantified by the dispersion parameter k of the offspring distribution. Quantifying heterogeneity is important as it affects other parameter estimates, it modulates the degree of unpredictability of an epidemic, and it needs to be accounted for in models of infection control. Aggregated data such as incidence time series are often not sufficiently informative to estimate k. Incorporating phylogenetic analysis can help to estimate k concurrently with other epidemiological parameters. We have developed an inference framework that uses particle Markov Chain Monte Carlo to estimate k and other epidemiological parameters using both incidence time series and the pathogen phylogeny. Using the framework to fit a modified compartmental transmission model that includes the parameter k to simulated data, we found that more accurate and less biased estimates of the reproductive number were obtained by combining epidemiological and phylogenetic analyses. However, k was most accurately estimated using pathogen phylogeny alone. Accurately estimating k was necessary for unbiased estimates of the reproductive number, but it did not affect the accuracy of reporting probability and epidemic start date estimates. We further demonstrated that inference was possible in the presence of phylogenetic uncertainty by sampling from the posterior distribution of phylogenies. Finally, we used the inference framework to estimate transmission parameters from epidemiological and genetic data collected during a poliovirus outbreak. Despite the large degree of phylogenetic uncertainty, we demonstrated that incorporating phylogenetic data in parameter inference improved the accuracy and precision of estimates.
Wymant CM, Hall M, Blanquart F, et al., 2017, Phylogenetics between and within hosts along the genome reveals transmission, dual infections, recombination and contamination, Publisher: JOHN WILEY & SONS LTD, Pages: 101-102
Wymant C, Hall M, Ratmann O, et al., 2017, PHYLOSCANNER: Inferring Transmission from Within‐ and Between-Host Pathogen Genetic Diversity
<jats:title>Abstract</jats:title><jats:p>A central feature of pathogen genomics is that different infectious particles (virions, bacterial cells, etc.) within an infected individual may be genetically distinct, with patterns of relatedness amongst infectious particles being the result of both within-host evolution and transmission from one host to the next. Here we present a new software tool, phyloscanner, which analyses pathogen diversity from multiple infected hosts. phyloscanner provides unprecedented resolution into the transmission process, allowing inference of the direction of transmission from sequence data alone. Multiply infected individuals are also identified, as they harbour subpopulations of infectious particles that are not connected by within-host evolution, except where recombinant types emerge. Low-level contamination is flagged and removed. We illustrate phyloscanner on both viral and bacterial pathogens, namely HIV-1 sequenced on Illumina and Roche 454 platforms, HCV sequenced with the Oxford Nanopore MinION platform, and <jats:italic>Streptococcus pneumoniae</jats:italic> with sequences from multiple colonies per individual. phyloscanner is available from <jats:underline><jats:ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="https://github.com/BDI-pathogens/phyloscanner">https://github.com/BDI-pathogens/phyloscanner</jats:ext-link></jats:underline>.</jats:p>
Ratmann O, Wymant C, Colijn C, et al., 2017, HIV-1 full-genome phylogenetics of generalized epidemics in sub-Saharan Africa: impact of missing nucleotide characters in next-generation sequences, Aids Research and Human Retroviruses, Vol: 33, Pages: 1083-1098, ISSN: 1931-8405
To characterize HIV-1 transmission dynamics in regions where the burden of HIV-1 is greatest, the “Phylogenetics and Networks for Generalised HIV Epidemics in Africa” consortium (PANGEA-HIV) is sequencing full-genome viral isolates from across sub-Saharan Africa. We report the first 3,985 PANGEA-HIV consensus sequences from four cohort sites (Rakai Community Cohort Study, n = 2,833; MRC/UVRI Uganda, n = 701; Mochudi Prevention Project, n = 359; Africa Health Research Institute Resistance Cohort, n = 92). Next-generation sequencing success rates varied: more than 80% of the viral genome from the gag to the nef genes could be determined for all sequences from South Africa, 75% of sequences from Mochudi, 60% of sequences from MRC/UVRI Uganda, and 22% of sequences from Rakai. Partial sequencing failure was primarily associated with low viral load, increased for amplicons closer to the 3′ end of the genome, was not associated with subtype diversity except HIV-1 subtype D, and remained significantly associated with sampling location after controlling for other factors. We assessed the impact of the missing data patterns in PANGEA-HIV sequences on phylogeny reconstruction in simulations. We found a threshold in terms of taxon sampling below which the patchy distribution of missing characters in next-generation sequences (NGS) has an excess negative impact on the accuracy of HIV-1 phylogeny reconstruction, which is attributable to tree reconstruction artifacts that accumulate when branches in viral trees are long. The large number of PANGEA-HIV sequences provides unprecedented opportunities for evaluating HIV-1 transmission dynamics across sub-Saharan Africa and identifying prevention opportunities. Molecular epidemiological analyses of these data must proceed cautiously because sequence sampling remains below the identified threshold and a considerable negative impact of missing characters on phyloge
Rutstein SE, Ananworanich J, Fidler S, et al., 2017, Clinical and public health implications of acute and early HIV detection and treatment: a scoping review, Journal of the International AIDS Society, Vol: 20, ISSN: 1758-2652
Introduction: The unchanged global HIV incidence may be related to ignoring acute HIV infection (AHI). This scoping review examines diagnostic, clinical, and public health implications of identifying and treating persons with AHI.Methods: We searched PubMed, in addition to hand-review of key journals identifying research pertaining to AHI detection and treatment. We focused on the relative contribution of AHI to transmission and the diagnostic, clinical, and public health implications. We prioritized research from low- and middle-income countries (LMICs) published in the last fifteen years.Results and Discussion: Extensive AHI research and limited routine AHI detection and treatment have begun in LMIC. Diagnostic challenges include ease-of-use, suitability for application and distribution in LMIC, and throughput for high-volume testing. Risk score algorithms have been used in LMIC to screen for AHI among individuals with behavioural and clinical characteristics more often associated with AHI. However, algorithms have not been implemented outside research settings. From a clinical perspective, there are substantial immunological and virological benefits to identifying and treating persons with AHI – evading the irreversible damage to host immune systems and seeding of viral reservoirs that occurs during untreated acute infection. The therapeutic benefits require rapid initiation of antiretrovirals, a logistical challenge in the absence of point-of-care testing. From a public health perspective, AHI diagnosis and treatment is critical to: decrease transmission via viral load reduction and behavioural interventions; improve pre-exposure prophylaxis outcomes by avoiding treatment initiation for HIV-seronegative persons with AHI; and, enhance partner services via notification for persons recently exposed or likely transmitting.Conclusions: There are undeniable clinical and public health benefits to AHI detection and treatment, but also substantial diagnostic and lo
Blanquart F, Wymant C, Cornelissen M, et al., 2017, Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe, PLoS Biology, Vol: 15, ISSN: 1544-9173
HIV-1 set-point viral load—the approximately stable value of viraemia in the first years of chronic infection—is a strong predictor of clinical outcome and is highly variable across infected individuals. To better understand HIV-1 pathogenesis and the evolution of the viral population, we must quantify the heritability of set-point viral load, which is the fraction of variation in this phenotype attributable to viral genetic variation. However, current estimates of heritability vary widely, from 6% to 59%. Here we used a dataset of 2,028 seroconverters infected between 1985 and 2013 from 5 European countries (Belgium, Switzerland, France, the Netherlands and the United Kingdom) and estimated the heritability of set-point viral load at 31% (CI 15%–43%). Specifically, heritability was measured using models of character evolution describing how viral load evolves on the phylogeny of whole-genome viral sequences. In contrast to previous studies, (i) we measured viral loads using standardized assays on a sample collected in a strict time window of 6 to 24 months after infection, from which the viral genome was also sequenced; (ii) we compared 2 models of character evolution, the classical “Brownian motion” model and another model (“Ornstein–Uhlenbeck”) that includes stabilising selection on viral load; (iii) we controlled for covariates, including age and sex, which may inflate estimates of heritability; and (iv) we developed a goodness of fit test based on the correlation of viral loads in cherries of the phylogenetic tree, showing that both models of character evolution fit the data well. An overall heritability of 31% (CI 15%–43%) is consistent with other studies based on regression of viral load in donor–recipient pairs. Thus, about a third of variation in HIV-1 virulence is attributable to viral genetic variation.
Mostowy RJ, Croucher NJ, De Maio N, et al., 2017, Pneumococcal capsule synthesis locus cps as evolutionary hotspot with potential to generate novel serotypes by recombination, Molecular Biology and Evolution, Vol: 34, Pages: 2537-2554, ISSN: 1537-1719
Diversity of the polysaccharide capsule in Streptococcus pneumoniae -- main surface antigen and the target of the currently used pneumococcal vaccines -- constitutes a major obstacle in eliminating pneumococcal disease. Such diversity is genetically encoded by almost 100 variants of the capsule biosynthesis locus, cps. However, the evolutionary dynamics of the capsule remains not fully understood. Here, using genetic data from 4,519 bacterial isolates, we found cps to be an evolutionary hotspot with elevated substitution and recombination rates. These rates were a consequence of relaxed purifying selection and positive, diversifying selection acting at this locus, supporting the hypothesis that the capsule has an increased potential to generate novel diversity compared to the rest of the genome. Diversifying selection was particularly evident in the region of wzd/wze genes, which are known to regulate capsule expression and hence the bacterium's ability to cause disease. Using a novel, capsule-centred approach, we analysed the evolutionary history of twelve major serogroups. Such analysis revealed their complex diversification scenarios, which were principally driven by recombination with other serogroups and other streptococci. Patterns of recombinational exchanges between serogroups could not be explained by serotype frequency alone, thus pointing to non-random associations between co-colonising serotypes. Finally, we discovered a previously unobserved mosaic serotype 39X, which was confirmed to carry a viable and structurally novel capsule. Adding to previous discoveries of other mosaic capsules in densely sampled collections, these results emphasise the strong adaptive potential of the bacterium by its ability to generate novel antigenic diversity by recombination.
Blanquart F, Lehtinen S, Fraser C, 2017, An evolutionary model to predict the frequency of antibiotic resistance under seasonal antibiotic use, and an application to <i>Streptococcus pneumoniae</i>, PROCEEDINGS OF THE ROYAL SOCIETY B-BIOLOGICAL SCIENCES, Vol: 284, ISSN: 0962-8452
- Author Web Link
- Cite
- Citations: 20
Fraser C, Li LM, 2017, Coalescent models for populations with time-varying population sizes and arbitrary offspring distributions
<jats:title>Abstract</jats:title><jats:p>The coalescent has been used to infer from gene genealogies the population dynamics of biological systems, such as the prevalence of an infectious disease. The offspring distribution affects the relationship between population dynamics and the genealogy, and for infectious diseases, the offspring distribution is often highly overdispersed. Here, we provide a general formula for the coalescent rate for populations with time-varying sizes and any offspring distribution. The formula is valid in the same large population limit as Kingman’s original derivation. By relating our derivation to existing formulations of the coalescent, we show that differences in the coalescent rate derived for many population models may be explained by differences in the offspring distribution. The coalescent derivations presented here could be used to quantify the overdispersion in the offspring distribution of infectious diseases, which is useful for accurate modelling disease outbreaks.</jats:p>
Dudas G, Carvalho LM, Bedford T, et al., 2017, Virus genomes reveal factors that spread and sustained the Ebola epidemic, Nature, Vol: 544, Pages: 309-315, ISSN: 0028-0836
The 2013–2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic ‘gravity’ model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.
Cobey S, Baskerville EB, Colijn C, et al., 2017, Host population structure and treatment frequency maintain balancing selection on drug resistance
<jats:title>Abstract</jats:title><jats:p>It is a truism that antimicrobial drugs select for resistance, but explaining pathogen- and population-specific variation in patterns of resistance remains an open problem. Like other common commensals, <jats:italic>Streptococcus pneumoniae</jats:italic> has demonstrated persistent coexistence of drug-sensitive and drug-resistant strains. Theoretically, this outcome is unlikely. We modeled the dynamics of competing strains of <jats:italic>S. pneumoniae</jats:italic> to investigate the impact of transmission dynamics and treatment-induced selective pressures on the probability of stable coexistence. We find that the outcome of competition is extremely sensitive to structure in the host population, although coexistence can arise from age-assortative transmission models with age-varying rates of antibiotic use. Moreover, we find that the selective pressure from antibiotics arises not so much from the rate of antibiotic use per se but from the frequency of treatment: frequent antibiotic therapy disproportionately impacts the fitness of sensitive strains. This same phenomenon explains why serotypes with longer durations of carriage tend to be more resistant. These dynamics may apply to other potentially pathogenic, microbial commensals and highlight how population structure, which is often omitted from models, can have a large impact.</jats:p>
Garske T, Cori A, Ariyarajah A, et al., 2017, Heterogeneities in the case fatality ratio in the West African Ebola outbreak 2013 – 2016, Philosophical Transactions of the Royal Society B: Biological Sciences, Vol: 372, ISSN: 1471-2970
The 2013–2016 Ebola outbreak in West Africa is the largest on record with 28 616 confirmed, probable and suspected cases and 11 310 deaths officially recorded by 10 June 2016, the true burden probably considerably higher. The case fatality ratio (CFR: proportion of cases that are fatal) is a key indicator of disease severity useful for gauging the appropriate public health response and for evaluating treatment benefits, if estimated accurately. We analysed individual-level clinical outcome data from Guinea, Liberia and Sierra Leone officially reported to the World Health Organization. The overall mean CFR was 62.9% (95% CI: 61.9% to 64.0%) among confirmed cases with recorded clinical outcomes. Age was the most important modifier of survival probabilities, but country, stage of the epidemic and whether patients were hospitalized also played roles. We developed a statistical analysis to detect outliers in CFR between districts of residence and treatment centres (TCs), adjusting for known factors influencing survival and identified eight districts and three TCs with a CFR significantly different from the average. From the current dataset, we cannot determine whether the observed variation in CFR seen by district or treatment centre reflects real differences in survival, related to the quality of care or other factors or was caused by differences in reporting practices or case ascertainment.
This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.