Derek M. Harkins

A novel method of consensus pan-chromosome assembly and large-scale comparative analysis reveal the highly flexible pan-genome of Acinetobacter baumannii

BackgroundInfections by pan-drug resistant Acinetobacter baumannii plague military and civilian healthcare systems. Previous A. baumannii pan-genomic studies used modest sample sizes of low diversity and comparisons to a single reference genome, limiting our understanding of gene order and content. A consensus representation of multiple genomes will provide a better framework for comparison. A large-scale comparative study will identify genomic determinants associated with their diversity and adaptation as a successful pathogen.ResultsWe determine draft-level genomic sequence of 50 diverse military isolates and conduct the largest bacterial pan-genome analysis of 249 genomes. The pan-genome of A. baumannii is open when the input genomes are normalized for diversity with 1867 core proteins and a paralog-collapsed pan-genome size of 11,694 proteins. We developed a novel graph-based algorithm and use it to assemble the first consensus pan-chromosome, identifying both the order and orientation of core genes and flexible genomic regions. Comparative genome analyses demonstrate the existence of novel resistance islands and isolates with increased numbers of resistance island insertions over time, from single insertions in the 1950s to triple insertions in 2011. Gene clusters responsible for carbon utilization, siderophore production, and pilus assembly demonstrate frequent gain or loss among isolates.ConclusionsThe highly variable and dynamic nature of the A. baumannii genome may be the result of its success in rapidly adapting to both abiotic and biotic environments through the gain and loss of gene clusters controlling fitness. Importantly, some archaic adaptation mechanisms appear to have reemerged among recent isolates.

TIGRFAMs and Genome Properties in 2013

TIGRFAMs, available online at http://www.jcvi.org/tigrfams is a database of protein family definitions. Each entry features a seed alignment of trusted representative sequences, a hidden Markov model (HMM) built from that alignment, cutoff scores that let automated annotation pipelines decide which proteins are members, and annotations for transfer onto member proteins. Most TIGRFAMs models are designated equivalog, meaning they assign a specific name to proteins conserved in function from a common ancestral sequence. Models describing more functionally heterogeneous families are designated subfamily or domain, and assign less specific but more widely applicable annotations. The Genome Properties database, available at http://www.jcvi.org/genome-properties, specifies how computed evidence, including TIGRFAMs HMM results, should be used to judge whether an enzymatic pathway, a protein complex or another type of molecular subsystem is encoded in a genome. TIGRFAMs and Genome Properties content are developed in concert because subsystems reconstruction for large numbers of genomes guides selection of seed alignment sequences and cutoff values during protein family construction. Both databases specialize heavily in bacterial and archaeal subsystems. At present, 4284 models appear in TIGRFAMs, while 628 systems are described by Genome Properties. Content derives both from subsystem discovery work and from biocuration of the scientific literature.

Complete Genome Sequence of the Multiresistant Taxonomic Outlier Pseudomonas aeruginosa PA7

Pseudomonas aeruginosa PA7 is a non-respiratory human isolate from Argentina that is multiresistant to antibiotics. We first sequenced gyrA, gyrB, parC, parE, ampC, ampR, and several housekeeping genes and found that PA7 is a taxonomic outlier. We report here the complete sequence of the 6,588,339 bp genome, which has only about 95% overall identity to other strains. PA7 has multiple novel genomic islands and a total of 51 occupied regions of genomic plasticity. These islands include antibiotic resistance genes, parts of transposons, prophages, and a pKLC102-related island. Several PA7 genes not present in PAO1 or PA14 are putative orthologues of other Pseudomonas spp. and Ralstonia spp. genes. PA7 appears to be closely related to the known taxonomic outlier DSM1128 (ATCC9027). PA7 lacks several virulence factors, notably the entire TTSS region corresponding to PA1690-PA1725 of PAO1. It has neither exoS nor exoU and lacks toxA, exoT, and exoY. PA7 is serotype O12 and pyoverdin type II. Preliminary proteomic studies indicate numerous differences with PAO1, some of which are probably a consequence of a frameshift mutation in the mvfR quorum sensing regulatory gene.

New Insights into Dissemination and Variation of the Health Care-Associated Pathogen Acinetobacter baumannii from Genomic Analysis

ABSTRACT Acinetobacter baumannii is a globally important nosocomial pathogen characterized by an increasing incidence of multidrug resistance. Routes of dissemination and gene flow among health care facilities are poorly resolved and are important for understanding the epidemiology of A. baumannii, minimizing disease transmission, and improving patient outcomes. We used whole-genome sequencing to assess diversity and genome dynamics in 49 isolates from one United States hospital system during one year from 2007 to 2008. Core single-nucleotide-variant-based phylogenetic analysis revealed multiple founder strains and multiple independent strains recovered from the same patient yet was insufficient to fully resolve strain relationships, where gene content and insertion sequence patterns added additional discriminatory power. Gene content comparisons illustrated extensive and redundant antibiotic resistance gene carriage and direct evidence of gene transfer, recombination, gene loss, and mutation. Evidence of barriers to gene flow among hospital components was not found, suggesting complex mixing of strains and a large reservoir of A. baumannii strains capable of colonizing patients. IMPORTANCE Genome sequencing was used to characterize multidrug-resistant Acinetobacter baumannii strains from one United States hospital system during a 1-year period to better understand how A. baumannii strains that cause infection are related to one another. Extensive variation in gene content was found, even among strains that were very closely related phylogenetically and epidemiologically. Several mechanisms contributed to this diversity, including transfer of mobile genetic elements, mobilization of insertion sequences, insertion sequence-mediated deletions, and genome-wide homologous recombination. Variation in gene content, however, lacked clear spatial or temporal patterns, suggesting a diverse pool of circulating strains with considerable interaction between strains and hospital locations. Widespread genetic variation among strains from the same hospital and even the same patient, particularly involving antibiotic resistance genes, reinforces the need for molecular diagnostic testing and genomic analysis to determine resistance profiles, rather than a reliance primarily on strain typing and antimicrobial resistance phenotypes for epidemiological studies. Genome sequencing was used to characterize multidrug-resistant Acinetobacter baumannii strains from one United States hospital system during a 1-year period to better understand how A. baumannii strains that cause infection are related to one another. Extensive variation in gene content was found, even among strains that were very closely related phylogenetically and epidemiologically. Several mechanisms contributed to this diversity, including transfer of mobile genetic elements, mobilization of insertion sequences, insertion sequence-mediated deletions, and genome-wide homologous recombination. Variation in gene content, however, lacked clear spatial or temporal patterns, suggesting a diverse pool of circulating strains with considerable interaction between strains and hospital locations. Widespread genetic variation among strains from the same hospital and even the same patient, particularly involving antibiotic resistance genes, reinforces the need for molecular diagnostic testing and genomic analysis to determine resistance profiles, rather than a reliance primarily on strain typing and antimicrobial resistance phenotypes for epidemiological studies.

Insights into the Environmental Resistance Gene Pool from the Genome Sequence of the Multidrug-Resistant Environmental Isolate Escherichia coli SMS-3-5

The increasing occurrence of multidrug-resistant pathogens of clinical and agricultural importance is a global public health concern. While antimicrobial use in human and veterinary medicine is known to contribute to the dissemination of antimicrobial resistance, the impact of microbial communities and mobile resistance genes from the environment in this process is not well understood. Isolated from an industrially polluted aquatic environment, Escherichia coli SMS-3-5 is resistant to a record number of antimicrobial compounds from all major classes, including two front-line fluoroquinolones (ciprofloxacin and moxifloxacin), and in many cases at record-high concentrations. To gain insights into antimicrobial resistance in environmental bacterial populations, the genome of E. coli SMS-3-5 was sequenced and compared to the genome sequences of other E. coli strains. In addition, selected genetic loci from E. coli SMS-3-5 predicted to be involved in antimicrobial resistance were phenotypically characterized. Using recombinant vector clones from shotgun sequencing libraries, resistance to tetracycline, streptomycin, and sulfonamide/trimethoprim was assigned to a single mosaic region on a 130-kb plasmid (pSMS35_130). The remaining plasmid backbone showed similarity to virulence plasmids from avian-pathogenic E. coli (APEC) strains. Individual resistance gene cassettes from pSMS35_130 are conserved among resistant bacterial isolates from multiple phylogenetic and geographic sources. Resistance to quinolones was assigned to several chromosomal loci, mostly encoding transport systems that are also present in susceptible E. coli isolates. Antimicrobial resistance in E. coli SMS-3-5 is therefore dependent both on determinants acquired from a mobile gene pool that is likely available to clinical and agricultural pathogens, as well, and on specifically adapted multidrug efflux systems. The association of antimicrobial resistance with APEC virulence genes on pSMS35_130 highlights the risk of promoting the spread of virulence through the extensive use of antibiotics.

Whole Genome Analysis of Leptospira licerasiae Provides Insight into Leptospiral Evolution and Pathogenicity

The whole genome analysis of two strains of the first intermediately pathogenic leptospiral species to be sequenced (Leptospira licerasiae strains VAR010 and MMD0835) provides insight into their pathogenic potential and deepens our understanding of leptospiral evolution. Comparative analysis of eight leptospiral genomes shows the existence of a core leptospiral genome comprising 1547 genes and 452 conserved genes restricted to infectious species (including L. licerasiae) that are likely to be pathogenicity-related. Comparisons of the functional content of the genomes suggests that L. licerasiae retains several proteins related to nitrogen, amino acid and carbohydrate metabolism which might help to explain why these Leptospira grow well in artificial media compared with pathogenic species. L. licerasiae strains VAR010T and MMD0835 possess two prophage elements. While one element is circular and shares homology with LE1 of L. biflexa, the second is cryptic and homologous to a previously identified but unnamed region in L. interrogans serovars Copenhageni and Lai. We also report a unique O-antigen locus in L. licerasiae comprised of a 6-gene cluster that is unexpectedly short compared with L. interrogans in which analogous regions may include >90 such genes. Sequence homology searches suggest that these genes were acquired by lateral gene transfer (LGT). Furthermore, seven putative genomic islands ranging in size from 5 to 36 kb are present also suggestive of antecedent LGT. How Leptospira become naturally competent remains to be determined, but considering the phylogenetic origins of the genes comprising the O-antigen cluster and other putative laterally transferred genes, L. licerasiae must be able to exchange genetic material with non-invasive environmental bacteria. The data presented here demonstrate that L. licerasiae is genetically more closely related to pathogenic than to saprophytic Leptospira and provide insight into the genomic bases for its infectiousness and its unique antigenic characteristics.

Pathogenomic inference of virulence-associated genes in Leptospira interrogans.

Leptospirosis is a globally important, neglected zoonotic infection caused by spirochetes of the genus Leptospira. Since genetic transformation remains technically limited for pathogenic Leptospira, a systems biology pathogenomic approach was used to infer leptospiral virulence genes by whole genome comparison of culture-attenuated Leptospira interrogans serovar Lai with its virulent, isogenic parent. Among the 11 pathogen-specific protein-coding genes in which non-synonymous mutations were found, a putative soluble adenylate cyclase with host cell cAMP-elevating activity, and two members of a previously unstudied ∼15 member paralogous gene family of unknown function were identified. This gene family was also uniquely found in the alpha-proteobacteria Bartonella bacilliformis and Bartonella australis that are geographically restricted to the Andes and Australia, respectively. How the pathogenic Leptospira and these two Bartonella species came to share this expanded gene family remains an evolutionary mystery. In vivo expression analyses demonstrated up-regulation of 10/11 Leptospira genes identified in the attenuation screen, and profound in vivo, tissue-specific up-regulation by members of the paralogous gene family, suggesting a direct role in virulence and host-pathogen interactions. The pathogenomic experimental design here is generalizable as a functional systems biology approach to studying bacterial pathogenesis and virulence and should encourage similar experimental studies of other pathogens.

Draft Genome Sequence of Turicibacter sanguinis PC909, Isolated from Human Feces

While the microbiota resident in the human gut is now known to provide a range of functions relevant to host health, many of the microbial members of the community have not yet been cultured or are represented by a limited number of isolates. We describe here the draft genome sequence of Turicibacter sanguinis PC909, isolated from a pooled healthy human fecal sample as part of the Australian Human Gut Microbiome Project.

Stool microbiota composition is associated with the prospective risk of Plasmodium falciparum infection.

BackgroundIn humans it is unknown if the composition of the gut microbiota alters the risk of Plasmodium falciparum infection or the risk of developing febrile malaria once P. falciparum infection is established. Here we collected stool samples from a cohort composed of 195 Malian children and adults just prior to an intense P. falciparum transmission season. We assayed these samples using massively parallel sequencing of the 16S ribosomal RNA gene to identify the composition of the gut bacterial communities in these individuals. During the ensuing 6-month P. falciparum transmission season we examined the relationship between the stool microbiota composition of individuals in this cohort and their prospective risk of both P. falciparum infection and febrile malaria.ResultsConsistent with prior studies, stool microbial diversity in the present cohort increased with age, although the overall microbiota profile was distinct from cohorts in other regions of Africa, Asia and North America. Age-adjusted Cox regression analysis revealed a significant association between microbiota composition and the prospective risk of P. falciparum infection; however, no relationship was observed between microbiota composition and the risk of developing febrile malaria once P. falciparum infection was established.ConclusionsThese findings underscore the diversity of gut microbiota across geographic regions, and suggest that strategic modulation of gut microbiota composition could decrease the risk of P. falciparum infection in malaria-endemic areas, potentially as an adjunct to partially effective malaria vaccines.

CharProtDB: a database of experimentally characterized protein annotations

CharProtDB (http://www.jcvi.org/charprotdb/) is a curated database of biochemically characterized proteins. It provides a source of direct rather than transitive assignments of function, designed to support automated annotation pipelines. The initial data set in CharProtDB was collected through manual literature curation over the years by analysts at the J. Craig Venter Institute (JCVI) [formerly The Institute of Genomic Research (TIGR)] as part of their prokaryotic genome sequencing projects. The CharProtDB has been expanded by import of selected records from publicly available protein collections whose biocuration indicated direct rather than homology-based assignment of function. Annotations in CharProtDB include gene name, symbol and various controlled vocabulary terms, including Gene Ontology terms, Enzyme Commission number and TransportDB accession. Each annotation is referenced with the source; ideally a journal reference, or, if imported and lacking one, the original database source.