iBsu1103: a new genome-scale metabolic model of Bacillus subtilis based on SEED annotations

Background  

Bacillus subtilis is an organism of interest because of its extensive industrial applications, its similarity to pathogenic organisms, and its role as the model organism for Gram-positive, sporulating bacteria. In this work, we introduce a new genome-scale metabolic model of B. subtilis 168 called iBsu1103. This new model is based on the annotated B. subtilis 168 genome generated by the SEED, one of the most up-to-date and accurate annotations of B. subtilis 168 available.     

Genome‐scale metabolic modeling reveals key features of a minimal gene set

Mesoplasma florum, a fast‐growing near‐minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein‐coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species‐specific constraints to produce the functional iJL208 genome‐scale model (GEM) of metabolism. Genome‐wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI‐syn3.0 and its parent JCVI‐syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping‐stone toward model‐driven whole‐genome engineering.     

Metabolic characterization of the genus Brucella. VI. Growth stimulation by i-erythritol compared with strain virulence for guinea pigs

Strains of Brucella abortus isolated 20 to 38 years ago and strains recently isolated were assessed for their virulence to guinea pigs and for their ability to grow in half-strength Tryptose Broth with and without i-erythritol. The recently isolated strains were virulent and i-erythritol enhanced their growth. The aged strains were avirulent and grew equally well in both media. Three of the recently isolated strains were subcultured serially every 24 hr alternately on Tryptose Agar slants and in half-strength Tryptose Broth without i-erythritol. After 8 to 13 such transfers, the growth of each strain was equivalent in both media. The subculture on which growth equivalence occurred was retested for virulence. None of the three strains had decreased in its virulence for guinea pigs. The conclusion was drawn that strain virulence for guinea pigs and growth enhancement by i-erythritol are independent characteristics.     

Evaluation of selective media for isolation and enumeration of vibrios from estuarine waters

We have compared the effectiveness of four media developed in the last years together with the medium GSTC(glucose-salt-tellurite-crystal violet), devised in our laboratory, for the recovery of vibrios from estuarine waters. In addition, a number of reference Vibrio and non-Vibrio strains have been tested for growth on the five media. TCBS and GSTC were the most selective media for reference strains of Vibrio spp. However, when the media were tested with samples of water from three different sites of an estuary, only TCBS was effective enough to recover vibrios inhibiting the growth of non-Vibrio populations. We also report here the usefulness of TCBS for isolation of the fish pathogen V. anguillarum, since a total of 81 strains isolated from diseased fish and water in various parts of the world grew on this medium. In conclusion, we consider the TCBS as the best medium to isolate Vibrio species pathogenic for humans and fish, and for recovery of vibrios from estuarine waters.     

Symbiosis Island Shuffling with Abundant Insertion Sequences in the Genomes of Extra-Slow-Growing Strains of Soybean Bradyrhizobia

Extra-slow-growing bradyrhizobia from root nodules of field-grown soybeans harbor abundant insertion sequences (ISs) and are termed highly reiterated sequence-possessing (HRS) strains. We analyzed the genome organization of HRS strains with the focus on IS distribution and symbiosis island structure. Using pulsed-field gel electrophoresis, we consistently detected several plasmids (0.07 to 0.4 Mb) in the HRS strains (NK5, NK6, USDA135, 2281, USDA123, and T2), whereas no plasmids were detected in the non-HRS strain USDA110. The chromosomes of the six HRS strains (9.7 to 10.7 Mb) were larger than that of USDA110 (9.1 Mb). Using MiSeq sequences of 6 HRS and 17 non-HRS strains mapped to the USDA110 genome, we found that the copy numbers of ISRj1, ISRj2, ISFK1, IS1632, ISB27, ISBj8, and IS1631 were markedly higher in HRS strains. Whole-genome sequencing showed that the HRS strain NK6 had four small plasmids (136 to 212 kb) and a large chromosome (9,780 kb). Strong colinearity was found between 7.4-Mb core regions of the NK6 and USDA110 chromosomes. USDA110 symbiosis islands corresponded mainly to five small regions (S1 to S5) within two variable regions, V1 (0.8 Mb) and V2 (1.6 Mb), of the NK6 chromosome. The USDA110 nif gene cluster (nifDKENXSBZHQW-fixBCX) was split into two regions, S2 and S3, where ISRj1-mediated rearrangement occurred between nifS and nifB. ISs were also scattered in NK6 core regions, and ISRj1 insertion often disrupted some genes important for survival and environmental responses. These results suggest that HRS strains of soybean bradyrhizobia were subjected to IS-mediated symbiosis island shuffling and core genome degradation.     

Comparative genomics and evolution of the HSP90 family of genes across all kingdoms of organisms

Background

Vibrio Parahaemolyticus is an aquatic, halophilic, Gram-negative bacterium, first discovered in 1950 in Japan during a food-poisoning outbreak. Infections resulting from consumption of V. Parahaemolyticus have increased globally in the last 10 years leading to the bacterium's classification as a newly emerging pathogen. In 1996 the first appearance of a pandemic V. Parahaemolyticus clone occurred, a new O3:K6 serotype strain that has now been identified worldwide as a major cause of seafood-borne gastroenteritis.

Results

We examined the sequenced genome of V. Parahaemolyticus RIMD2210633, an O3:K6 serotype strain isolated in Japan in 1996, by bioinformatic analyses to uncover genomic islands (GIs) that may play a role in the emergence and pathogenesis of pandemic strains. We identified 7 regions ranging in size from 10 kb to 81 kb that had the characteristics of GIs such as aberrant base composition compared to the core genome, presence of phage-like integrases, flanked by direct repeats and the absence of these regions from closely related species. Molecular analysis of worldwide clinical isolates of V. Parahaemolyticus recovered over the last 33 years demonstrated that a 24 kb region named V. Parahaemolyticus island-1 (VPaI-1) encompassing ORFs VP0380 to VP0403 is only present in new O3:K6 and related strains recovered after 1995. We investigated the presence of 3 additional regions, VPaI-4 (VP2131 to VP2144), VPaI-5 (VP2900 to VP2910) and VPaI-6 (VPA1254 to VPA1270) by PCR assays and Southern blot analyses among the same set of V. Parahaemolyticus isolates. These 3 VPaI regions also gave similar distribution patterns amongst the 41 strains examined.

Conclusion

The 4 VPaI regions examined may represent DNA acquired by the pandemic group of V. Parahaemolyticus isolates that increased their fitness either in the aquatic environment or in their ability to infect humans.     

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: C_L1, C_L2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (C_L1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages C_L1, C_L2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages C_L1, C_L2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.     

Canthaxanthin production with modified Mucor circinelloides strains

Canthaxanthin is a natural diketo derivative of β-carotene primarily used by the food and feed industries. Mucor circinelloides is a β-carotene-accumulating zygomycete fungus and one of the model organisms to study the carotenoid biosynthesis in fungi. In this study, the β-carotene ketolase gene (crtW) of the marine bacterium Paracoccus sp. N81106 fused with fungal promoter and terminator regions was integrated into the M. circinelloides genome to construct stable canthaxanthin-producing strains. Different transformation methods including polyethylene glycol-mediated transformation with linear DNA fragments, restriction enzyme-mediated integration and Agrobacterium tumefaciens-mediated transformation were tested to integrate the crtW gene into the Mucor genome. Mitotic stability, site of integration and copy number of the transferred genes were analysed in the transformants, and several stable strains containing the crtW gene in high copy number were isolated. Carotenoid composition of selected transformants and effect of culturing conditions, such as temperature, carbon sources and application of certain additives in the culturing media, on their carotenoid content were analysed. Canthaxanthin-producing transformants were able to survive at higher growth temperature than the untransformed strain, maybe due to the effect of canthaxanthin on the membrane fluidity and integrity. With the application of glucose, trehalose, dihydroxyacetone and l-aspartic acid as sole carbon sources in minimal medium, the crtW-expressing M. circinelloides strain, MS12+pCA8lf/1, produced more than 200 μg/g (dry mass) of canthaxanthin.     

Comparative Genomic Analysis of Sulfurospirillum cavolei MES Reconstructed from the Metagenome of an Electrosynthetic Microbiome

Sulfurospirillum spp. play an important role in sulfur and nitrogen cycling, and contain metabolic versatility that enables reduction of a wide range of electron acceptors, including thiosulfate, tetrathionate, polysulfide, nitrate, and nitrite. Here we describe the assembly of a Sulfurospirillum genome obtained from the metagenome of an electrosynthetic microbiome. The ubiquity and persistence of this organism in microbial electrosynthesis systems suggest it plays an important role in reactor stability and performance. Understanding why this organism is present and elucidating its genetic repertoire provide a genomic and ecological foundation for future studies where Sulfurospirillum are found, especially in electrode-associated communities. Metabolic comparisons and in-depth analysis of unique genes revealed potential ecological niche-specific capabilities within the Sulfurospirillum genus. The functional similarities common to all genomes, i.e., core genome, and unique gene clusters found only in a single genome were identified. Based upon 16S rRNA gene phylogenetic analysis and average nucleotide identity, the Sulfurospirillum draft genome was found to be most closely related to Sulfurospirillum cavolei. Characterization of the draft genome described herein provides pathway-specific details of the metabolic significance of the newly described Sulfurospirillum cavolei MES and, importantly, yields insight to the ecology of the genus as a whole. Comparison of eleven sequenced Sulfurospirillum genomes revealed a total of 6246 gene clusters in the pan-genome. Of the total gene clusters, 18.5% were shared among all eleven genomes and 50% were unique to a single genome. While most Sulfurospirillum spp. reduce nitrate to ammonium, five of the eleven Sulfurospirillum strains encode for a nitrous oxide reductase (nos) cluster with an atypical nitrous-oxide reductase, suggesting a utility for this genus in reduction of the nitrous oxide, and as a potential sink for this potent greenhouse gas.     

Genetic and Phenotypic Comparison of Facultative Methylotrophy between Methylobacterium extorquens Strains PA1 and AM1

Methylobacterium extorquens AM1, a strain serendipitously isolated half a century ago, has become the best-characterized model system for the study of aerobic methylotrophy (the ability to grow on reduced single-carbon compounds). However, with 5 replicons and 174 insertion sequence (IS) elements in the genome as well as a long history of domestication in the laboratory, genetic and genomic analysis of M. extorquens AM1 face several challenges. On the contrary, a recently isolated strain - M. extorquens PA1- is closely related to M. extorquens AM1 (100% 16S rRNA identity) and contains a streamlined genome with a single replicon and only 20 IS elements. With the exception of the methylamine dehydrogenase encoding gene cluster (mau), genes known to be involved in methylotrophy are well conserved between M. extorquens AM1 and M. extorquens PA1. In this paper we report four primary findings regarding methylotrophy in PA1. First, with a few notable exceptions, the repertoire of methylotrophy genes between PA1 and AM1 is extremely similar. Second, PA1 grows faster with higher yields compared to AM1 on C₁ and multi-C substrates in minimal media, but AM1 grows faster in rich medium. Third, deletion mutants in PA1 throughout methylotrophy modules have the same C₁ growth phenotypes observed in AM1. Finally, the precision of our growth assays revealed several unexpected growth phenotypes for various knockout mutants that serve as leads for future work in understanding their basis and generality across Methylobacterium strains.     

Selection of phosphate-solubilizing diazotrophic Herbaspirillum and Burkholderia strains and their effect on rice crop yield and nutrient uptake

Background and Aims

Plant growth-promoting bacteria, mainly diazotrophs and phosphate solubilizers, can reduce the use of chemical fertilizers for rice crops. Here, diazotrophic bacteria isolated from rice were screened for their ability to solubilize inorganic P (P_i) in vitro and in association with rice plants cultivated in pots.

Methods

Forty-nine isolates were tested for the ability to solubilize P_i on NBRIP and GL agar plate media and seven selected strains were further evaluated in NBRIP liquid medium. Three of these strains were inoculated in rice plants grown in soil pots containing ¹⁵N-labeled fertilizer and two sources of P: tricalcium phosphate (TCP) or simple superphosphate (SSP). The dry matter, yield, N, P, and the ¹⁵N content accumulated in plant tissues were measured at 135 days after planting.

Results

Seven strains belonging to the genera Herbaspirillum and Burkholderia formed a halo of solubilized P_i on agar plates. The Burkholderia strains showed peak soluble P (around 200 mg P L^?1) on the fifth day when grown in NBRIP liquid medium for 14 days. Inoculation of Herbaspirillum strains (H18, ZA15) and a Burkholderia vietaminensis strain (AR114) increased rice grain yield from 33 to 47 % with TCP and 18 to 44 % with TSS, respectively. The bacterial inoculation led to enhanced N-use efficiency of the ¹⁵N-labeled fertilizer.

Conclusion

These results suggest that the selection and use of P-solubilizing diazotrophic bacteria are a good strategy to promote P solubilization and/or N use efficiency in rice plants.     

Drug resistance marker-aided genome shuffling to improve acetic acid tolerance in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

Acetic acid existing in a culture medium is one of the most limiting constraints in yeast growth and viability during ethanol fermentation. To improve acetic acid tolerance in Saccharomyces cerevisiae strains, a drug resistance marker-aided genome shuffling approach with higher screen efficiency of shuffled mutants was developed in this work. Through two rounds of genome shuffling of ultraviolet mutants derived from the original strain 308, we obtained a shuffled strain YZ2, which shows significantly faster growth and higher cell viability under acetic acid stress. Ethanol production of YZ2 (within 60 h) was 21.6% higher than that of 308 when 0.5% (v/v) acetic acid was added to fermentation medium. Membrane integrity, higher in vivo activity of the H⁺-ATPase, and lower oxidative damage after acetic acid treatment are the possible reasons for the acetic acid-tolerance phenotype of YZ2. These results indicated that this novel genome shuffling approach is powerful to rapidly improve the complex traits of industrial yeast strains.     

Osmotic Stress Leads to Decreased Intracellular pH of Listeria monocytogenes as Determined by Fluorescence Ratio-Imaging Microscopy

Intracellular pH (pH_i) of Listeria monocytogenes was determined after exposure to NaCl or sorbitol in liquid and solid media (agar). Both compounds decreased pH_i, and recovery on solid medium was impaired compared to that in liquid medium. N,N′-dicyclohexylcarbodiimide abolished pH_i recovery, and lowering a_w with glycerol showed no effect on pH_i.     

Frequency of the Insertion Sequence IS4Bsu1 among Bacillus subtilis Strains Isolated from Fermented Soybean Foods in Southeast Asia

Among 45 Bacillus subtilis strains isolated from non-salted types of fermented soybeans produced in several Southeast Asian countries, 20 had the insertion sequence IS4Bsu1 in the chromosome. In contrast, none of 49 B. subtilis strains of non-food origin contained IS4Bsu1. Frequent occurrence of this mobile DNA element in the soybean-fermenting B. subtilis would reflect the fact that few strains flourish on soybeans and thereby contribute to soybean fermentation.     

Diverse bacteria isolated from microtherm oil-production water

In total, 435 pure bacterial strains were isolated from microtherm oil-production water from the Karamay Oilfield, Xinjiang, China, by using four media: oil-production water medium (Cai medium), oil-production water supplemented with mineral salt medium (CW medium), oil-production water supplemented with yeast extract medium (CY medium), and blood agar medium (X medium). The bacterial isolates were affiliated with 61 phylogenetic groups that belong to 32 genera in the phyla Actinobacteria, Firmicutes, and Proteobacteria. Except for the Rhizobium, Dietzia, and Pseudomonas strains that were isolated using all the four media, using different media led to the isolation of bacteria with different functions. Similarly, nonheme diiron alkane monooxygenase genes (alkB/alkM) also clustered according to the isolation medium. Among the bacterial strains, more than 24 % of the isolates could use n-hexadecane as the sole carbon source for growth. For the first time, the alkane-degrading ability and alkB/alkM were detected in Rhizobium, Rhodobacter, Trichococcus, Micrococcus, Enterococcus, and Bavariicoccus strains, and the alkM gene was detected in Firmicutes strains.     

Filamentous Bacteriophages of Vibrio parahaemolyticus as a Possible Clue to Genetic Transmission

We have previously reported the isolation and characterization of two filamentous bacteriophages of Vibrio parahaemolyticus, designated Vf12 and Vf33. In this study, to understand the potential of these phages as tools for genetic transmission, we investigated the gene structures of replicative-form (RF) DNAs of their genomes and the distribution of these DNAs on chromosomal and extrachromosomal DNAs. The 7,965-bp nucleotide sequences of Vf12 and Vf33 were determined. An analysis of the overall gene structures revealed that Vf12 and Vf33 had conserved regions and distinctive regions. The gene organization of their conserved regions was similar to that of CTX phage of Vibrio cholerae and coliphage Ff of Escherichia coli, while their distinctive regions were characteristic of Vf12 and Vf33 phage genomes. Southern blot hybridization testing revealed that the filamentous phage genomes integrated into chromosomal DNA of V. parahaemolyticus at the distinctive region of the phage genome and were also distributed on some plasmids of V. parahaemolyticus and total cellular DNAs of one Vibrio damsela and one nonagglutinable Vibrio strain tested. These results strongly suggest the possibilities of genetic interaction among the bacteriophage Vf12 and Vf33 genomes and chromosomal and plasmid-borne DNAs of V. parahaemolyticus strains and of genetic transmission among strains through these filamentous phages.     

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

Background

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

Results

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

Conclusions

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-720) contains supplementary material, which is available to authorized users.