Transcriptome analysis revealed growth phase-associated changes of a centenarian-originated probiotic Bifidobacterium animalis subsp. lactis A6
BMC Microbiology volume 22, Article number: 61 (2022)
The physiology and application characteristics of probiotics are closely associated with the growth phase. Bifidobacterium animalis subsp. lactis A6 is a promising probiotic strain isolated from the feces of a healthy centenarian in China. In this study, RNA-seq was carried out to investigate the metabolic mechanism between the exponential and the stationary phase in B. lactis A6.
Differential expression analysis showed that a total of 815 genes were significantly changed in the stationary phase compared to the exponential phase, which consisted of 399 up-regulated and 416 down-regulated genes. The results showed that the transport and metabolism of cellobiose, xylooligosaccharides and raffinose were enhanced at the stationary phase, which expanded carbon source utilizing profile to confront with glucose consumption. Meanwhile, genes involved in cysteine-cystathionine-cycle (CCC) pathway, glutamate dehydrogenase, branched-chain amino acids (BCAAs) biosynthesis, and Clp protease were all up-regulated in the stationary phase, which may enhance the acid tolerance of B. lactis A6 during stationary phase. Acid tolerance assay indicated that the survival rate of stationary phase cells was 51.07% after treatment by pH 3.0 for 2h, which was 730-fold higher than that of 0.07% with log phase cells. In addition, peptidoglycan biosynthesis was significantly repressed, which is comparable with the decreased growth rate during the stationary phase. Remarkably, a putative gene cluster encoding Tad pili was up-regulated by 6.5 to 12.1-fold, which is consistent with the significantly increased adhesion rate to mucin from 2.38% to 4.90% during the transition from the exponential phase to the stationary phase.
This study reported growth phase-associated changes of B. lactis A6 during fermentation, including expanded carbon source utilizing profile, enhanced acid tolerance, and up-regulated Tad pili gene cluster responsible for bacterial adhesion in the stationary phase. These findings provide a novel insight into the growth phase associated characteristics in B. lactis A6 and provide valuable information for further application in the food industry.
Probiotics are defined as “live microorganisms which when administered in adequate amounts confer a health benefit on the host” . Some probiotics confer a number of health-promoting benefits to human host, such as competitive exclusion of pathogenic bacteria, modulation of immune system and enhancement of epithelial barrier function [2,3,4]. Therefore, probiotics, mainly including lactobacilli and bifidobacteria, have been widely applied as food and dietary supplements . However, the physiology and application characteristics of probiotics are closely associated with the growth phase. After spray drying, the viability of Lactobacillus rhamnosus GG was approximately 20 times higher at stationary phase than that of cells from the exponential phase . Meanwhile, the maximum adhesion ability of L. rhamnosus GG to Caco-2 cells at the early stationary phase was ten times higher than that of cells at the exponential phase . Therefore, it is noteworthy to investigate the growth phase-associated metabolic mechanisms in probiotics.
Omics techniques provide an effective strategy to reveal growth phase-associated metabolic mechanisms. Transcriptomic analysis revealed a shift from the cell division gene expression at the lag phase to a carbohydrate metabolism-related gene repertoire at the exponential growth phase in Bifidobacterium bifidum PRL2010 . Meanwhile, transcriptomic analysis of Lactobacillus casei Zhang showed that genes involved in carbohydrate metabolism, inorganic ion transport, and chaperones were highly expressed at the stationary phase, whereas genes related to nucleotide transport and metabolism, energy production and conversion were predominantly expressed during the exponential phase . Furthermore, combined transcriptomic and proteomic analyses demonstrated that the shift from glucose fermentation to galactose utilization and the transition from homolactic to mixed acid fermentation were observed from the exponential to the stationary growth phase in L. rhamnosus GG .
Bifidobacterium animalis subsp. lactis A6 (B. lactis A6) was isolated from a healthy centenarian in the Bama County of the Guangxi Zhuang Autonomous Region in China , which is famous for having a population with a high life-expectancy. Previous study showed that B. lactis A6 had high acid resistance to low pH and could improve obesity in mice [12, 13]. Unlike DNA microarray method which is constrained by the certain levels of background signal and the limited dynamic detection range, RNA-Seq, as the Next-Generation Sequencing technology, provides higher efficiency and sensitivity for microbial transcriptomic analysis . In this study, RNA-seq was performed to investigate the growth phase-related changes between the exponential and stationary phase in B. lactis A6. Our results will provide new insights for the rational application of B. lactis A6 as probiotics.
Growth of B. lactis A6 strain
In order to investigate the changes of growth situation, B. lactis A6 was cultured in 15 ml de Man-Rogosa-Sharp broth supplemented with 0.05% (w/v) L-cysteine (MRSc) for 24 h. As shown in Fig. 1, the OD600 value reached 5.4 at the end-point of fermentation at 24 h, and the pH was reduced to 4.6 from an initial value of 6.5. The bacteria entered into the exponential growth phase at time point T0 with OD600 ~ 0.6 and the stationary growth phase at time point T1 with OD600 ~ 4.0.
Global gene expression profiles in the stationary phase compared to the exponential phase
The number of clean reads and the reads mapping to the reference genome are shown in Table S1 and S2. After removing low quality reads (a single sequencing read contains more than 50% bases with Phred Quality Score ≤ 20) and adaptor sequences, more than 19.6 million high-quality reads per sample were generated. The clean reads were aligned to the whole reference genome sequence, and more than 99% of the clean reads for each sample were mapped to the genome. Correlation of gene expression level between three independent biological replicates was shown by the Pearson’s correlation coefficient, which was more than 0.872 for each group (Fig. S1). These results indicated the obtained transcriptome data were suitable for further analysis. Differential expression analysis showed that a total of 815 differential expression genes (DEGs) consisted of 399 up-regulated and 416 down-regulated genes in the stationary phase (A6_WD) compared to the exponential phase (A6_log) (Fig. 2A). The Veen diagram showed 1651 co-expressed genes (FPKM value > 1) both in the stationary and log phase, 18 specifically expressed genes (SEGs) in the stationary phase, and 1 SEG in the log phase (Fig. 2B). Among these co-expressed genes, ~ 60 top expressed genes are summarized in Fig. 2C, which are mainly involved in the core metabolic pathways, including transcription and translation (ribosomal protein), energy production (ATP synthase and glycolysis), and redox reaction (oxidoreductase). The SEGs in each growth phase are listed in Table S3, however, most of these genes encoded hypothetical proteins. The KEGG pathway enrichment analysis of A6_WD compared to A6_log was summarized separately in Fig. 2D. It revealed that the up-regulated genes are mainly participated in the branched-chain amino acids biosynthesis, ABC transporters and metabolic pathways, while down-regulated genes are associated with ribosome, secondary metabolites biosynthesis, and carbon metabolism (glycolysis pathway). Moreover, a list of predicted promoters of constitutively or specifically expressed genes in both growth phases was summarized in the Supplemental profile. Based on the KEGG pathway analysis, the DEGs involved in carbohydrate uptake and metabolism, oligopeptides uptake, peptidoglycan biosynthesis, acid tolerance resistance and adhesion to mucin were discussed in detail as follows.
Carbon source utilization profile was expanded during the stationary phase
In this study, genes BAA6_0517-0518 and BAA6_1585-1587 encoding two ABC-type transporters for xylooligosaccharides (XOS) and raffinose transport were up-regulated over 2.0-fold during the stationary phase, however, genes BAA6_0488-0491, BAA6_1555, and BAA6_1548-1549 encoding three ABC-type transporters for galactooligosaccharides (GOS), maltotriose, and pullulan transport were down-regulated over 1.8-fold. Interestingly, BAA6_1415 encoding a putative cellobiose ABC transporter permease protein CebG and BAA6_1416 encoding cellobiose phosphorylase CbpA were found simultaneously up-regulated over 3.6-fold in the stationary phase. Exogenous cellobiose is transported into the bacterial cell by CebG, and then catalyzed by CbpA to yield one molecule of glucose-1-phosphate and glucose . Moreover, genes BAA6_0481 and BAA6_0050 encoding two putative MFS transporters for the uptake of GOS and gentiobiose were up-regulated 6.2 and 3.6-fold, respectively (Table S4 and Fig. 3). A phenomenon of selectively expressed genes involved in uptake of exogenous carbon sources was observed during transition from the exponential phase to the stationary phase in B. lactis A6.
For glucose transport and metabolism, gene BAA6_1092 encoding a putative glucose uptake permease was down-regulated 6.1-fold during the stationary phase. Meanwhile, eighteen genes involved in glucose metabolism called “the bifid shunt” were also down-regulated 1.5~3.6-fold during the stationary phase (Table S4 and Fig. 3). Among them, gene BAA6_0968 encoding a bifunctional fructose-6-phosphate/xylulose-5-phosphate phosphoketolase (F6PPK/Xfp) was decreased 1.7-fold, which catalyzes the conversion of xylulose 5-phosphate or fructose 6-phosphate to acetyl phosphate . Furthermore, gene pyk, ldh2, pfl, and ackA encoding pyruvate kinase, lactate dehydrogenase, formate acetyltransferase and acetate kinase were down-regulated 1.5, 2.6, 3.6 and 1.7-fold, respectively, which are involved in the production of pyruvate, L-lactate, formate and acetate . Remarkably, genes araA, araB, and araD for L-arabinose metabolism to xylulose-5-P and gene xylB for xylulose metabolism to xylulose-5-phosphate were all up-regulated over 2.0-fold at the stationary phase . These results showed that the glycolysis pathway was slowed down during the stationary phase with glucose consumption, but B. lactis A6 tends to increase gene expression for selective utilization of exogeneous carbon sources, such as XOS, cellobiose, and L-arabinose.
Oligopeptides uptake was enhanced to adapt to changes in growth
Genes aspC, metE, pheA, argFGH, proC, and glnA involved in aspartate, methionine, tyrosine, phenylalanine, arginine, proline, and glutamine biosynthesis were down-regulated 1.6~16.4-fold during the stationary phase (Table S4). Meanwhile, genes BAA6_1269-1270, BAA6_0680-0683, and BAA6_1082-1083 encoding three pairs of ABC-type amino acid transporters for aspartate, glutamate and methionine uptake were also down-regulated 2.0~4.5-fold. However, the opp operon, BAA6_0566-0570, encoding an ABC oligopeptide transporter OppABCDF was upregulated 2.0~5.2-fold, which are involved in the uptake of oligopeptides for bacteria . Simultaneously, genes dppD and dppC encoding two components of dipeptide transporter DppABCDF complex were increased 8.4 and 5.8-fold, respectively. DppABCDF is specialized in dipeptide transportation to support bacterial growth in Escherichia coli . Remarkably, genes BAA6_0186 and BAA6_0230 encoding two dipeptidases were all up-regulated over 1.7-fold. Dipeptidase catalyzes the hydrolysis of dipeptides to produce free amino acids for bacterial growth . Therefore, B. lactis A6 tends to utilize exogenous amino acids by oligopeptides uptake, rather than de novo synthesis of single amino acid during the stationary phase.
Peptidoglycan biosynthesis were repressed in the stationary phase
Gene BAA6_0562 encoding glutamine-fructose-6-phosphate transaminase was down-regulated 2.2-fold during the stationary phase, which catalyzes the first step of peptidoglycan biosynthesis from fructose-6P (Fig. S2). The formed glucosamine-6-phosphate was then catalyzed by phosphoglucosamine mutase encoded by BAA6_1325 to glucosamine-1-phosphate (GlcN-1P). The gene glmU encoding a bifunctional protein GlmU was also down-regulated 3.1-fold, which catalyzes the biosynthesis of UDP-N-acetylglucosamine from GlcN-1P . Additionally, murC, murD, and murF were down-regulated over 1.5-fold, which are involved in the addition of L-alanine, D-glutamate, and D-alanyl-D-alanine to UDP-N-acetylmuramate to form the skeletal structure of a DAP-type peptidoglycan. Moreover, gene BAA6_0083 encoding transpeptidase (DD-TPase) and murJ encoding a putative peptidoglycan lipid II flippase were also down-regulated over 1.6-fold. DD-TPase was reported to synthesize cross-linked peptidoglycan from lipid intermediates, and then the lipid-linked peptidoglycan precursors were transported from the inner cytoplasmic membrane to the outer leaflet by flippase during cell wall formation [23, 24]. Furthermore, gene BAA6_0084 encoding a bacterial cell division membrane protein and genes BAA6_0550, BAA6_1198, and BAA6_1201 encoding cell division protein FtsX, FtsQ, and FtsW were all down-regulated around 2.0~3.6-fold during the stationary phase (Table S4). Taken all together, genes involved in peptidoglycan biosynthesis, a key cell wall component in bifidobacteria, and cell division were generally repressed during the stationary phase in B. lactis A6, reflecting a retarded growth situation at the stationary growth phase.
Acid tolerance response
Gene BAA6_0572 encoding cystathionine beta-synthase (CysK) and BAA6_0947 encoding cystathionine beta-lyase (MetC) were upregulated 7.1 and 2.8-fold, respectively. (Table S4). CysK catalyzes the synthesis of cystathionine from homocysteine and serine, and MetC catalyzes the cleavage of cystathionine to homocysteine, pyruvate and NH3, which constitute the cysteine-cystathionine-cycle (CCC) in response to acid stress . Meanwhile, gene serB encoding phosphoserine phosphatase for the biosynthesis of serine was also up-regulated 1.8-fold during the stationary phase. Moreover, gene BAA6_0010 encoding a NADP-specific glutamate dehydrogenase (GdhA) was up-regulated 3.2-fold, which catalyzes glutamate to 2-oxoglutarate and NH3. In contrast, gene BAA6_0708 and BAA6_1220 encoding glutamine synthetase 1 and 2 were respectively down-regulated 7.4- and 2.2-fold, which are involved in catalyzing the biosynthesis of glutamine from glutamate and NH3. In addition, gene livK encoding a substrate-binding protein for the transport of branched-chain amino acids was up-regulated 3.1-fold during the stationary phase. Meanwhile, BAA6_0283 encoding acetolactate synthase and BAA6_0146 encoding ketol-acid reductoisomerase were also up-regulated 2.6- and 4.1-fold, respectively. Acetolactate synthase catalyzes the biosynthesis of acetolactate and 2-aceto-2-hydroxybutanoate from pyruvate, which are further catalyzed by ketol-acid reductoisomerase for valine, leucine and isoleucine biosynthesis . Meanwhile, BAA6_0675 encoding ATP-binding subunit ClpE and BAA6_1072~1073 encoding two Clp protease proteolytic subunits ClpP1 and ClpP2 were also up-regulated 1.6, 1.6 and 3.3-fold, respectively (Table S4). Thus, genes involved in cysteine-cystathionine-cycle (CCC) pathway, glutamate dehydrogenase, branched-chain amino acids (BCAAs) biosynthesis, and Clp protease were all up-regulated in the stationary phase, which may contribute to enhance acid tolerance of B. lactis A6 during fermentation.
To investigate the general acid tolerance of B. lactis A6 under the log and stationary growth phase, an average of ~108 CFU/mL bacterial cells of each growth phase were treated with acidic pH 4.0, 3.0, 2.0 for 2 h, respectively. As shown in Fig. 4, no significant difference was observed in survival rate between the two growth phases after treatment with pH 4.0 for 2 h. However, the survival rate of log phase cells was dramatically decreased to 0.07% after treatment with pH 3.0 for 2 h, significantly lower than that of 51.07% with stationary phase cells. No viable cells were detectable in log phase cells after treatment with pH 2.0 for 2 h, while ~5 log CFU/mL decrease of viable cells was observed in stationary phase cells from 8.86 to 3.99 log CFU/mL and the survival rate was 0.002%. These results suggested B. lactis A6 cells in stationary phase are more acid-tolerant than the log phase.
Adherence to mucin was enhanced during the stationary phase
In this study, a putative gene cluster of BAA6_0205~0212 involved in type IVb tight adherence (Tad) pili biosynthesis was found up-regulated 7.0~12.4-fold during the stationary phase compared to the log phase (Table S4). Among them, genes flp, tadE and tadF encoding a fimbrial protein prepilin and two pseudopilins were upregulated 7.5, 7.7 and 8.0-fold, respectively, which are participated in prepilin proteins biosynthesis. Meanwhile, gene tadV encoding a prepilin peptidase was up-regulated 1.5-fold in the stationary phase. Prepilin peptidase catalyzes the formation of mature pilin from prepilin proteins . In addition, tadZ and tadABC encoding four Flp pilus assembly proteins involved in pilus assembly and localization were also up-regulated 10.4, 7.0, 12.4 and 10.6-fold, respectively. A gene cluster - “tad2003” locus (Bbr_0132~0138 and Bbr_0901) encoding Tad pili has been reported in Bifidobacterium breve UCC2003 (GenBank accession number NC_020517.1), which are responsible for colonization in murine gut . The amino acid sequence alignment of Tad pili between B. breve UCC2003 and B. lactis A6 showed an amino acid sequence identity among 34.2%-59.3% (Table S5 and Fig. 5). Tad pili is an important adhesive structure identified in bifidobacteria , and the adhesion of B. lactis A6 to mucin showed that the adhesion rate was 4.90% for stationary cells, significantly higher than that of 2.38% for exponential cells (Fig. 6). These findings further confirmed the adhesion to mucin was expressed in a growth phase-dependent manner in B. lactis A6 strain, and the enhanced expression of Tad pili may facilitate its colonization and maintenance in the intestine of host receptor.
Probiotics are widely accepted to be beneficial for the maintenance of the gut homeostasis, improvement of the immune system, and amelioration of various metabolic disease [30, 31]. In recent years, the consumption of probiotics to promote health has grown rapidly worldwide and become an independent industry . Notably, the physiology and application characteristics of probiotics are closely associated with the growth phase. Many studies have focused on the overall growth phase-associated changes during fermentation based on transcriptome and/or proteome technology among L. casei, L. rhamnosus, and B. longum [9, 10, 33]. For the species of B. lactis, transcriptomic techniques have been applied to reveal the oligosaccharides uptake and metabolism in B. lactis Bl-04 , the xylo-oligosaccharides utilization in B. lactis BB12 , and the tetracycline response in B. lactis strains Bl-04 and HN019 . However, to the best of our knowledge, no transcriptomic study about the growth phase-associated changes has been reported in this species. Andersen and colleagues identified the genetic loci involved in the uptake and catabolism of 11 potential prebiotics oligosaccharides based on transcriptional analysis in B. lactis Bl-04. The results showed that B. lactis Bl-04 exhibits extensive capabilities to utilize indigestible carbohydrates in the human gastrointestinal tract. For cellobiose metabolism, genes Balac_1569~1570 encoding two permease components and Balac_1572 encoding oligosaccharide-binding protein of an ABC transporter responsible for cellobiose uptake were found significantly up-regulated when cellobiose was used as carbon source compared with glucose. However, in the present study, we observed two genes BAA6_1415 encoding cellobiose ABC transporter permease protein CebG and BAA6_1416 encoding cellobiose phosphorylase CbpA were increased 3.6 and 7.7-fold in the stationary phase compared with the log phase, respectively. B. lactis A6 tends to increase gene expression for utilization of exogeneous cellobiose with glucose consumption in the stationary growth phase. The strain of B. lactis A6 was isolated from the feces of a Bama centenarian, whose diet was rich in plant-based fibers . Plant derived carbohydrates, such as cellulose, could be catalyzed to produce soluble cellobiose by cellobiohydrolases in the human gut . Therefore, the up-regulation of cellobiose transport and utilization may contribute to the adaptation of B. lactis A6 in the intestine of healthy centenarians.
Peptidoglycan is the major component of Gram-positive bacterial cell wall and is essential for the maintenance of bacterial integrity and shape . In this study, most of genes involved in peptidoglycan biosynthesis and cell division were down-regulated in the stationary phase compared with the log phase, which suggested that the bacterial growth was generally repressed during this stage. Laakso et al. has been reported the growth phase-related changes of L. rhamnosus GG when grew in industrial-type whey medium, and the results showed that two genes murC and murD involved in peptidoglycan biosynthesis were significantly down-regulated both at the transcriptomic and proteomic level in the stationary growth phase compared to the mid-log phase . This result is consistent with what has been observed in our study. Interestingly, genes involved in oligopeptides uptake and dipeptides metabolism were found up-regulated in the stationary phase. In contrast, genes involved in de novo synthesis of amino acid, such as aspartate, methionine, tyrosine, phenylalanine, arginine, proline, and glutamine were highly down-regulated. It has been reported that lactic acid bacteria could lower the synthesis of DNA and proteins and shift physiological state to the stationary phase due to carbohydrate starvation conditions during fermentation . At this stage, an increased catabolism of carbon sources like peptides and free amino acids was observed for energy production . Therefore, the up-regulated oligopeptides uptake and dipeptides metabolism may produce more free amino acids for energy production in bacterial cells of B. lactis A6 and then elaborate a potential survival strategy to adapt to glucose consumption during growth.
For acid tolerance response, B. lactis A6 exhibited various strategies to regulate the pH reduction during growth. NH3 could neutralize intracellular excess H+ to respond acid tolerance in bacteria. Three metabolic pathways have been reported with NH3 production: cysteine-cystathionine-cycle (CCC) pathway and deamination of glutamate and BCAAs. In B. longum BBMN68, three genes of malY1, metC3, and cysK involved in CCC pathway were up-regulated 9.25-, 2.71-, and 2.84-fold after acid treatment, respectively . Glutamate dehydrogenase (GDH) catalyzes the reversible oxidative deamination of glutamate to α-ketoglutarate and NH3. Under low pH conditions, this reaction towards to produce more NH3 . In addition, the overexpressed BCAAs biosynthesis has been reported in an acid-pH-resistant mutant B. longum 8809dpH after growth at an acidic condition. A higher intracellular NH4+ concentration was observed to respond acid treatment, which are positively associated with the deamination of BCAAs . In this study, genes involved in the CCC pathway (cysK and metC), glutamate dehydrogenase (gdhA) and BCAAs biosynthesis (ilvB and ilvC) were all up-regulated over 2.6-fold in the stationary phase compared to the log phase, which may contribute to NH3 production for acid tolerance response in B. lactis A6. Moreover, Clp protease has been reported for participating general stress responses. Clp proteases could recognizes and degrades irreversibly damaged proteins in response to acid stress . The up-regulated Clp proteases for degradation of damaged proteins could facilitate amino acids recycling in B. lactis A6 under acid stress. Furthermore, the general acid tolerance between the log phase and stationary phase cells has been verified by in vitro experiment, the results showed that the stationary cells are more acid-tolerant than the log cells in B. lactis A6.
Colonization in the mammalian gastrointestinal tract by microbes is believed to play an essential part to confer health benefits . A Tad pili encoding gene cluster - “tad2003” locus was identified in B. breve UCC2003, which was proved to be essential for colonization and persistence in murine gut as a conserved host-colonization factor . The “tad2003” locus encodes eight proteins for the biosynthesis of a mature Tad pilus. TadE and TadF are participated in the prepilin precursor biosynthesis, and TadA is an ATPase to provide energy for Flp-pilus assembly. TadB and TadC are thought to channel the energy of ATP hydrolysis into Flp-pilus polymerization or serve as an inner membrane scaffold for the assembly of the pilus biogenesis apparatus. Then the Tad apparatus is directed to the cell poles by TadZ [27, 44]. In this study, a putative Tad pili encoding gene cluster BAA6_0205 ~ 0212 was highly expressed during the stationary phase, which exhibited similar gene cluster structure and higher amino acid identity of 34.2%-59.3% with “tad2003” locus. In addition, Tad pili could bind to carbohydrate moieties present in glycoproteins or glycolipids receptors of mucus layer, which covered on the surface of intestinal epithelial cells [45, 46]. Thus, the enhanced expression of Tad pili in B. lactis A6 strain during the stationary phase may facilitate its colonization and maintenance in the intestine of the host. To verify the hypothesis, the adhesion rate of B. lactis A6 to mucin was performed in vitro with the log and stationary phase bacterial cells, respectively. A significantly increased adhesion rate to mucin from 2.38% to 4.90% was observed during the transition from the exponential phase to the stationary phase.
Probiotic bacteria harvested at different phases of growth has been shown to cause different application characteristics and mucosal responses in human. In this study, growth phase-associated changes were investigated in B. lactis A6 based on RNA-seq analysis. The bacterium expanded the carbon source utilizing profile to transport cellobiose, XOS and raffinose in the stationary phase and tended to metabolize L-arabinose and xylose with reduced glycolysis in the stationary phase. Meanwhile, genes involved in acid tolerance response were upregulated, and in vitro experiment proved that the stationary phase cells are more acid-tolerant than the log phase cells. Remarkably, a putative gene cluster encoding Tad pili was significantly up-regulated at the stationary phase, accompanied with an enhanced bacterial adhesion to mucin at this stage. These findings provide a novel insight into the growth phase associated characteristics in B. lactis A6 and provide valuable information for further application of B. lactis A6 in the food industry.
Bacterial strains and culture conditions
B. lactis A6 was isolated from the feces of a centenarian in Bama, Guangxi, China and collected from the China General Microbiological Culture Collection Center (CGMCC) with CGMCC No. 9273. A -80°C glycerol stock of B. lactis A6 was anaerobically (100% N2) incubated at 37°C in MRSc broth for 24 h, followed by 1% (v/v) inoculation to a 15 mL MRSc for overnight. The basal MRSc broth contained the following sources (per liter of deionized water): 10 g of beef extract, 5 g of yeast extract, 10 g of tryptone peptone, 20 g of glucose, 2 g of triamine citrate, 2 g of K2HPO4, 0.5 g of MgSO4, 0.25 g of MnSO4, 5 g of sodium acetate, 1 mL of Tween-80 and 0.5 g of L-cysteine. The pH and optical density at 600 nm (OD600) were recorded at 2-h intervals.
Overnight culture of B. lactis A6 strain was anaerobically (100% N2) inoculated (1% v/v) into a 15 mL fresh MRSc broth and incubated at 37°C. Bacterial cells were harvested in the exponential phase (OD600 = 0.6-0.8) and the stationary phase (OD600 = 4.0-4.2) by centrifugation at 4°C, 8000 × g for 5 min. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA) and purified using the TURBO DNA-free Kit (Ambion, Austin, TX, USA) according to the manufacturer’s instructions. RNA concentration was measured using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA).
Sequencing libraries were generated using NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB, USA). Index codes were added to attribute sequences to each sample. The library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA) to select cDNA fragments of 150~200 bp in length. Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37°C for 15 min followed by 95°C for 5 min. The PCR was performed using Phusion High-Fidelity DNA polymerase (NEB) with Universal PCR primers and Index (X) primer. Finally, PCR products were purified by AMPure XP system (Beckman Coulter) and the library quality was assessed on the Agilent Bioanalyzer 2100 system (Agilent Technologies). The library preparations were sequenced on an Illumina Hiseq platform and paired-end reads were generated.
Transcriptomic data processing
Raw data of fastq format were firstly processed through in-house perl scripts. Clean reads were obtained by removing reads containing adapter, reads containing ploy-N or reads with low quality from raw data. All clean reads were aligned to the genome of B. lactis A6 with GenBank Accession No. NZ_CP010433.1 by Bowtie2-2.2.3 . HTSeq v0.6.1 was then used to count the numbers of reads mapped to each gene and calculate the expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced (FPKM) value . Differential expression analysis of samples at the stationary growth phase compared with exponential phase was performed using the DESeq R package (V1.18.0) . Genes with an adjusted P-value < 0.05 and |Fold change| ≥ 1.5 were assigned as differentially expressed. The metabolic pathways involved in these DEGs were further analyzed by KEGG enrichment analysis .
Acid tolerance test
The acid tolerance of B. lactis A6 was performed according to Waddington et al. with some modifications . Briefly, the bacteria were anaerobically (100% N2) incubated at 37°C in MRSc broth. Cells were harvested at the log growth phase (OD600 = 0.6-0.8) and the stationary phase (OD600 = 4.0-4.2). After centrifugation at 8000 × g, 22°C for 5 min, cell pellets were washed twice with 0.85% NaCl solution and adjusted the initial bacterial numbers to ~108 CFU/mL. Aliquots of 1 mL cell suspension were collected and centrifugated at 8000 × g, 22°C for 5 min. The cell pellets were then suspended in 1 mL fresh MRSc broth at pH values of 4.0, 3.0, or 2.0 (adjusted by lactic acid). After anaerobic incubation at 37°C for 2 h, the bacterial number (CFU/mL) was determined by plating 10-fold serial dilutions on MRSc agar and incubating anaerobically at 37°C for 40 h. The bacterial cells untreated with acid were included as a control. The survival rate was calculated as follows: Survival rate (%) = [(CFU/mL) acid treatment / (CFU/mL) control] × 100. The results were obtained by three independent experiments.
Bacterial adhesion to mucin
The amino acid sequence alignment of Tad pili between B. breve UCC2003 and B. lactis A6 was performed with SnapGene software (from Insightful Science, version 5.3) in a local alignment (Smith-Waterman) model. The adherence of bacteria to mucin was performed as described by Xiong et al. . Briefly, Nunc MaxiSorp 96-well microplate (Thermo Fisher Scientific, Denmark) was coated with 100 μL of mucin (from porcine stomach; Sigma-Aldrich, USA) and incubated for 16 h at 4°C at a concentration of 2.5 pmol per well. The wells were then washed twice with PBS and incubated with 2% (w/v) bovine serum albumin for 2 h at 37°C. The wells were subsequently washed three times with PBS. The bacteria collected at the exponential and stationary phase were washed twice with PBS and resuspended to a final OD600 equivalent to ~2-3 × 108 CFU/mL. Then 200 μL of bacterial suspension were added to coated 96-well plates and incubated at 37°C for 1 h. The unattached bacteria were removed by washing the wells for 3 times with PBS, and then treated with 200 μL 0.05% (v/v) Triton X-100 and incubated at 37°C for 30 min. The adhesion ratios (%) were calculated by comparing the bacterial counts after adhesion to the number of cells originally added to the plate wells. The results were obtained by three independent experiments.
Data were analyzed using GraphPad Prism 8 software for Windows (GraphPad Software, Inc., La Jolla, CA, USA). Results were presented as the mean value ± standard deviation. An unpaired Student t-test was used to calculate P values when two groups were compared.
Availability of data and materials
All the raw data of transcriptome sequencing could be downloaded from GEO with the Accession No. GSE173957.
- B. lactis A6:
Bifidobacterium animalis subsp. lactis A6
- L. rhamnosus GG:
Lactobacillus rhamnosus GG
- L. casei Zhang:
Lactobacillus casei Zhang
- B. breve UCC2003:
Bifidobacterium breve UCC2003
de Man-Rogosa-Sharp broth with L-cysteine
Differential expression genes
Kyoto Encyclopedia of Genes and Genomes
Branched-chain amino acids
- Tad pili:
Type IVb tight adherence pili
China General Microbiological Culture Collection Center
Expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced
Specifically expressed genes
Lactic acid bacteria
Hill C, Guarner F, Reid G, Gibson GR, Merenstein DJ, Pot B, et al. The International Scientific Association for Probiotics and Prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat. Rev. Gastroenterol Hepatol. 2014;11(8):506–14.
Bron PA, Kleerebezem M, Brummer RJ, Cani PD, Mercenier A, Macdonald TT, et al. Can probiotics modulate human disease by impacting intestinal barrier function? Br. J. Nutr. 2017;117(01):93–107.
Mörkl S, Butler MI, Holl A, Cyran JF, Dinan TG. Probiotics and the microbiota-gut-brain axis: focus on psychiatry. Curr Nutr Rep. 2020;9(8):171–82.
Sanders ME, Merenstein DJ, Reid G, Gibson GR, Rastall RA. Probiotics and prebiotics in intestinal health and disease: from biology to the clinic. Nat. Rev. Gastroenterol. Hepatol. 2019;16(Suppl. 1):605–16.
Rivera-Espinoza Y, Gallardo-Navarro Y. Non-dairy probiotic products. Food Microbiol. 2010;27(1):1–11.
Broeckx G, Kiekens S, Jokicevic K, Byl E, Henkens T, Vandenheuvel D, et al. Effects of initial cell concentration, growth phase, and process parameters on the viability of Lactobacillus rhamnosus GG after spray drying. Dry Technol. 2020;38(11):1474–92.
Deepika G, Green RJ, Frazier RA, Charalampopoulos D. Effect of growth time on the surface and adhesion properties of Lactobacillus rhamnosus GG. J Appl Microbiol. 2009;107(4):1230–40.
Turroni F, Foroni E, Montanini B, Viappiani A, Strati F, Duranti S, et al. Global genome transcription profiling of Bifidobacterium bifidum PRL2010 under in vitro conditions and identification of reference genes for quantitative real-time PCR. Appl Environ Microbiol. 2011;77(24):8578–87.
Wang J, Zhang W, Zhi Z, Wei A, Bao Q, Yong Z, et al. Gene expression profile of probiotic Lactobacillus casei Zhang during the late stage of milk fermentation. Food Control. 2012;25(1):321–7.
Laakso K, Koskenniemi K, Koponen J, Kankainen M, Surakka A, Salusjärvi T, et al. Growth phase-associated changes in the proteome and transcriptome of Lactobacillus rhamnosus GG in industrial-type whey medium. Microb Biotechnol. 2011;4(6):746–66.
Zhao L, Xu W, Ibrahim SA, Jin J, Feng J, Jiang J, et al. Effects of age and region on fecal microflora in elderly subjects living in Bama, Guangxi. China. Curr Microbiol. 2011;62(1):64–70.
Sun E, Zhao L, Ren F, Liu S, Zhang M, Guo HY. Complete genome sequence of Bifidobacterium animalis subsp. lactis A6, a probiotic strain with high acid resistance ability. J Biotechnol. 2015;200:8–9.
Huo Y, Lu X, Wang X, Wang X, Chen L, Guo H, et al. Bifidobacterium animalis subsp. lactis A6 alleviates obesity associated with promoting mitochondrial biogenesis and function of adipose tissue in mice. Molecules. 2020;25(7):1490.
MäDer U, Nicolas P, Richard H, Bessières P, Aymerich S. Comprehensive identification and quantification of microbial transcriptomes by genome-wide unbiased methods. Curr. Opin. Biotechnol. 2011;22(1):32–41.
Andersen JM, Barrangou R, Abou Hachem M, Lahtinen SJ, Goh YJ, Sevensson B, et al. Transcriptional analysis of oligosaccharide utilization by Bifidobacterium lactis Bl-04. BMC Genomics. 2013;14:312.
Glenn K, Smith KS. Allosteric regulation of Lactobacillus plantarum xylulose 5-phosphate/fructose 6-phosphate phosphoketolase (Xfp). J. Bacteriol. 2015;197(7):1157–63.
Fushinobu S. Unique sugar metabolic pathways of bifidobacteria. Biosci. Biotechnol. Biochem. 2010;74(12):2374–84.
Pokusaeva K, Fitzgerald GF, van Sinderen D. Carbohydrate metabolism in Bifidobacteria. Genes Nutr. 2011;6(3):285–306.
Higgins CF. ABC transporters: from microorganisms to man. Annu. Rev. Cell Biol. 1992;8:67–113.
Abouhamad WN, Manson M, Gibson MM, Higgins CF. Peptide transport and chemotaxis in Escherichia coli and Salmonella typhimurium: characterization of the dipeptide permease (Dpp) and the dipeptide-binding protein. Mol Microbiol. 2010;5(5):1035–47.
An H, Douillard FP, Wang G, Zhai Z, Yang J, Song S, et al. Integrated transcriptomic and proteomic analysis of the bile stress response in a centenarian-originated probiotic Bifidobacterium longum BBMN68. Mol. Cell Proteomics. 2014;13(10):2558–72.
Mengin-Lecreulx D, Heijenoort J. Copurification of glucosamine-1-phosphate acetyltransferase and N-acetylglucosamine-1-phosphate uridyltransferase activities of Escherichia coli: characterization of the glmU gene product as a bifunctional enzyme catalyzing two subsequent steps in the pathway for UDP-N-acetylglucosamine synthesis. J. Bacteriol. 1994;176(18):5788–95.
Sung MT, Lai YT, Huang CY, Chou LY, Shih HW, Cheng WC, et al. Crystal structure of the membrane-bound bifunctional transglycosylase PBP1b from Escherichia coli. Proc. Natl. Acad. Sci. USA. 2009;106(22):8824–9.
Sham LT, Butler EK, Lebar M, Kahne D, Bernhardt TG, Ruiz N. MurJ is the flippase of lipid-linked precursors for peptidoglycan biogenesis. Science. 2014;345(6193):220.
Jin J, Bing Z, Guo H, Cui J, Jiang L, Song S, et al. Mechanism analysis of acid tolerance response of Bifidobacterium longum subsp. longum BBMN 68 by gene expression profile using RNA-Sequencing. PLoS One. 2012;7(12):e50777.
Huo Y, Zhan Y, Wang Q, Li S, Yang S, Nomura CT, et al. Acetolactate synthase (AlsS) in Bacillus licheniformis WX-02: enzymatic properties and efficient functions for acetoin/butanediol and l-valine biosynthesis. Bioproc. Biosyst. Eng. 2018;41(1):87–96.
Tomich M, Planet PJ, Figurshi DH. The tad locus: postcards from the widespread colonization island. Nat. Rev. Microbiol. 2007;5(5):363–75.
O'Connell Motherway M, Zomer A, Leahy SC, Sinead C, Reunanen J, Bottacini F, et al. Functional genome analysis of Bifidobacterium breve UCC2003 reveals type IVb tight adherence (Tad) pili as an essential and conserved host-colonization factor. Proc. Natl. Acad. Sci. USA. 2011;108(27):11217–22.
Grimm V, Westermann C, Riedel CU. Bifidobacteria-host interactions—an update on colonisation factors. Biomed Res Int. 2014;960826.
Zhao Y, Zeng Y, Zeng D, Wang H, Zhou M, Sun N, et al. Probiotics and microRNA: their roles in the host-microbe interactions. Front Microbiol. 2021;604462.
Varsha KK, Maheshwari AP, Nampoothiri KM. Accomplishment of probiotics in human health pertaining to immunoregulation and disease control. Clin Nutr ESPEN. 2021;44:26–37.
Lu K, Dong SW, Wu XY, Jin RM, Chen HB. Probiotics in Cancer. Front. Oncol. 2021;11:638148.
Mozzetti V, Grattepanche F, Moine D, Berger B, Rezzonico E, Arigoni F, et al. Transcriptome analysis and physiology of Bifidobacterium longum NCC2705 cells under continuous culture conditions. Benef. Microbes. 2012;3(4):261–72.
Gilad O, Jacobsen S, Stuer-Lauridsen B, Pedersen MB, Garrigues C, Svensson B. Combined transcriptome and proteome analysis of Bifidobacterium animalis subsp. lactis BB-12 grown on xylo-oligosaccharides and a model of their utilization. Appl. Environ. Microbiol. 2010;76(21):7285–91.
Morovic W, Roos P, Zabel B, Hidalgo-Cantabrana C, Kiefer A, Barrangou R. Transcriptional and functional analysis of Bifidobacterium animalis subsp. lactis exposure to tetracycline. Appl. Environ. Microbiol. 2018;84(23):e01999–18.
Wang F, Yu T, Huang GH, Cai D, Liang XL, Sun HY, et al. Gut microbiota community and Its assembly associated with age and diet in Chinese centenarians. J Microbiol Biotechn. 2015;25(8):1195–204.
Cantarel BL, Lombard V, Henrissat B. Complex carbohydrate utilization by the healthy human microbiome. PLoS One. 2012;7(6):e28742.
Regulski K, Courtin P, Meyrand M, Claes IJ, Lebeer S, Vanderleyden J, et al. Analysis of the peptidoglycan hydrolase complement of Lactobacillus casei and characterization of the major γ-D-glutamyl-L-lysyl- endopeptidase. PLoS One. 2012;7(2):e32301.
Papadimitriou K, Alegría A, Bron PA, de Angelis M, Gobbetti M, Kleerebezem M, et al. Stress physiology of lactic acid bacteria. Microbiol. Mol. Biol. Rev. 2016;80(3):837–90.
Ganesan B, Stuart MR, Weimer BC. Carbohydrate starvation causes a metabolically active but nonculturable state in Lactococcus lactis. Appl. Environ. Microbiol. 2007;73(8):2498–512.
Smith HQ, Li C, Stanley CA, Smith TJ. Glutamate dehydrogenase, a complex enzyme at a crucial metabolic branch point. Neurochem. Res. 2019;44(1):117–32.
Sánchez B, Champomier-Vergès MC, Anglade P, Baraige F, Reyes-Gavilan CGD, Margolles A, et al. Proteomic analysis of global changes in protein expression during bile salt exposure of Bifidobacterium longum NCIMB 8809. J Bacteriol. 2005;187(16):5799–808.
Shanahan F. The colonic microflora and probiotic therapy in health and disease. Curr. Opin. Gastroenterol. 2011;27(1):61–5.
Xu Q, Christen B, Chiu HJ, Jaroszewski L, Klock HE, Knuth MW, et al. Structure of the pilus assembly protein TadZ from Eubacterium rectale: implications for polar localization. Mol. Microbiol. 2012;83(4):712–27.
Ventura M, Turroni F, Motherway MO, Macsharry J, Sinderen D. Host-microbe interactions that facilitate gut colonization by commensal bifidobacteria. Trends Microbiol. 2012;20(10):467–76.
Paone P, Cani PD. Mucus barrier, mucins and gut microbiota: The expected slimy partners? Gut. 2020;69(12):2232–43.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat. Methods. 2012;9(4):357–9.
Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;2:166–9.
Anders S. Analysing RNA-Seq data with the DESeq package. Mol. Biol. 2010;43:1–17.
Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93.
Waddington L, Cyr T, Hefford M, Hansen LT, Kalmokoff M. Understanding the acid tolerance response of bifidobacteria. J Appl. Microbiol. 2010;108(4):1408–20.
Xiong Y, Zhai Z, Lei Y, Xiao B, Hao Y. A novel major pilin subunit protein FimM is involved in adhesion of Bifidobacterium longum BBMN68 to intestinal epithelial cells. Front. Microbiol. 2020;11:590435.
This work was supported by the Key Laboratory of Functional Dairy, Co-constructed by Ministry of Education and Beijing Municipality, and College of Food Science and Nutritional Engineering, China Agricultural University, Beijing, China.
This work was supported by the National Key Research and Development Program of China (2018YFC1604303).
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, H., An, J., Fan, C. et al. Transcriptome analysis revealed growth phase-associated changes of a centenarian-originated probiotic Bifidobacterium animalis subsp. lactis A6. BMC Microbiol 22, 61 (2022). https://doi.org/10.1186/s12866-022-02474-5