Crucial Role for Insertion Sequence Elements in Lactobacillus helveticus Evolution as Revealed by Interstrain Genomic Comparison

Lactobacillus helveticus is a versatile dairy bacterium found to possess heterogeneous genotypes depending on the ecosystem from which it was isolated. The recently published genome sequence showed the remarkable flexibility of its structure, demonstrated by a substantial level of insertion sequence (IS) element expansion in association with massive gene decay. To assess this diversity and examine the level of genome plasticity within the L. helveticus species, an array-based comparative genome hybridization (aCGH) experiment was designed in which 10 strains were analyzed. The aCGH experiment revealed 16 clusters of open reading frames (ORFs) flanked by IS elements. Four of these ORFs are associated with restriction/modification which may have played a role in accelerated evolution of strains in a commercially intensive ecosystem undoubtedly challenged through successive phage attack. Furthermore, analysis of the IS-flanked clusters demonstrated that the most frequently encountered ISs were also those most abundant in the genome (IS1201, ISL2, ISLhe1, ISLhe2, ISLhe65, and ISLhe63). These findings contribute to the overall viewpoint of the versatile character of IS elements and the role they may play in bacterial genome plasticity.Lactobacillus helveticus is a gram-positive, homofermentative lactic acid bacterium which is widely used in the manufacture of cheeses, such as Swiss cheese and some Cheddar-type cheeses (22, 25). It is also commonly used in the production of different types of Italian cheeses, such as Parmigiano Reggiano (18) and Grana Padano, where it contributes to the formation of specific flavor compounds (42).Phylogenetic analysis of ribosomal protein sequences derived from lactobacilli and streptococci classified L. helveticus in the same group along with both gastrointestinal (GI) tract and dairy-specific species (14). Comparative analysis of the 16S rRNA of L. helveticus DPC4571 revealed 98.4% identity with Lactobacillus acidophilus NCFM and indicated that this probiotic strain was closely related to strain DPC4571, despite the different environments these two lactobacilli inhabit (4). The results of genomic analysis of L. helveticus suggested that two major events have occurred in the diversification process of L. helveticus from a common ancestor with L. acidophilus, selective gene loss and acquisition of a large number of insertion sequence (IS) elements (4). IS elements are DNA sequences capable of independent transposition within and between bacterial genomes (31). Their capacity for independent mobility demonstrates the parasitic nature of these elements (11); however, they can also be regarded as having a positive influence, as they assist in promoting genetic variation (1). Thus, even though the primal character of these elements remains unclear in that they may be considered simply as selfish DNA elements, their impact on the architecture of microbial genomes is undeniable. It has already been demonstrated that IS-related mutations occur in Escherichia coli (44), Lactococcus lactis (10), Mycobacterium tuberculosis (33), and Francisella tularensis (39). Their active role was also demonstrated in the evolution of Paracoccus methylutens DM12 plasmids (3). Early bioinformatic analysis of the L. helveticus DPC4571 genome sequence resulted in identification of IS-associated truncations in genes associated with cellobiose transport, acetaldehyde dehydrogenase and diacetyl reductase (6). Considering the extraordinary abundance of IS elements in the L. helveticus DPC4571 chromosome (213 in total), it is noteworthy that very few open reading frames (ORFs) are directly affected by their presence. Presumably, the vast majority of insertion events proved detrimental to some aspect of the strain''s competitiveness and so were not selected in the ensuing population. We believe that the phenomenonal abundance of IS elements in L. helveticus makes it a very suitable system in which to study the role of IS elements in the evolution of bacterial genomes, particularly in ecosystems which impose challenging selective pressures.The level of chromosomal synteny that exists between L. helveticus DPC4571 and L. acidophilus NCFM is surprising, especially since the latter strain contains only 17 IS elements, and this observation highlighted the need for further studies of mobile genetic elements in the L. helveticus species. In order to address this issue, we employed DNA microarray technology to compare the overall genetic complement and specific genes associated with IS elements in different strains of L. helveticus. The use of comparative whole-genome array-based comparative genome hybridization (aCGH) has already been successfully applied to the identification of genetic differences within many closely related microorganisms. For example, large genomic deletions were identified among pathogenic Mycobacterium avium subsp. avium and Mycobacterium avium subsp. paratuberculosis (49) and Tropheryma whipplei strains (27). In addition, the absence of five Streptomyces coelicolor genomic islands were reported in Streptomyces lividans (24), and differences in gene content were detected in other species, Salmonella enterica (38) and Xylella fastidiosa (26). In this work we compared the genomes of nine strains of L. helveticus which were isolated from the dairy environment.     

Comparative genomics of lactic acid bacteria reveals a niche-specific gene set

Background  

The recently sequenced genome of Lactobacillus helveticus DPC4571 [1] revealed a dairy organism with significant homology (75% of genes are homologous) to a probiotic bacteria Lb. acidophilus NCFM [2]. This led us to hypothesise that a group of genes could be determined which could define an organism's niche.     

Characterization of the Genome of the Dairy Lactobacillus helveticus Bacteriophage ΦAQ113

The complete genomic sequence of the dairy Lactobacillus helveticus bacteriophage ΦAQ113 was determined. Phage ΦAQ113 is a Myoviridae bacteriophage with an isometric capsid and a contractile tail. The final assembled consensus sequence revealed a linear, circularly permuted, double-stranded DNA genome with a size of 36,566 bp and a G+C content of 37%. Fifty-six open reading frames (ORFs) were predicted, and a putative function was assigned to approximately 90% of them. The ΦAQ113 genome shows functionally related genes clustered together in a genome structure composed of modules for DNA replication/regulation, DNA packaging, head and tail morphogenesis, cell lysis, and lysogeny. The identification of genes involved in the establishment of lysogeny indicates that it may have originated as a temperate phage, even if it was isolated from natural cheese whey starters as a virulent phage, because it is able to propagate in a sensitive host strain. Additionally, we discovered that the ΦAQ113 phage genome is closely related to Lactobacillus gasseri phage KC5a and Lactobacillus johnsonii phage Lj771 genomes. The phylogenetic similarities between L. helveticus phage ΦAQ113 and two phages that belong to gut species confirm a possible common ancestral origin and support the increasing consideration of L. helveticus as a health-promoting organism.     

Phenotypic and Genotypic Analysis of Amino Acid Auxotrophy in Lactobacillus helveticus CNRZ 32

The conversion of amino acids into volatile and nonvolatile compounds by lactic acid bacteria in cheese is thought to represent the rate-limiting step in the development of mature flavor and aroma. Because amino acid breakdown by microbes often entails the reversible action of enzymes involved in biosynthetic pathways, our group investigated the genetics of amino acid biosynthesis in Lactobacillus helveticus CNRZ 32, a commercial cheese flavor adjunct that reduces bitterness and intensifies flavor notes. Most lactic acid bacteria are auxotrophic for several amino acids, and L. helveticus CNRZ 32 requires 14 amino acids. The reconstruction of amino acid biosynthetic pathways from a draft-quality genome sequence for L. helveticus CNRZ 32 revealed that amino acid auxotrophy in this species was due primarily to gene absence rather than point mutations, insertions, or small deletions, with good agreement between gene content and phenotypic amino acid requirements. One exception involved the phenotypic requirement for Asp (or Asn), which genome predictions suggested could be alleviated by citrate catabolism. This prediction was confirmed by the growth of L. helveticus CNRZ 32 after the addition of citrate to a chemically defined medium that lacked Asp and Asn. Genome analysis also predicted that L. helveticus CNRZ 32 possessed ornithine decarboxylase activity and would therefore catalyze the conversion of ornithine to putrescine, a volatile biogenic amine. However, experiments to confirm ornithine decarboxylase activity in L. helveticus CNRZ 32 by the use of several methods were unsuccessful, which indicated that this bacterium likely does not contribute to putrescine production in cheese.     

IS elements in Aliivibrio salmonicida LFI1238: Occurrence,variability and impact on adaptability

Insertion sequence (IS) elements are short, self-replicating DNA sequences that are capable of efficiently spreading over the host genome. Possessing varied integration specificity IS elements are capable of the irreversible inactivation of genes, which diversifies the pool of intact genetic determinants in host populations. In the current study, we performed a complex analysis of IS elements (Vsa IS) in the previously sequenced genome of Aliivibrio salmonicida LFI1238 and proposed a model of the spread of the Vsa IS elements over the genome of this microorganism. Along with the prediction of the integration sites for Vsa IS elements, the current study provides an overview of the properties of A. salmonicida IS elements, as well as information regarding their occurrence in different bacterial classes. An analysis of individual alleles of the IS elements has allowed us to depict a history of the accumulation of mutations and to describe distinctive microevolution lines for actively transposing Vsa IS elements in the genome of A. salmonicida LFI1238. Our results demonstrate the high importance of the dead end microevolution of actively transposing Vsa IS elements for the inactivation of genes in A. salmonicida LFI1238.     

Comparative High-Density Microarray Analysis of Gene Expression during Growth of Lactobacillus helveticus in Milk versus Rich Culture Medium

Lactobacillus helveticus CNRZ32 is used by the dairy industry to modulate cheese flavor. The compilation of a draft genome sequence for this strain allowed us to identify and completely sequence 168 genes potentially important for the growth of this organism in milk or for cheese flavor development. The primary aim of this study was to investigate the expression of these genes during growth in milk and MRS medium by using microarrays. Oligonucleotide probes against each of the completely sequenced genes were compiled on maskless photolithography-based DNA microarrays. Additionally, the entire draft genome sequence was used to produce tiled microarrays in which noninterrupted sequence contigs were covered by consecutive 24-mer probes and associated mismatch probe sets. Total RNA isolated from cells grown in skim milk or in MRS to mid-log phase was used as a template to synthesize cDNA, followed by Cy3 labeling and hybridization. An analysis of data from annotated gene probes identified 42 genes that were upregulated during the growth of CNRZ32 in milk (P < 0.05), and 25 of these genes showed upregulation after applying Bonferroni's adjustment. The tiled microarrays identified numerous additional genes that were upregulated in milk versus MRS. Collectively, array data showed the growth of CNRZ32 in milk-induced genes encoding cell-envelope proteinases, oligopeptide transporters, and endopeptidases as well as enzymes for lactose and cysteine pathways, de novo synthesis, and/or salvage pathways for purines and pyrimidines and other functions. Genes for a hypothetical phosphoserine utilization pathway were also differentially expressed. Preliminary experiments indicate that cheese-derived, phosphoserine-containing peptides increase growth rates of CNRZ32 in a chemically defined medium. These results suggest that phosphoserine is used as an energy source during the growth of L. helveticus CNRZ32.     

Complete genome sequence of Lactobacillus helveticus H10

Lactobacillus helveticus strain H10 was isolated from traditional fermented milk in Tibet, China. We sequenced the whole genome of strain H10 and compared it to the published genome sequence of Lactobacillus helveticus DPC4571.     

Complete Genome Sequence of Lactobacillus helveticus R0052, a Commercial Probiotic Strain

Lactobacillus helveticus R0052 is a commercially available strain that is widely used in probiotic preparations. The genome sequence consisted of 2,129,425 bases. Comparative analysis showed that it was unique among L. helveticus strains in that it contained genes encoding mucus-binding proteins similar to those found in Lactobacillus acidophilus.     

Population dynamics of lactobacilli in Grana cheese

A survey on the presence, microbial diversity, and population dynamics of lactobacilli in Grana cheese is presented. Evolution of thermophilic rod lactic acid bacteria within the first two days from cheese making and during ripening was different according to different bacterial groups, which were selectively enumerated and identified by molecular methods. Species-specific microbial counts indicated prevalence ofLactobacillus helveticus in both the whey starter and the cheese at moulding, and ofLactobacillus delbrueckii subsp.lactis in cheese after two months of ripening. In more advanced ripening, a decrease of total thermophilic lactobacilli and an increase of mesophilic lactobacilli (mostly belonging toLactobacillus casei/paracasei andLactobacillus rhamnosus) was observed. PCR fingerprinting of lactobacilli, which was performed by PCR-fingerprinting, indicated a marked microbial heterogeneity within theLactobacillus spp. populations, which enabled strain (or group)-specific fingerprints to be observed.     

ISLpl1 Is a Functional IS30-Related Insertion Element in Lactobacillus plantarum That Is Also Found in Other Lactic Acid Bacteria

We describe the first functional insertion sequence (IS) element in Lactobacillus plantarum. ISLpl1, an IS30-related element, was found on the pLp3 plasmid in strain FB335. By selection of spontaneous mutants able to grow in the presence of uracil, it was demonstrated that the IS had transposed into the uracil phosphoribosyltransferase-encoding gene upp on the FB335 chromosome. The plasmid-carried IS element was also sequenced, and a second potential IS element was found: ISLpl2, an IS150-related element adjacent to ISLpl1. When Southern hybridization was used, the copy number and genome (plasmid versus chromosome) distribution data revealed different numbers and patterns of ISLpl1-related sequences in different L. plantarum strains as well as in Pediococcus strains. The ISLpl1 pattern changed over many generations of the strain L. plantarum NCIMB 1406. This finding strongly supports our hypothesis that ISLpl1 is a mobile element in L. plantarum. Database analysis revealed five quasi-identical ISLpl1 elements in Lactobacillus, Pediococcus, and Oenococcus strains. Three of these elements may be cryptic IS, since point mutations or 1-nucleotide deletions were found in their transposase-encoding genes. In some cases, ISLpl1 was linked to genes involved in cold shock adaptation, bacteriocin production, sugar utilization, or antibiotic resistance. ISLpl1 is transferred among lactic acid bacteria (LAB) and may play a role in LAB genome plasticity and adaptation to their environment.     

Peptidoglycan Hydrolases as Species-Specific Markers to Differentiate Lactobacillus helveticus from Lactobacillus gallinarum and Other Closely Related Homofermentative Lactobacilli

We propose a new method that allows accurate discrimination of Lactobacillus helveticus from other closely related homofermentative lactobacilli, especially Lactobacillus gallinarum. This method is based on the amplification by PCR of two peptidoglycan hydrolytic genes, Lhv_0190 and Lhv_0191. These genes are ubiquitous and show high homology at the intra-species level. The PCR method gave two specific PCR products, of 542 and 747 bp, for 25 L. helveticus strains coming from various sources. For L. gallinarum, two amplicons were obtained, the specific 542 bp amplicon and another one with a size greater than 1,500 bp. No specific PCR products were obtained for 12 other closely related species of lactobacilli, including the L. acidophilus complex, L. delbrueckii, and L. ultunensis. The developed PCR method provided rapid, precise, and easy identification of L. helveticus. Moreover, it enabled differentiation between the two closely phylogenetically related species L. helveticus and L. gallinarum.     

In vitro Effects of Prebiotics and Synbiotics on Apis cerana Gut Microbiota

This study aimed to investigate in vitro effects of the selected prebiotics alone, and in combination with two potential probiotic Lactobacillus strains on the microbial composition of Apis cerana gut microbiota and acid production. Four prebiotics, inulin, fructo-oligosaccharides, xylo-oligosaccharides, and isomalto-oligosaccharides were chosen, and glucose served as the carbon source. Supplementation of this four prebiotics increased numbers of Bifidobacterium and lactic acid bacteria while decreasing the pH value of in vitro fermentation broth inoculated with A. cerana gut microbiota compared to glucose. Then, two potential probiotics derived from A. cerana gut at different dosages, Lactobacillus helveticus KM7 and Limosilactobacillus reuteri LP4 were added with isomalto-oligosaccharides in fermentation broth inoculated with A. cerana gut microbiota, respectively. The most pronounced impact was observed with isomalto-oligosaccharides. Compared to isomalto-oligosaccharides alone, the combination of isomalto-oligosaccharides with both lactobacilli strains induced the growth of Bifidobacterium, LAB, and total bacteria and reduced the proliferation of Enterococcus and fungi. Consistent with these results, the altered metabolic activity was observed as lowered pH in in vitro culture of gut microbiota supplemented with isomalto-oligosaccharides and lactobacilli strains. The symbiotic impact varied with the types and concentration of Lactobacillus strains and fermentation time. The more effective ability was observed with IMO combined with L. helveticus KM7. These results suggested that isomalto-oligosaccharides could be a potential prebiotic and symbiotic with certain lactobacilli strains on A. cerana gut microbiota.     

Lactobacillus paracasei Comparative Genomics: Towards Species Pan-Genome Definition and Exploitation of Diversity

Lactobacillus paracasei is a member of the normal human and animal gut microbiota and is used extensively in the food industry in starter cultures for dairy products or as probiotics. With the development of low-cost, high-throughput sequencing techniques it has become feasible to sequence many different strains of one species and to determine its “pan-genome”. We have sequenced the genomes of 34 different L. paracasei strains, and performed a comparative genomics analysis. We analysed genome synteny and content, focussing on the pan-genome, core genome and variable genome. Each genome was shown to contain around 2800–3100 protein-coding genes, and comparative analysis identified over 4200 ortholog groups that comprise the pan-genome of this species, of which about 1800 ortholog groups make up the conserved core. Several factors previously associated with host-microbe interactions such as pili, cell-envelope proteinase, hydrolases p40 and p75 or the capacity to produce short branched-chain fatty acids (bkd operon) are part of the L. paracasei core genome present in all analysed strains. The variome consists mainly of hypothetical proteins, phages, plasmids, transposon/conjugative elements, and known functions such as sugar metabolism, cell-surface proteins, transporters, CRISPR-associated proteins, and EPS biosynthesis proteins. An enormous variety and variability of sugar utilization gene cassettes were identified, with each strain harbouring between 25–53 cassettes, reflecting the high adaptability of L. paracasei to different niches. A phylogenomic tree was constructed based on total genome contents, and together with an analysis of horizontal gene transfer events we conclude that evolution of these L. paracasei strains is complex and not always related to niche adaptation. The results of this genome content comparison was used, together with high-throughput growth experiments on various carbohydrates, to perform gene-trait matching analysis, in order to link the distribution pattern of a specific phenotype to the presence/absence of specific sets of genes.     

Host specific diversity in Lactobacillus johnsonii as evidenced by a major chromosomal inversion and phage resistance mechanisms

Genetic diversity and genomic rearrangements are a driving force in bacterial evolution and niche adaptation. We sequenced and annotated the genome of Lactobacillus johnsonii DPC6026, a strain isolated from the porcine intestinal tract. Although the genome of DPC6026 is similar in size (1.97 mbp) and GC content (34.8%) to the sequenced human isolate L. johnsonii NCC 533, a large symmetrical inversion of approximately 750 kb differentiated the two strains. Comparative analysis among 12 other strains of L. johnsonii including 8 porcine, 3 human and 1 poultry isolate indicated that the genome architecture found in DPC6026 is more common within the species than that of NCC 533. Furthermore a number of unique features were annotated in DPC6026, some of which are likely to have been acquired by horizontal gene transfer (HGT) and contribute to protection against phage infection. A putative type III restriction-modification system was identified, as were novel Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) elements. Interestingly, these particular elements are not widely distributed among L. johnsonii strains. Taken together these data suggest intra-species genomic rearrangements and significant genetic diversity within the L. johnsonii species and indicate towards a host-specific divergence of L. johnsonii strains with respect to genome inversion and phage exposure.     

Heterogeneity of Putative Surface Layer Proteins in Lactobacillus helveticus

The S-layer-encoding genes of 21 Lactobacillus helveticus strains were characterized. Phylogenetic analysis based on the identified S-layer genes revealed two main clusters, one which includes a sequence similar to that of the slpH1 gene of L. helveticus CNRZ 892 and a second cluster which includes genes similar to that of prtY. These results were further confirmed by Southern blot hybridization. This study demonstrates S-layer gene variability in the species L. helveticus.     

Detection and Genomic Characterization of Motility in Lactobacillus curvatus: Confirmation of Motility in a Species outside the Lactobacillus salivarius Clade

Lactobacillus is the largest genus within the lactic acid bacteria (LAB), with almost 180 species currently identified. Motility has been reported for at least 13 Lactobacillus species, all belonging to the Lactobacillus salivarius clade. Motility in lactobacilli is poorly characterized. It probably confers competitive advantages, such as superior nutrient acquisition and niche colonization, but it could also play an important role in innate immune system activation through flagellin–Toll-like receptor 5 (TLR5) interaction. We now report strong evidence of motility in a species outside the L. salivarius clade, Lactobacillus curvatus (strain NRIC 0822). The motility of L. curvatus NRIC 0822 was revealed by phase-contrast microscopy and soft-agar motility assays. Strain NRIC 0822 was motile at temperatures between 15°C and 37°C, with a range of different carbohydrates, and under varying atmospheric conditions. We sequenced the L. curvatus NRIC 0822 genome, which revealed that the motility genes are organized in a single operon and that the products are very similar (>98.5% amino acid similarity over >11,000 amino acids) to those encoded by the motility operon of Lactobacillus acidipiscis KCTC 13900 (shown for the first time to be motile also). Moreover, the presence of a large number of mobile genetic elements within and flanking the motility operon of L. curvatus suggests recent horizontal transfer between members of two distinct Lactobacillus clades: L. acidipiscis in the L. salivarius clade and L. curvatus in the L. sakei clade. This study provides novel phenotypic, genetic, and phylogenetic insights into flagellum-mediated motility in lactobacilli.     

Insertion sequence-caused large-scale rearrangements in the genome of Escherichia coli

A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 × 10⁻⁴ insertions per genome per generation and the rate of IS homologous recombination is 4.5 × 10⁻⁵ recombinations per genome per generation. These events are mostly contributed by the IS elements IS1, IS2, IS5 and IS186. Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS186 target-sites, 5′-GGGG(N6/N7)CCCC-3′. We also detected 48 long deletions not involving IS elements.     

Extensive intra-phylotype diversity in lactobacilli and bifidobacteria from the honeybee gut

Background

In the honeybee Apis mellifera, the bacterial gut community is consistently colonized by eight distinct phylotypes of bacteria. Managed bee colonies are of considerable economic interest and it is therefore important to elucidate the diversity and role of this microbiota in the honeybee. In this study, we have sequenced the genomes of eleven strains of lactobacilli and bifidobacteria isolated from the honey crop of the honeybee A. mellifera.

Results

Single gene phylogenies confirmed that the isolated strains represent the diversity of lactobacilli and bifidobacteria in the gut, as previously identified by 16S rRNA gene sequencing. Core genome phylogenies of the lactobacilli and bifidobacteria further indicated extensive divergence between strains classified as the same phylotype. Phylotype-specific protein families included unique surface proteins. Within phylotypes, we found a remarkably high level of gene content diversity. Carbohydrate metabolism and transport functions contributed up to 45% of the accessory genes, with some genomes having a higher content of genes encoding phosphotransferase systems for the uptake of carbohydrates than any previously sequenced genome. These genes were often located in highly variable genomic segments that also contained genes for enzymes involved in the degradation and modification of sugar residues. Strain-specific gene clusters for the biosynthesis of exopolysaccharides were identified in two phylotypes. The dynamics of these segments contrasted with low recombination frequencies and conserved gene order structures for the core genes. Hits for CRISPR spacers were almost exclusively found within phylotypes, suggesting that the phylotypes are associated with distinct phage populations.

Conclusions

The honeybee gut microbiota has been described as consisting of a modest number of phylotypes; however, the genomes sequenced in the current study demonstrated a very high level of gene content diversity within all three described phylotypes of lactobacilli and bifidobacteria, particularly in terms of metabolic functions and surface structures, where many features were strain-specific. Together, these results indicate niche differentiation within phylotypes, suggesting that the honeybee gut microbiota is more complex than previously thought.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1476-6) contains supplementary material, which is available to authorized users.