Comparative genomic analysis of the gut bacterium Bifidobacterium longum reveals loci susceptible to deletion during pure culture growth
© Lee et al; licensee BioMed Central Ltd. 2008
Received: 16 November 2007
Accepted: 27 May 2008
Published: 27 May 2008
Bifidobacteria are frequently proposed to be associated with good intestinal health primarily because of their overriding dominance in the feces of breast fed infants. However, clinical feeding studies with exogenous bifidobacteria show they don't remain in the intestine, suggesting they may lose competitive fitness when grown outside the gut.
To further the understanding of genetic attenuation that may be occurring in bifidobacteria cultures, we obtained the complete genome sequence of an intestinal isolate, Bifidobacterium longum DJO10A that was minimally cultured in the laboratory, and compared it to that of a culture collection strain, B. longum NCC2705. This comparison revealed colinear genomes that exhibited high sequence identity, except for the presence of 17 unique DNA regions in strain DJO10A and six in strain NCC2705. While the majority of these unique regions encoded proteins of diverse function, eight from the DJO10A genome and one from NCC2705, encoded gene clusters predicted to be involved in diverse traits pertinent to the human intestinal environment, specifically oligosaccharide and polyol utilization, arsenic resistance and lantibiotic production. Seven of these unique regions were suggested by a base deviation index analysis to have been precisely deleted from strain NCC2705 and this is substantiated by a DNA remnant from within one of the regions still remaining in the genome of NCC2705 at the same locus. This targeted loss of genomic regions was experimentally validated when growth of the intestinal B. longum in the laboratory for 1,000 generations resulted in two large deletions, one in a lantibiotic encoding region, analogous to a predicted deletion event for NCC2705. A simulated fecal growth study showed a significant reduced competitive ability of this deletion strain against Clostridium difficile and E. coli. The deleted region was between two IS30 elements which were experimentally demonstrated to be hyperactive within the genome. The other deleted region bordered a novel class of mobile elements, termed mobile integrase cassettes (MIC) substantiating the likely role of these elements in genome deletion events.
Deletion of genomic regions, often facilitated by mobile elements, allows bifidobacteria to adapt to fermentation environments in a very rapid manner (2 genome deletions per 1,000 generations) and the concomitant loss of possible competitive abilities in the gut.
Recent molecular studies into the microbial diversity of the human intestine reveal a much greater diversity than previously recognized and very little is currently known of the contribution of individual groups to the human organism . One numerically dominant group of microbes, the bifidobacteria, is often suggested to be associated with good intestinal health given their overriding dominance in the feces of breast fed infants . This phenomenon led to their discovery in 1899 by the pediatrician Henri Tissier and his subsequent use of these bacteria for the treatment of infantile diarrhea . The proposed beneficial effect of bifidobacteria is further supported by the decrease of these bacteria in geriatric individuals and the concomitant increase of other microbial groups, most notably clostridia and E. coli [4–6]. This has led to the growing worldwide interest of including bifidobacteria in foods specifically for their potential intestinal health benefits . However, clinical feeding studies with bifidobacteria show that while the strains can be detected in subject's feces during feeding trials, they are rapidly lost upon cessation of the studies pointing to a possible loss of competitive fitness of the strains for competition within the human intestinal environment [7–9]. This may be due to attenuation of the strains, as the fermentation environment is very different to the buffered and anaerobic environment of the human colon.
To further the understanding of genetic attenuation that may occur in bifidobacteria, the complete genomic sequence of a numerically dominant human intestinal isolate of Bifidobacterium longum, that was grown for less than 20 generations in laboratory media, was deciphered. This was compared to an available genome sequence of a culture collection strain, B. longum NCC2705  and analysis of the functional attributes of its unique sequences have contributed to a better understanding of attenuation that can occur in these bacteria in a fermentation environment.
Results and Discussion
Genomic sequencing of a minimally cultured B. longum strain
General characteristics of the B. longum DJO10A genome
Overall characteristics of the genomes of B. longum strains DJO10A and NCC2705.
Size of chromosome (bp)
Overall G+C %
Number of plasmids
2 (10 and 3.6 kb)
1 (3.6 kb)
Average gene length (bp)
Gene density (genes/kb)
Gene coding percentage (%)
Gene G+C %
Strain-specific unique genes
Number of unique regionsa
Number of genes in unique regions
Number of genes in prophage region
RNAs and Repeat Sequences
Mobile Elements b
Mobile integrase cassette (MIC)
Organization of the origin and terminus of replication
An oriC and terC were found in identical locations in the genome of strain DJO10A and the updated genome sequence of strain NCC2705 (Additional file 4). These regions are extremely highly conserved in both genomes (> 99.9% identity) and consist of three oriC clusters and a terC, which is consistent with the predicted replication regions from other bacterial genomes . However, the location of the observed oriC region in both genomes is slightly different from the predicted location based on genome asymmetry, a feature that has previously been seen in the Helicobacter pylori 26695 genome [16, 17]. As well as the multiple oriC clusters, there are 7 different types of DnaA boxes, consistent with the majority of sequenced genomes and are proposed to be involved in controlling initiation of chromosome replication .
Restriction and modification (R-M) systems
The protective role that R-M systems impart on bacteria has been compared to the immune system of higher organisms . The presence of these systems in numerous bacteria demonstrates their important role for bacterial survival in nature. Both of the B. longum genomes encode type I and two type II R-M systems that are highly conserved (Additional file 5). They also contain an Mrr system that is predicted to restrict methylated DNA (usually Hha II or Pst I methylated DNA) that is 100% conserved between both strains (Additional file 5A). The clustering of Mrr with the type I R-M system is similar to E. coli K12 (GenBank U00096). The low identity (40%) between the HsdS proteins in the two strains likely reflects the independent evolution of this type I R-M system in these strains following their evolutionary divergence, as these systems evolve by changing their specificity components (HsdS) to enable it to recognize different sequences. This is substantiated by the existence of an hsdS gene that was inactivated by an IS256 insertion event and both parts of this disrupted gene exhibit much higher conservation, suggesting the insertion event occurred before their evolutionary divergence (Additional file 5A). Upstream from this locus in strain DJO10A there is another restriction gene, McrA (restricts DNA methylated by Hpa II or Sss I), that is not present in NCC2705. The conserved type II R-M systems in both strains are isoschizomers of Sau 3AI and Eco RII which restrict DNA at very frequently occurring sites (Additional file 5B and 5C). This, together with the range of restriction systems present, may be a factor in limiting the incursion of foreign DNA into these bacteria thus explaining the very low electroporation frequencies reported for bifidobacteria to date.
Unique genome regions in the B. longum strains
One unique region in each genome corresponds to a prophage. The prophage in strain NCC2705, which is truncated, appears to be a longtime resident of the genome as it does not correspond with a Base Deviation Index (BDI) peak (Fig. 2A), as this analysis predicts recent horizontal gene transfer (HGT) events. This appears to have been replaced in the genome of strain DJO10A with a different prophage, that is complete and inducible  and corresponds with a significant BDI peak substantiating this recent HGT event. The other five unique regions in strain NCC2705 contain largely hypothetical genes or genes of diverse functions without any significant gene clusters. However one of these regions (unique region 4') does encode putative xylan degradation genes, which is a function predicted to be important for competition in the large intestine. As this region corresponds to a BDI peak, it suggests it may be a recent acquisition by this strain and its evolution in the large intestine would provide the selective pressure for acquiring this unique region. Of the other 16 unique regions in the strain DJO10A, eight encode significant gene clusters involved in functions predicted to be important for competition in the large intestine, specifically oligosaccharide and polyol utilization, arsenic resistance and lantibiotic production.
Oligosaccharide and polyol utilization
Genome attenuation of B. longum in fermentation environments
Given the large number of unique DNA regions in the genome of strain DJO10A, that are predicted to have been lost from strain NCC2705, it suggests that deletion of DNA regions that are not useful may reflect the rapid adaptation of B. longum to a new and very different environment than exists in the human large gut. This would suggest an elevated mutation frequency. A comparative nucleotide substitution analysis between strains DJO10A and NCC2705 shows the majority of genes are highly conserved (Additional file 8), which is to be expected with two closely related strains. However, analysis of the 52 least conserved genes (listed as 'positive selection' in Additional file 8) indicates that of the mutations that can be attributed to one strain or the other (frameshifts, deletions, insertions and stop mutations), 11 are from strain NCC2705 and three from strain DJO10A (Additional file 9). Further substantiation of an increased mutation frequency in strain NCC2705 comes from comparing genes encoding surface protein homologs between the two strains. A search of the DJO10A genome for LPXTG motifs, which is a signature of one class of cell surface anchoring proteins found four potential proteins and SignalP analysis of these proteins (BLD1468, BLD1511, BLD1637 and BLD1638) confirmed the presence of a signal sequence in each case (Additional file 10). In addition, BLASTP analysis of these four proteins showed that they are very similar to other known surface proteins containing the LPXTG motif. The NCC2705 showed three of these gene homologs (BLD1468, BLD1637 and BLD1638), and had a predicted protein exhibiting 99% amino acid identity to BLD1511, but was missing the LPXTG motif due to an ISL3 insertion in the 3' end of the gene. This further highlights the rapid evolutionary status of bifidobacteria when they are removed from the human gut into pure-culture fermentation environments.
IS30 'jumping' in the B. longum genome
Adaptation of B. longum DJO10A to a pure-culture environment
To test the hypothesis that the switch from a variable and complex environment like the gut to a relatively stable and simplified, pure-culture one, results in hyper IS30 activity and rapid DNA loss of regions that are no longer beneficial to the new environment, strain DJO10A was grown in a typical laboratory medium without pH control for ~1,000 generations. Isolated colonies were then screened for seven unique regions encoding functions predicted to be useful for survival in the human gut. One of these regions (no. 12) involved in the lantibiotic production was found to be missing from 40% of the isolates (Additional file 11) substantiating this hypothesis. Analysis of this adapted strain, DJO10A-JH1, shows the deletion extended over the full lantibiotic region very similar to strain NCC2705 (Fig. 6A). It is further noted using Pulsed Field Gel Electrophoresis (PFGE) that the 39.9 kb Xba I band containing this region is missing from strain DJO10A-JH1 (Fig. 6B). The loss of the complete lantibiotic gene cluster from 40% of the culture was intriguing as the cluster also encodes the immunity gene to protect cells from the lantibiotic activity. However, analysis of lantibiotic production by strain DJO10A showed that none occurred during growth in broth media, and a solid surface such as agar, was needed for production (Fig. 6C) similar to streptin production from Streptococcus pyogenes . The loss of the complete lantibiotic gene cluster renders strain DJO10A-JH1 sensitive to this pronase-E sensitive lantibiotic, which is also active against a wide spectrum of bacteria (Fig. 6C). Interestingly, the lantibiotic genome region that was deleted during the adaptation of strain DJO10A to the pure-culture environment was located between two IS30 elements, suggesting its role in genome deletion events.
It was also noted that the pure-culture adapted strain, DJO10A-JH1, was also missing a 140.7 kb Xba I band (Fig. 6B). It is intriguing that this band contains one of the four MIC elements, suggesting it may have been involved. PCR analysis of the loci immediately bordering this MIC element revealed the deletion extended between 10 and 50 kb directly downstream from this element substantiating its likely role in this deletion event. This further substantiated the rapid loss of DNA, and the prominent role of mobile elements, during evolutionary adaptation by these bacteria.
Southern hybridization of strains DJO10A and DJO10A-JH1 substantiate the IS30 'jumping' during growth in a pure-culture environment, while the positions of the other IS elements (IS21, IS256 and ISL3) remained stable (Fig. 7B). This IS30 hyperactivity in B. longum genomes may play an important role in deletion events and genome reduction during adaptation to new environments.
Competitive 'fitness' of the adapted B. longum strainDJO1A-JH1
The rapid genome reduction experienced by B. longum DJO10A during pure-culture growth in fermentation conditions suggested that the genomic regions lost may have been important for competition in the intestine. To test if this in vitro adaptation affected the 'fitness' of the strain, a simulated fecal competitive approach was developed. Bifidobacteria are frequently proposed to successfully compete against members of the clostridia and the enterobacteriae in the intestinal environment. A member of both of these bacterial groups was selected to test the relative competitive abilities of B. longum DJO10A and its in vitro adapted derivative, strain DJO10A-JH1. To ensure that the selected competitor strains were not attenuated in any way, new isolates were obtained from fresh feces by plating on selective media and speciated using a sequence analysis of the 16S rRNA gene. This resulted in the isolation of Clostridium difficile DJOcd1 and E. coli DJOec1, which were minimally cultured prior to undergoing fecal competitive experiments with the B. longum strains. An in vitro growth rate analysis established that E. coli DJOec1 had the fastest growth rate, followed by C. difficile DJOcd1, B. longum DJO10A-JH1 and B. longum DJO10A (Additional file 12). The noticeable increased growth rate of B. longum DJO10A-JH1 compared to strain DJO10A substantiated that the genome reduction of strain DJO10A-JH1 favored the in vitro growth environment.
While the simulated fecal competition studies suggested that the lantibiotic encoding genome region was important for competition in the human intestine, in vivo studies in an intestinal model would be necessary to verify this hypothesis.
This study compares the genomic sequences of two strains of B. longum and suggests that bifidobacteria can rapidly loose genome regions during pure-culture growth that may be important for intestinal survival. This genomic prediction was experimentally validated during pure-culture growth of strain DJO10A and the genome reduction was shown to reduce competitive 'fitness' in a simulated fecal environment. The rapid loss of genomic regions that may be important for intestinal competition may compromise the ability of exogenous bifidobacteria to re-colonize human intestines.
Bacterial strains and growth conditions
Bifidobacterium longum strain DJO10A was isolated from a healthy young adult's feces  and B. animalis subsp. lactis BB12 was obtained from Chr. Hansen. B. animalis subsp. lactis strains S1, S2, and S14 are genetically distinct isolates from fermented foods (J.-H. Lee. and D.J.O'Sullivan, unpublished). Clostridium difficile DJOcd1 was isolated from fresh feces by plating on Clostridium difficile Selective Agar (BD Diagnostics) and speciated using a sequence analysis of its 16S rRNA gene. E. coli DJOec1 was obtained from fresh feces by plating on MacConkey agar (Difco) and speciated using a sequence analysis of its 16S rRNA gene. E. coli K12 was obtained from the American Type Culture Collection (ATCC). Bifidobacteria were cultivated at 37°C in MRS (Difco) supplemented with 0.05% L-cysteine·HCl (Sigma), Bifidobacteria Low-Iron Medium (BLIM)  or Bifidobacteria Fermentation Medium (BFM) (2% proteose peptone, 0.15% K2HPO4, 0.15% MgSO4·7H2O, 0.5% D-glucose) under anaerobic conditions using either the BBL Anaerobic system (BBL) or the Bactron II Anaerobic/Environmental Chamber (Sheldon Manufacturing).
Genome sequencing and assembly
Whole-genome shotgun sequencing was carried out at the US Department of Energy Joint Genome Institute (JGI). Sequences were assembled into 227 contigs using the Phred/Phrep/Consed software and the sequence coverage was 9.2-fold. Gap closure and genome sequence finishing was carried out at Fidelity Systems using ThermoFidelase-Fimer direct genome sequencing technology . Shotgun reads with and without IS30 elements covering A5, A6 and A7 loci were identified and assembled separately. The presence and location of long repeated sequences in genomic DNA samples were verified by direct genomic sequencing of the unique/repeat junctions. The resolution of the most complex high GC-rich repeats was achieved by sequencing of PCR products amplified with a hybrid TopoTaq DNA polymerase with increased strand displacement capacity.
Annotation of all open reading frames (ORFs) was carried out using Glimmer, GeneMark, JGI annotation tools and GAMOLA , before manual checking of all predicted genes. A comparative analysis of the two B. longum genomes was conducted using MUMmer3, ACT4 and ClustalX. The origin of replication and terminus were predicted using OriLoc . Codon usage was analyzed using the General Codon Usage Analysis (GCUA) program . The base-deviation index (BDI) was performed by scaled χ2 analysis of Artemis8. To predict gene functions, the two conserved protein domain databases of GAMOLA and InterProScan were used. COG functional categories were used for functional classification of all genes in both B. longum genome sequences.
General sequencing was conducted using a Big-Dye terminator and ABI Prism 3730xl Auto sequencer (Applied Biosystems). All PCR primers are listed in Additional file 13. For Southern blot analysis of unique region 12, a 646 bp probe from the lanM gene was obtained using PCR with LANT-F and LANT-R primers. Probes for IS elements were also PCR amplified. Probes were DIG-labeled and hybridized with digested genomic DNA according to the manufacturer's instructions (Roche). Pulsed field gel electrophoresis of Xba I-digested B. longum genomes was performed using a CHEF-DR III Variable Angle Pulsed Field Electrophoresis System according to manufacturer's instructions (Bio-Rad).
Identification of gene homologs between the two B. longum genomes
Comparative nucleotide substitution analysis by Nei and Gojobori's algorithm  was used to identify gene homologs. The predicted genes of both genome sequences were compared using the local BlastN program in the NCBI toolkit and 1,590 aligned genes were used for the nucleotide substitution analysis by Nei's unweighted method I . According to the ratio of dN:dS, all matched genes were categorized into three groups, highly conserved (< 0.035), normal, and positive selection (> 1).
Minimal inhibitory concentration of arsenic
To determine the minimal inhibitory concentration of arsenic, BLIM was supplemented with different concentrations of sodium arsenite (AsO2-, 1 to 100 mM) and sodium arsenate (AsO3-, 1 to 500 mM). Freshly grown cultures were sub-inoculated into the arsenite/arsenate media and incubated anaerobically at 37°C for 48 h.
Adaptation of B. longum DJO10A to in vitro fermentation conditions
B. longum DJO10A was grown in BFM continuously up to ~1,000 generations. The culture was then serially diluted and plated on BFM agar. Ten colonies were randomly selected for analysis.
Mapping the deletions in strain DJO10A-JH1
To find the precise location of the deletion of the lantibiotic operon in the B. longum DJO10A-JH1 genome, PCR was used to test for several genes within the lantibiotic operon. The two primers F3 (position 1,974,570–1,974,587 bp) and R3 (position 1,996,024–1,996,005 bp) were used to amplify a ~1.8 kb region spanning the deletion and sequencing located the precise borders (Figure 6). To map the position of the deletion in the 140.7 kb Xba I fragment, primers MIC-F1 (position, 1,539,767–1,539,768) and MIC-R1 (position, 1,542,535–1,542,553) were used to amplify the upstream region of MIC III and primers MIC-F2 (position, 1,543,406–1,543,424) and MIC-R2 (position, 1,545,713–1,545,732) were used to amplify the downstream region.
Bioassay for lantibiotic activity
B. longum DJO10A was inoculated into the center of an MRS agar plate and incubated anaerobically at 37°C for 2 days. After incubation, molten 0.5% top agar of the same medium containing 1% of an indicator strain was overlaid on the plates prior to incubation.
Simulated fecal competitive analysis of bifidobacteria
To access the competitive 'fitness' of the wild-type B. longum DJO10A compared to its in vitro adapted derivative strain DJO10A-JH1, a simulated in vitro fecal system was developed. Triplicate experiments for each strain were used. Each experiment was conducted in 10 g sterilized feces in an anaerobic chamber, to which 0.38 g Reinforced Clostridial Medium (RCM) and 0.02 g mucin (Porcine gastric type III) was added. The two competitor bacteria were added to all tubes at calculated concentrations of 1.2 × 107 cfu/g for E. coli DJOec1 and 5.1 × 107 for Clostridium difficile DJOcd1. B. longum DJO10A was added to three tubes at a calculated concentration of 4.0 × 107 cfu/g and strain DJO10A-JH1 to the other three tubes at 4.4 cfu/g. Standard viable plate counts were used to calculate all bacterial concentrations. After thorough mixing in an anaerobic environment, the tubes were left at 37°C for 3 days, whereby the entire fecal samples were homogenized in 90 ml peptone water to conduct an accurate serial plate count analysis.
Nucleotide sequence accession number
The Sequence and annotation data have been deposited in GenBank under the accession number CP000605.
This study was supported by the Minnesota Agricultural Experiment Station, Dairy Management Inc., the Midwest Dairy Research Center, and the U.S. Department of Energy (including SBIR grants DE-FG02-98ER82577 and DE-FG02-00ER83009 to SK). The other members of the Lactic Acid Bacteria Genome Consortium, especially T. Klaenhammer, L. McKay, J. Broadbent, J. Steele, B. Hutkins and F. Breidt are acknowledged for their help in initiating the genome sequencing effort. We thank W. Xu and the University of Minnesota Supercomputing Institute for computational analysis support.
- Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312: 1355-1359. 10.1126/science.1124234.PubMedPubMed CentralView Article
- Yoshioka H, Iseki K, Fujita K: Development and differences of intestinal flora in the neonatal period in breast-fed and bottle-fed infants. Pediatrics. 1983, 72: 317-321.PubMed
- Tissier H: Traitement des infections intestinales par la methode de la flore bacterienne de l'intestin. Crit Rev Soc Biol. 1906, 60: 359-361.
- Mitsuoka T, Hayakawa K: The fecal flora in man. I. Composition of the fecal flora of various age groups. Zentralbl Bakteriol [Orig A]. 1973, 223: 333-342.
- Hopkins MJ, Sharp R, Macfarlane GT: Age and disease related changes in intestinal bacterial populations assessed by cell culture, 16S rRNA abundance, and community cellular fatty acid profiles. Gut. 2001, 48: 198-205. 10.1136/gut.48.2.198.PubMedPubMed CentralView Article
- Ishibashi N, Yaeshima T, Hayasawa H: Bifidobacteria: their significance in human intestinal health. Mal J Nutr. 1997, 3: 149-159.
- O'Sullivan DJ: Primary Sources of Probiotic Cultures. Probiotics in food safety and human health. Edited by: Goktepe I, Juneja VK, Ahmedna M. 2006, Boca Raton: Taylor & Francis/CRC Press, 91-107.
- Fukushima Y, Kawata Y, Hara H, Terada A, Mitsuoka T: Effect of a probiotic formula on intestinal immunoglobulin A production in healthy children. Int J Food Microbiol. 1998, 42: 39-44. 10.1016/S0168-1605(98)00056-7.PubMedView Article
- Su P, Henriksson A, Tandianus JE, Park JH, Foong F, Dunn NW: Detection and quantification of Bifidobacterium lactis LAFTI B94 in human faecal samples from a consumption trial. FEMS Microbiol Lett. 2005, 244: 99-103. 10.1016/j.femsle.2005.01.022.PubMedView Article
- Schell MA, Karmirantzou M, Snel B, Vilanova D, Berger B, Pessi G, Zwahlen MC, Desiere F, Bork P, Delley M, Pridmore RD, Arigoni F: The genome sequence of Bifidobacterium longum reflects its adaptation to the human gastrointestinal tract. Proc Natl Acad Sci USA. 2002, 99: 14422-14427. 10.1073/pnas.212527599.PubMedPubMed CentralView Article
- O'Sullivan DJ: Bifidobacteria and Siderophores Produced Thereby and Methods of Use. US Patent No. 6,746,672. 2004
- Lee JH, O'Sullivan DJ: Sequence Analysis of Two Cryptic Plasmids from Bifidobacterium longum DJO10A and Construction of a Shuttle Cloning Vector. Appl Environ Microbiol. 2006, 72: 527-535. 10.1128/AEM.72.1.527-535.2006.PubMedPubMed CentralView Article
- Ventura M, Lee JH, Canchaya C, Zink R, Leahy S, Moreno-Munoz JA, O'Connell-Motherway M, Higgins D, Fitzgerald GF, O'Sullivan DJ, van Sinderen D: Prophage-like elements in bifidobacteria: insights from genomics, transcription, integration, distribution, and phylogenetic analysis. Appl Environ Microbiol. 2005, 71: 8692-8705. 10.1128/AEM.71.12.8692-8705.2005.PubMedPubMed CentralView Article
- Skouloubris S, Ribas de Pouplana L, De Reuse H, Hendrickson TL: A noncognate aminoacyl-tRNA synthetase that may resolve a missing link in protein evolution. Proc Natl Acad Sci USA. 2003, 100: 11297-11302. 10.1073/pnas.1932482100.PubMedPubMed CentralView Article
- Min B, Pelaschier JT, Graham DE, Tumbula-Hansen D, Soll D: Transfer RNA-dependent amino acid biosynthesis: an essential route to asparagine formation. Proc Natl Acad Sci USA. 2002, 99: 2678-2683. 10.1073/pnas.012027399.PubMedPubMed CentralView Article
- Mackiewicz P, Zakrzewska-Czerwinska J, Zawilak A, Dudek MR, Cebrat S: Where does bacterial replication start? Rules for predicting the oriC region. Nucleic Acids Res. 2004, 32: 3781-3791. 10.1093/nar/gkh699.PubMedPubMed CentralView Article
- Zawilak A, Cebrat S, Mackiewicz P, Krol-Hulewicz A, Jakimowicz D, Messer W, Gosciniak G, Zakrzewska-Czerwinska J: Identification of a putative chromosomal replication origin from Helicobacter pylori and its interaction with the initiator protein DnaA. Nucleic Acids Res. 2001, 29: 2251-2259. 10.1093/nar/29.11.2251.PubMedPubMed CentralView Article
- Price C, Bickle TA: A possible role for DNA restriction in bacterial evolution. Microbiol Sci. 1986, 3: 296-299.PubMed
- Nilsson AI, Koskiniemi S, Eriksson S, Kugelberg E, Hinton JC, Andersson DI: Bacterial genome size reduction by experimental evolution. Proc Natl Acad Sci USA. 2005, 102: 12112-12116. 10.1073/pnas.0503654102.PubMedPubMed CentralView Article
- Fleischmann RD, Alland D, Eisen JA, Carpenter L, White O, Peterson J, DeBoy R, Dodson R, Gwinn M, Haft D, Hickey E, Kolonay JF, Nelson WC, Umayam LA, Ermolaeva M, Salzberg SL, Delcher A, Utterback T, Weidman J, Khouri H, Gill J, Mikula A, Bishai W, Jacobs WR, Venter JC, Fraser CM: Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains. J Bacteriol. 2002, 184: 5479-5490. 10.1128/JB.184.19.5479-5490.2002.PubMedPubMed CentralView Article
- Berger B, Pridmore RD, Barretto C, Delmas-Julien F, Schreiber K, Arigoni F, Brussow H: Similarity and differences in the Lactobacillus acidophilus group identified by polyphasic analysis and comparative genomics. J Bacteriol. 2007, 189: 1311-1321. 10.1128/JB.01393-06.PubMedPubMed CentralView Article
- Molenaar D, Bringel F, Schuren FH, de Vos WM, Siezen RJ, Kleerebezem M: Exploring Lactobacillus plantarum genome diversity by using microarrays. J Bacteriol. 2005, 187: 6119-6127. 10.1128/JB.187.17.6119-6127.2005.PubMedPubMed CentralView Article
- Tatusov RL, Galperin MY, Natale DA, Koonin EV: The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000, 28: 33-36. 10.1093/nar/28.1.33.PubMedPubMed CentralView Article
- Gostner A, Blaut M, Schaffer V, Kozianowski G, Theis S, Klingeberg M, Dombrowski Y, Martin D, Ehrhardt S, Taras D, Schwiertz A, Kleessen B, Luhrs H, Schauber J, Dorbath D, Menzel T, Scheppach W: Effect of isomalt consumption on faecal microflora and colonic metabolism in healthy volunteers. Br J Nutr. 2006, 95: 40-50. 10.1079/BJN20051589.PubMedView Article
- Ratnaike RN: Acute and chronic arsenic toxicity. Postgrad Med J. 2003, 79: 391-396. 10.1136/pmj.79.933.391.PubMedPubMed CentralView Article
- Olasz F, Kiss J, Konig P, Buzas Z, Stalder R, Arber W: Target specificity of insertion element IS30. Mol Microbiol. 1998, 28: 691-704. 10.1046/j.1365-2958.1998.00824.x.PubMedView Article
- Wescombe PA, Tagg JR: Purification and characterization of streptin, a type A1 lantibiotic produced by Streptococcus pyogenes. Appl Environ Microbiol. 2003, 69: 2737-2747. 10.1128/AEM.69.5.2737-2747.2003.PubMedPubMed CentralView Article
- Islam A: Iron reversible inhibition by bifidobacteria and microbial diversity of the hman intestine. MS thesis. 2006, University of Minnesota, Department of Food Science and Nutrition
- Slesarev AI, Mezhevaya KV, Makarova KS, Polushin NN, Shcherbinina OV, Shakhova VV, Belova GI, Aravind L, Natale DA, Rogozin IB, Tatusov RL, Wolf YI, Stetter KO, Malykh AG, Koonin EV, Kozyavkin SA: The complete genome of hyperthermophile Methanopyrus kandleri AV19 and monophyly of archaeal methanogens. Proc Natl Acad Sci USA. 2002, 99: 4644-4649. 10.1073/pnas.032671499.PubMedPubMed CentralView Article
- Altermann E, Klaenhammer TR: GAMOLA: a new local solution for sequence annotation and analyzing draft and finished prokaryotic genomes. OMICS. 2003, 7: 161-169. 10.1089/153623103322246557.PubMedView Article
- Frank AC, Lobry JR: Oriloc: prediction of replication boundaries in unannotated bacterial chromosomes. Bioinformatics. 2000, 16: 560-561. 10.1093/bioinformatics/16.6.560.PubMedView Article
- McInerney JO: GCUA: general codon usage analysis. Bioinformatics. 1998, 14: 372-373. 10.1093/bioinformatics/14.4.372.PubMedView Article
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMed
- Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, Azevedo V, Bertero MG, Bessieres P, Bolotin A, Borchert S, Borriss R, Boursier L, Brans A, Braun M, Brignell SC, Bron S, Brouillet S, Bruschi CV, Caldwell B, Capuano V, Carter NM, Choi SK, Codani JJ, Connerton IF, Danchin A: The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature. 1997, 390: 249-256. 10.1038/36786.PubMedView Article
- Xu J, Bjursell MK, Himrod J, Deng S, Carmichael LK, Chiang HC, Hooper LV, Gordon JI: A genomic view of the human-Bacteroides thetaiotaomicron symbiosis. Science. 2003, 299: 2074-2076. 10.1126/science.1080029.PubMedView Article
- Makarova K, Slesarev A, Wolf Y, Sorokin A, Mirkin B, Koonin E, Pavlov A, Pavlova N, Karamychev V, Polouchine N, Shakhova V, Grigoriev I, Lou Y, Rohksar D, Lucas S, Huang K, Goodstein DM, Hawkins T, Plengvidhya V, Welker D, Hughes J, Goh Y, Benson A, Baldwin K, Lee JH, Diaz-Muniz I, Dosti B, Smeianov V, Wechter W, Barabote R, Lorca G, Altermann E, Barrangou R, Ganesan B, Xie Y, Rawsthorne H, Tamir D, Parker C, Breidt F, Broadbent J, Hutkins R, O'Sullivan D, Steele J, Unlu G, Saier M, Klaenhammer T, Richardson P, Kozyavkin S, Weimer B, Mills D: Comparative genomics of the lactic acid bacteria. Proc Natl Acad Sci USA. 103: 15611-15616. 10.1073/pnas.0607117103.
- Kleerebezem M, Boekhorst J, van Kranenburg R, Molenaar D, Kuipers OP, Leer R, Tarchini R, Peters SA, Sandbrink HM, Fiers MW, Stiekema W, Lankhorst RM, Bron PA, Hoffer SM, Groot MN, Kerkhoven R, de Vries M, Ursing B, de Vos WM, Siezen RJ: Complete genome sequence of Lactobacillus plantarum WCFS1. Proc Natl Acad Sci USA. 2003, 100: 1990-1995. 10.1073/pnas.0337704100.PubMedPubMed CentralView Article
- Pridmore RD, Berger B, Desiere F, Vilanova D, Barretto C, Pittet AC, Zwahlen MC, Rouvet M, Altermann E, Barrangou R, Mollet B, Mercenier A, Klaenhammer T, Arigoni F, Schell MA: The genome sequence of the probiotic intestinal bacterium Lactobacillus johnsonii NCC 533. Proc Natl Acad Sci USA. 2004, 101: 2512-2517. 10.1073/pnas.0307327101.PubMedPubMed CentralView Article
- Sofia HJ, Burland V, Daniels DL, Plunkett G, Blattner FR: Analysis of the Escherichia coli genome. V. DNA sequence of the region from 76.0 to 81.5 minutes. Nucleic Acids Res. 1994, 22: 2576-2586. 10.1093/nar/22.13.2576.PubMedPubMed CentralView Article
- van Kranenburg R, Golic N, Bongers R, Leer RJ, de Vos WM, Siezen RJ, Kleerebezem M: Functional analysis of three plasmids from Lactobacillus plantarum. Appl Environ Microbiol. 2005, 71: 1223-1230. 10.1128/AEM.71.3.1223-1230.2005.PubMedPubMed CentralView Article