Evolution of the U.S. Biological Select Agent Rathayibacter toxicus

Rathayibacter toxicus is a toxin-producing species found in Australia and is often fatal to grazing animals. The threat of introduction of the species into the United States led to its inclusion in the Federal Select Agent Program, which makes R. toxicus a highly regulated species. This work provides novel insights into the evolution of R. toxicus. R. toxicus is the only species in the genus to have acquired a CRISPR adaptive immune system to protect against bacteriophages. Results suggest that coexistence with the bacteriophage NCPPB3778 led to the massive shrinkage of the R. toxicus genome, species divergence, and the maintenance of low genetic diversity in extant bacterial groups. This work contributes to an understanding of the evolution and ecology of an agriculturally important species of bacteria. ABSTRACT Rathayibacter toxicus is a species of Gram-positive, corynetoxin-producing bacteria that causes annual ryegrass toxicity, a disease often fatal to grazing animals. A phylogenomic approach was employed to model the evolution of R. toxicus to explain the low genetic diversity observed among isolates collected during a 30-year period of sampling in three regions of Australia, gain insight into the taxonomy of Rathayibacter, and provide a framework for studying these bacteria. Analyses of a data set of more than 100 sequenced Rathayibacter genomes indicated that Rathayibacter forms nine species-level groups. R. toxicus is the most genetically distant, and evidence suggested that this species experienced a dramatic event in its evolution. Its genome is significantly reduced in size but is colinear to those of sister species. Moreover, R. toxicus has low intergroup genomic diversity and almost no intragroup genomic diversity between ecologically separated isolates. R. toxicus is the only species of the genus that encodes a clustered regularly interspaced short palindromic repeat (CRISPR) locus and that is known to host a bacteriophage parasite. The spacers, which represent a chronological history of infections, were characterized for information on past events. We propose a three-stage process that emphasizes the importance of the bacteriophage and CRISPR in the genome reduction and low genetic diversity of the R. toxicus species. IMPORTANCE Rathayibacter toxicus is a toxin-producing species found in Australia and is often fatal to grazing animals. The threat of introduction of the species into the United States led to its inclusion in the Federal Select Agent Program, which makes R. toxicus a highly regulated species. This work provides novel insights into the evolution of R. toxicus. R. toxicus is the only species in the genus to have acquired a CRISPR adaptive immune system to protect against bacteriophages. Results suggest that coexistence with the bacteriophage NCPPB3778 led to the massive shrinkage of the R. toxicus genome, species divergence, and the maintenance of low genetic diversity in extant bacterial groups. This work contributes to an understanding of the evolution and ecology of an agriculturally important species of bacteria.

[1]  R. Irizarry ggplot2 , 2019, Introduction to Data Science.

[2]  D. Luster,et al.  Rathayibacter agropyri (non O'Gara 1916) comb. nov., nom. rev., isolated from western wheatgrass (Pascopyrum smithii). , 2018, International journal of systematic and evolutionary microbiology.

[3]  Y. Ye,et al.  Reconstituting the History of Cronobacter Evolution Driven by Differentiated CRISPR Activity , 2018, Applied and Environmental Microbiology.

[4]  Peter C. Fineran,et al.  CRISPR-Cas-Mediated Phage Resistance Enhances Horizontal Gene Transfer by Transduction , 2018, mBio.

[5]  Elizabeth A. Savory,et al.  Evolutionary transitions between beneficial and phytopathogenic Rhodococcus challenge disease management , 2017, eLife.

[6]  W. Schneider,et al.  Complete Genome Sequence of Rathayibacter toxicus Phage NCPPB3778 , 2017, Genome Announcements.

[7]  C. Fennessey,et al.  Whole genome sequence of two Rathayibacter toxicus strains reveals a tunicamycin biosynthetic cluster similar to Streptomyces chartreusis , 2017, PloS one.

[8]  S. Subbotin,et al.  Rathayibacter toxicus, Other Rathayibacter Species Inducing Bacterial Head Blight of Grasses, and the Potential for Livestock Poisonings. , 2017, Phytopathology.

[9]  A. Rojas,et al.  A non-canonical mismatch repair pathway in prokaryotes , 2017, Nature Communications.

[10]  Niklaus J Grünwald,et al.  vcfr: a package to manipulate and visualize variant call format data in R , 2017, Molecular ecology resources.

[11]  Rolf Backofen,et al.  Characterizing leader sequences of CRISPR loci , 2016, Bioinform..

[12]  Javier F. Tabima,et al.  Gall-ID: tools for genotyping gall-causing phytopathogenic bacteria , 2016, PeerJ.

[13]  K. Datsenko,et al.  Highly efficient primed spacer acquisition from targets destroyed by the Escherichia coli type I-E CRISPR-Cas interfering complex , 2016, Proceedings of the National Academy of Sciences.

[14]  Peter C. Fineran,et al.  Inactivation of CRISPR-Cas systems by anti-CRISPR proteins in diverse bacterial species , 2016, Nature Microbiology.

[15]  Sanzhen Liu,et al.  Emergence of a New Population of Rathayibacter toxicus: An Ecologically Complex, Geographically Isolated Bacterium , 2016, PloS one.

[16]  A. Buckling,et al.  The diversity-generating benefits of a prokaryotic adaptive immune system , 2016, Nature.

[17]  Peer Bork,et al.  Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees , 2016, Nucleic Acids Res..

[18]  Christine L. Sun,et al.  Metagenomic reconstructions of bacterial CRISPR loci constrain population histories , 2015, The ISME Journal.

[19]  E. Charpentier,et al.  Adaptation in CRISPR-Cas Systems. , 2016, Molecular cell.

[20]  James Mallet,et al.  What Is Speciation? , 2016, PLoS genetics.

[21]  Hon Wai Leong,et al.  GI-SVM: A sensitive method for predicting genomic islands based on unannotated sequence of a single genome , 2016, J. Bioinform. Comput. Biol..

[22]  Rotem Sorek,et al.  CRISPR–Cas adaptation: insights into the mechanism of action , 2016, Nature Reviews Microbiology.

[23]  Dongwan D. Kang,et al.  Genome-wide selective sweeps and gene-specific sweeps in natural bacterial populations , 2016, The ISME Journal.

[24]  Sita J. Saunders,et al.  An updated evolutionary classification of CRISPR–Cas systems , 2015, Nature Reviews Microbiology.

[25]  S. Donaldson,et al.  Selective Sweeps and Parallel Pathoadaptation Drive Pseudomonas aeruginosa Evolution in the Cystic Fibrosis Lung , 2015, mBio.

[26]  Mike Boots,et al.  Parasite Exposure Drives Selective Evolution of Constitutive versus Inducible Defense , 2015, Current Biology.

[27]  Asaf Levy,et al.  CRISPR adaptation biases explain preference for acquisition of foreign DNA , 2015, Nature.

[28]  Eugene V Koonin,et al.  No evidence of inhibition of horizontal gene transfer by CRISPR–Cas on evolutionary timescales , 2015, The ISME Journal.

[29]  T. Lumley,et al.  gplots: Various R Programming Tools for Plotting Data , 2015 .

[30]  Jos Boekhorst,et al.  Degenerate target sites mediate rapid primed CRISPR adaptation , 2014, Proceedings of the National Academy of Sciences.

[31]  B. Koskella,et al.  Bacteria–phage coevolution as a driver of ecological and evolutionary processes in microbial communities , 2014, FEMS microbiology reviews.

[32]  Keith A. Jolley,et al.  A Reference Pan-Genome Approach to Comparative Bacterial Genomics: Identification of Novel Epidemiological Markers in Pathogenic Campylobacter , 2014, PloS one.

[33]  Ying Zhang,et al.  Pan-genome analyses identify lineage- and niche-specific markers of evolution and adaptation in Epsilonproteobacteria , 2014, Front. Microbiol..

[34]  Matthew Fraser,et al.  InterProScan 5: genome-scale protein function classification , 2014, Bioinform..

[35]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[36]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[37]  B. Contreras-Moreira,et al.  GET_HOMOLOGUES, a Versatile Software Package for Scalable and Robust Microbial Pangenome Analysis , 2013, Applied and Environmental Microbiology.

[38]  Aaron A. Klammer,et al.  Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data , 2013, Nature Methods.

[39]  Chris M. Brown,et al.  CRISPRTarget , 2013, RNA biology.

[40]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[41]  A. Aertsen,et al.  Phage–host interactions during pseudolysogeny , 2013, Bacteriophage.

[42]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[43]  K. Sneppen,et al.  Spatial Structure and Lamarckian Adaptation Explain Extreme Genetic Diversity at CRISPR Locus , 2012, mBio.

[44]  Liam J. Revell,et al.  phytools: an R package for phylogenetic comparative biology (and other things) , 2012 .

[45]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[46]  Otto X. Cordero,et al.  Population Genomics of Early Events in the Ecological Differentiation of Bacteria , 2012, Science.

[47]  Wayne M. Getz,et al.  Persisting Viral Sequences Shape Microbial CRISPR-based Immunity , 2012, PLoS Comput. Biol..

[48]  U. Qimron,et al.  Proteins and DNA elements essential for the CRISPR adaptation process in Escherichia coli , 2012, Nucleic acids research.

[49]  N. Moran,et al.  Extreme genome reduction in symbiotic bacteria , 2011, Nature Reviews Microbiology.

[50]  Tal Pupko,et al.  GLOOME: gain loss mapping engine , 2010, Bioinform..

[51]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[52]  N. Perna,et al.  progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement , 2010, PloS one.

[53]  L. Marraffini,et al.  CRISPR interference: RNA-directed adaptive immunity in bacteria and archaea , 2010, Nature Reviews Genetics.

[54]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[55]  Howard Ochman,et al.  Deletional Bias across the Three Domains of Life , 2009, Genome biology and evolution.

[56]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

[57]  C. Fraser,et al.  The Bacterial Species Challenge: Making Sense of Genetic and Ecological Diversity , 2009, Science.

[58]  N. L. Held,et al.  Viral biogeography revealed by signatures in Sulfolobus islandicus genomes. , 2009, Environmental microbiology.

[59]  L. Marraffini,et al.  CRISPR Interference Limits Horizontal Gene Transfer in Staphylococci by Targeting DNA , 2008, Science.

[60]  David R. Riley,et al.  Comparative genomics: the bacterial pan-genome. , 2008, Current opinion in microbiology.

[61]  Philippe Horvath,et al.  Diversity, Activity, and Evolution of CRISPR Loci in Streptococcus thermophilus , 2007, Journal of bacteriology.

[62]  J. Banfield,et al.  Rapidly evolving CRISPRs implicated in acquired resistance of microorganisms to viruses. , 2007, Environmental microbiology.

[63]  M. Glickman,et al.  Bacterial DNA repair by non-homologous end joining , 2007, Nature Reviews Microbiology.

[64]  B. Tsvetanova,et al.  Biosynthesis of the Tunicamycins: A Review , 2007, The Journal of Antibiotics.

[65]  Ibtissem Grissa,et al.  CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats , 2007, Nucleic Acids Res..

[66]  R. Barrangou,et al.  CRISPR Provides Acquired Resistance Against Viruses in Prokaryotes , 2007, Science.

[67]  D. Cahill,et al.  Development and application of polymerase chain reaction-based assays for Rathayibacter toxicus and a bacteriophage associated with annual ryegrass (Lolium rigidum) toxicity , 2007 .

[68]  J. R. Lobry,et al.  SeqinR 1.0-2: A Contributed Package to the R Project for Statistical Computing Devoted to Biological Sequences Retrieval and Analysis , 2007 .

[69]  P. Vandamme,et al.  DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. , 2007, International journal of systematic and evolutionary microbiology.

[70]  Robert S. Harris,et al.  Improved pairwise alignment of genomic dna , 2007 .

[71]  V. Kunin,et al.  Evolutionary conservation of sequence and secondary structures in CRISPR repeats , 2007, Genome Biology.

[72]  E. Postnikova,et al.  Genetic Characterization and Diversity of Rathayibacter toxicus. , 2006, Phytopathology.

[73]  K. Konstantinidis,et al.  Genomic insights that advance the species definition for prokaryotes. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[74]  Pavel A. Pevzner,et al.  De novo identification of repeat families in large genomes , 2005, ISMB.

[75]  E. Postnikova,et al.  Xylella fastidiosa subspecies: X. fastidiosa subsp. [correction] fastidiosa [correction] subsp. nov., X. fastidiosa subsp. multiplex subsp. nov., and X. fastidiosa subsp. pauca subsp. nov. , 2004, Systematic and applied microbiology.

[76]  Korbinian Strimmer,et al.  APE: Analyses of Phylogenetics and Evolution in R language , 2004, Bioinform..

[77]  Anders Gorm Pedersen,et al.  RevTrans: multiple alignment of coding DNA from aligned amino acid sequences , 2003, Nucleic Acids Res..

[78]  L. Hurst The Ka/Ks ratio: diagnosing the form of sequence evolution. , 2002, Trends in genetics : TIG.

[79]  F. Cohan What are bacterial species? , 2002, Annual review of microbiology.

[80]  I. Longden,et al.  EMBOSS: the European Molecular Biology Open Software Suite. , 2000, Trends in genetics : TIG.

[81]  R. Hendrix,et al.  Evolutionary relationships among diverse bacteriophages and prophages: all the world's a phage. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[82]  A. Kerr,et al.  Association of bacteriophage particles with toxin production by Clavibacter toxicus, the causal agent of annual ryegrass toxicity , 1993 .

[83]  J. Peterson,et al.  Inhibition of glycosylation by corynetoxin, the causative agent of annual ryegrass toxicity: a comparison with tunicamycin. , 1983, Chemico-biological interactions.

[84]  P. Kloot The Genus Lolium in Australia , 1983 .

[85]  C. Kado,et al.  Selective media for isolation of Agrobacterium, Corynebacterium, Erwinia, Pseudomonas, and Xanthomonas. , 1970, Phytopathology.

[86]  G. Bertani,et al.  STUDIES ON LYSOGENESIS I , 1951, Journal of bacteriology.