- Short genome report
- Open Access
Permanent draft genome sequence of the probiotic strain Propionibacterium freudenreichii CIRM-BIA 129 (ITG P20)
Standards in Genomic Sciencesvolume 11, Article number: 6 (2016)
Propionibacterium freudenreichii belongs to the class Actinobacteria (Gram positive with a high GC content). This “Generally Recognized As Safe” (GRAS) species is traditionally used as (i) a starter for Swiss-type cheeses where it is responsible for holes and aroma production, (ii) a vitamin B12 and propionic acid producer in white biotechnologies, and (iii) a probiotic for use in humans and animals because of its bifidogenic and anti-inflammatory properties. Until now, only strain CIRM-BIA1T had been sequenced, annotated and become publicly available. Strain CIRM-BIA129 (commercially available as ITG P20) has considerable anti-inflammatory potential. Its gene content was compared to that of CIRM-BIA1 T. This strain contains 2384 genes including 1 ribosomal operon, 45 tRNA and 30 pseudogenes.
Propionibacterium freudenreichii belongs to the class of Actinobacteria (Gram positive bacteria with a high GC content). This ‘Generally Recognized As Safe’ species is traditionally used as (i) a starter for Swiss-type cheeses where it is responsible for holes and aroma production (ii) a vitamin B12 and propionic acid producer in white biotechnologies, and (iii) a probiotic for use in humans and animals because of its bifidogenic properties (enhancing intestinal transit). A recent screening of 23 strains belonging to this species revealed the considerable anti-inflammatory properties of strain CIRM-BIA129 (ITG-P20) . Its anti-inflammatory potential is superior to that of the previously sequenced CIRM-BIA1 T strain. CIRM-BIA129 (ITG-P20) is currently used to make cheese at an industrial scale and indeed is consumed in large quantities because one gram of Swiss-type cheese contains at least 109 CFU. The genetic basis for the anti-inflammatory properties of P. freudenreichii are still poorly understood. Sequencing CIRM-BIA129 (ITG-P20) would enable investigation of the genomic determinants responsible for the important anti-inflammatory properties of this strain. To our knowledge, CIRM-BIA129 is the second strain in this species to be sequenced, annotated and made publicly available.
Classification and features
P. freudenreichii CIRM-BIA129 (ITG P20) cells are Gram-positive, microaerophilic, pleiomorphic (coccoid to rod shape forming ‘Chinese characters’) bacillae (1.0–1.5 μm × 0.5–0.8 μm wide) forming creamy-white colonies on YEL agar plates (Table 1). Cells grow at the bottom of liquid medium tubes (to escape oxygen) and tend to clot in liquid culture at the beginning of the stationary phase. Transmission electron microscopy pictures of liquid-grown cultures revealed a thick cell wall (Fig. 1) made of peptidoglycan. No exopolysaccharides were observed at the surface of the bacteria, unlike what is seen in numerous strains of P. freudenreichii , including the type strain CIRM-BIA1 T (Fig. 2).
A BLASTn of the 16S sequence of type strain CIRM-BIA1 T against the contigs of CIRM-BIA129 (ITG-P20) confirmed the affiliation of this strain to the species (100 % identity, 100 % coverage).
Representative genomic 16S rRNA sequences of the strains were compared with those of other type strains belonging to Actinobacteria present in the Ribosomal Database Project. (Figure 2) shows the phylogenetic tree.
CIRM-BIA129 (ITG P20) was found to be able to utilize glycerol, erythritol, L-arabinose, adonitol, galactose, D-glucose, D-fructose, D-mannose, inositol, arbutine, esculine, lactose, lactate and gluconate as carbon sources according to the gallery-API results.
CIRM-BIA129 (ITG P20) was shown to grow with up to 1 M NaCl in a chemically defined medium in presence of osmoprotectant (G. Jan, personal communication).
Genome sequencing information
Genome project history
Propionibacterium freudenreichii CIRM-BIA 129 (ITG P20) genome was sequenced to obtain information regarding mechanism(s) or molecule(s) responsible for anti-inflammatory properties of the strain. Project information and associated MIGS are shown in Table 2.
Growth conditions and genomic DNA preparation
The P. freudenreichii strain CIRM-BIA129 (ITG-P20), isolated from Emmental cheese by Actalia Dairy Products (Institut Technique du Gruyère, Actalia, Rennes, France), was provided by the CIRM-BIA Biological Resource Centre (Centre International de Ressources Microbiennes-Bactéries d’Intérêt Alimentaire, INRA, Rennes, France). It was cultivated at 30 °C without shaking in YEL . Growth was monitored spectrophotometrically at 650 nm, as well as by counting colony-forming units in YEL medium containing 1.5 % agar.
A cell pellet (equivalent to 2 × 1010 CFU) was obtained by centrifugation for 10 min at 5000 × g from a one-day exponentional phase culture of CIRM-BIA129 (ITG P20). DNA was extracted using the Blood & Cell Culture DNA Midi Kit (Qiagen) according to the manufacturer’s recommendations with the following modifications. Briefly, complete bacterial lysis was obtained by adding 220 mg of lysozyme powder (Qbiogene) to 3.5 ml of B1 buffer (Qiagen) followed by 2.5 h of incubation at 37 °C. High molecular weight genomic DNA was purified by gravity flow and anion exchange chromatography, eluted in 5 ml QF buffer (Qiagen) and precipitated with 3.5 ml of isopropanol. DNA was collected by centrifugation for 10 min at 4 °C and 15000 g and then air dried. DNA was resuspended in 100 μL TE 1X buffer at pH8 (Sigma) for two hours at 55 °C. DNA integrity was checked on an 0.8 % agarose gel. The OD260/280nm was 1.9.
Genome sequencing and assembly
The genome of the CIRM-BIA129 (ITG-P20) strain was obtained using a whole-genome strategy based on Illumina 36 nucleotide paired-end sequencing (Illumina genome analyser II) (Table 2). After filtration and trimming, the number of reads was downsized to 50 million of 34.8 bases mean read length, representing a 681-fold genome coverage. The reads were assembled using Velvet version 1.1.03  with a k-mer size of 31. The assembly resulted in 59 scaffolds of a length greater than 1 kb, with an N50 of 123 kb.
Automatic and manual annotations were carried out with the AGMIAL platform , using P. freudenreichii CIRM-BIA 1 T as the reference . Sub-cellular localization was predicted using SurfG+ , a software specific to sub-cellular prediction in Gram-positive bacteria. A comparison of the predicted proteome of CIRM-BIA129 with that of the sequenced type strain CIRM-BIA1 was performed with Koriblast software (Korilog, France) with its default parameters. Genes present in CIRM-BIA129 and not in the type strain CIRM-BIA1 were searched by tblastn of the CIRM-BIA129 predicted proteome against the complete genome of CIRM-BIA1. No hit proteins were specific to CIRM-BIA129. In order to search for genomic islands specific to CIRM-BIA129, a blastp of the CIRM-BIA129 proteome against the CIRM-BIA1 proteome was performed. Proteins filtered for a ‘hit coverage’ of less than 75 % and those corresponding to locus-tags colocalized along the chromosome were attributed to genomic islands. To search for probable pseudogenes in CIRM-BIA1, a blastp analysis of CIRM-BIA1 protein against CIRM-BIA129 proteins filtered on a ‘hit coverage’ of less than 75 % was performed. Among them, the corresponding genes were declared as pseudogenes in the event of two adjacent genes having the same ‘best hit’. S-layer homology domains inside the proteins were sought from the SLH.hmm pattern using hmmer with default parameters.
The genome was found to be composed of 111 contigs for a total size of 2,588,969 bp (67.3 % GC content). The contig length was encompassed between 158 bp and 179,656 bp with an average at 23,324 bp. The N50 is 46,655 bp. One ribosomal operon (containing one 16S rRNA gene and one 23S rRNA gene) was present, on contig 24, A total of 2354 complete genes were predicted, 2,307 of them protein-coding genes and 45 tRNA encoding genes. Thirty genes were pseudogenes. According to the COG results, 1,236 protein coding genes were assigned to a putative function, the remainder being annotated as hypothetical proteins. The properties and statistics of the genome are summarized in (Tables 3 and 4).
Insights into the genome sequence
Surface proteins encoding genes
In a previous research study , the authors observed that the removal of surface-associated proteins led to a marked decrease in anti-inflammatory properties. Numerous studies dealing with the anti-inflammatory properties of food or probiotic bacteria have also identified the S-layer proteins or other surface compounds responsible for this trait [8, 9].
The S-layer proteins were found to be paracrystalline mono-layered assemblies of proteins which coat the surface of bacteria. S-layer proteins were associated with the cell wall via an SLH domain, with a cell wall polymer serving as the anchoring structure. This SLH domain comprised about 40–50 amino-acids, and could be found as one or more copies in the protein. Proteins other than S-layer proteins are known to be anchored by SLH domains to the cell wall of bacteria and are called S-layer associated proteins . As P. freudenreichii CIRM-BIA129 is known for its anti-inflammatory properties [7, 11], it was necessary to determine the presence of protein sequences containing SLH domains. This search led to the identification of eight genes: slpB (PFCIRM129_00700) that contained five SLH domains; slpE, (PFCIRM129_05460) containing four SLH domains; slpF (PFCIRM129_01545), slpG (PFCIRM129_09890), slpA (PFCIRM129_09350) and inlA (PFCIRM129_12235) that contained three SLH domains at the C-terminal part, and slh2 (PFCIRM129_03800) and slpD (PFCIRM129_11775), containing two SLH domains. Interestingly, CIRM-BIA129 (ITG-P20) did not present a typical, thick S-layer at its surface, so these genes most probably code for SLAPs.
Comparison with the reference strain
A tblastn analysis revealed 77 new genes (no hit) that were not present in the type strain CIRM-BIA1. Most of them were of unknown function except for four, which encoded two transporters PFCIRM129_07475 (pseudogene) and PFCIRM129_08355, an amidohydrolase PFCIRM129_03910 and a transcription factor PFCIRM129_03595. The blastp results revealed three genomic islands that were present in CIRM-BIA129 but absent from CIRM-BIA1: (i) from PFCIRM129_09365 to 09660 encompassed by conjugal transfer protein gene PFCIRM_09215 and conjugative transfer gene complex protein PFCIRM129_09665, probably corresponding to an integrative conjugal element, (ii) from _10870 to _10910 comprising several relaxase genes and an exopolyphosphatase gene PFCIRM129_10900, and (iii) from PFCIRM129_10830 to 10905 containing three alpha, beta hydrolase genes (PFCIRM129_10860, 10890, 10895).
A blastp analysis revealed a gene encoding a ribulokinase PFCIRM129_07785 which was pseudogenized in CIRM-BIA129 but appeared to be functional in CIRM-BIA1.
By contrast, the glucokinase gene PFREUD_00950 involved in gluconate degradation was pseudogenized in CIRM-BIA1 but complete in CIRM-BIA129. This difference may explain the ability of CIRM-BIA129 to degrade gluconate. In the same way, inlA (containing three SLH domains, see above) was pseudogenized in CIRM-BIA1 but complete in CIRM-BIA129. This suggests a better ability to interact with eukaryotic cells, as internalin A was described as a bacterial adhesin.
The genome of CIRM-BIA129 revealed new genes that had never previously been described in the species. Some of them, which probably encode surface exposed proteins, may be of considerable importance to the adaptation of the bacterium to the intestinal tract. Its ability to degrade gluconate may enable it to survive in the intestine where this sugar is abundant. Genes including the SLH domain may be candidates for the immune properties of the strain because SLH domains enable the protein to be anchored in the cell wall.
yeast extract lactate
Foligne B, Breton J, Mater D, Jan G. Tracking the microbiome functionality: focus on Propionibacterium species. Gut. 2013;62:1227–8.
Malik AC, Reinbold GW, Vedamuthu ER. Evaluation of the taxonomy of the Propionibacterium. Can J Microbiol. 1968;11:1185–91.
Zerbino DR, Birney E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
Bryson K, Loux V, Bossy R, Nicolas P, Chaillou S, van de Guchte M, et al. AGMIAL: implementing an annotation strategy for prokaryote genomes as a distributed system. Nucleic Acids Res. 2006;34:3533–45.
Falentin H, Deutsch S-M, Jan G, Loux V, Thierry A, Parayre S, et al. The Complete Genome of Propionibacterium freudenreichii CIRM-BIA1T, a Hardy Actinobacterium with Food and Probiotic Applications. PLoS ONE. 2010;5:e11748.
Barinov A, Loux V, Hammani A, Nicolas P, Langella P, Ehrlich D, et al. Prediction of surface exposed proteins in Streptococcus pyogenes, with a potential application to other Gram-positive bacteria. Proteomics. 2009;9:61–73.
Foligne B, Deutsch SM, Breton J, Cousin FJ, Dewulf J, Samson M, et al. Promising immunomodulatory effects of selected strains of dairy propionibacteria evidenced in vitro and in vivo. Appl Environ Microbiol. 2010;76:8259–64.
Konstantinov SR, Smidt H, de Vos WM, Bruijns SC, Singh SK, Valence F, et al. S layer protein A of Lactobacillus acidophilus NCFM regulates immature dendritic cell and T cell functions. Proc Natl Acad Sci U S A. 2008;105:19474–8.
Taverniti V, Stuknyte M, Minuzzo M, Arioli S, De Noni I, Scabiosi C, et al. S-layer protein mediates the stimulatory effect of Lactobacillus helveticus MIMLh5 on innate immunity. Appl Environ Microbiol. 2013;79:1221–31.
Johnson B, Selle K, O’Flaherty S, Goh YJ, Klaenhammer T. Identification of extracellular surface-layer associated proteins in Lactobacillus acidophilus NCFM. Microbiology. 2013;159:2269–82.
Le Marechal C, Peton V, Ple C, Vroland C, Jardin J, Briard-Bion V, et al. Surface proteins of Propionibacterium freudenreichii are involved in its anti-inflammatory properties. J Proteomics. 2015;113:447–61.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci U S A. 2004;101:11030–5.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011;28:2731–9.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci. 1990;87:4576–9.
Stackebrandt E, Rainey FA, Ward-Rainey NL. Proposal for a New Hierarchic Classification System, Actinobacteria classis nov. Int J Syst Bacteriol. 1997;47:479–91.
Whitman WB, Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Ludwig W, et al. Bergey’s Manual of Systematic Bacteriology: Volume 5a: The Actinobacteria. 2nd ed. Springer; 2012.
Zhi X-Y, Li W-J, Stackebrandt E. An update of the structure and 16S rRNA gene sequence-based definition of higher ranks of the class Actinobacteria, with the proposal of two new suborders and four new families and emended descriptions of the existing higher taxa. Int J Syst Evol Microbiol. 2009;59:589–608.
Charfreitag O, Stackebrandt E. Inter- and Intrageneric Relationships of the Genus Propionibacterium as Determined by 16S rRNA Sequences. Microbiology. 1989;135:2065–70.
Judicial Commission of the International Committee on Systematics of Prokaryotes. Necessary corrections to the Approved Lists of Bacterial Names according to Rule 40d (formerly Rule 46). Opinion 86. Int J Syst Evol Microbiol. 2008;58:1975.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
We would like to thank Actalia for supplying the CIRM-BIA129 (ITG P20) strain. Work on sequencing was funded by INRA through the “AIP-Bioressource” call for projects.
Authors do not disclose any competing interests.
HF extracted the DNA, annotated the genome, supervised the project and wrote the paper. SMD annotated surface proteins and provided pictures. VL cured the database and deposited the genome in the EMBL database. AH, JB and JMA assembled the genome. SP phenotyped the strain. VC stored the strain and checked its identity. VB sequenced the strain. GJ and YLL participated to the writing and provided funding for the research. All authors read and approved the final manuscript.