Building the repertoire of dispensable chromosome regions in Bacillus subtilis entails major refinement of cognate large-scale metabolic model

The nonessential regions in bacterial chromosomes are ill-defined due to incomplete functional information. Here, we establish a comprehensive repertoire of the genome regions that are dispensable for growth of Bacillus subtilis in a variety of media conditions. In complex medium, we attempted deletion of 157 individual regions ranging in size from 2 to 159 kb. A total of 146 deletions were successful in complex medium, whereas the remaining regions were subdivided to identify new essential genes (4) and coessential gene sets (7). Overall, our repertoire covers ∼76% of the genome. We screened for viability of mutant strains in rich defined medium and glucose minimal media. Experimental observations were compared with predictions by the iBsu1103 model, revealing discrepancies that led to numerous model changes, including the large-scale application of model reconciliation techniques. We ultimately produced the iBsu1103V2 model and generated predictions of metabolites that could restore the growth of unviable strains. These predictions were experimentally tested and demonstrated to be correct for 27 strains, validating the refinements made to the model. The iBsu1103V2 model has improved considerably at predicting loss of viability, and many insights gained from the model revisions have been integrated into the Model SEED to improve reconstruction of other microbial models.

[1]  M. Itaya,et al.  Bottom-up genome assembly using the Bacillus subtilis genome vector , 2008, Nature Methods.

[2]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

[3]  Nicola Zamboni,et al.  Genome engineering reveals large dispensable regions in Bacillus subtilis. , 2003, Molecular biology and evolution.

[4]  C. Walsh,et al.  Investigation of anticapsin biosynthesis reveals a four-enzyme pathway to tetrahydrotyrosine in Bacillus subtilis. , 2010, Biochemistry.

[5]  Vinay Satish Kumar,et al.  GrowMatch: An Automated Method for Reconciling In Silico/In Vivo Growth Predictions , 2009, PLoS Comput. Biol..

[6]  Satoru Miyano,et al.  Prediction of Transcriptional Terminators in Bacillus subtilis and Related Species , 2005, PLoS Comput. Biol..

[7]  S. Kanaya,et al.  Enhanced Recombinant Protein Productivity by Genome Reduction in Bacillus subtilis , 2008, DNA research : an international journal for rapid publication of reports on genes and genomes.

[8]  A. van Loon,et al.  Regulation of Riboflavin Biosynthesis inBacillus subtilis Is Affected by the Activity of the Flavokinase/Flavin Adenine Dinucleotide Synthetase Encoded byribC , 1998, Journal of bacteriology.

[9]  J. Hoch,et al.  Characterization of the gene for a protein kinase which phosphorylates the sporulation-regulatory proteins Spo0A and Spo0F of Bacillus subtilis , 1989, Journal of bacteriology.

[10]  S. Altman,et al.  Reaction in vitro of some mutants of RNase P with wild-type and temperature-sensitive substrates. , 1989, Journal of molecular biology.

[11]  Rick L. Stevens,et al.  Building the blueprint of life , 2010, Biotechnology journal.

[12]  N. Pace,et al.  In vitro selection for altered divalent metal specificity in the RNase P RNA. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Ali R. Zomorrodi,et al.  Genome-scale gene/reaction essentiality and synthetic lethality analysis , 2009, Molecular systems biology.

[14]  G. Wagner,et al.  EVOLUTION AND DETECTION OF GENETIC ROBUSTNESS , 2003 .

[15]  S. Ehrlich,et al.  Essential Bacillus subtilis genes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[16]  C. Anagnostopoulos,et al.  REQUIREMENTS FOR TRANSFORMATION IN BACILLUS SUBTILIS , 1961, Journal of bacteriology.

[17]  Guy Plunkett,et al.  Engineering a reduced Escherichia coli genome. , 2002, Genome research.

[18]  P. Noirot,et al.  A new mutation delivery system for genome‐scale approaches in Bacillus subtilis , 2002, Molecular microbiology.

[19]  Y. Fujita,et al.  Specific recognition of the Bacillus subtilis gnt cis‐acting catabolite‐responsive element by a protein complex formed between CcpA and seryl‐phosphorylated HPr , 1995, Molecular microbiology.

[20]  Matthew D. Jankowski,et al.  Group contribution method for thermodynamic analysis of complex metabolic networks. , 2008, Biophysical journal.

[21]  B. Palsson,et al.  Genome-scale Reconstruction of Metabolic Network in Bacillus subtilis Based on High-throughput Phenotyping and Gene Essentiality Data* , 2007, Journal of Biological Chemistry.

[22]  Jun Hyoung Lee,et al.  Metabolic engineering of a reduced-genome strain of Escherichia coli for L-threonine production , 2009, Microbial cell factories.

[23]  D. Henner,et al.  The organization and nucleotide sequence of the Bacillus subtilis hisH, tyrA and aroE genes. , 1986, Gene.

[24]  Naotake Ogasawara,et al.  Combined Effect of Improved Cell Yield and Increased Specific Productivity Enhances Recombinant Enzyme Production in Genome-Reduced Bacillus subtilis Strain MGB874 , 2011, Applied and Environmental Microbiology.

[25]  U. Sauer,et al.  Large-scale in vivo flux analysis shows rigidity and suboptimal performance of Bacillus subtilis metabolism , 2005, Nature Genetics.

[26]  Jun Hyoung Lee,et al.  Minimization of the Escherichia coli genome using a Tn5-targeted Cre/loxP excision system , 2002, Nature Biotechnology.

[27]  Vinay Satish Kumar,et al.  Optimization based automated curation of metabolic reconstructions , 2007, BMC Bioinformatics.

[28]  Hiroshi Mizoguchi,et al.  Cell size and nucleoid organization of engineered Escherichia coli cells with a reduced genome , 2004, Molecular microbiology.

[29]  Rick L Stevens,et al.  iBsu1103: a new genome-scale metabolic model of Bacillus subtilis based on SEED annotations , 2009, Genome Biology.

[30]  Rick L. Stevens,et al.  High-throughput generation, optimization and analysis of genome-scale metabolic models , 2010, Nature Biotechnology.

[31]  R. Hartmann,et al.  Analysis of RNase P Protein (rnpA) Expression in Bacillus subtilis Utilizing Strains with Suppressible rnpA Expression , 2006, Journal of bacteriology.

[32]  Jörg Stülke,et al.  Connecting parts with processes: SubtiWiki and SubtiPathways integrate gene and pathway annotation for Bacillus subtilis. , 2010, Microbiology.

[33]  Matthew D. Jankowski,et al.  Genome-scale thermodynamic analysis of Escherichia coli metabolism. , 2006, Biophysical journal.

[34]  Naryttza N. Diaz,et al.  The Subsystems Approach to Genome Annotation and its Use in the Project to Annotate 1000 Genomes , 2005, Nucleic acids research.

[35]  Sven Panke,et al.  Microbial systems engineering: First successes and the way ahead , 2010, BioEssays : news and reviews in molecular, cellular and developmental biology.

[36]  Kenta Nakai,et al.  DBTBS: a database of transcriptional regulation in Bacillus subtilis containing upstream intergenic conservation information , 2007, Nucleic Acids Res..

[37]  M. Itaya,et al.  Combining two genomes in one cell: stable cloning of the Synechocystis PCC6803 genome in the Bacillus subtilis 168 genome. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Insertion of Unmarked DNA Sequences in Multiple Loci of the Bacillus subtilis 168 Genome: an Efficient Selection Method , 2005, Bioscience, biotechnology, and biochemistry.

[39]  Costas D. Maranas,et al.  Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data , 2010, BMC Systems Biology.

[40]  Luis Serrano,et al.  Synthetic biology: promises and challenges , 2007, Molecular systems biology.

[41]  E. Raineri,et al.  Evolvability and hierarchy in rewired bacterial gene networks , 2008, Nature.

[42]  A. Danchin,et al.  From a consortium sequence to a unified sequence: the Bacillus subtilis 168 reference genome a decade later , 2009, Microbiology.

[43]  S. Gerdes,et al.  Discovery and Characterization of HemQ , 2010, The Journal of Biological Chemistry.

[44]  B. Palsson,et al.  Transcriptional regulation in constraints-based metabolic models of Escherichia coli Covert , 2002 .

[45]  H. Jenkinson Altered arrangement of proteins in the spore coat of a germination mutant of Bacillus subtilis. , 1983, Journal of general microbiology.

[46]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[47]  S Rozen,et al.  Primer3 on the WWW for general users and for biologist programmers. , 2000, Methods in molecular biology.

[48]  J. Stelling,et al.  Robustness of Cellular Functions , 2004, Cell.

[49]  Andrew R. Joyce,et al.  Experimental and Computational Assessment of Conditionally Essential Genes in Escherichia coli , 2006, Journal of bacteriology.

[50]  G. Hong,et al.  Nucleic Acids Research , 2015, Nucleic Acids Research.

[51]  Priscilla E. M. Purnick,et al.  The second wave of synthetic biology: from modules to systems , 2009, Nature Reviews Molecular Cell Biology.

[52]  F. Blattner,et al.  Emergent Properties of Reduced-Genome Escherichia coli , 2006, Science.