Evaluation of phylogenetic reconstruction methods using bacterial whole genomes: a simulation based study

Background: Phylogenetic reconstruction is a necessary first step in many analyses which use whole genome sequence data from bacterial populations. There are many available methods to infer phylogenies, and these have various advantages and disadvantages, but few unbiased comparisons of the range of approaches have been made. Methods: We simulated data from a defined "true tree" using a realistic evolutionary model. We built phylogenies from this data using a range of methods, and compared reconstructed trees to the true tree using two measures, noting the computational time needed for different phylogenetic reconstructions. We also used real data from Streptococcus pneumoniae alignments to compare individual core gene trees to a core genome tree. Results: We found that, as expected, maximum likelihood trees from good quality alignments were the most accurate, but also the most computationally intensive. Using less accurate phylogenetic reconstruction methods, we were able to obtain results of comparable accuracy; we found that approximate results can rapidly be obtained using genetic distance based methods. In real data we found that highly conserved core genes, such as those involved in translation, gave an inaccurate tree topology, whereas genes involved in recombination events gave inaccurate branch lengths. We also show a tree-of-trees, relating the results of different phylogenetic reconstructions to each other. Conclusions: We recommend three approaches, depending on requirements for accuracy and computational time. Quicker approaches that do not perform full maximum likelihood optimisation may be useful for many analyses requiring a phylogeny, as generating a high quality input alignment is likely to be the major limiting factor of accurate tree topology. We have publicly released our simulated data and code to enable further comparisons

Dominant Role of Nucleotide Substitution in the Diversification of Serotype 3 Pneumococci over Decades and during a Single Infection

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex’s entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that upregulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype

Genomic analysis and comparison of two gonorrhoea outbreaks

© 2016 Didelot et al.Gonorrhea is a sexually transmitted disease causing growing concern, with a substantial increase in reported incidence over the past few years in the United Kingdom and rising levels of resistance to a wide range of antibiotics. Understanding its epidemiology is therefore of major biomedical importance, not only on a population scale but also at the level of direct transmission. However, the molecular typing techniques traditionally used for gonorrhea infections do not provide sufficient resolution to investigate such fine-scale patterns. Here we sequenced the genomes of 237 isolates from two local collections of isolates from Sheffield and London, each of which was resolved into a single type using traditional methods. The two data sets were selected to have different epidemiological properties: the Sheffield data were collected over 6 years from a predominantly heterosexual population, whereas the London data were gathered within half a year and strongly associated with men who have sex with men. Based on contact tracing information between individuals in Sheffield, we found that transmission is associated with a median time to most recent common ancestor of 3.4 months, with an upper bound of 8 months, which we used as a criterion to identify likely transmission links in both data sets. In London, we found that transmission happened predominantly between individuals of similar age, sexual orientation, and location and also with the same HIV serostatus, which may reflect serosorting and associated risk behaviors. Comparison of the two data sets suggests that the London epidemic involved about ten times more cases than the Sheffield outbreak. IMPORTANCE: The recent increases in gonorrhea incidence and antibiotic resistance are cause for public health concern. Successful intervention requires a better understanding of transmission patterns, which is not uncovered by traditional molecular epidemiology techniques. Here we studied two outbreaks that took place in Sheffield and London, United Kingdom. We show that whole-genome sequencing provides the resolution to investigate direct gonorrhea transmission between infected individuals. Combining genome sequencing with rich epidemiological information about infected individuals reveals the importance of several transmission routes and risk factors, which can be used to design better control measures

Microevolution of extensively drug-resistant tuberculosis in Russia.

Extensively drug-resistant (XDR) tuberculosis (TB), which is resistant to both first- and second-line antibiotics, is an escalating problem, particularly in the Russian Federation. Molecular fingerprinting of 2348 Mycobacterium tuberculosis isolates collected in Samara Oblast, Russia, revealed that 72%belonged to the Beijing lineage, a genotype associated with enhanced acquisition of drug resistance and increased virulence. Whole-genome sequencing of 34 Samaran isolates, plus 25 isolates representing global M. tuberculosis complex diversity, revealed that Beijing isolates originating in Eastern Europe formed a monophyletic group. Homoplasic polymorphisms within this clade were almost invariably associated with antibiotic resistance, indicating that the evolution of this population is primarily driven by drug therapy. Resistance genotypes showed a strong correlation with drug susceptibility phenotypes. A novel homoplasic mutation in rpoC, found only in isolates carrying a common rpoB rifampicin-resistance mutation, may play a role in fitness compensation. Most multidrug-resistant (MDR) isolates also had mutations in the promoter of a virulence gene, eis, which increase its expression and confer kanamycin resistance. Kanamycin therapy may thus select for mutants with increased virulence, helping preserve bacterial fitness and promoting transmission of drug-resistant TB strains. The East European clade was dominated by two MDR clusters, each disseminated across Samara. Polymorphisms conferring fluoroquinolone resistance were independently acquired multiple times within each cluster, indicating that XDR TB is currently not widely transmitted. © 2012 by Cold Spring Harbor Laboratory Press

Recombination in Streptococcus pneumoniae Lineages Increase with Carriage Duration and Size of the Polysaccharide Capsule.

Streptococcus pneumoniae causes a high burden of invasive pneumococcal disease (IPD) globally, especially in children from resource-poor settings. Like many bacteria, the pneumococcus can import DNA from other strains or even species by transformation and homologous recombination, which has allowed the pneumococcus to evade clinical interventions such as antibiotics and pneumococcal conjugate vaccines (PCVs). Pneumococci are enclosed in a complex polysaccharide capsule that determines the serotype; the capsule varies in size and is associated with properties including carriage prevalence and virulence. We determined and quantified the association between capsule and recombination events using genomic data from a diverse collection of serotypes sampled in Malawi. We determined both the amount of variation introduced by recombination relative to mutation (the relative rate) and how many individual recombination events occur per isolate (the frequency). Using univariate analyses, we found an association between both recombination measures and multiple factors associated with the capsule, including duration and prevalence of carriage. Because many capsular factors are correlated, we used multivariate analysis to correct for collinearity. Capsule size and carriage duration remained positively associated with recombination, although with a reduced P value, and this effect may be mediated through some unassayed additional property associated with larger capsules. This work describes an important impact of serotype on recombination that has been previously overlooked. While the details of how this effect is achieved remain to be determined, it may have important consequences for the serotype-specific response to vaccines and other interventions. IMPORTANCE: The capsule determines >90 different pneumococcal serotypes, which vary in capsule size, virulence, duration, and prevalence of carriage. Current serotype-specific vaccines elicit anticapsule antibodies. Pneumococcus can take up exogenous DNA by transformation and insert it into its chromosome by homologous recombination. This mechanism has disseminated drug resistance and generated vaccine escape variants. It is hence crucial to pneumococcal evolutionary response to interventions, but there has been no systematic study quantifying whether serotypes vary in recombination and whether this is associated with serotype-specific properties such as capsule size or carriage duration. Larger capsules could physically inhibit DNA uptake, or given the longer carriage duration for larger capsules, this may promote recombination. We find that recombination varies among capsules and is associated with capsule size, carriage duration, and carriage prevalence and negatively associated with invasiveness. The consequence of this work is that serotypes with different capsules may respond differently to selective pressures like vaccines

High-throughput, quantitative analyses of genetic interactions in E. coli.

Large-scale genetic interaction studies provide the basis for defining gene function and pathway architecture. Recent advances in the ability to generate double mutants en masse in Saccharomyces cerevisiae have dramatically accelerated the acquisition of genetic interaction information and the biological inferences that follow. Here we describe a method based on F factor-driven conjugation, which allows for high-throughput generation of double mutants in Escherichia coli. This method, termed genetic interaction analysis technology for E. coli (GIANT-coli), permits us to systematically generate and array double-mutant cells on solid media in high-density arrays. We show that colony size provides a robust and quantitative output of cellular fitness and that GIANT-coli can recapitulate known synthetic interactions and identify previously unidentified negative (synthetic sickness or lethality) and positive (suppressive or epistatic) relationships. Finally, we describe a complementary strategy for genome-wide suppressor-mutant identification. Together, these methods permit rapid, large-scale genetic interaction studies in E. coli

Designing a physical activity parenting course : parental views on recruitment, content and delivery

Background Many children do not engage in sufficient levels of physical activity (PA) and spend too much time screen-viewing (SV). High levels of SV (e.g. watching TV, playing video games and surfing the internet) and low levels of PA have been associated with adverse health outcomes. Parenting courses may hold promise as an intervention medium to change children’s PA and SV. The current study was formative work conducted to design a new parenting programme to increase children’s PA and reduce their SV. Specifically, we focussed on interest in a course, desired content and delivery style, barriers and facilitators to participation and opinions on control group provision. Methods In-depth telephone interviews were conducted with thirty two parents (29 female) of 6–8 year olds. Data were analysed thematically. An anonymous online survey was also completed by 750 parents of 6–8 year old children and descriptive statistics calculated. Results Interview participants were interested in a parenting course because they wanted general parenting advice and ideas to help their children be physically active. Parents indicated that they would benefit from knowing how to quantify their child’s PA and SV levels. Parents wanted practical ideas of alternatives to SV. Most parents would be unable to attend unless childcare was provided. Schools were perceived to be a trusted source of information about parenting courses and the optimal recruitment location. In terms of delivery style, the majority of parents stated they would prefer a group-based approach that provided opportunities for peer learning and support with professional input. Survey participants reported the timing of classes and the provision of childcare were essential factors that would affect participation. In terms of designing an intervention, the most preferred control group option was the opportunity to attend the same course at a later date. Conclusions Parents are interested in PA/SV parenting courses but the provision of child care is essential for attendance. Recruitment is likely to be facilitated via trusted sources. Parents want practical advice on how to overcome barriers and suggest advice is provided in a mutually supportive group experience with expert input

Access to primary care for socio-economically disadvantaged older people in rural areas: a qualitative study

Objective: We aim to explore the barriers to accessing primary care for socio-economically disadvantaged older people in rural areas. Methods: Using a community recruitment strategy, fifteen people over 65 years, living in a rural area, and receiving financial support were recruited for semi-structured interviews. Four focus groups were held with rural health professionals. Interviews and focus groups were audio-recorded and transcribed. Thematic analysis was used to identify barriers to primary care access. Findings: Older people’s experience can be understood within the context of a patient perceived set of unwritten rules or social contract – an individual is careful not to bother the doctor in return for additional goodwill when they become unwell. However, most found it difficult to access primary care due to engaged telephone lines, availability of appointments, interactions with receptionists; breaching their perceived social contract. This left some feeling unwelcome, worthless or marginalised, especially those with high expectations of the social contract or limited resources, skills and/or desire to adapt to service changes Health professionals’ described how rising demands and expectations coupled with service constraints had necessitated service development, such as fewer home visits, more telephone consultations, triaging calls and modifying the appointment system. Conclusion: Multiple barriers to accessing primary care exist for this group. As primary care is re-organised to reduce costs, commissioners and practitioners must not lose sight of the perceived social contract and models of care that form the basis of how many older people interact with the service

Novel AroA from Pseudomonas putida Confers Tobacco Plant with High Tolerance to Glyphosate

Glyphosate is a non-selective broad-spectrum herbicide that inhibits 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS, also designated as AroA), a key enzyme in the aromatic amino acid biosynthesis pathway in microorganisms and plants. Previously, we reported that a novel AroA (PpAroA1) from Pseudomonas putida had high tolerance to glyphosate, with little homology to class I or class II glyphosate-tolerant AroA. In this study, the coding sequence of PpAroA1 was optimized for tobacco. For maturation of the enzyme in chloroplast, a chloroplast transit peptide coding sequence was fused in frame with the optimized aroA gene (PparoA1optimized) at the 5′ end. The PparoA1optimized gene was introduced into the tobacco (Nicotiana tabacum L. cv. W38) genome via Agrobacterium-mediated transformation. The transformed explants were first screened in shoot induction medium containing kanamycin. Then glyphosate tolerance was assayed in putative transgenic plants and its T1 progeny. Our results show that the PpAroA1 from Pseudomonas putida can efficiently confer tobacco plants with high glyphosate tolerance. Transgenic tobacco overexpressing the PparoA1optimized gene exhibit high tolerance to glyphosate, which suggest that the novel PpAroA1 is a new and good candidate applied in transgenic crops with glyphosate tolerance in future

Tipping the Balance: Robustness of Tip Cell Selection, Migration and Fusion in Angiogenesis

Vascular abnormalities contribute to many diseases such as cancer and diabetic retinopathy. In angiogenesis new blood vessels, headed by a migrating tip cell, sprout from pre-existing vessels in response to signals, e.g., vascular endothelial growth factor (VEGF). Tip cells meet and fuse (anastomosis) to form blood-flow supporting loops. Tip cell selection is achieved by Dll4-Notch mediated lateral inhibition resulting, under normal conditions, in an interleaved arrangement of tip and non-migrating stalk cells. Previously, we showed that the increased VEGF levels found in many diseases can cause the delayed negative feedback of lateral inhibition to produce abnormal oscillations of tip/stalk cell fates. Here we describe the development and implementation of a novel physics-based hierarchical agent model, tightly coupled to in vivo data, to explore the system dynamics as perpetual lateral inhibition combines with tip cell migration and fusion. We explore the tipping point between normal and abnormal sprouting as VEGF increases. A novel filopodia-adhesion driven migration mechanism is presented and validated against in vivo data. Due to the unique feature of ongoing lateral inhibition, ‘stabilised’ tip/stalk cell patterns show sensitivity to the formation of new cell-cell junctions during fusion: we predict cell fates can reverse. The fusing tip cells become inhibited and neighbouring stalk cells flip fate, recursively providing new tip cells. Junction size emerges as a key factor in establishing a stable tip/stalk pattern. Cell-cell junctions elongate as tip cells migrate, which is shown to provide positive feedback to lateral inhibition, causing it to be more susceptible to pathological oscillations. Importantly, down-regulation of the migratory pathway alone is shown to be sufficient to rescue the sprouting system from oscillation and restore stability. Thus we suggest the use of migration inhibitors as therapeutic agents for vascular normalisation in cancer