Search

Pathogens of theMycobacterium tuberculosiscomplex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termedMycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome ofM. tuberculosisis at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Purpose. Differentiation of the Mycobacterium tuberculosis complex (MTBc) from non-tuberculous mycobacteria (NTM) is important for tuberculosis diagnosis and is a prerequisite for reliable phenotypic drug-resistance testing. We evaluated the performance of the rapid MPT64 antigen identification test for the detection of Mycobacterium africanum lineage 5 (MAF L5). Methodology. Smear-positive tuberculosis patients' sputa were included prospectively. Culture was performed on Lowenstein-Jensen medium and, when positive, the MPT64 test and the classical para-nitro benzoic acid susceptibility and heat-labile catalase (PNB/catalase) identification tests were performed. The MPT64 test was repeated 14 days after an initially negative first testing. Direct spoligotyping was performed for MTBc lineage determination. Results. In total, 333 isolates were tested for all of the methods. Three hundred and twenty-two (96.7 %) were pure MTBc, by agreement between spoligotyping and PNB/catalase, and 11 were NTM or a mixture of MTBc/NTM. The MPT64 test conducted on day zero of culture-positivity correctly identified most of the pure MTBc isolates (93.2 %, 300/322), but it failed to detect 24% of the L5 isolates (18/75) versus 2% (4/202) of the L4 ones [OR=15.6 (5.3-45.8), P

Coscolla et al., Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum . Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum , and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally., M.C. is supported by the Ramón y Cajal programme from the Ministerio de Ciencia, Innovación y Universidades. This work was supported by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) (research award to M.C.), Ministerio de Ciencia, Innovación y Universidades (grant number RTI2018-094399-A-I00 to M.C.) and Consellería de Educació de la Generalitat Valenciana (grant number SEJI/2019/011 to M.C.), the Swiss National Science Foundation (grants 310030_188888, IZRJZ3_164171, IZLSZ3_170834 and CRSII5_177163 to S. G.), the European Research Council (883582-ECOEVODRTB to S. G.) and Wellcome (grant number 097134/Z/11/Z to D. Y.-M).

Neisseria gonorrhoeae is an urgent public health threat due to rapidly increasing incidence and antibiotic resistance. In contrast with the trend of increasing resistance, clinical isolates that have reverted to susceptibility regularly appear, prompting questions about which pressures compete with antibiotics to shape gonococcal evolution. Here, we used genome-wide association to identify loss-of-function (LOF) mutations in the efflux pump mtrCDE operon as a mechanism of increased antibiotic susceptibility and demonstrate that these mutations are overrepresented in cervical relative to urethral isolates. This enrichment holds true for LOF mutations in another efflux pump, farAB, and in urogenitally-adapted versus typical N. meningitidis, providing evidence for a model in which expression of these pumps in the female urogenital tract incurs a fitness cost for pathogenic Neisseria. Overall, our findings highlight the impact of integrating microbial population genomics with host metadata and demonstrate how host environmental pressures can lead to increased antibiotic susceptibility., Antibiotic resistance in Neisseria gonorrhoeae is rising, yet sometimes strains emerge that have reverted to susceptibility. Here, the authors find that selective pressures from the host may influence susceptibility through loss-of-function mutations in genes that encode for efflux pumps.

Pathogens of theMycobacterium tuberculosiscomplex (MTBC) are considered monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate the different MTBC lineages (L), especially L5 and L6 (traditionally termedMycobacterium africanum), from each other. However, genome variability and gene content especially of L5 and L6 strains have not been fully explored and may be potentially important for pathobiology and current approaches for genomic analysis of MTBC isolates, including transmission studies.We compared the genomes of 358 L5 clinical isolates (including 3 completed genomes and 355 Illumina WGS (whole genome sequenced) isolates) to the L5 complete genomes and H37Rv, and identified multiple genes differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sublineage into L5.3.1 and L5.3.2. These gene content differences had a small knock on effect on transmission cluster estimation, with clustering rates influenced by the selection of reference genome, and with potential over-estimation of recent transmission when using H37Rv as the reference genome.Our data show that the use of H37Rv as reference genome results in missing SNPs in genes unique for L5 strains. This potentially leads to an underestimation of the diversity present in the genome of L5 strains and in turn affects the transmission clustering rates. As such, a full capture of the gene diversity, especially for high resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most WGS data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome ofM. tuberculosisis at least several kilobases larger than previously thought, implying a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.Data summarySequence data for the Illumina dataset are available at European Genome-phenome Archive (EGA;https://www.ebi.ac.uk/ega/) under the study accession numbers PRJEB38317 and PRJEB38656. Individual runs accession numbers are indicated in Table S8.PacBio raw reads for the L5 Benin genome are available on the ENA accession SAME3170744. The assembled L5 Benin genome is available on NCBI with accession PRJNA641267. To ensure naming conventions of the genes in the three L5 genomes can be followed, we have uploaded these annotated GFF files to figshare athttps://doi.org/10.6084/m9.figshare.12911849.v1.Custom python scripts used in this analysis can be found athttps://github.com/conmeehan/pathophy.

Neisseria gonorrhoeaeis an urgent public health threat due to rapidly increasing incidence and antibiotic resistance. In contrast with the trend of increasing resistance, clinical isolates that have reverted to susceptibility regularly appear, prompting questions about which pressures compete with antibiotics to shape gonococcal evolution. Here, we used genome-wide association on the largest collection ofN. gonorrhoeaeisolates to date (n=4852) to identify loss-of-function (LOF) mutations in the efflux pumpmtrCDEoperon as a mechanism of increased antibiotic susceptibility and demonstrate that these mutations are overrepresented in cervical isolates relative to urethral isolates (odds ratio (OR) = 3.74, 95% CI [1.98-6.70]). In support of a model in which pump expression incurs a fitness cost in this niche, cervical isolates were also enriched relative to urethral isolates in LOF mutations in themtrCDEactivatormtrA(OR = 8.60, 95% CI [4.96-14.57]) and infarA,a subunit of the FarAB efflux pump (OR = 6.25, 95% CI [3.90-9.83]). In total, approximately 2 in 5 cervical isolates (42.6%) contained a LOF mutation in either the efflux pump componentsmtrCorfarAor the activatormtrA.Our findings extend beyondN. gonorrhoeaeto otherNeisseria:mtrCLOF mutations are rare (N. meningitidisin a collection of 14,798 genomes but enriched in a heterosexual urethritis-associated lineage (8.6%, p = 9.90×10−5), indicating that efflux pump downregulation contributes broadly to the adaptation of pathogenicNeisseriato the female urogenital tract. Overall, our findings highlight the impact of integrating microbial population genomics with host metadata and demonstrate how host environmental pressures can lead to increased antibiotic susceptibility.

An antimicrobial activity screen of Burkholderia gladioli BCC0238, a clinical isolate from a cystic fibrosis patient, led to the discovery of gladiolin, a novel macrolide antibiotic with potent activity against Mycobacterium tuberculosis H37Rv. Gladiolin is structurally related to etnangien, a highly unstable antibiotic from Sorangium cellulosum that is also active against Mycobacteria. Like etnangien, gladiolin was found to inhibit RNA polymerase, a validated drug target in M. tuberculosis. However, gladiolin lacks the highly labile hexaene moiety of etnangien and was thus found to possess significantly increased chemical stability. Moreover, gladiolin displayed low mammalian cytotoxicity and good activity against several M. tuberculosis clinical isolates, including four that are resistant to isoniazid and one that is resistant to both isoniazid and rifampicin. Overall, these data suggest that gladiolin may represent a useful starting point for the development of novel drugs to tackle multidrug-resistant tuberculosis. The B. gladioli BCC0238 genome was sequenced using Single Molecule Real Time (SMRT) technology. This resulted in four contiguous sequences: two large circular chromosomes and two smaller putative plasmids. Analysis of the chromosome sequences identified 49 putative specialized metabolite biosynthetic gene clusters. One such gene cluster, located on the smaller of the two chromosomes, encodes a trans-acyltransferase (trans-AT) polyketide synthase (PKS) multienzyme that was hypothesized to assemble gladiolin. Insertional inactivation of a gene in this cluster encoding one of the PKS subunits abrogated gladiolin production, confirming that the gene cluster is responsible for biosynthesis of the antibiotic. Comparison of the PKSs responsible for the assembly of gladiolin and etnangien showed that they possess a remarkably similar architecture, obfuscating the biosynthetic mechanisms responsible for most of the structural differences between the two metabolites.

Summary Whole-genome sequencing can be used to support or refute suspected links between hospital water systems and Legionnaires’ disease cases. However, caveats regarding the interpretation of genomic data from Legionella pneumophila are described that should be considered in future investigations., Background. Legionnaires’ disease is an important cause of hospital-acquired pneumonia and is caused by infection with the bacterium Legionella. Because current typing methods often fail to resolve the infection source in possible nosocomial cases, we aimed to determine whether whole-genome sequencing (WGS) could be used to support or refute suspected links between cases and hospitals. We focused on cases involving a major nosocomial-associated strain, L. pneumophila sequence type (ST) 1. Methods. WGS data from 229 L. pneumophila ST1 isolates were analyzed, including 99 isolates from the water systems of 17 hospitals and 42 clinical isolates from patients with confirmed or suspected hospital-acquired infections, as well as isolates obtained from or associated with community-acquired sources of Legionnaires’ disease. Results. Phylogenetic analysis demonstrated that all hospitals from which multiple isolates were obtained have been colonized by 1 or more distinct ST1 populations. However, deep sampling of 1 hospital also revealed the existence of substantial diversity and ward-specific microevolution within the population. Across all hospitals, suspected links with cases were supported with WGS, although the degree of support was dependent on the depth of environmental sampling and available contextual information. Finally, phylogeographic analysis revealed that hospitals have been seeded with L. pneumophila via both local and international spread of ST1. Conclusions. WGS can be used to support or refute suspected links between hospitals and Legionnaires’ disease cases. However, deep hospital sampling is frequently required due to the potential coexistence of multiple populations, existence of substantial diversity, and similarity of hospital isolates to local populations.

14 páginas, 6 figuras, Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants., This work was funded by projects of the European Research Council (638553-TB-ACCELERATE) and Ministerio de Economía y Competitividad (Spanish government) research grant SAF2016-77346-R (to I.C.); BFU2014-58656-R and BFU2017-89594-R from Ministerio de Economía y Competitividad (Spanish government) and PROMETEO/2016/122 from Generalitat Valenciana (to F.G.-C); and the Swiss National Science Foundation (grants 310030_166687, IZRJZ3_164171, IZLSZ3_170834, and CRSII5_177163), the European Research Council (309540-EVODRTB), and SystemsX.ch (to S.G.). Á.C.-O. is the recipient of an FPU fellowship from Ministerio de Ciencia, Innovación y Universidades FPU13/00913 (Spanish government). L.S.-B. was funded by Wellcome grant number 098051.

Burkholderia is a multi-talented genus of Gram-negative bacteria, which in recent years has become increasingly recognised as a promising source of bioactive natural products. Metabolite profiling of Burkholderia gladioli BCC0238 showed that it produces the asymmetric lipopeptidiolide antibiotic icosalide A1, originally isolated from a fungus. Comparative bioinformatics analysis of several genome-sequenced B. gladioli isolates identified a gene encoding a nonribosomal peptide synthase (NRPS) with an unusual architecture that was predicted to be responsible for icosalide biosynthesis. Inactivation of this gene in B. gladioli BCC0238 abolished icosalide production. PCR analysis and sequencing of total DNA from the original fungal icosalide A1 producer revealed it has a B. gladioli strain associated with it that harbours an NRPS with an identical architecture to that responsible for icosalide A1 assembly in B. gladioli BCC0238. Sequence analysis of the icosalide NRPS indicated that it contains two chain-initiating condensation (CI) domains. One of these is appended to the N-terminus of module 1 - a common architecture for NRPSs involved in lipopeptide assembly. The other is embedded in module 3, immediately downstream of a putative chain-elongating condensation domain. Analysis of the reactions catalysed by a tridomain construct from module 3 of the NRPS using intact protein mass spectrometry showed that the embedded CI domain initiates assembly of a second lipopeptide chain, providing key insights into the mechanism for asymmetric diolide assembly. Open Access article ispartof: Chemical Science vol:10 issue:21 ispartof: location:England status: published

The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease., Claire Chewapreecha et al. combine 753 newly sequenced Thai Burkholderia pseudomallei isolates with 258 Australian isolates to identify genes associated with either clinical or environmental strains. They find 47 genes, 12 of which replicate in both locations, that may provide clues to the strategy used by this microbe to adapt to survive in wide range of ecological niches, including human hosts.

Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans. The disease is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Currently, the major hurdles facing disease control are incomplete understandings of both the mode of transmission and environmental reservoirs of M. ulcerans. As decades of spasmodic environmental sampling surveys have not brought us much closer to overcoming these hurdles, the Buruli ulcer research community has recently switched to using comparative genomics. The significance of our research is in how we used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our approach shows that it might be possible to use bacterial population genomics to assess the impact of health interventions, providing valuable feedback for managers of disease control programs in areas where health surveillance infrastructure is poor., Buruli ulcer is a neglected tropical disease of skin and subcutaneous tissue caused by infection with the pathogen Mycobacterium ulcerans. Many critical issues for disease control, such as understanding the mode of transmission and identifying source reservoirs of M. ulcerans, are still largely unknown. Here, we used genomics to reconstruct in detail the evolutionary trajectory and dynamics of M. ulcerans populations at a central African scale and at smaller geographical village scales. Whole-genome sequencing (WGS) data were analyzed from 179 M. ulcerans strains isolated from all Buruli ulcer foci in the Democratic Republic of the Congo, The Republic of Congo, and Angola that have ever yielded positive M. ulcerans cultures. We used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our phylogeographic analysis revealed one almost exclusively predominant sublineage of M. ulcerans that arose in Central Africa and proliferated in its different regions of endemicity during the Age of Discovery. We observed how the best sampled endemic hot spot, the Songololo territory, became an area of endemicity while the region was being colonized by Belgium (1880s). We furthermore identified temporal parallels between the observed past population fluxes of M. ulcerans from the Songololo territory and the timing of health policy changes toward control of the Buruli ulcer epidemic in that region. These findings suggest that an intervention based on detecting and treating human cases in an area of endemicity might be sufficient to break disease transmission chains, irrespective of other reservoirs of the bacterium. IMPORTANCE Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans. The disease is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Currently, the major hurdles facing disease control are incomplete understandings of both the mode of transmission and environmental reservoirs of M. ulcerans. As decades of spasmodic environmental sampling surveys have not brought us much closer to overcoming these hurdles, the Buruli ulcer research community has recently switched to using comparative genomics. The significance of our research is in how we used both temporal associations and the study of the mycobacterial demographic history to estimate the contribution of humans as a reservoir in Buruli ulcer transmission. Our approach shows that it might be possible to use bacterial population genomics to assess the impact of health interventions, providing valuable feedback for managers of disease control programs in areas where health surveillance infrastructure is poor.

[Background]: Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen., [Methods]: Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes., [Results]: Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p, [Conclusion]: Our data indicate potential zoonotic transmission of bTB in Ghana and hence calls for intensified public education on bTB, especially among risk groups., This work was supported by the Wellcome Trust Intermediate Fellowship awarded to DYM (Grant Number 097134/Z/11/Z).

The ability to DNA fingerprint Mycobacterium bovis (M. bovis) isolates helped to define the role of wildlife in the persistnace of bovine tuberculosis in New Zealand. DNA fingerprinting results currently help to guide wildlife control measures and also aid in tracing the source of infections that result from movement of livestock. During the last 5 years we have developed the ability to distinguish New Zealand (NZ) M. bovis isolates by comparing the sequences of whole genome sequenced (WGS) M. bovis samples. WGS provides much higher resolution than our other established typing methods and greatly improved the definition of the regional localization of NZ M. bovis types. Three outbreak investigations are described and results demonstrate how WGS analysis has led to the confirmation of epidemiological sourcing of infection, to better definition of the source of new infections by ruling out other possible sources, and has revealed probable wildlife infection in an area considered to be free of infected wildlife. The routine use of WGS analyses for sourcing new M. bovis infections will be an important component of the strategy employed to eradicate bovine TB from NZ livestock and wildlife.

Sequence-based typing (SBT), analogous to multilocus sequence typing (MLST), is the current “gold standard” typing method for investigation of legionellosis outbreaks caused by Legionella pneumophila . However, as common sequence types (STs) cause many infections, some investigations remain unresolved. In this study, various whole-genome sequencing (WGS)-based methods were evaluated according to published guidelines, including (i) a single nucleotide polymorphism (SNP)-based method, (ii) extended MLST using different numbers of genes, (iii) determination of gene presence or absence, and (iv) a kmer-based method. L. pneumophila serogroup 1 isolates ( n = 106) from the standard “typing panel,” previously used by the European Society for Clinical Microbiology Study Group on Legionella Infections (ESGLI), were tested together with another 229 isolates. Over 98% of isolates were considered typeable using the SNP- and kmer-based methods. Percentages of isolates with complete extended MLST profiles ranged from 99.1% (50 genes) to 86.8% (1,455 genes), while only 41.5% produced a full profile with the gene presence/absence scheme. Replicates demonstrated that all methods offer 100% reproducibility. Indices of discrimination range from 0.972 (ribosomal MLST) to 0.999 (SNP based), and all values were higher than that achieved with SBT (0.940). Epidemiological concordance is generally inversely related to discriminatory power. We propose that an extended MLST scheme with ∼50 genes provides optimal epidemiological concordance while substantially improving the discrimination offered by SBT and can be used as part of a hierarchical typing scheme that should maintain backwards compatibility and increase discrimination where necessary. This analysis will be useful for the ESGLI to design a scheme that has the potential to become the new gold standard typing method for L. pneumophila .

Global spread of aggressive mycobacteria Many mycobacteria, in addition to those causing leprosy and tuberculosis, are capable of infecting humans. Some can be particularly dangerous in patients suffering from immunosuppression or chronic disease, such as cystic fibrosis (CF). Bryant et al. observed clusters of near-identical isolates of drug-resistant Mycobacterium abscessus in patients reporting to CF clinics. The similarity of the isolates suggests transmission between patients, rather than environmental acquisition. Although this bacterium is renowned for its environmental resilience, the mechanism for its long-distance transmission among the global CF patient community remains a puzzle. Science , this issue p. 751

Background Although Mycoplasma genitalium is a common sexually transmitted pathogen causing clinically distinct diseases both in male and females, few genomes have been sequenced up to now, due mainly to its fastidious nature and slow growth. Hence, we lack a robust phylogenetic framework to provide insights into the population structure of the species. Currently our understanding of the nature and diversity of M. genitalium relies on molecular tests targeting specific genes or regions of the genome and knowledge is limited by a general under-testing internationally. This is set against a background of drug resistance whereby M. genitalium has developed resistance to mainly all therapeutic antimicrobials. Results We sequenced 28 genomes of Mycoplasma genitalium from temporally (1980–2010) and geographically (Europe, Japan, Australia) diverse sources. All the strain showed essentially the same genomic content without any accessory regions found. However, we identified extensive recombination across their genomes with a total of 25 regions showing heightened levels of SNP density. These regions include the MgPar loci, associated with host interactions, as well as other genes that could also be involved in this role. Using these data, we generated a robust phylogeny which shows that there are two main clades with differing degrees of genomic variability. SNPs found in region V of 23S rRNA and parC were consistent with azithromycin/erythromycin and fluoroquinolone resistances, respectively, and with their phenotypic MIC data. Conclusions The sequence data here generated is essential for designing rational approaches to type and track Mycoplasma genitalium as antibiotic resistance increases. It represents a first approach to its population genetics to better appreciate the role of this organism as a sexually transmitted pathogen. Electronic supplementary material The online version of this article (10.1186/s12864-017-4399-6) contains supplementary material, which is available to authorized users.

Mycobacterium africanum(Maf) causes up to half of human tuberculosis in West Africa, but little is known on this pathogen. We compared the genomes of 253Mafclinical isolates from Ghana, including both L5 and L6. We found that the genomic diversity of L6 was higher than in L5, and the selection pressures differed between both groups. Regulatory proteins appeared to evolve neutrally in L5 but under purifying selection in L6. Conversely, human T cell epitopes were under purifying selection in L5, but under positive selection in L6. Although only 10% of the T cell epitopes were variable, mutations were mostly lineage-specific. Our findings indicate thatMafL5 and L6 are genomically distinct, possibly reflecting different ecological niches.

Antimicrobial resistance (AMR) is one of the major threats to human and animal health worldwide, yet few high-throughput tools exist to analyse and predict the resistance of a bacterial isolate from sequencing data. Here we present a new tool, ARIBA, that identifies AMR-associated genes and single nucleotide polymorphisms directly from short reads, and generates detailed and customisable output. The accuracy and advantages of ARIBA over other tools are demonstrated on three datasets from Gram-positive and Gram-negative bacteria, with ARIBA outperforming existing methods. ARIBA is available at https://github.com/sanger-pathogens/ariba.

A major outbreak of whooping cough, or pertussis, occurred in 2012 in the United Kingdom (UK), with nearly 10 000 laboratory-confirmed cases and 14 infant deaths attributed to pertussis. A worldwide resurgence of pertussis has been linked to switch to the use of acellular pertussis vaccines and the evolution of Bordetella pertussis away from vaccine-mediated immunity. We have conducted genomic analyses of multiple strains from the UK outbreak. We show that the UK outbreak was polyclonal in nature, caused by multiple distinct but closely related strains. Importantly, we demonstrate that acellular vaccine antigen-encoding genes are evolving at higher rates than other surface protein-encoding genes. This was true even prior to the introduction of pertussis vaccines but has become more pronounced since the introduction of the current acellular vaccines. The fast evolution of vaccine antigen-encoding genes has serious consequences for the ability of current vaccines to continue to control pertussis.

Iceland is free of the major infectious diseases of horses. However, in 2010 an epidemic of respiratory disease of unknown cause spread through the country’s native horse population of 77,000. Microbiological investigations ruled out known viral agents but identified the opportunistic pathogen Streptococcus equi subsp. zooepidemicus (S. zooepidemicus) in diseased animals. We sequenced the genomes of 257 isolates of S. zooepidemicus to differentiate epidemic from endemic strains. We found that although multiple endemic clones of S. zooepidemicus were present, one particular clone, sequence type 209 (ST209), was likely to have been responsible for the epidemic. Concurrent with the epidemic, ST209 was also recovered from a human case of septicemia, highlighting the pathogenic potential of this strain. Epidemiological investigation revealed that the incursion of this strain into one training yard during February 2010 provided a nidus for the infection of multiple horses that then transmitted the strain to farms throughout Iceland. This study represents the first time that whole-genome sequencing has been used to investigate an epidemic on a national scale to identify the likely causative agent and the link to an associated zoonotic infection. Our data highlight the importance of national biosecurity to protect vulnerable populations of animals and also demonstrate the potential impact of S. zooepidemicus transmission to other animals, including humans. IMPORTANCE An epidemic of respiratory disease affected almost the entire native Icelandic horse population of 77,000 animals in 2010, resulting in a self-imposed ban on the export of horses and significant economic costs to associated industries. Although the speed of transmission suggested that a viral pathogen was responsible, only the presence of the opportunistic pathogen Streptococcus zooepidemicus was consistent with the observed clinical signs. We applied genomic sequencing to differentiate epidemic from endemic strains and to shed light on the rapid transmission of the epidemic strain throughout Iceland. We further highlight the ability of epidemic and endemic strains of S. zooepidemicus to infect other animals, including humans. This study represents the first time that whole-genome sequencing has been used to elucidate an outbreak on a national scale and identify the likely causative agent., Wellcome Trust 098051

This work was supported by the Animal Health Trust. The zoonotic bacterium Streptococcus equi subsp. zooepidemicus (S. zooepidemicus) is a diverse, opportunistic pathogen that can cause mastitis in dairy sheep and goats. We used multilocus sequence typing (MLST) to define the genetic diversity of 60 isolates of S. zooepidemicus, which were recovered from sheep and goats in Spain between 2003 and 2010. We identify a novel clonal complex based on sequence type (ST), ST-236, which accounted for 39 of the 60 isolates. A representative ST-236 strain, S. zooepidemicus strain C7 (SzC7), was sequenced and interrogated for the presence of novel nutritional uptake or utilisation systems, the acquisition of which have previously been shown to be important for environmental adaptation in other streptococcal pathogens. A novel phosphoenolpyruvate sugar phosphotransferase system (PTS), which enabled the utilisation of arbutin, was identified. Functionality of the PTS was confirmed following deletion of the PTS from SzC7. Arbutin is found in multiple animal foodstuffs and we propose that the ability to utilise arbutin may have conferred a selective advantage to strains infecting animals, the diet of which contains this sugar. Postprint

Legionella pneumophila is an environmental bacterium and the causative agent of Legionnaires’ disease. Previous genomic studies have shown that recombination accounts for a high proportion (>96%) of diversity within several major disease-associated sequence types (STs) of L. pneumophila. This suggests that recombination represents a potentially important force shaping adaptation and virulence. Despite this, little is known about the biological effects of recombination in L. pneumophila, particularly with regards to homologous recombination (whereby genes are replaced with alternative allelic variants). Using newly available population genomic data, we have disentangled events arising from homologous and non-homologous recombination in six major disease-associated STs of L. pneumophila (subsp. pneumophila), and subsequently performed a detailed characterisation of the dynamics and impact of homologous recombination. We identified genomic “hotspots” of homologous recombination that include regions containing outer membrane proteins, the lipopolysaccharide (LPS) region and Dot/Icm effectors, which provide interesting clues to the selection pressures faced by L. pneumophila. Inference of the origin of the recombined regions showed that isolates have most frequently imported DNA from isolates belonging to their own clade, but also occasionally from other major clades of the same subspecies. This supports the hypothesis that the possibility for horizontal exchange of new adaptations between major clades of the subspecies may have been a critical factor in the recent emergence of several clinically important STs from diverse genomic backgrounds. However, acquisition of recombined regions from another subspecies, L. pneumophila subsp. fraseri, was rarely observed, suggesting the existence of a recombination barrier and/or the possibility of ongoing speciation between the two subspecies. Finally, we suggest that multi-fragment recombination may occur in L. pneumophila, whereby multiple non-contiguous segments that originate from the same molecule of donor DNA are imported into a recipient genome during a single episode of recombination., Author summary Legionella pneumophila is an environmental bacterium that causes Legionnaires’ disease, a serious and potentially fatal pneumonia. Previous studies have shown that members of this species undergo a process called recombination, whereby DNA is imported from another bacterial cell into the recipient genome. The imported DNA can either replace an equivalent segment of the genome (homologous recombination) or can comprise novel genes that are new to the recipient genome (non-homologous recombination). Whilst recombination plays an undoubtedly important role in L. pneumophila evolution, accounting for more than 96% of the diversity observed within some lineages, little is known about its biological impact. In this study, we performed a detailed characterisation of the dynamics and effect of homologous recombination on L. pneumophila evolution in six clinically important lineages of L. pneumophila. We identified “hotspot” regions of the genome in which an excess of homologous recombination events was observed, which provided important clues to the selection pressures faced by L. pneumophila. By determining the donors of the recombined regions, we also revealed that recombination has occurred most frequently between isolates from the same clade, but also occurred between isolates from different major clades. This demonstrates the possibility of new adaptations arising in one lineage and being transferred to another distantly related lineage, which we predict has been an important factor in the emergence of several major disease-causing strains from diverse genomic backgrounds.

Background Chlamydia abortus (formerly Chlamydophila abortus) is an economically important livestock pathogen, causing ovine enzootic abortion (OEA), and can also cause zoonotic infections in humans affecting pregnancy outcome. Large-scale genomic studies on other chlamydial species are giving insights into the biology of these organisms but have not yet been performed on C. abortus. Our aim was to investigate a broad collection of European isolates of C. abortus, using next generation sequencing methods, looking at diversity, geographic distribution and genome dynamics. Results Whole genome sequencing was performed on our collection of 57 C. abortus isolates originating primarily from the UK, Germany, France and Greece, but also from Tunisia, Namibia and the USA. Phylogenetic analysis of a total of 64 genomes shows a deep structural division within the C. abortus species with a major clade displaying limited diversity, in addition to a branch carrying two more distantly related Greek isolates, LLG and POS. Within the major clade, seven further phylogenetic groups can be identified, demonstrating geographical associations. The number of variable nucleotide positions across the sampled isolates is significantly lower than those published for C. trachomatis and C. psittaci. No recombination was identified within C. abortus, and no plasmid was found. Analysis of pseudogenes showed lineage specific loss of some functions, notably with several Pmp and TMH/Inc proteins predicted to be inactivated in many of the isolates studied. Conclusions The diversity within C. abortus appears to be much lower compared to other species within the genus. There are strong geographical signatures within the phylogeny, indicating clonal expansion within areas of limited livestock transport. No recombination has been identified within this species, showing that different species of Chlamydia may demonstrate different evolutionary dynamics, and that the genome of C. abortus is highly stable. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3657-y) contains supplementary material, which is available to authorized users.

An epidemic of post-surgical wound infections, caused by a non-tuberculous mycobacterium, has been on-going in Brazil. It has been unclear whether one or multiple lineages are responsible and whether their wide geographical distribution across Brazil is due to spread from a single point source or is the result of human-mediated transmission. 188 isolates, collected from nine Brazilian states, were whole genome sequenced and analysed using phylogenetic and comparative genomic approaches. The isolates from Brazil formed a single clade, which was estimated to have emerged in 2003. We observed temporal and geographic structure within the lineage that enabled us to infer the movement of sub-lineages across Brazil. The genome size of the Brazilian lineage was reduced relative to most strains in the three subspecies of Mycobacterium abscessus and contained a novel plasmid, pMAB02, in addition to the previously described pMAB01 plasmid. One lineage, which emerged just prior to the initial outbreak, is responsible for the epidemic of post-surgical wound infections in Brazil. Phylogenetic analysis indicates that multiple transmission events led to its spread. The presence of a novel plasmid and the reduced genome size suggest that the lineage has undergone adaptation to the surgical niche.

Background Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is an important livestock disease raising public health and economic concerns around the world. In New Zealand, a number of wildlife species are implicated in the spread and persistence of bTB in cattle populations, most notably the brushtail possum (Trichosurus vulpecula). Whole Genome Sequenced (WGS) M. bovis isolates sourced from infected cattle and wildlife across New Zealand were analysed. Bayesian phylogenetic analyses were conducted to estimate the substitution rate of the sampled population and investigate the role of wildlife. In addition, the utility of WGS was examined with a view to these methods being incorporated into routine bTB surveillance. Results A high rate of exchange was evident between the sampled wildlife and cattle populations but directional estimates of inter-species transmission were sensitive to the sampling strategy employed. A relatively high substitution rate was estimated, this, in combination with a strong spatial signature and a good agreement to previous typing methods, acts to endorse WGS as a typing tool. Conclusions In agreement with the current knowledge of bTB in New Zealand, transmission of M. bovis between cattle and wildlife was evident. Without direction, these estimates are less informative but taken in conjunction with the low prevalence of bTB in New Zealand’s cattle population it is likely that, currently, wildlife populations are acting as the main bTB reservoir. Wildlife should therefore continue to be targeted if bTB is to be eradicated from New Zealand. WGS will be a considerable aid to bTB eradication by greatly improving the discriminatory power of molecular typing data. The substitution rates estimated here will be an important part of epidemiological investigations using WGS data. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3569-x) contains supplementary material, which is available to authorized users.

The environmental bacterium Burkholderia pseudomallei causes an estimated 165,000 cases of human melioidosis per year worldwide and is also classified as a biothreat agent. We used whole genome sequences of 469 B. pseudomallei isolates from 30 countries collected over 79 years to explore its geographic transmission. Our data point to Australia as an early reservoir, with transmission to Southeast Asia followed by onward transmission to South Asia and East Asia. Repeated reintroductions were observed within the Malay Peninsula and between countries bordered by the Mekong River. Our data support an African origin of the Central and South American isolates with introduction of B. pseudomallei into the Americas between 1650 and 1850, providing a temporal link with the slave trade. We also identified geographically distinct genes/variants in Australasian or Southeast Asian isolates alone, with virulence-associated genes being among those over-represented. This provides a potential explanation for clinical manifestations of melioidosis that are geographically restricted. Postprint

The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail. These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.

The molecular mechanisms determining the transmissibility and prevalence of drug-resistant tuberculosis in a population were investigated through whole-genome sequencing of 1,000 prospectively obtained patient isolates from Russia. Two-thirds belonged to the Beijing lineage, which was dominated by two homogeneous clades. Multidrug-resistant (MDR) genotypes were found in 48% of isolates overall and in 87% of the major clades. The most common rpoB mutation was associated with fitness-compensatory mutations in rpoA or rpoC, and a new intragenic compensatory substitution was identified. The proportion of MDR cases with extensively drug-resistant (XDR) tuberculosis was 16% overall, with 65% of MDR isolates harboring eis mutations, selected by kanamycin therapy, which may drive the expansion of strains with enhanced virulence. The combination of drug resistance and compensatory mutations displayed by the major clades confers clinical resistance without compromising fitness and transmissibility, showing that, in addition to weaknesses in the tuberculosis control program, biological factors drive the persistence and spread of MDR and XDR tuberculosis in Russia and beyond.

Prediction of antibiotic resistance phenotypes from whole genome sequencing data by machine learning methods has been proposed as a promising platform for the development of sequence-based diagnostics. However, there has been no systematic evaluation of factors that may influence performance of such models, how they might apply to and vary across clinical populations, and what the implications might be in the clinical setting. Here, we performed a meta-analysis of seven large Neisseria gonorrhoeae datasets, as well as Klebsiella pneumoniae and Acinetobacter baumannii datasets, with whole genome sequence data and antibiotic susceptibility phenotypes using set covering machine classification, random forest classification, and random forest regression models to predict resistance phenotypes from genotype. We demonstrate how model performance varies by drug, dataset, resistance metric, and species, reflecting the complexities of generating clinically relevant conclusions from machine learning-derived models. Our findings underscore the importance of incorporating relevant biological and epidemiological knowledge into model design and assessment and suggest that doing so can inform tailored modeling for individual drugs, pathogens, and clinical populations. We further suggest that continued comprehensive sampling and incorporation of up-to-date whole genome sequence data, resistance phenotypes, and treatment outcome data into model training will be crucial to the clinical utility and sustainability of machine learning-based molecular diagnostics., Author summary Machine learning-based prediction of antibiotic resistance from bacterial genome sequences represents a promising tool to rapidly determine the antibiotic susceptibility profile of clinical isolates and reduce the morbidity and mortality resulting from inappropriate and ineffective treatment. However, while there has been much focus on demonstrating the diagnostic potential of these modeling approaches, there has been little assessment of potential caveats and prerequisites associated with implementing predictive models of drug resistance in the clinical setting. Our results highlight significant biological and technical challenges facing the application of machine learning-based prediction of antibiotic resistance as a diagnostic tool. By outlining specific factors affecting model performance, our findings provide a framework for future work on modeling drug resistance and underscore the necessity of continued comprehensive sampling and reporting of treatment outcome data for building reliable and sustainable diagnostics.

Comparison of the genome sequences of three Pseudomonas fluorescens strains reveals a heterogeneity reminiscent of a species complex rather than a single species, Background Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

The native horse population of Iceland has remained free of major infectious diseases. Between May and July 2010 an epidemic of respiratory disease swept through the population. Initial microbiological investigations ruled out known equine viral agents as the cause of the infections, but identified the opportunistic pathogenStreptococcus zooepidemicusas being frequently isolated from diseased animals. This diverse bacterial species has a broad host range and is usually regarded as a commensal of horses. By genome sequencingS. zooepidemicusrecovered from horses during the epidemic we show that although multiple clones ofS. zooepidemicuswere present in the population, one particular clone, ST209, was responsible for the epidemic. Concurrent with the epidemic, ST209 caused zoonotic infections, highlighting the pathogenic potential of this clone. Phylogenetic analysis suggests that the original ST209 strain entered Iceland in late 2008 or early 2009. Epidemiological investigation revealed that the incursion of this strain into a training yard that utilized a submerged treadmill between the 5thand 19thof February 2010 was a critical trigger for the ensuing epidemic of disease, provided a nidus for the infection of multiple horses, and subsequent distribution of these animals to multiple sites in Iceland.

Legionella pneumophila is an environmental bacterium and the leading cause of Legionnaires’ disease. Just five sequence types (ST), from more than 2000 currently described, cause nearly half of disease cases in northwest Europe. Here, we report the sequence and analyses of 364 L. pneumophila genomes, including 337 from the five disease-associated STs and 27 representative of the species diversity. Phylogenetic analyses revealed that the five STs have independent origins within a highly diverse species. The number of de novo mutations is extremely low with maximum pairwise single-nucleotide polymorphisms (SNPs) ranging from 19 (ST47) to 127 (ST1), which suggests emergences within the last century. Isolates sampled geographically far apart differ by only a few SNPs, demonstrating rapid dissemination. These five STs have been recombining recently, leading to a shared pool of allelic variants potentially contributing to their increased disease propensity. The oldest clone, ST1, has spread globally; between 1940 and 2000, four new clones have emerged in Europe, which show long-distance, rapid dispersal. That a large proportion of clinical cases is caused by recently emerged and internationally dispersed clones, linked by convergent evolution, is surprising for an environmental bacterium traditionally considered to be an opportunistic pathogen. To simultaneously explain recent emergence, rapid spread and increased disease association, we hypothesize that these STs have adapted to new man-made environmental niches, which may be linked by human infection and transmission.

Background Mycobacterium avium subspecies paratuberculosis (Map) is an infectious enteric pathogen that causes Johne’s disease in livestock. Determining genetic diversity is prerequisite to understanding the epidemiology and biology of Map. We performed the first whole genome sequencing (WGS) of 141 global Map isolates that encompass the main molecular strain types currently reported. We investigated the phylogeny of the Map strains, the diversity of the genome and the limitations of commonly used genotyping methods. Results Single nucleotide polymorphism (SNP) and phylogenetic analyses confirmed two major lineages concordant with the former Type S and Type C designations. The Type I and Type III strain groups are subtypes of Type S, and Type B strains are a subtype of Type C and not restricted to Bison species. We found that the genome-wide SNPs detected provided greater resolution between isolates than currently employed genotyping methods. Furthermore, the SNP used for IS1311 typing is not informative, as it is likely to have occurred after Type S and C strains diverged and does not assign all strains to the correct lineage. Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) differentiates Type S from Type C but provides limited resolution between isolates within these lineages and the polymorphisms detected do not necessarily accurately reflect the phylogenetic relationships between strains. WGS of passaged strains and coalescent analysis of the collection revealed a very high level of genetic stability, with the substitution rate estimated to be less than 0.5 SNPs per genome per year. Conclusions This study clarifies the phylogenetic relationships between the previously described Map strain groups, and highlights the limitations of current genotyping techniques. Map isolates exhibit restricted genetic diversity and a substitution rate consistent with a monomorphic pathogen. WGS provides the ultimate level of resolution for differentiation between strains. However, WGS alone will not be sufficient for tracing and tracking Map infections, yet importantly it can provide a phylogenetic context for affirming epidemiological connections. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2234-5) contains supplementary material, which is available to authorized users.

Host adaptation is a key factor contributing to the emergence of new bacterial, viral and parasitic pathogens. Many pathogens are considered promiscuous because they cause disease across a range of host species, while others are host-adapted, infecting particular hosts1. Host adaptation can potentially progress to host restriction where the pathogen is strictly limited to a single host species and is frequently associated with more severe symptoms. Host-adapted and host-restricted bacterial clades evolve from within a broader host-promiscuous species and sometimes target different niches within their specialist hosts, such as adapting from a mucosal to a systemic lifestyle. Genome degradation, marked by gene inactivation and deletion, is a key feature of host adaptation, although the triggers initiating genome degradation are not well understood. Here, we show that a chronic systemic non-typhoidal Salmonella infection in an immunocompromised human patient resulted in genome degradation targeting genes that are expendable for a systemic lifestyle. We present a genome-based investigation of a recurrent blood-borne Salmonella enterica serotype Enteritidis (S. Enteritidis) infection covering 15 years in an interleukin (IL)-12 β-1 receptor-deficient individual that developed into an asymptomatic chronic infection. The infecting S. Enteritidis harbored a mutation in the mismatch repair gene mutS that accelerated the genomic mutation rate. Phylogenetic analysis and phenotyping of multiple patient isolates provides evidence for a remarkable level of within-host evolution that parallels genome changes present in successful host-restricted bacterial pathogens but never before observed on this timescale. Our analysis identifies common pathways of host adaptation and demonstrates the role that immunocompromised individuals can play in this process.

Objectives The global emergence of OXA-48-producing Klebsiella pneumoniae clones is a significant threat to public health. We used WGS and phylogenetic analysis of Spanish isolates to investigate the population structure of blaOXA-48-like-expressing K. pneumoniae ST11 and ST405 and to determine the distribution of resistance genes and plasmids encoding blaOXA-48-like carbapenemases. Methods SNPs identified in whole-genome sequences were used to reconstruct phylogenetic trees, identify resistance determinants and de novo assemble the genomes of 105 blaOXA-48-like-expressing K. pneumoniae isolates. Results Genome variation was generally lower in outbreak-associated isolates compared with those associated with sporadic infections. The relatively limited variation observed within the outbreak-associated isolates was on average 7-10 SNPs per outbreak. Of 24 isolates from suspected sporadic infections, 7 were very closely related to isolates causing hospital outbreaks and 17 were more diverse and therefore probably true sporadic cases. On average, 14 resistance genes were identified per isolate. The 17 ST405 isolates from sporadic cases of infection had four distinct resistance gene profiles, while the resistance gene profile differed in all ST11 isolates from sporadic cases. Sequence analysis of 94 IncL/M plasmids carrying blaOXA-48-like genes revealed an average of two SNP differences, indicating a conserved plasmid clade. Conclusions Whole-genome sequence analysis enabled the discrimination of outbreak and sporadic isolates. Significant inter-regional spread within Spain of highly related isolates was evident for both ST11 and ST405 K. pneumoniae. IncL/M plasmids carrying blaOXA-48-like carbapenemase genes were highly conserved geographically and across the outbreaks, sporadic cases and clones.

Sequencing of the bla KPC -positive strain Proteus mirabilis AOUC-001 was performed using both the MiSeq and PacBio RS II platforms and yielded a single molecule of 4,272,433 bp, representing the complete chromosome. Genome analysis showed the presence of several acquired resistance determinants, including two copies of bla KPC-2 carried on a fragment of a KPC-producing plasmid previously described in Klebsiella pneumoniae .

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

Summary The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations. PaperClip, Graphical Abstract Highlights ► Whole-genome sequences of the Tasmanian devil and two distant cancer subclones ► The Tasmanian devil cancer lineage originated recently in a female devil ► The devil cancer genome is relatively stable despite ongoing evolution ► Clonal divergence and geographic spread elucidated through patterns of mutation, Whole-genome sequences of the Tasmanian devil and two devil cancer subclones suggest that the cancer first arose from a female devil and that the clone has subsequently genetically diverged during its spread across Tasmania.

Objective/background Mycobacterium africanum (MAF) remains an important TB causing pathogen in West Africa; however, little is known about its population structure and actual diversity which may have implications for diagnostics and vaccines. We carried out comparative genomics analysis of candidate Mycobacterium tuberculosis (MTB) and MAF using whole genome sequencing. Methods Clinical MTB complex strains (n = 187) comprising L4 (n = 22), L5 (n = 126), and L6 (n = 39) isolated over 8 years from Ghana were whole genome sequenced. The reads were mapped onto a reference genome for phylogenetic and functional genomics analysis. A maximum likelihood tree with 100 bootstraps was constructed from the single nucleotide polymorphisms (SNPs) found using RAxML and clustered with hierBAPS. A total of 147 (18 L4, 36 L6, and 93 L5) of the genomes were de novo assembled and annotated for comparative pangenome analysis using Roary. Results The population structure of MAF revealed at least five clusters of L5 as compared to three for L6. We also identified a group of three multi-drug-resistants (MDRs) within a single cluster of L5 strains from Southern Ghana isolated in 2013. Among the global collection of MTB complex, there were four Ghana-specific L5 clusters of which one (L5.1.1) had traits of clonal expansion. From the 5947 pan genes extracted from the collection, 3215 (54.1%) were core to all the 147 genomes whereas 719 (12.1%) were found in single genomes. Most of the variable genes were PE-PGRS/PPE (1,281) duplicates of other genes (431). The genome degradation was more pronounced in Lineages 4 and 6 as compared to Lineage 5. We identified the absence of some unique genes among specific lineages and/or clades with possible clinical implications. For example, mpt64 and mlaD encoding respectively an immunogenic protein and a mammalian cell entry protein were missing from all L6 genomes. In addition, all L5 strains had an amino acid substitution I43 N within the mpt64 gene. Analysis of SNPs within some genes encoding proteins for substrate metabolism, ion transport and secretory systems showed higher proportion of SNPs among L6 compared to L5 and L4. We also identified a number of lineage/sublineage specific SNPs and indels that may be utilized in rapid PCR based genotyping of MTB complex. Conclusion This work emphasizes on the possibility that the mpt64-based rapid diagnostic kit would not be effective in MAF endemic settings. More mutations in ESAT-6 secretory system of MAF compared to MTB sensu stricto can affect efficacy of ESAT-6-based vaccines in the future.

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. Conclusions 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.

[This corrects the article DOI: 10.1371/journal.pgen.1006855.].

The origin of the eukaryotic genetic apparatus is thought to be central to understanding the evolution of the eukaryotic cell. Disagreement about the source of the relevant genes has spawned competing hypotheses for the origins of the eukaryote nuclear lineage. The iconic rooted 3-domains tree of life shows eukaryotes and archaebacteria as separate groups that share a common ancestor to the exclusion of eubacteria. By contrast, the eocyte hypothesis has eukaryotes originating within the archaebacteria and sharing a common ancestor with a particular group called the Crenarchaeota or eocytes. Here, we have investigated the relative support for each hypothesis from analysis of 53 genes spanning the 3 domains, including essential components of the eukaryotic nucleic acid replication, transcription, and translation apparatus. As an important component of our analysis, we investigated the fit between model and data with respect to composition. Compositional heterogeneity is a pervasive problem for reconstruction of ancient relationships, which, if ignored, can produce an incorrect tree with strong support. To mitigate its effects, we used phylogenetic models that allow for changing nucleotide or amino acid compositions over the tree and data. Our analyses favor a topology that supports the eocyte hypothesis rather than archaebacterial monophyly and the 3-domains tree of life.

Yersinia enterocolitica is a common cause of food-borne gastroenteritis worldwide. Recent work defining the phylogeny of the genus Yersinia subdivided Y. enterocolitica into six distinct phylogroups. Here, we provide detailed analyses of the evolutionary processes leading to the emergence of these phylogroups. The dominant phylogroups isolated from human infections, PG3-5, show very little diversity at the sequence level, but do present marked patterns of gain and loss of functions, including those involved in pathogenicity and metabolism, including the acquisition of phylogroup-specific O-antigen loci. We tracked gene flow across the species in the core and accessory genome, and show that the non-pathogenic PG1 strains act as a reservoir for diversity, frequently acting as donors in recombination events. Analysis of the core and accessory genome also suggested that the different Y. enterocolitica phylogroups may be ecologically separated, in contrast to the long-held belief of common shared ecological niches across the Y. enterocolitica species.

© 2015 Elsevier Ltd. Background: Shigellosis is an acute, severe bacterial colitis that, in high-income countries, is typically associated with travel to high-risk regions (Africa, Asia, and Latin America). Since the 1970s, shigellosis has also been reported as a sexually transmitted infection in men who have sex with men (MSM), in whom transmission is an important component of shigellosis epidemiology in high-income nations. We aimed to use sophisticated subtyping and international sampling to determine factors driving shigellosis emergence in MSM linked to an outbreak in the UK. Methods: We did a large-scale, cross-sectional genomic epidemiological study of shigellosis cases collected from 29 countries between December, 1995, and June 8, 2014. Focusing on an ongoing epidemic in the UK, we collected and whole-genome sequenced clinical isolates of Shigella flexneri serotype 3a from high-risk and low-risk regions, including cases associated with travel and sex between men. We examined relationships between geographical, demographic, and clinical patient data with the isolate antimicrobial susceptibility, genetic data, and inferred evolutionary relationships. Findings: We obtained 331 clinical isolates of S flexneri serotype 3a, including 275 from low-risk regions (44 from individuals who travelled to high-risk regions), 52 from high-risk regions, and four outgroup samples (ie, closely related, but genetically distinct isolates used to determine the root of the phylogenetic tree). We identified a recently emerged lineage of S flexneri 3a that has spread intercontinentally in less than 20 years throughout regions traditionally at low risk for shigellosis via sexual transmission in MSM. The lineage had acquired multiple antimicrobial resistance determinants, and prevailing sublineages were strongly associated with resistance to the macrolide azithromycin. Eight (4%) of 206 isolates from the MSM-associated lineage were obtained from patients who had previously provided an isolate; these serial isolations indicated atypical infection patterns (eg, reinfection). Interpretation: We identified transmission-facilitating behaviours and atypical course(s) of infection as precipitating factors in shigellosis-affected MSM. The intercontinental spread of antimicrobial-resistant shigella through established transmission routes emphasises the need for new approaches to tackle the public health challenge of sexually transmitted infections in MSM. Funding: Wellcome Trust (grant number 098051).

The assembly of DNA sequence data is undergoing a renaissance thanks to emerging technologies capable of producing reads tens of kilobases long. Assembling complete bacterial and small eukaryotic genomes is now possible, but the final step of circularizing sequences remains unsolved. Here we present Circlator, the first tool to automate assembly circularization and produce accurate linear representations of circular sequences. Using Pacific Biosciences and Oxford Nanopore data, Circlator correctly circularized 26 of 27 circularizable sequences, comprising 11 chromosomes and 12 plasmids from bacteria, the apicoplast and mitochondrion of Plasmodium falciparum and a human mitochondrion. Circlator is available at http://sanger-pathogens.github.io/circlator/. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0849-0) contains supplementary material, which is available to authorized users.

Bordetella pertussis is the causative agent of whooping cough, a highly contagious, acute respiratory illness that has seen resurgence despite the use of vaccines. We present the complete genome sequence of a clinical strain of B. pertussis , D420, which is representative of a currently circulating clade of this pathogen.