BACKGROUND
lactobacillus rhamnosus is a facultatively heterofermentative lactic acid bacterium  frequently encountered in many dairy products, where it can be added as a probiotic microorganism or can be naturally present arising from raw milk.
lab may play different roles in cheese manufacture: some species participate in the fermentation process and contribute to acid production acting as starter lab , whereas others, called non-starter lab , are mainly implicated in the maturation process.
in particular, l. rhamnosus plays a significant role during ripening, leading to the formation of flavor  <cit>  and, for this reason, members of this species are generally recognized as nslab.
it is noteworthy that nslab generally have a high tolerance to hostile environments, such as those with high salt concentration, low moisture,  <dig> - <dig>  ph values, low temperatures and deficiency of nutrients  <cit> .
moreover, several studies have reported that in long-ripened cheese varieties, nslab populations dominate during aging after slab decline due to autolysis  <cit> .
increasing by about four to five orders of magnitude within a few months, nslab can have a major impact in determining curd maturation and final characteristics of cheese  <cit> .
in particular, l. rhamnosus has been shown to become dominant within nslab population in several cheeses, including parmigiano reggiano   <cit> .
it persists throughout the whole time of pr cheese ripening  and this implies its capacity to adapt to changing environmental conditions  <cit> .
notably, different l. rhamnosus strains have been detected in relation to specific steps of cheese ripening, suggesting that these strains may have specific metabolic activities, which could account for the adaptation to the changing microenvironment of cheese ripening  <cit> .
up to date, the molecular mechanisms allowing the growth of l. rhamnosus in cheese are still poorly understood.
nslab development during ripening can be attributed to their ability to use the major nutrient sources available in ripened cheese, that is lactose-free.
the nutrient sources available include milk components modified by technological treatment  and starter lab development, starter lab metabolites and cell lysis products.
potential substrate for microbial growth are represented by small peptides or amino acids  <cit> , citrate, lactate, and free fatty acids  <cit> .
additionally, sugars and phospholipids, nucleic acids and peptides can be released in the cheese matrix when slab autolysis starts to occur  <cit> .
these compounds represent carbon sources that could yield intracellular pyruvate  or be converted into different metabolites.
to investigate the metabolic pathways occurring in l. rhamnosus during cheese ripening, bove and colleagues  <cit>  recently compared the proteomic profiles of 10 l. rhamnosus strains grown in mrs and a cheese-like medium .
differently from mrs, which is the standard laboratory medium for lactobacilli  <cit>  and is considered to be a rich substrate with glucose as the primary carbon source for microbial growth, cb is an experimental medium formulated with 20-month-ripened pr cheese  <cit> , that tries to mimic the nutritional composition of pr cheese during ripening.
pr cheese, and thus cb, is considered to be a substrate poor in carbohydrates and characterized by the absence of milk sugar, lactose.
during the curd acidification step of pr cheese production, the conversion of lactose into lactic acid is the main biochemical process that occurs.
lactose is completely depleted within  <dig> to 48 h  <cit> .
the composition of cb, prepared as the protocol of neviani et al.
<cit> , is the following : proteins,  <dig> ; lactose,  <dig> ; glucose,  <dig> ; galactose,  <dig> ; lactic acid,  <dig> ; nacl,  <dig> ; sodium citrate,  <dig> .
according to the findings of bove et al.
<cit>  compared to the cultivation in mrs, the differentially expressed proteins under cheese-like conditions were mainly linked to protein biosynthesis and catabolism, nucleotide and carbohydrate metabolisms, citrate metabolism, cell wall and exopolysaccharide biosynthesis, cell regulation, oxidation/reduction processes, and stress response  <cit> .
notably, l. rhamnosus produced lactic acid as the primary end product when growing in mrs, whereas in cb low levels of lactic acid together with high levels of acetic acid were detected for all strains.
despite this common trend, the authors also observed strain-specific physiological responses, suggesting a strain variability in the adaptation to changing environmental conditions in accordance with genetic polymorphism studies  <cit> .
among all strains, l. rhamnosus pr <dig>  isolated from pr 4-month-ripened cheese, has shown the greatest ability to growth in cb reaching after 48 h the highest cell numbers with one of the lowest levels of lactic acid and one of the highest levels of acetic acid.
for this reason, in the present work, we focused our attention only on this strain, with the aim to identify the genes that could concur to explain its growth ability in cb and its acid acetic production.
the physiological adaptation of l. rhamnosus pr <dig> in cb was evaluated using a transcriptomic approach, based on cdna-amplified fragment length polymorphism  and quantitative real-time reverse transcription-pcr .
cdna-aflp is one of the most robust and sensitive transcriptomic technologies for genome-wide expression studies, with the main advantage of not requiring any prior knowledge of gene sequences while allowing the detection of lowly expressed genes through transcript amplification  <cit> .
using this approach, we identified a set of genes resulted over-expressed in cb compared to mrs, potentially involved in alternative metabolic pathways.
interesting genes were searched in other nslab and slab genomes with the aim to explore their diversity.
overall, the results described in this work highlight mechanisms of adaptation leading to the production of acetic acid coupled with atp generation, that could support the l. rhamnosus growth in cheese during ripening.
methods
bacterial growth conditions
l. rhamnosus pr <dig> was isolated from parmigiano reggiano  at 4 months of ripening on cheese based medium  <cit>  plate counts and identified by 16s rdna gene sequencing  <cit>  and species-specific pcr  <cit> .
the strain was cultivated in mrs broth  or cheese broth  at 30°c, under anaerobiosis, for  <dig> or 48 h, respectively.
cb, a culture medium that mimics raw-milk long-ripened cheese, was prepared according to the modified protocol described by bove et al.
<cit> .
rna extraction and cdna synthesis
the growth of l. rhamnosus pr <dig> in mrs and cb broth was monitored by measuring optical density  at 600 nm.
about  <dig> cells at the top of logarithmic phase were harvested, and total rna, stabilized with rnaprotect bacteria reagent , was isolated using rneasy protect bacteria mini kit .
three independent biological experiment were made.
rna was quantified using a nanodrop nd- <dig> spectrophotometer  and visualized by formaldehyde agarose gel electrophoresis according to standard procedures.
all rnas were of sufficient quantity  and high quality .
after a step of mrna enrichment and polyadenylation of rna transcripts, cdna was synthesized by reverse transcription  using a biotinylated oligo , following the protocol reported by bove et al.
<cit> .
for qpcr, single-stranded cdna was prepared in a final volume of 20 μl from total rna using an oligo  primer and superscript™ ii reverse transcriptase , according to the manufacturer’s instructions.
briefly, 12-μl reaction mixtures containing 500 ng of oligo  primer, 2 μg total rna and 10 nmol dntp mix in depc-treated h2o were heated to 65°c for 5 min, added with 4 μl of 5x first-strand buffer  and 200 nmol dtt, and then incubated at 42°c for 2 min.
rt reactions were started by the addition of  <dig> u of enzyme, incubated at 42°c for 50 min and inactivated by heating at 70°c for 15 min.
rt step was carried out in duplicate.
cdna-aflp
cdna-aflp analysis was carried out as described by bove et al.
<cit> .
the protocol is based on the production of cdna-aflp fragments that are detected using infrared dye  detection technology and the odyssey infrared imaging system.
briefly, after cdna synthesis, a double digestion was carried out with ecori and msei restriction enzymes and fragments were captured with the aid of streptavidin-coated magnetic beads.
digested cdna fragments were subsequently ligated with adaptors to allow selective amplification with ecori primers labeled with an infrared dye , and unlabeled msei-n .
three primer combinations were used to selectively amplify the expressed genes: dy-ecori-ac/msei-at, dy-ecori-at/msei-ac and dy-ecori-at/msei-at  <cit> .
ligators and primers used are reported in table  <dig>  separation of cdna-aflp fragments was carried out in a polyacrylamide gel and visualized by odissey  at 700 nm.
primer sets were designed as reported by bove et al.
<cit> .
cdna-aflp fragment isolation, re-amplification and sequencing
transcript-derived fragments  of interest were cut from polyacrylamide gels as reported by vuylsteke et al.
<cit> , resuspended in  <dig> μl of distilled water and subsequently re-amplified using the re-amplification and selective pcr primers ecori-ac/msei-at, ecori-at/msei-ac and ecori-at/msei-at  according to the origin of cdna-aflp fragments.
amplification reactions were performed in a final volume of  <dig> μl containing  <dig> μl of resuspended dna fragment, 25 mm mgcl <dig>  10x pcr buffer,  <dig> μm ecori-n primer,  <dig> μm msei-n primer, 5 mm dntps,  <dig>  μl of amplitaq  <dig> dna polymerase  and  <dig> μl of  <dig> gc enhancer .
pcr consisted of: i) 30 s of denaturation step at 94°c, 30 s of annealing step at 65°c  and 1 min of extension step at 72°c for  <dig> total cycles; ii) 10 s of denaturation step at 94°c, 30 s of annealing step at 56°c and 1 min of extension step at 72°c  for  <dig> total cycles.
these cycles were preceded by a common denaturation step of 2 min at 94°c and followed by a final 10-min extension at 72°c, and were carried out in a mastercycler ep gradient s thermal cycler .
amplified products were checked on a 1% agarose gel with a 100-bp marker  and subsequently purified using the wizard® sv gel and pcr clean-up system according to the manufacturer’s instructions .
amplified fragments were then cloned in e.coli using the pgem-t easy vector system kit , and plasmids from selected clones were purified using pureyield miniprep system kit  referring to the producer’s manual.
cloned fragments were finally sequenced by eurofins mwg operon using primers m <dig> and sequences were analysed by blast alignment  <cit> .
tdf sequences were deposited in the ddbj database under the accession numbers ab <dig> to ab <dig>
qpcr and data processing
qpcr was carried out using the lightcycler sybr green system  as previously described  <cit> .
briefly, 1 μl of cdna template was used in each reaction along with 4 μl of sybr green pcr master mix  and  <dig> pmol of the appropriate gene-specific primers in a final volume of 20 μl.
the following cycle profile was used: 10 min at 95°c,  <dig> repeats of 15 s at 95°c, 25 s at 58°c for spxb, ulae and 16s rdna genes or 55°c for xfp, 72°c for 20 s  and an additional 5-s incubation step at 81°c for fluorescence acquisition.
oligonucleotide sequence information and detailed primer-specific conditions are given in table  <dig>  two technical replicates were done for each combination of cdna and primer pair.
to assess background and residual dna contamination, a no-template control  and a no-reverse transcription control  were performed for each target.
dna contamination was considered to be negligible when the difference in cq  between the sample and the respective nort was above 5 cycles.
product detection and pcr specificity were checked post-amplification by examining the dissociation curves.
pcr amplicons were resolved by 2% agarose gel electrophoresis to verify the expected size.
to evaluate repeatability and reproducibility of the qpcr assay, intra- and inter-assay coefficients of variation  were assessed.
the intra-assay cv was from  <dig>  to  <dig> % whereas the inter-assay cv ranged from  <dig>  to  <dig> %.
amplification efficiency was calculated from the slope of standard curves generated with two-fold serial dilutions of the same cdna sample, as e =  <dig>  relative expression of target genes was determined using the ΔΔct method after pfaffl correction  <cit> .
16s rdna was used as a reference gene.
a
b
aprimer sets were designed based on the sequences of cdna-aflp fragments.
primers for 16s rdna gene were designed as reported by giraffa et al.
<cit> .
btarget gene expression was calculated relative to 16s rdna as a reference gene using the efficiency-corrected ΔΔct method  <cit> .
the relative expression ratios in cb compared to mrs are shown.
in silico analysis
tdf sequences were annotated using blast search.
pathway assignment was performed according to cog   <cit>  functional categories and kegg   <cit>  pathway database.
gene synteny across nslab and slab genomes was explored through the web server synttax  <cit> .
genome mining for promoter and terminator elements was performed using pepper toolbox  <cit> .
translated protein sequences were subjected to pfam motif analysis  <cit> .
protein alignments were performed using clustalw <dig>  <cit>  and used for phylogenetic tree construction at the interactive tree of life  <cit> .
multisequence amino acid alignments were represented using clc-bio sequence viewer  <cit> .
RESULTS
cdna-aflp analysis
in this study, the cdna-aflp technique  <cit>  was applied to profile the transcriptome of a l. rhamnosus strain grown in conditions mimicking cheese ripening.
despite it is not widely used in bacteria, cdna-aflp can be considered an ideal system for genome-wide expression analysis, mainly for the detection of lowly expressed genes.
three primer combinations were used to selectively amplify the genes expressed by l. rhamnosus pr <dig> in cb and mrs, allowing to generate different cdna-aflp profiles with a fragment size ranging from  <dig> to 500 bp .
a total of  <dig> and  <dig> tdfs were detected in mrs and cb, respectively.
in order to investigate the main adaptations of l. rhamnosus to the pr cheese environment, we focused on tdfs over-expressed in cb.
identification of tdfs over-expressed in cb
twenty tdfs strongly over-expressed by l. rhamnosus in cb compared to mrs were extracted from gel and used as templates for re-amplification by pcr.
direct sequencing of the pcr products failed probably because the sequences were contaminated with co-migrating bands or were too small .
products from bands were then cloned and pcr amplicons from  <dig> individual colonies were sequenced.
reliable  sequences were obtained for  <dig> of the  <dig> tdfs.
each sequence was identified by similarity search using the blast program against the genbank non-redundant  public sequence database.
as shown in table  <dig>   <dig> transcripts  showed a significant similarity to sequences with known function,  <dig> transcripts  were closely related to l. casei plasmid sequences,  <dig> transcripts  were annotated as hypothetical protein-coding sequences, and  <dig> transcripts  were identified as 5s rrna.
l. rhamnosus
pr <dig> over-expressed in cb compared to mrs
a
b
c
d
cog
e
awhen available, ec numbers assigned to the putative enzymatic reactions are provided.
bthe column indicates the microorganism of the best hit from blastx search.
cmax identity and e-value from the best hit of blastx search are provided.
dpathway assignment was performed according to cog functional categories and kegg pathway database.
ee, amino acid transport and metabolism; f, nucleotide transport and metabolism; g, carbohydrate transport and metabolism; m, cell wall/membrane/envelope biogenesis; r, general function prediction only.
it is known that plasmids often carry genes that might be essential for survival under harsh conditions, encoding important traits, such as enzymes involved in secondary metabolic pathways  <cit> .
plasmids are known to be a source of lab genetic and phenotypic diversity which occasionally confers adaptive advantages to host strains  <cit> .
however, further studies are clearly needed to better explore the role of plasmid sequences in the l. rhamnosus adaptation to the cheese ripening environment.
to validate the cdna-aflp expression profiles,  <dig> genes, encoding pyruvate oxidase , l-xylulose 5-phosphate 3-epimerase , and xylulose-5-phosphate phosphoketolase  were selected for qpcr.
the relative mrna abundances were normalized by that of the commonly used reference gene 16s rdna, and expressed as a ratio of cb to mrs levels.
amplification efficiency for all assays ranged between  <dig> and 105%.
confirming the reliability of cdna-aflp results, all transcripts were more abundant in cb, with expression ratios over 5-fold .
to investigate a possible role for these genes in allowing l. rhamnosus growth in cheese during ripening, in silico analyses were carried out.
spxb
in silico analysis of tdf no.
<dig> , encoding  <dig> amino acid residues, revealed the highest identity in amino acid sequence  with a pyruvate oxidase  from l. rhamnosus gg .
lower levels of identity were observed for spxb of other members of l. casei group .
blastx search also returned a number of pyruvate oxidases of other nslab, such as l. curvatus , l. buchneri , l. brevis , l. plantarum  and l. pentosus , as well as of non-lactobacillus bacteria.
spxb is an enzyme involved in the pyruvate metabolism pathway.
lab can metabolize pyruvate into lactate by lactate dehydrogenase  or into acetate via pyruvate formate lyase , phosphotransacetylase  and acetate kinase , or via pyruvate oxidase  pathway  <cit> .
in the latter, pyruvate is oxidized with the production of hydrogen peroxide and acetyl phosphate, followed by acetate production and atp generation via ack .
to our knowledge, this is the first evidence of activation of the pox pathway in l. rhamnosus.
on the contrary, pox activity has been extensively described to date in l. plantarum and involved with acetate production in its survival during the stationary phase of aerobic growth  <cit> .
in particular, accumulation of acetate instead of lactate is thought to play a role in ensuring the ph homeostasis with an overall beneficial effect for the cell  <cit> .
the additional atp generated via ack has been shown to enhance the biomass production  <cit> .
interestingly, lorquet et al.
<cit>  showed that in the late stationary phase, when the production of acetate stopped, an od decrease resulting from lytic processes occurred.
the hypothesis is that in the absence of atp production, protons can no longer be extruded by atpases with a consequent dissipation of the proton motive force, which has been shown to be one of the mechanisms triggering autolysis of gram-positive bacteria.
interestingly, high levels of acetic acid and low levels of lactic acid have been recently observed in l. rhamnosus strains grown in cb under the same conditions of our study  <cit>  furthermore, by a proteomic approach, bove et al.
<cit>  showed an increase in expression of pta and ack, which are involved in the synthesis of acetic acid in a branch of the pyruvate metabolism other than pox pathway , during l. rhamnosus growth in cb compared to mrs. highlighting a possible alternative route of degradation of pyruvate to acetate , our transcriptomic results seem to complement data from proteomics, strengthening the hypothesis that l. rhamnosus can utilize pyruvate as a growth substrate during cheese ripening.
pyruvate is an intracellular metabolite that could be produced through different metabolic routes using the carbon sources present in cheese .
moreover, pyruvate can be released in the cheese matrix with starter lysis.
liu et al.
<cit>  showed that the activity of pox in l. plantarum could be related to the catabolism of l-serine.
according to the authors, l-serine is deaminated via a serine dehydratase into pyruvate, which is subsequently converted into acetate by the pox enzyme  <cit> .
pyruvate conversion by pox has been recently supposed also in l. casei <cit> .
looking at the closely related taxonomic groups of these heterofermentative lactobacilli and considering that they share similar ecological niches, including ripened cheese, it appears interesting to further investigate in the future the role of pyruvate metabolism during cheese ripening.
notably,  <dig> genes encoding putative pyruvate oxidases are harbored in the completely sequenced genomes of l. rhamnosus gg and l. casei atcc  <dig>  whereas  <dig> and  <dig> pox genes were retrieved in the genome sequences of l. buchneri cd <dig> and l. plantarum wcfs <dig>  respectively.
goffin et al.
<cit>  reported that among the predicted pox genes encoded in the l. plantarum lp <dig> genome, only poxb and poxf appeared to be involved in the generation of acetate from lactate during the stationary phase of aerobic growth.
interestingly, poxb and poxf genes shared  <dig> and 61% amino acid similarity with tdf  <dig>  respectively.
to date, only one gene potentially encoding for pyruvate oxidase has been located in the complete genome sequences of the slab l. helveticus r <dig> and l. delbrueckii subsp.
bulgaricus atcc  <dig>  the pyruvate oxidase gene of l. rhamnosus gg with the highest homology to tdf  <dig> is flanked by genes whose order and transcriptional orientation are partially shared with l. casei atcc  <dig> but not with l. buchneri cd <dig>  l. plantarum wcfs <dig>  l. helveticus r <dig>  l. delbrueckii subsp.
bulgaricus atcc  <dig> and l. brevis atcc  <dig> .
in particular, spxb locus in l. rhamnosus and l. casei genomes is preceded by three genes encoding putative hydroxymethylglutaryl-coa synthase, hydroxymethylglutaryl-coa reductase and acetyl-coa acetyltransferase.
these enzymes are known to be involved in the mevalonate pathway, routing acetyl-coa towards isoprenoid biosynthesis.
however, whether these proteins are actually expressed in l. rhamnosus and play a role in deviating the flow of acetyl-coa from the acetate production via pta and ack during cheese ripening still remain to be determined.
according to pepper, spxb gene from l. rhamnosus gg was predicted to be monocistronically transcribed.
phylogenetic tree showed a clear segregation of putative pyruvate oxidases from l. casei group .
as expected, a subgroup was represented by pox proteins from the slab l. helveticus, l. delbrueckii subsp.
bulgaricus and l. delbrueckii subsp.
lactis.
l. plantarum and l. pentosus homologues clustered together and close to l. buchneri.
multiple sequence alignment of tdf  <dig> and pyruvate oxidase protein sequences from several nslab and slab is shown in additional file 1: figure s1a.
ulae
blastx analysis of tdf no.
<dig> , putatively encoding  <dig> amino acid residues, showed the maximum identity  to a protein annotated as l-xylulose 5-phosphate 3-epimerase  from l. rhamnosus gg .
eighty-four percent of identity was exhibited to the same putative protein from other l. casei group members .
homologues were also found in nslab known to play a role in flavor generation and other ripening processes: l. suebicus , l. coryniformis  and carnobacterium maltaromaticum .
ulae is an epimerase involved with other enzymes  in the production of d-xylulose 5-phosphate  <cit> , an intermediate in the pentose phosphate pathway.
according to synttax, regions up and downstream of ulae gene from l. rhamnosus gg shared a conserved gene order with l. casei atcc  <dig>  whereas no synteny was found in l. buchneri cd <dig>  l. plantarum wcfs <dig>  l. helveticus r <dig>  l. delbrueckii subsp.
bulgaricus atcc  <dig> and l. brevis atcc  <dig> genomes .
according to pepper analysis of l. rhamnosus gg genome, a potential terminator stem-loop structure was identified 82 bp downstream from the arad gene stop codon.
no putative promoters were predicted up to 5000 bp upstream of ulae gene.
interestingly, the upstream lgg_ <dig> gene was annotated as a transcriptional regulator, belonging to deor family.
phylogenetic analysis of l-xylulose 5-phosphate 3-epimerase homologues revealed that ulae predicted protein from l. rhamnosus clustered close to the putative enzymes from other l. casei group members and l. coryniformis .
multiple sequence alignment of tdf  <dig> and homologs from several nslab is shown in additional file 1: figure s1b.
xfp
tdf no.
<dig>  displayed the highest identity  in amino acid sequence with a putative phosphoketolase  from l. rhamnosus gg .
percentages of identity > 95% were found with other l. casei group members .
blastx search also revealed a significant match to a predicted xylulose-5-phosphate phosphoketolase from l. coryniformis .
interestingly, lower levels of identity were obtained with slab, such as l. delbrueckii subsp.
bulgaricus , l. delbrueckii subsp.
lactis  and l. helveticus .
facultatively heterofermentative lab, like l. rhamnosus, degrade hexoses via the embden-meyerhoff-parnas pathway and pentoses via the phosphoketolase pathway .
xylulose 5-phosphate phosphoketolase is the central enzyme of pkp.
in the presence of inorganic phosphate this enzyme converts xylulose 5-phosphate into glyceraldehyde 3-phosphate and acetylphosphate   <cit> .
recently, mcleod et al.
<cit>  studied the transcriptome response of l. sakei during growth on ribose, demonstrating that the ribose uptake and catabolic machinery are highly regulated and closely linked with the catabolism of nucleotides.
it is known that ribonucleosides are source of ribose as a fermentable carbohydrate in environments where free carbohydrates are lacking.
for example, in the meat, a rich environment but carbohydrate-poor substrate for microorganisms, the ability of l. sakei to use nucleosides offers a competitive advantage  <cit> .
nucleosides represent a potential energy source also in the cheese environment, where microbial autolysis occurs, releasing ribose- and desoxyribose-containing nucleic acids  <cit> .
notably, it has been observed that ribose released after lysis of slab decreased steadily in parallel with the growth of facultatively heterofermentative lactobacilli, strongly suggesting that these bacteria used ribose as a growth substrate  <cit> .
the over-expression of xfp mrna levels in l. rhamnosus grown in cb, as found in our study, seems to support this hypothesis.
moreover, our findings are in agreement with the proteomic data of bove and colleagues  <cit> , who observed an increase in expression level of ribose-5-phosphate isomerase  after l. rhamnosus growth in cb compared to mrs. this enzyme acts in a step upstream of xfp in the pathway that leads from ribose 5-phosphate  to the production of acetate, catalyzing the conversion of r5p to ribulose 5-phosphate .
according to pfam search, tdf 40-deduced  <dig> amino acid sequence contains a portion of the xfp c-terminal domain .
the genetic organization and location of xfp gene on l. rhamnosus gg and l. casei atcc  <dig> chromosomes were shown to be highly similar .
in particular, xfp genes are preceded by a divergently transcribed orf, encoding a major facilitator superfamily transporter, and are followed by several genes predicted to encode components of abc transporter and pts systems for sugar uptake.
according to pepper, no high-scoring promoter consensus sequences were identified in the 5000-bp upstream region of xfp gene in l. rhamnosus gg.
a rho-independent transcriptional terminator was predicted to occur 43 bp downstream from the xfp stop codon.
phylogram showed that xfp proteins from l. casei group made a separate cluster, close to the putative enzyme from l. coryniformis .
analogously, different clusters were observed for the slab l. helveticus, l. delbrueckii subsp.
lactis and l. delbrueckii subsp.
bulgaricus.
additional file 1: figure s1c displays a multiple sequence alignment of tdf  <dig> and putative phosphoketolases from several slab and nslab.
CONCLUSIONS
in this study, we applied a transcriptomic approach, based on cdna-aflp and qpcr, to investigate the physiological adaptation of l. rhamnosus to the cheese environment.
l. rhamnosus is known to be one of the few nslab species able to survive and grow during long ripening of sseveral cheeses.
in particular, the strain l. rhamnosus pr <dig>  isolated from 4-month-ripened pr cheese, has previously shown a great ability to growth in cb coupled with high levels of production of acetic acid.
by comparing the gene expression profiles of l. rhamnosus pr <dig> in cb respect to mrs, we identified among others as over-expressed in cb, genes linked to the conversion of pyruvate to acetate as well as to the pathway of ribose degradation.
notably, the activation of pox pathway in l. rhamnosus has never been observed before.
pyruvate is a intracellular metabolite that could be produced by different metabolism using the carbon source present in cheese and can be released in the cheese matrix with the starter lysis.
similarly the ribonucleosides release with starter lysis could be carriers of ribose that represents a fermentable carbohydrate in an environments such cheese where carbohydrates are lacking.
both pyruvate degradation and ribose catabolism induce a metabolite flux toward acetate, coupled with atp production via acetate kinase.
taking into account these consideration, and in agreement with previous findings  <cit>  we assume that l. rhamnosus when growing in media poor in carbohydrates, such as cb, arguably uses different metabolic pathways to produce energy.
notably, the transcriptomic approach employed in this study evidenced the over-expression in cb of enzymes other than those identified through proteomics by bove et al.
<cit> , acting at different steps or in different branches of the ribose and pyruvate utilization pathways.
this discrepancy, probably owing to issues of technique sensitivity and resolution, highlighted the need to integrate transcriptomic and proteomic data in order to get a view as complete as possible of the l. rhamnosus metabolic adaptations during cheese ripening.
since, to our knowledge, this is the first study that showed the activation of pox pathway in l. rhamnosus, further work will be directed to investigate more in depth the role of the pyruvate metabolism in the growth of this specie in cheese.
abbreviations
lab: lactic acid bacteria; slab: starter lactic acid bacteria; nslab: non-starter lactic acid bacteria; l.: lactobacillus; cdna-aflp: cdna-amplified fragment length polymorphism; qpcr: quantitative real-time reverse transcription-pcr; cb: cheese broth; tdfs: transcript-derived fragments; spxb/pox: pyruvate oxidase; ulae: l-xylulose 5-phosphate 3-epimerase; xfp: xylulose-5-phosphate phosphoketolase.
competing interests
the authors declare that they have no competing interests.
author’s contributions
cl: conceived the study, participated in its coordination and drafted the manuscript.
st: carried out genetic and bioinformatic analysis and helped to draft the manuscript.
am: carried out genetic analysis, participated in data collection and their interpretation.
es: carried out genetic analysis, participated in data collection and their interpretation.
en: critically revised the paper.
pb: critically revised the paper.
mg: conceived the study, supervised the research work and critically revised the paper.
all authors read and approved the final manuscript.
supplementary material
additional file 1: figure s1
multiple sequence alignment of spxb , ulae  and xfp  sequences from diverse lactobacilli.
conservation plots and consensus sequences are shown at the bottom.
protein alignments were performed and represented using clc-bio sequence viewer  <cit> .
reference organisms: l. rhamnosus gg, l. casei atcc  <dig>  l. paracasei subsp.
paracasei atcc  <dig>  l. zeae , l. buchneri cd <dig>  l. plantarum wcfs <dig>  l. helveticus r <dig>  l. delbrueckii subsp.
lactis dsm  <dig>  l. delbrueckii subsp.
bulgaricus atcc  <dig>  l. curvatus crl  <dig>  l. brevis atcc  <dig>  l. pentosus kca <dig>  l. coryniformis .
click here for file
acknowledgments
the authors are grateful to dr. claudio giorgio bove for technical assistance.

 lactobacillus rhamnosus is a facultatively heterofermentative lactic acid bacterium frequently encountered in many dairy products, where it can be added as a probiotic microorganism or can be naturally present arising from raw milk. lab may play different roles in cheese manufacture: some species participate in the fermentation process and contribute to acid production acting as starter lab, whereas others, called non-starter lab, are mainly implicated in the maturation process. in particular, l. rhamnosus plays a significant role during ripening, leading to the formation of flavor and, for this reason, members of this species are generally recognized as nslab. it is noteworthy that nslab generally have a high tolerance to hostile environments, such as those with high salt concentration, low moisture, - ph values, low temperatures and deficiency of nutrients. moreover, several studies have reported that in long-ripened cheese varieties, nslab populations dominate during aging after slab decline due to autolysis. increasing by about four to five orders of magnitude within a few months, nslab can have a major impact in determining curd maturation and final characteristics of cheese. in particular, l. rhamnosus has been shown to become dominant within nslab population in several cheeses, including parmigiano reggiano. it persists throughout the whole time of pr cheese ripening and this implies its capacity to adapt to changing environmental conditions. notably, different l. rhamnosus strains have been detected in relation to specific steps of cheese ripening, suggesting that these strains may have specific metabolic activities, which could account for the adaptation to the changing microenvironment of cheese ripening. up to date, the molecular mechanisms allowing the growth of l. rhamnosus in cheese are still poorly understood. nslab development during ripening can be attributed to their ability to use the major nutrient sources available in ripened cheese, that is lactose-free. the nutrient sources available include milk components modified by technological treatment and starter lab development, starter lab metabolites and cell lysis products. potential substrate for microbial growth are represented by small peptides or amino acids, citrate, lactate, and free fatty acids. additionally, sugars and phospholipids, nucleic acids and peptides can be released in the cheese matrix when slab autolysis starts to occur. these compounds represent carbon sources that could yield intracellular pyruvate or be converted into different metabolites. to investigate the metabolic pathways occurring in l. rhamnosus during cheese ripening, bove and colleagues recently compared the proteomic profiles of 10 l. rhamnosus strains grown in mrs and a cheese-like medium. differently from mrs, which is the standard laboratory medium for lactobacilli and is considered to be a rich substrate with glucose as the primary carbon source for microbial growth, cb is an experimental medium formulated with 20-month-ripened pr cheese, that tries to mimic the nutritional composition of pr cheese during ripening. pr cheese, and thus cb, is considered to be a substrate poor in carbohydrates and characterized by the absence of milk sugar, lactose. during the curd acidification step of pr cheese production, the conversion of lactose into lactic acid is the main biochemical process that occurs. lactose is completely depleted within to 48 h. the composition of cb, prepared as the protocol of neviani et al., is the following: proteins,; lactose,; glucose,; galactose,; lactic acid,; nacl,; sodium citrate,. according to the findings of bove et al. compared to the cultivation in mrs, the differentially expressed proteins under cheese-like conditions were mainly linked to protein biosynthesis and catabolism, nucleotide and carbohydrate metabolisms, citrate metabolism, cell wall and exopolysaccharide biosynthesis, cell regulation, oxidation/reduction processes, and stress response. notably, l. rhamnosus produced lactic acid as the primary end product when growing in mrs, whereas in cb low levels of lactic acid together with high levels of acetic acid were detected for all strains. despite this common trend, the authors also observed strain-specific physiological responses, suggesting a strain variability in the adaptation to changing environmental conditions in accordance with genetic polymorphism studies. among all strains, l. rhamnosus pr isolated from pr 4-month-ripened cheese, has shown the greatest ability to growth in cb reaching after 48 h the highest cell numbers with one of the lowest levels of lactic acid and one of the highest levels of acetic acid. for this reason, in the present work, we focused our attention only on this strain, with the aim to identify the genes that could concur to explain its growth ability in cb and its acid acetic production. the physiological adaptation of l. rhamnosus pr in cb was evaluated using a transcriptomic approach, based on cdna-amplified fragment length polymorphism and quantitative real-time reverse transcription-pcr. cdna-aflp is one of the most robust and sensitive transcriptomic technologies for genome-wide expression studies, with the main advantage of not requiring any prior knowledge of gene sequences while allowing the detection of lowly expressed genes through transcript amplification. using this approach, we identified a set of genes resulted over-expressed in cb compared to mrs, potentially involved in alternative metabolic pathways. interesting genes were searched in other nslab and slab genomes with the aim to explore their diversity. overall, the results described in this work highlight mechanisms of adaptation leading to the production of acetic acid coupled with atp generation, that could support the l. rhamnosus growth in cheese during ripening. methods bacterial growth conditions l. rhamnosus pr was isolated from parmigiano reggiano at 4 months of ripening on cheese based medium plate counts and identified by 16s rdna gene sequencing and species-specific pcr. the strain was cultivated in mrs broth or cheese broth at 30°c, under anaerobiosis, for or 48 h, respectively. cb, a culture medium that mimics raw-milk long-ripened cheese, was prepared according to the modified protocol described by bove et al.. rna extraction and cdna synthesis the growth of l. rhamnosus pr in mrs and cb broth was monitored by measuring optical density at 600 nm. about cells at the top of logarithmic phase were harvested, and total rna, stabilized with rnaprotect bacteria reagent, was isolated using rneasy protect bacteria mini kit. three independent biological experiment were made. rna was quantified using a nanodrop nd- spectrophotometer and visualized by formaldehyde agarose gel electrophoresis according to standard procedures. all rnas were of sufficient quantity and high quality. after a step of mrna enrichment and polyadenylation of rna transcripts, cdna was synthesized by reverse transcription using a biotinylated oligo, following the protocol reported by bove et al.. for qpcr, single-stranded cdna was prepared in a final volume of 20 μl from total rna using an oligo primer and superscript™ ii reverse transcriptase, according to the manufacturer’s instructions. briefly, 12-μl reaction mixtures containing 500 ng of oligo primer, 2 μg total rna and 10 nmol dntp mix in depc-treated h2o were heated to 65°c for 5 min, added with 4 μl of 5x first-strand buffer and 200 nmol dtt, and then incubated at 42°c for 2 min. rt reactions were started by the addition of u of enzyme, incubated at 42°c for 50 min and inactivated by heating at 70°c for 15 min. rt step was carried out in duplicate. cdna-aflp cdna-aflp analysis was carried out as described by bove et al.. the protocol is based on the production of cdna-aflp fragments that are detected using infrared dye detection technology and the odyssey infrared imaging system. briefly, after cdna synthesis, a double digestion was carried out with ecori and msei restriction enzymes and fragments were captured with the aid of streptavidin-coated magnetic beads. digested cdna fragments were subsequently ligated with adaptors to allow selective amplification with ecori primers labeled with an infrared dye, and unlabeled msei-n. three primer combinations were used to selectively amplify the expressed genes: dy-ecori-ac/msei-at, dy-ecori-at/msei-ac and dy-ecori-at/msei-at. ligators and primers used are reported in table  separation of cdna-aflp fragments was carried out in a polyacrylamide gel and visualized by odissey at 700 nm. primer sets were designed as reported by bove et al.. cdna-aflp fragment isolation, re-amplification and sequencing transcript-derived fragments of interest were cut from polyacrylamide gels as reported by vuylsteke et al., resuspended in μl of distilled water and subsequently re-amplified using the re-amplification and selective pcr primers ecori-ac/msei-at, ecori-at/msei-ac and ecori-at/msei-at according to the origin of cdna-aflp fragments. amplification reactions were performed in a final volume of μl containing μl of resuspended dna fragment, 25 mm mgcl 10x pcr buffer, μm ecori-n primer, μm msei-n primer, 5 mm dntps, μl of amplitaq dna polymerase and μl of gc enhancer. pcr consisted of: i) 30 s of denaturation step at 94°c, 30 s of annealing step at 65°c and 1 min of extension step at 72°c for total cycles; ii) 10 s of denaturation step at 94°c, 30 s of annealing step at 56°c and 1 min of extension step at 72°c for total cycles. these cycles were preceded by a common denaturation step of 2 min at 94°c and followed by a final 10-min extension at 72°c, and were carried out in a mastercycler ep gradient s thermal cycler. amplified products were checked on a 1% agarose gel with a 100-bp marker and subsequently purified using the wizard® sv gel and pcr clean-up system according to the manufacturer’s instructions. amplified fragments were then cloned in e.coli using the pgem-t easy vector system kit, and plasmids from selected clones were purified using pureyield miniprep system kit referring to the producer’s manual. cloned fragments were finally sequenced by eurofins mwg operon using primers m and sequences were analysed by blast alignment. tdf sequences were deposited in the ddbj database under the accession numbers ab to ab qpcr and data processing qpcr was carried out using the lightcycler sybr green system as previously described. briefly, 1 μl of cdna template was used in each reaction along with 4 μl of sybr green pcr master mix and pmol of the appropriate gene-specific primers in a final volume of 20 μl. the following cycle profile was used: 10 min at 95°c, repeats of 15 s at 95°c, 25 s at 58°c for spxb, ulae and 16s rdna genes or 55°c for xfp, 72°c for 20 s and an additional 5-s incubation step at 81°c for fluorescence acquisition. oligonucleotide sequence information and detailed primer-specific conditions are given in table  two technical replicates were done for each combination of cdna and primer pair. to assess background and residual dna contamination, a no-template control and a no-reverse transcription control were performed for each target. dna contamination was considered to be negligible when the difference in cq between the sample and the respective nort was above 5 cycles. product detection and pcr specificity were checked post-amplification by examining the dissociation curves. pcr amplicons were resolved by 2% agarose gel electrophoresis to verify the expected size. to evaluate repeatability and reproducibility of the qpcr assay, intra- and inter-assay coefficients of variation were assessed. the intra-assay cv was from to % whereas the inter-assay cv ranged from to %. amplification efficiency was calculated from the slope of standard curves generated with two-fold serial dilutions of the same cdna sample, as e =  relative expression of target genes was determined using the ΔΔct method after pfaffl correction. 16s rdna was used as a reference gene. a b aprimer sets were designed based on the sequences of cdna-aflp fragments. primers for 16s rdna gene were designed as reported by giraffa et al.. btarget gene expression was calculated relative to 16s rdna as a reference gene using the efficiency-corrected ΔΔct method. the relative expression ratios in cb compared to mrs are shown. in silico analysis tdf sequences were annotated using blast search. pathway assignment was performed according to cog functional categories and kegg pathway database. gene synteny across nslab and slab genomes was explored through the web server synttax. genome mining for promoter and terminator elements was performed using pepper toolbox. translated protein sequences were subjected to pfam motif analysis. protein alignments were performed using clustalw and used for phylogenetic tree construction at the interactive tree of life. multisequence amino acid alignments were represented using clc-bio sequence viewer.  cdna-aflp analysis in this study, the cdna-aflp technique was applied to profile the transcriptome of a l. rhamnosus strain grown in conditions mimicking cheese ripening. despite it is not widely used in bacteria, cdna-aflp can be considered an ideal system for genome-wide expression analysis, mainly for the detection of lowly expressed genes. three primer combinations were used to selectively amplify the genes expressed by l. rhamnosus pr in cb and mrs, allowing to generate different cdna-aflp profiles with a fragment size ranging from to 500 bp. a total of and tdfs were detected in mrs and cb, respectively. in order to investigate the main adaptations of l. rhamnosus to the pr cheese environment, we focused on tdfs over-expressed in cb. identification of tdfs over-expressed in cb twenty tdfs strongly over-expressed by l. rhamnosus in cb compared to mrs were extracted from gel and used as templates for re-amplification by pcr. direct sequencing of the pcr products failed probably because the sequences were contaminated with co-migrating bands or were too small. products from bands were then cloned and pcr amplicons from individual colonies were sequenced. reliable sequences were obtained for of the tdfs. each sequence was identified by similarity search using the blast program against the genbank non-redundant public sequence database. as shown in table  transcripts showed a significant similarity to sequences with known function, transcripts were closely related to l. casei plasmid sequences, transcripts were annotated as hypothetical protein-coding sequences, and transcripts were identified as 5s rrna. l. rhamnosus pr over-expressed in cb compared to mrs a b c d cog e awhen available, ec numbers assigned to the putative enzymatic reactions are provided. bthe column indicates the microorganism of the best hit from blastx search. cmax identity and e-value from the best hit of blastx search are provided. dpathway assignment was performed according to cog functional categories and kegg pathway database. ee, amino acid transport and metabolism; f, nucleotide transport and metabolism; g, carbohydrate transport and metabolism; m, cell wall/membrane/envelope biogenesis; r, general function prediction only. it is known that plasmids often carry genes that might be essential for survival under harsh conditions, encoding important traits, such as enzymes involved in secondary metabolic pathways. plasmids are known to be a source of lab genetic and phenotypic diversity which occasionally confers adaptive advantages to host strains. however, further studies are clearly needed to better explore the role of plasmid sequences in the l. rhamnosus adaptation to the cheese ripening environment. to validate the cdna-aflp expression profiles, genes, encoding pyruvate oxidase, l-xylulose 5-phosphate 3-epimerase, and xylulose-5-phosphate phosphoketolase were selected for qpcr. the relative mrna abundances were normalized by that of the commonly used reference gene 16s rdna, and expressed as a ratio of cb to mrs levels. amplification efficiency for all assays ranged between and 105%. confirming the reliability of cdna-aflp results, all transcripts were more abundant in cb, with expression ratios over 5-fold. to investigate a possible role for these genes in allowing l. rhamnosus growth in cheese during ripening, in silico analyses were carried out. spxb in silico analysis of tdf no., encoding amino acid residues, revealed the highest identity in amino acid sequence with a pyruvate oxidase from l. rhamnosus gg. lower levels of identity were observed for spxb of other members of l. casei group. blastx search also returned a number of pyruvate oxidases of other nslab, such as l. curvatus, l. buchneri, l. brevis, l. plantarum and l. pentosus, as well as of non-lactobacillus bacteria. spxb is an enzyme involved in the pyruvate metabolism pathway. lab can metabolize pyruvate into lactate by lactate dehydrogenase or into acetate via pyruvate formate lyase, phosphotransacetylase and acetate kinase, or via pyruvate oxidase pathway. in the latter, pyruvate is oxidized with the production of hydrogen peroxide and acetyl phosphate, followed by acetate production and atp generation via ack. to our knowledge, this is the first evidence of activation of the pox pathway in l. rhamnosus. on the contrary, pox activity has been extensively described to date in l. plantarum and involved with acetate production in its survival during the stationary phase of aerobic growth. in particular, accumulation of acetate instead of lactate is thought to play a role in ensuring the ph homeostasis with an overall beneficial effect for the cell. the additional atp generated via ack has been shown to enhance the biomass production. interestingly, lorquet et al. showed that in the late stationary phase, when the production of acetate stopped, an od decrease resulting from lytic processes occurred. the hypothesis is that in the absence of atp production, protons can no longer be extruded by atpases with a consequent dissipation of the proton motive force, which has been shown to be one of the mechanisms triggering autolysis of gram-positive bacteria. interestingly, high levels of acetic acid and low levels of lactic acid have been recently observed in l. rhamnosus strains grown in cb under the same conditions of our study furthermore, by a proteomic approach, bove et al. showed an increase in expression of pta and ack, which are involved in the synthesis of acetic acid in a branch of the pyruvate metabolism other than pox pathway, during l. rhamnosus growth in cb compared to mrs. highlighting a possible alternative route of degradation of pyruvate to acetate, our transcriptomic results seem to complement data from proteomics, strengthening the hypothesis that l. rhamnosus can utilize pyruvate as a growth substrate during cheese ripening. pyruvate is an intracellular metabolite that could be produced through different metabolic routes using the carbon sources present in cheese. moreover, pyruvate can be released in the cheese matrix with starter lysis. liu et al. showed that the activity of pox in l. plantarum could be related to the catabolism of l-serine. according to the authors, l-serine is deaminated via a serine dehydratase into pyruvate, which is subsequently converted into acetate by the pox enzyme. pyruvate conversion by pox has been recently supposed also in l. casei. looking at the closely related taxonomic groups of these heterofermentative lactobacilli and considering that they share similar ecological niches, including ripened cheese, it appears interesting to further investigate in the future the role of pyruvate metabolism during cheese ripening. notably, genes encoding putative pyruvate oxidases are harbored in the completely sequenced genomes of l. rhamnosus gg and l. casei atcc whereas and pox genes were retrieved in the genome sequences of l. buchneri cd and l. plantarum wcfs respectively. goffin et al. reported that among the predicted pox genes encoded in the l. plantarum lp genome, only poxb and poxf appeared to be involved in the generation of acetate from lactate during the stationary phase of aerobic growth. interestingly, poxb and poxf genes shared and 61% amino acid similarity with tdf respectively. to date, only one gene potentially encoding for pyruvate oxidase has been located in the complete genome sequences of the slab l. helveticus r and l. delbrueckii subsp. bulgaricus atcc the pyruvate oxidase gene of l. rhamnosus gg with the highest homology to tdf is flanked by genes whose order and transcriptional orientation are partially shared with l. casei atcc but not with l. buchneri cd l. plantarum wcfs l. helveticus r l. delbrueckii subsp. bulgaricus atcc and l. brevis atcc. in particular, spxb locus in l. rhamnosus and l. casei genomes is preceded by three genes encoding putative hydroxymethylglutaryl-coa synthase, hydroxymethylglutaryl-coa reductase and acetyl-coa acetyltransferase. these enzymes are known to be involved in the mevalonate pathway, routing acetyl-coa towards isoprenoid biosynthesis. however, whether these proteins are actually expressed in l. rhamnosus and play a role in deviating the flow of acetyl-coa from the acetate production via pta and ack during cheese ripening still remain to be determined. according to pepper, spxb gene from l. rhamnosus gg was predicted to be monocistronically transcribed. phylogenetic tree showed a clear segregation of putative pyruvate oxidases from l. casei group. as expected, a subgroup was represented by pox proteins from the slab l. helveticus, l. delbrueckii subsp. bulgaricus and l. delbrueckii subsp. lactis. l. plantarum and l. pentosus homologues clustered together and close to l. buchneri. multiple sequence alignment of tdf and pyruvate oxidase protein sequences from several nslab and slab is shown in additional file 1: figure s1a. ulae blastx analysis of tdf no., putatively encoding amino acid residues, showed the maximum identity to a protein annotated as l-xylulose 5-phosphate 3-epimerase from l. rhamnosus gg. eighty-four percent of identity was exhibited to the same putative protein from other l. casei group members. homologues were also found in nslab known to play a role in flavor generation and other ripening processes: l. suebicus, l. coryniformis and carnobacterium maltaromaticum. ulae is an epimerase involved with other enzymes in the production of d-xylulose 5-phosphate, an intermediate in the pentose phosphate pathway. according to synttax, regions up and downstream of ulae gene from l. rhamnosus gg shared a conserved gene order with l. casei atcc whereas no synteny was found in l. buchneri cd l. plantarum wcfs l. helveticus r l. delbrueckii subsp. bulgaricus atcc and l. brevis atcc genomes. according to pepper analysis of l. rhamnosus gg genome, a potential terminator stem-loop structure was identified 82 bp downstream from the arad gene stop codon. no putative promoters were predicted up to 5000 bp upstream of ulae gene. interestingly, the upstream lgg_ gene was annotated as a transcriptional regulator, belonging to deor family. phylogenetic analysis of l-xylulose 5-phosphate 3-epimerase homologues revealed that ulae predicted protein from l. rhamnosus clustered close to the putative enzymes from other l. casei group members and l. coryniformis. multiple sequence alignment of tdf and homologs from several nslab is shown in additional file 1: figure s1b. xfp tdf no. displayed the highest identity in amino acid sequence with a putative phosphoketolase from l. rhamnosus gg. percentages of identity > 95% were found with other l. casei group members. blastx search also revealed a significant match to a predicted xylulose-5-phosphate phosphoketolase from l. coryniformis. interestingly, lower levels of identity were obtained with slab, such as l. delbrueckii subsp. bulgaricus, l. delbrueckii subsp. lactis and l. helveticus. facultatively heterofermentative lab, like l. rhamnosus, degrade hexoses via the embden-meyerhoff-parnas pathway and pentoses via the phosphoketolase pathway. xylulose 5-phosphate phosphoketolase is the central enzyme of pkp. in the presence of inorganic phosphate this enzyme converts xylulose 5-phosphate into glyceraldehyde 3-phosphate and acetylphosphate. recently, mcleod et al. studied the transcriptome response of l. sakei during growth on ribose, demonstrating that the ribose uptake and catabolic machinery are highly regulated and closely linked with the catabolism of nucleotides. it is known that ribonucleosides are source of ribose as a fermentable carbohydrate in environments where free carbohydrates are lacking. for example, in the meat, a rich environment but carbohydrate-poor substrate for microorganisms, the ability of l. sakei to use nucleosides offers a competitive advantage. nucleosides represent a potential energy source also in the cheese environment, where microbial autolysis occurs, releasing ribose- and desoxyribose-containing nucleic acids. notably, it has been observed that ribose released after lysis of slab decreased steadily in parallel with the growth of facultatively heterofermentative lactobacilli, strongly suggesting that these bacteria used ribose as a growth substrate. the over-expression of xfp mrna levels in l. rhamnosus grown in cb, as found in our study, seems to support this hypothesis. moreover, our findings are in agreement with the proteomic data of bove and colleagues, who observed an increase in expression level of ribose-5-phosphate isomerase after l. rhamnosus growth in cb compared to mrs. this enzyme acts in a step upstream of xfp in the pathway that leads from ribose 5-phosphate to the production of acetate, catalyzing the conversion of r5p to ribulose 5-phosphate. according to pfam search, tdf 40-deduced amino acid sequence contains a portion of the xfp c-terminal domain. the genetic organization and location of xfp gene on l. rhamnosus gg and l. casei atcc chromosomes were shown to be highly similar. in particular, xfp genes are preceded by a divergently transcribed orf, encoding a major facilitator superfamily transporter, and are followed by several genes predicted to encode components of abc transporter and pts systems for sugar uptake. according to pepper, no high-scoring promoter consensus sequences were identified in the 5000-bp upstream region of xfp gene in l. rhamnosus gg. a rho-independent transcriptional terminator was predicted to occur 43 bp downstream from the xfp stop codon. phylogram showed that xfp proteins from l. casei group made a separate cluster, close to the putative enzyme from l. coryniformis. analogously, different clusters were observed for the slab l. helveticus, l. delbrueckii subsp. lactis and l. delbrueckii subsp. bulgaricus. additional file 1: figure s1c displays a multiple sequence alignment of tdf and putative phosphoketolases from several slab and nslab.  in this study, we applied a transcriptomic approach, based on cdna-aflp and qpcr, to investigate the physiological adaptation of l. rhamnosus to the cheese environment. l. rhamnosus is known to be one of the few nslab species able to survive and grow during long ripening of sseveral cheeses. in particular, the strain l. rhamnosus pr isolated from 4-month-ripened pr cheese, has previously shown a great ability to growth in cb coupled with high levels of production of acetic acid. by comparing the gene expression profiles of l. rhamnosus pr in cb respect to mrs, we identified among others as over-expressed in cb, genes linked to the conversion of pyruvate to acetate as well as to the pathway of ribose degradation. notably, the activation of pox pathway in l. rhamnosus has never been observed before. pyruvate is a intracellular metabolite that could be produced by different metabolism using the carbon source present in cheese and can be released in the cheese matrix with the starter lysis. similarly the ribonucleosides release with starter lysis could be carriers of ribose that represents a fermentable carbohydrate in an environments such cheese where carbohydrates are lacking. both pyruvate degradation and ribose catabolism induce a metabolite flux toward acetate, coupled with atp production via acetate kinase. taking into account these consideration, and in agreement with previous findings we assume that l. rhamnosus when growing in media poor in carbohydrates, such as cb, arguably uses different metabolic pathways to produce energy. notably, the transcriptomic approach employed in this study evidenced the over-expression in cb of enzymes other than those identified through proteomics by bove et al., acting at different steps or in different branches of the ribose and pyruvate utilization pathways. this discrepancy, probably owing to issues of technique sensitivity and resolution, highlighted the need to integrate transcriptomic and proteomic data in order to get a view as complete as possible of the l. rhamnosus metabolic adaptations during cheese ripening. since, to our knowledge, this is the first study that showed the activation of pox pathway in l. rhamnosus, further work will be directed to investigate more in depth the role of the pyruvate metabolism in the growth of this specie in cheese. abbreviations lab: lactic acid bacteria; slab: starter lactic acid bacteria; nslab: non-starter lactic acid bacteria; l.: lactobacillus; cdna-aflp: cdna-amplified fragment length polymorphism; qpcr: quantitative real-time reverse transcription-pcr; cb: cheese broth; tdfs: transcript-derived fragments; spxb/pox: pyruvate oxidase; ulae: l-xylulose 5-phosphate 3-epimerase; xfp: xylulose-5-phosphate phosphoketolase. competing interests the authors declare that they have no competing interests. author’s contributions cl: conceived the study, participated in its coordination and drafted the manuscript. st: carried out genetic and bioinformatic analysis and helped to draft the manuscript. am: carried out genetic analysis, participated in data collection and their interpretation. es: carried out genetic analysis, participated in data collection and their interpretation. en: critically revised the paper. pb: critically revised the paper. mg: conceived the study, supervised the research work and critically revised the paper. all authors read and approved the final manuscript. supplementary material additional file 1: figure s1 multiple sequence alignment of spxb, ulae and xfp sequences from diverse lactobacilli. conservation plots and consensus sequences are shown at the bottom. protein alignments were performed and represented using clc-bio sequence viewer. reference organisms: l. rhamnosus gg, l. casei atcc l. paracasei subsp. paracasei atcc l. zeae, l. buchneri cd l. plantarum wcfs l. helveticus r l. delbrueckii subsp. lactis dsm l. delbrueckii subsp. bulgaricus atcc l. curvatus crl l. brevis atcc l. pentosus kca l. coryniformis. click here for file acknowledgments the authors are grateful to dr. claudio giorgio bove for technical assistance.
the aim of the present study was to identify genes potentially contributing to the growth ability of l. rhamnosus pr in cheese-like medium using a transcriptomic approach, based on cdna-amplified fragment length polymorphism and quantitative real-time reverse transcription-pcr. lactobacillus rhamnosuscheesecdna-amplified fragment length polymorphismquantitative real-time reverse transcription-pcr pyruvate is a intracellular metabolite that can be produced through different metabolic routes starting from the carbon sources present in cheese, and can be released in the cheese matrix with the starter lysis. the cdna-aflp results were validated on selected regulated genes by qpcr. in long-ripened cheeses it persists throughout the whole time of ripening due to its capacity to adapt to changing environmental conditions. the versatile adaptability of l. rhamnosus to different ecosystems has been associated with the capacity to use non-conventional energy sources, regulating different metabolic pathways. however, the molecular mechanisms allowing the growth of l. rhamnosus in the cheese dairy environment are still poorly understood. using three primer combinations, a total of and transcript-derived fragments were obtained for l. rhamnosus pr grown in commercial mrs medium and cb, respectively. taking into account these considerations, we suggest that the energy produced through these pathways may concur to explain the great ability of l. rhamnosus pr to grow on cb. by a transcriptomic approach we identified a set of genes involved in alternative metabolic pathways in l. rhamnosus that could be responsible for l. rhamnosus growth in cheese during ripening. both pyruvate degradation and ribose catabolism induce a metabolite flux toward acetate, coupled with atp production via acetate kinase. lactobacillus rhamnosus is a non-starter lactic acid bacterium that plays a significant role during cheese ripening, leading to the formation of flavor. similarly the ribonucleosides released with starter lysis could deliver ribose that represents a fermentable carbohydrate in environments, such as cheese, where free carbohydrates are lacking. in order to investigate the main adaptations to growth in a cheese-mimicking system, we focused on transcripts over-expressed in cb with respect to mrs. it is worth noting the presence of transcripts involved in the degradation of pyruvate and ribose.
 lactobacillus rhamnosus is a non-starter lactic acid bacterium that plays a significant role during cheese ripening, leading to the formation of flavor. in long-ripened cheeses it persists throughout the whole time of ripening due to its capacity to adapt to changing environmental conditions. the versatile adaptability of l. rhamnosus to different ecosystems has been associated with the capacity to use non-conventional energy sources, regulating different metabolic pathways. however, the molecular mechanisms allowing the growth of l. rhamnosus in the cheese dairy environment are still poorly understood. the aim of the present study was to identify genes potentially contributing to the growth ability of l. rhamnosus pr in cheese-like medium using a transcriptomic approach, based on cdna-amplified fragment length polymorphism and quantitative real-time reverse transcription-pcr.  using three primer combinations, a total of and transcript-derived fragments were obtained for l. rhamnosus pr grown in commercial mrs medium and cb, respectively. the cdna-aflp results were validated on selected regulated genes by qpcr. in order to investigate the main adaptations to growth in a cheese-mimicking system, we focused on transcripts over-expressed in cb with respect to mrs. it is worth noting the presence of transcripts involved in the degradation of pyruvate and ribose. pyruvate is a intracellular metabolite that can be produced through different metabolic routes starting from the carbon sources present in cheese, and can be released in the cheese matrix with the starter lysis. similarly the ribonucleosides released with starter lysis could deliver ribose that represents a fermentable carbohydrate in environments, such as cheese, where free carbohydrates are lacking. both pyruvate degradation and ribose catabolism induce a metabolite flux toward acetate, coupled with atp production via acetate kinase. taking into account these considerations, we suggest that the energy produced through these pathways may concur to explain the great ability of l. rhamnosus pr to grow on cb.  by a transcriptomic approach we identified a set of genes involved in alternative metabolic pathways in l. rhamnosus that could be responsible for l. rhamnosus growth in cheese during ripening. lactobacillus rhamnosuscheesecdna-amplified fragment length polymorphismquantitative real-time reverse transcription-pcr 