Original Research ARTICLE
Intra-species Genomic and Physiological Variability Impact Stress Resistance in Strains of Probiotic Potential
- 1Division of Gastroenterology and Hepatology, Department of Medicine, Microbiome Core Facility, Center for Gastrointestinal Biology and Disease, School of Medicine, University of North Carolina, Chapel Hill, NC, United States
- 2Department of Chemistry, College of Arts and Sciences, University of North Carolina, Chapel Hill, NC, United States
- 3Research Computing, University of North Carolina, Chapel Hill, NC, United States
Large-scale microbiome studies have established that most of the diversity contained in the gastrointestinal tract is represented at the strain level; however, exhaustive genomic and physiological characterization of human isolates is still lacking. With increased use of probiotics as interventions for gastrointestinal disorders, genomic and functional characterization of novel microorganisms becomes essential. In this study, we explored the impact of strain-level genomic variability on bacterial physiology of two novel human Lactobacillus rhamnosus strains (AMC143 and AMC010) of probiotic potential in relation to stress resistance. The strains showed differences with known probiotic strains (L. rhamnosus GG, Lc705, and HN001) at the genomic level, including nucleotide polymorphisms, mutations in non-coding regulatory regions, and rearrangements of genomic architecture. Transcriptomics analysis revealed that gene expression profiles differed between strains when exposed to simulated gastrointestinal stresses, suggesting the presence of unique regulatory systems in each strain. In vitro physiological assays to test resistance to conditions mimicking the gut environment (acid, alkali, and bile stress) showed that growth of L. rhamnosus AMC143 was inhibited upon exposure to alkaline pH, while AMC010 and control strain LGG were unaffected. AMC143 also showed a significant survival advantage compared to the other strains upon bile exposure. Reverse transcription qPCR targeting the bile salt hydrolase gene (bsh) revealed that AMC143 expressed bsh poorly (a consequence of a deletion in the bsh promoter and truncation of bsh gene in AMC143), while AMC010 had significantly higher expression levels than AMC143 or LGG. Insertional inactivation of the bsh gene in AMC010 suggested that bsh could be detrimental to bacterial survival during bile stress. Together, these findings show that coupling of classical microbiology with functional genomics methods for the characterization of bacterial strains is critical for the development of novel probiotics, as variability between strains can dramatically alter bacterial physiology and functionality.
The gastrointestinal tract is host to one of the densest and most diverse microbial communities on the planet (Gill et al., 2006; Li et al., 2012). The composition and function of the gut microbiota is critical to maintenance of host gastrointestinal health (Jones, 2016). In fact, diseases including diabetes (Tilg and Moschen, 2014), obesity (Carding et al., 2015), and colorectal cancer (Sobhani et al., 2011; Borges-Canha et al., 2015) have been shown to have gut dysbioses associated with progression. Conversely, modulation of the gut microbiota can alleviate or even eliminate disorders like Clostridium difficile infections (Petrof et al., 2013; Allegretti et al., 2014), inflammatory bowel disease (D'Haens et al., 2014), and lactose intolerance (Azcarate-Peril et al., 2017).
Intra-species genetic polymorphisms constitute the majority of the diversity within the microbiota of the human gut (Greenblum et al., 2015; Zhang and Zhao, 2016). Coupling classical microbiology approaches with next generation sequencing provides an opportunity to study physiological characteristics of individual microbial strains, and to identify unique genomic elements associated with those phenotypes. Moreover, advances in cultivation technologies have improved isolation of novel microorganisms from human subjects, provided the ability to study physiology of difficult to grow microbes (Faith et al., 2010), and ultimately develop advanced microbiota-derived treatments for gastrointestinal diseases (Forster and Lawley, 2015).
Use of probiotics, live microbes that convey a benefit to their host when administered in adequate amounts, as interventions for gastrointestinal disorders has gained increasing support (Grover et al., 2012; Vandenplas et al., 2015); however, intra-species variations can impact their functionality as effective probiotics (Chapman et al., 2012). Isolation and characterization of novel intestinal organisms are important steps toward the development of new and improved probiotics. Moreover, high throughput sequencing technology allow for cost-effective whole genome sequencing of bacterial isolates, providing a wealth of genomic data. Nevertheless, identification and characterization of novel probiotics also require careful examination of microbial physiology (Papadimitriou et al., 2015).
One essential physiological characteristic of probiotics is their ability to survive the environmental conditions of the gastrointestinal tract. Gut microorganisms have evolved highly conserved mechanisms for tolerance to gastrointestinal stresses, which include variations in pH and the antimicrobial effects of bile (Azcarate-Peril et al., 2004; Bruno-Bárcena et al., 2004). These mechanisms are controlled by well-regulated genetic systems (Azcarate-Peril et al., 2004), often with rapid response times to acute stresses, allowing for immediate resistance to environmental changes (Anderson et al., 2010; Bove et al., 2013). Bacterial responses to environmental stress occur in a highly regulated manner. First, stress-specific genes aim to aleviate the immediate stress (maintenance of cytoplasmic pH homeostasis, transport/degradation of toxic compounds, and others), followed by universal stress response systems aimed to repair DNA, protein, membrane, and cell wall damage.
If exposed to an acidic environment, H+ pump of ATPases (F0F1-ATPases) maintain pH homeostasis within the cytoplasm by pumping excess H+ ions back into the environment (Wang et al., 2013; Cusumano and Caparon, 2015). Other mechanisms involved in resistance to acid stress in Lactobacillus are arginine deiminases (ArcABC) (Fulde et al., 2014) that synthesize ammonia from arginine and free H+ within the cytoplasm (Casiano-Colón and Marquis, 1988) and amino acid decarboxylases coupled with amino acid/biogenic amines antiporter systems (Azcarate-Peril et al., 2004; Papadimitriou et al., 2016). Upon surviving the acidic conditions in the stomach and passing into the duodenum, the pH of the environment shifts from highly acidic to alkaline. Alkaline Shock Proteins (Asp) are encoded by lactic acid bacteria including Lactobacillus rhamnosus (Arnold et al., 2017), and have been implicated in resistance to sudden pH shifts (Kuroda et al., 1995; Seetharaaman, 2008; Anderson et al., 2010).
Entry into the duodenum not only involves a change in pH, but also exposes microorganisms to bile salts, which act as detergents, causing cell damage and cytotoxicity (Andreichin, 1980; Begley et al., 2005). Lactic acid bacteria including Lactobacillus and Bifidobacterium species encode bile salt hydrolases (bsh), which have been shown in some cases to provide resistance to bile toxicity (Grill et al., 2000; Patel et al., 2010; Lin et al., 2014). Lactobacillus acidophilus and Bifidobacterium encode multiple bsh genes, which exhibit substrate specificity for different conjugated bile salt targets (McAuliffe et al., 2005; Jarocki and Targonski, 2013). As different bile salts have varying toxicity, individual bsh genes may or may not provide survival advantages against bile (Fang et al., 2009). In addition to providing survival advantages in the gastrointestinal environment, bile salt hydrolases encoded by lactic acid bacteria (Jarocki and Targonski, 2013; Jarocki et al., 2014; Pithva et al., 2014) have been implicated in modulation of host cholesterol and lipid metabolism (Patel et al., 2010; Joyce et al., 2014), and as a mechanism for microbe-host signaling (Ridlon et al., 2006).
Universal bacterial stress-response systems encode mechanisms to cope with the deleterious consequences of exposure to environmental stresses, including DNA damage, protein misfolding, and loss of cell wall/membrane integrity. Bacterial DNA damage is repaired by MutS and RecA-like recombinases (Lee and Pi, 2010; Rossi et al., 2016; Calderini et al., 2017; Overbeck et al., 2017). Protein misfolding is reduced and repaired by molecular chaperones including heat-shock proteins (Ruiz et al., 2013; Calderini et al., 2017). Cell wall integrity of gram positive bacteria is maintained through lipoteichoic acid (LTA), exopolysaccharide, and fatty acid metabolism/synthesis genes (Koskenniemi et al., 2011). Membrane proteins including metalloproteases provide additional stress resistance in lactic acid bacteria (Bove et al., 2012). Lactobacillus also encodes two component regulatory systems (2CRS) that promote rapid responses to environmental stresses. These systems usually consist of a histidine kinase and a response regulator gene that work to sense and react to changes in the environment (Yu et al., 2014; Monedero et al., 2017), facilitating stress resistance (Morel-Deville et al., 1998; Alcántara et al., 2011). Together, universal stress-response genes work with specific stress resistance systems to ensure bacterial survival to gastrointestinal conditions.
Douillard et al. (2013) in their analysis of 100 L. rhamnosus genomes, showed that strains had unique genomic elements coupled with physiological traits that were influenced by the strain origin. In particular, bile resistance was higher in L. rhamnosus isolates derived from the gastrointestinal tract as opposed to dairy, oral, or vaginal isolates. In this study we analyzed genomic, transcriptomic, and physiological characteristics of two novel human L. rhamnosus strains of probiotic potential for which our group recently generated genomic sequence information (Thompson et al., 2015; Arnold et al., 2017), and compared them with an established probiotic strain of human origin (LGG). Lactobacillus rhamnosus GG (LGG) is a well-characterized probiotic lactic acid bacteria isolated from a healthy human host in 1983. The strain provides a number of benefits to its host including immunomodulation (Lebeer et al., 2012; Segers and Lebeer, 2014), pathogen exclusion (De Keersmaecker et al., 2006; Makras et al., 2006), and modulation of host gene expression (Kankainen et al., 2009; Yan et al., 2013; Segers and Lebeer, 2014). LGG sets a precedent for L. rhamnosus as a functional probiotic as well as for the isolation and identification of novel probiotics from healthy human hosts. Additionally, the strains were compared with L. rhamnosus of dairy origin (Lc705, HN001; Ceapa et al., 2015) to demonstrate the impact of intra-species genetic polymorphisms on tolerance to gastrointestinal stress.
Materials and Methods
Strains and Culture Media
Bacterial strains used in this study are presented in Table 1. Strains were propagated in MRS broth (Pronadisa, Madrid) at 37°C without agitation, or on MRS agar plates containing 1.5% agar. For growth assays, MRS supplemented with 0.1–1.0% Oxgall was used to test resistance to bile. MRS was titrated with HCl–pH 4 to test growth in acidic conditions, and titrated with NaOH to pH 8 to test growth under alkaline stress.
Bacterial Growth Assays
Lactobacillus rhamnosus strains were grown statically for 16 h in MRS broth at 37°C. Freshly harvested cells were washed with MRS without glucose and diluted 1:100 in MRS (pH 6.6) containing either bile (0, 0.1, 0.3, 0.5, and 1.0% w/v oxgall), or adjusted to pH 4 (HCl) or pH 8 (NaOH). 200 μl of bacterial suspensions were placed into each well of a 96 well plate, sealed and placed into Tecan Infinite 200 Pro spectrophotometer (Tecan, Switzerland), where they were grown for 24 h at 37°C. Measurements of optical density at 600 nm (OD600nm) were taken every 15 min during this 24-h period. Maximum specific growth rates (μmax) were then calculated for each treatment type and plotted in Origin2016 (OriginLab, Northampton, MA).
Bacterial Survival Assays
Bacterial cells grown to mid-log growth phase (OD600nm = ~0.6) were centrifuged, washed twice with sterile PBS and diluted 1:10,000 in simulated gastric juice (0.5% pepsin, 0.5% NaCl, pH3), simulated intestinal juice (0.3% pancreatin, 0.5% NaCl, pH8), or bile (0.5% Oxgall in MRS + 0% Glucose). Cell suspensions were incubated for 2 h at 37°C and then plated onto MRS agar plates by WASP Spiral Plater (Don Whitley Scientific, West Yorkshire, UK) at 1 and 2 h post inoculation. Plates were incubated at 37°C for 48 h. Survival was measured using ProtoCOL colony counter (Synbiosis, Frederick, MD) to count colonies present on MRS plates of treated and untreated cells. Three biological replicates with three technical replicates were included in survival experiments. Pairwise comparisons using a two-tailed Student's t-test were performed to determine statistically significant differences between treatments.
AMC010, AMC143, and LGG were grown to early log phase (OD600nm = ~0.4) in MRS prior to harvesting. Cells were then washed twice with sterile PBS and re-suspended in simulated gastric juice, simulated intestinal juice, or bile (0.5% w/v Oxgall in sterile PBS) and incubated for 10 min at 37°C. Three biological replicates were included in expression experiments. After incubation, cells were centrifuged, flash frozen in RNAlater (Thermo Fisher Scientific, Waltham, MA), and stored in −80°C until RNA isolation was performed. RNA isolation was performed using the Qiagen RNeasy PowerMicrobiome kit (Qiagen, Valencia, CA) as directed. RNA was eluted in 50 μl RNase free water and quantified using the 2200TapeStation (Agilent Technologies, Santa Clara, CA).
Real Time Quantitative PCR
Five ng of total RNA isolated from bacterial cultures exposed to bile (0.5% w/v Oxgall in sterile PBS for 30 min at 37°C) were reverse transcribed using qScript cDNA SuperMix (Quanta BioSciences, Gaithersburg, MD). To generate a standard curve, genomic DNA from AMC010 was amplified using primer set BSHqPCR-F (TTGGCGCTGACGACTTGC), BSHqPCR-R (AATCTTGACGCCTTGACC) (this study) and purified via gel extraction in 1.5% agarose gel. The purified PCR product was serially diluted and used in RT–qPCR experiments. The qPCR master mix included 10 μl of Power SYBR® Green 2x PCR Master Mix (Applied Biosystems, Foster City, CA), 2 μl of each primer (BSHqPCR-F, BSHqPCR-R) (1 μM stock), 1 μl PCR grade water, and 5 μl of template cDNA (1 ng/μl). The reaction was run for 40 cycles of melting (95°C) 15 s, annealing/extension (65°C) 45 s followed by SYBR detection was carried out on the 7,500 Fast Real Time PCR System (Applied Biosystems, Foster City, CA) after 10 min denaturation step at 95°C. Samples and standards were run in triplicate.
Ribosomal RNA was depleted from total RNA using the Ribo-Zero Gold Bacterial rRNA Removal Reagent (Epidemiology Kit) (Illumina, San Diego, CA) according to manufacturer's instructions. Briefly, the rRNA-specific magnetic beads were removed from storage buffer and mixed with 500 ng of total sample RNA. Subsequently, rRNA removal solution was added and samples were incubated for 10 min at 65°C. Finally, samples were placed on magnetic stand for 15 min at 22°C and mRNA was removed and immediately processed with TruSeq Stranded mRNA HT kit (Illumina, San Diego, CA) according to manufacturer's instructions. Briefly, RNA was mixed with Fragment-Prime mix and incubated at 94°C for 8 min. Samples were immediately subject to first strand and second strand cDNA synthesis reactions, respectively, followed by 3′ end repair, adenylation and adapter ligation. After adapter ligation, the libraries were enriched by PCR using the following thermal cycling conditions: 98°C for 30 s followed by 15 cycles of 98°C for 10 s, 60°C for 30 s and 72°C for 30 s. Final extension step of 70°C for 5 min was carried out following the last cycle. After enrichment, libraries were purified with Beckman Coulter magnetic beads (Brea, CA), washed with 80% ethanol and eluted in Tris pH 8.5. Following enrichment, cDNA was barcoded for multiplexing via PCR, using dual-index barcodes [index 1(i7) and index 2(i5)] (Illumina, San Diego, CA) in a combination unique to each sample. Final ds cDNA was purified with Beckman Coulter magnetic beads (Brea, CA). Library concentrations and quality were measured via TapeStation2200 (Agilent Technologies, Santa Clara, CA). Barcoded libraries were pooled at equimolar concentrations and sequenced on the Illumina HiSeq platform (Illumina, San Diego, CA).
mRNA Sequencing Data Analysis
Sequencing output from the Illumina HiSeq platform was converted to FASTQ format and demultiplexed using Illumina BclFastq 22.214.171.124 (Illumina, San Diego, CA). Quality control was performed via FastQC on both raw and processed sequencing reads (Babrahm Institute, Cambridge, UK). Demultiplexed FASTQ sequence files were uploaded to Geneious software (Biomatters, New Zealand). Sequencing reads from each treatment were mapped against LGG genome in Geneious software using “Geneious for RNA Seq” mapper, a minimum mapping quality of 30 (99.9% confidence), allowing gaps with a maximum of 10% per read, 20% maximum mismatches per read, with word length of 20 bases. Expression levels were calculated in Geneious and compared between treatment types. Genes identified as significantly differentially regulated between treatments (p ≤ 0.05) were further filtered to include only genes with 2-fold expression differences. Genes identified as significantly differentially regulated at ≥2-fold in at least one treatment type for each isolate were plotted as heat maps in OriginLab software (Origin Lab, Northampton, MA), and compared between strains to show differences.
Preparation of Electrocompetent Cells of L. rhamnosus
Electrocompetent L. rhamnosus cells were generated using a combination of previously described methods (Kim et al., 2005; Welker et al., 2015), with minor adjustments. Briefly, an overnight culture of L. rhamnosus AMC010 was diluted to 106 cells/ml in pre-warmed MRS broth supplemented with 2% glycine and was incubated overnight at 37°C without agitation. After incubation, 5 ml of culture was inoculated in 100 ml of freshly prepared, pre-warmed MRS medium supplemented with 2% glycine. The culture was incubated at 37°C without agitation until it reached and OD600nm of 0.2. Ampicillin was then added to the culture at a final concentration of 10 μg/ml, and incubation at 37°C continued until the culture reached an OD600nm of 0.4. Cells were then harvested by centrifugation at room temperature (10 min at 6,000 g), washed twice at room temperature with electroporation buffer (0.5 M sucrose, 7 mM potassium phosphate pH7.4, 1 mM MgCl2), resuspended in 1 ml of the same buffer, and placed on ice. The electrocompetent cells were used immediately for electroporation.
Generation of the bsh Insertional Mutant Strain
Insertional inactivation of the bsh gene in L. rhamnosus AMC010 was done as described (Welker et al., 2015). Briefly, primers BshFHindIII (TATTAAGCTTTTCAGACAGAGGCGGCTTTGC) and BshREcoRI (AATAGAATTCAATCTTGACGCCTTGACCAC) were designed to amplify a 652 bp region of the bsh gene in AMC010. The primers included EcoR1 and HindIII restriction sites for cloning into the pFAJ-5301 vector (Lebeer et al., 2012). The resulting vector (pFAJ-BSHi) was used to transform AMC010 by electroporation using an Bio-Rad Gene Pulser (peak voltage, 1.7 kV; capacitance, 25 μF; resistance, 200 Ω with time constant of 2–4 ms). Electroporated AMC010 cells were allowed to recover for 24 h in MRS at 37°C without agitation, and subsequently plated on MRS containing erythromycin at a final concentration of 2 μg/ml. Individual erythromycin resistant colonies were selected after 48 h of growth at 37°C and sub-cultured in MRS containing 2 μg/ml erythromycin overnight without agitation at 37°C. Disruption of the bsh genes was verified by PCR amplification using primers Bsh1753F (ATTGCCTGACCTAGATGCAGG) and Bsh1753R (AACACCGGCGACAGGTCCATC). Genomic integration of the erythromycin resistance cassette was also verified by PCR amplification with the primers Ery625F (CTACTTAATCTGATAAGTGAGC) and Ery625R (TCAGCACAGTTCATTATCAACC). AMC010::bsh mutants were cultivated in MRS broth containing 2 μg/ml erythromycin and stored at −80°C in 15% glycerol.
Comparative Analysis of the L. rhamnosus Stress Genomic Complement
The genomes of L. rhamnosus AMC 143 and AMC010 isolated from infant stools (Thompson et al., 2015; Arnold et al., 2017) were aligned to the genome of the well-characterized probiotic L. rhamnosus GG (Segers and Lebeer, 2014) for comparison (Figure 1A). Our study focused on stress genes and systems encoded by AMC143 and AMC010 in comparison with other L. rhamnosus strains of dairy and human origin.
Figure 1. Comparative genomic analysis of Lactobacillus rhamnosus isolates. (A) Genome alignment of L. rhamnosus strains. Strains AMC43, AMC010, Lc705, and HN001 were compared to LGG using BRIG genome alignment software. Relevant stress-response genes are annotated. Comparison of relevant stress-response genes in L. rhamnosus, (B) ABC-F ATPase, (C) arcD, (D) bsh were aligned and compared in Geneious 10.1.3 software. The bottom figure in each graph represents single nucleotide polymorphisms within the compared gene.
When exposed to environmental pH fluctuations, proton translocating F0F1-ATPases are the primary genes responsible for maintaining cytoplasmic pH homeostasis (Azcarate-Peril et al., 2004). Each strain in this study encoded ATPase genes with high nucleotide identity to one another (>98%), however the genomic architecture surrounding the ABC-F Family ATPase varied between strains (Figure 1B). In addition to ATPases, cytoplasmic pH homeostasis can be maintained by amino acid decarboxylases and their corresponding antiporters (Azcarate-Peril et al., 2004). Both AMC010 and AMC143 encoded the arginine/ornithine antiporter arcD gene; however, the 12 kb region that included this unique arcD gene was absent in LGG and HN001 (Figure 1C). This region encompassed arcD as well as an NADH oxidase, two diguanylate cyclases, cellulose synthase, glycosyl transferase, glycosyl hydrolase and a hypothetical protein. An ornithine decarboxylase was identified in all strains with over 97% nucleotide identity with no architectural deviations between strains. Another amino acid decarboxylase and antiporter system, the glutamate decarboxylase/glutamate gamma-aminobutyrate antiporter (gadC) has been shown to reduce the impact of acid stress in Lactobacillus (Azcarate-Peril et al., 2004); however, the gadC gene, as well as the corresponding decarboxylase gene were absent in both AMC010 and AMC143. Finally, an uncharacterized amino acid decarboxylase gene was present in LGG (LGG_RS12751) with 97.2% nucleotide identity to its homolog in AMC010, but this gene was absent in AMC143.
Similarly to acid exposure, cytoplasmic pH homeostasis is critical to maintain when cells are exposed to alkaline stress. Sodium-proton antiporters and potassium-proton antiporters maintain cytoplasmic pH in alkaline extracellular environments (Lee et al., 2011; Nyanga-Koumou et al., 2012). The strains in this study encoded 3 sodium-proton antiporters (98.2, 98, and 97.1% nucleotide identity). LGG encoded a unique sodium-proton antiporter absent in AMC010 and AMC143. Each strain also encoded a single potassium transporter, with 97.3% nucleotide identity between strains. Alkaline shock proteins (Asp) have been implicated in providing protection from damage caused by increases in extracellular pH change (Kuroda et al., 1995; Seetharaaman, 2008; Anderson et al., 2010). Each of the strains in this study encoded 4 asp genes including asp23, each with minimum of 98.3% nucleotide identity between strains.
The anti-microbial/detergent properties of bile salts require cells to respond upon exposure in order to resist cytotoxic effects. Each strain in this study encoded a single copy of a bile salt hydrolase (bsh). This gene has been implicated in bile-stress tolerance in other lactobacilli (McAuliffe et al., 2005). The sequence identity of the bsh gene exceeded 98.2% between strains, however AMC143 contained a deletion of 208 base pairs upstream of the bsh gene, eliminating a putative promoter, ribosome binding site, and inducing a 38 bp truncation of the bsh gene itself (Figures 1D, 5B). Aside from the deletion in AMC143, the bsh upstream non-coding region has 93% nucleotide identity between strains. Each strain also encoded cell wall synthesis/repair pathways including DltABCD (Koskenniemi et al., 2011).
In addition to specific stress-response mechanisms, universal stress response systems were present in AMC143 and AMC010. Genes associated with DNA repair were highly conserved, including radA, recN, recO, recU, radC, recA, and mutS, each of which had over 96.5% nucleotide identity between strains. Additionally, each strain encoded a suite of highly conserved molecular chaperones, which mitigate stress-associated protein misfolding. These genes included dnaK, groEL, groES, HSP20, HSP33, dnaI, and clpB, each exhibiting 98.2% nucleotide identity or higher between strains. Finally, bacteria membrane integrity is often compromised as a result of environmental stress, allowing metal ions otherwise excluded from the cytoplasm to permeate through the membrane, causing further DNA damage and cell death (Pratviel, 2012; Qiao and Ma, 2013). One mechanism that bacteria have to survive declines in membrane integrity is to actively pump cytotoxic metal ions out of their cytoplasm utilizing metal translocating ATPases (Chien et al., 2013). Each strain encoded three different non-specific metal transporting ATPases each with 97.6% nucleotide identity between strains, as well as specific ATPases for transport of copper (97.0% nucleotide identity between strains) and magnesium (97.5% nucleotide identity between strains).
Gene Expression Analysis under Stress Conditions
Differential transcription profiles were generated for AMC010 and AMC143 from mRNA sequencing data obtained from cells exposed to acid (pH 3 adjusted with HCl), alkaline (pH8 adjusted with NaOH), or bile (0.5% w/v oxgall in PBS) stress, and compared to untreated cells. A total 216 genes were identified as differentially regulated (>2-fold change, p < 0.05) for AMC010 in response to all treatments, while 79 genes were identified in AMC143. Figure 2 shows a heatmap indicating genes that were up or down regulated in each strain. We mapped our strains to L. rhamnosus GG to provide gene homolog references. Our data confirmed differential regulation of universal stress-response genes and specific stress response genes, although most specific genes were differentially regulated below either the statistical or fold change cutoffs (Supplementary Table 1).
Figure 2. Transcriptomics analysis. mRNA sequencing data from each treatment (pH3, pH8, 0.5%, w/v Bile) was mapped to LGG genome and compared to untreated cells to calculate Log2 differential expression ratio and p-value. The 50 genes with the highest differential expression ratio (FDR corrected p ≤ 0.05) from each strain were categorized by function (cell division, DNA/protein repair, transcription regulation, and molecular transport) and plotted in OriginLab software. Highly differentially expressed hypothetical genes and genes with functions that did not fit into a defined category were also included. Green text highlights genes differentially regulated in both AMC010 and AMC143.
Exposure to pH 3 increased expression of a copper-translocating P-type ATPase homologous to LGG_RS08675 by approximately 4-fold in both strains. F0F1-ATPases were not differentially regulated between strains, however a metal transporting ATPase homologous to LGG_RS00670 was up-regulated approximately 1.5-fold in both AMC010 and AMC143 upon acid exposure. Amino acid decarboxylases were not differentially regulated between strains. All treatments induced expression of a gene encoding a small heat shock protein Hsp20 homologous to LGG_RS13420 in both AMC143 and AMC010. Exposure of AMC010 to acid induced expression of an additional uncharacterized Hsp20-like molecular chaperone (LGG_RS03170, 11.3-fold), the padR transcription regulator (LGG_RS03330, 7.5-fold), and two hypothetical proteins (uncharacterized, Lactobacillus-specific LGG_RS11075, 9.2-fold, and fliK-like LGG_RS13395, 6-fold). Conversely, a gene cluster including 5 genes involved in purine biosynthesis (LGG_RS08715, LGG_RS08720, LGG_RS08725, LGG_RS08730, and LGG_RS08735) was inhibited in AMC010 across all conditions.
Exposure to alkaline conditions induced expression of the alkaline-shock protein Asp23/Gls24 family envelope stress response protein homologous to LGG_RS01130 in AMC010 by 2.3-fold. This gene also showed a non-statistically significant 2-fold induction in AMC143. Other alkaline shock proteins homologous to LGG_RS01125 and LGG_RS08100 were induced by at least 2-fold in both strains, while the putative alkaline shock protein homologous to LGG_RS01130 was repressed by approximately 2-fold under the same condition, though p-values for each of these genes fell below the significance cutoffs set for this study. Alkaline stress also induced the lacI transcription regulator (homologous to LGG_RS01785) in both AMC010 (6.1-fold) and AMC143 (12.1-fold), the Hsp20-like molecular chaperone (LGG_RS03170, 19.7-fold in AMC143), which was also induced by exposure to acid in AMC010, and uncharacterized hypothetical proteins homologous to LGG_RS00635 (7-fold) and LGG_RS01080 (5.7-fold) in AMC010. Additionally, exposure to pH 8 resulted in decreased expression of the translation initiation factor IF-3 (LGG_RS08400) and a ribosomal-processing cysteine protease (Prp) homologous to LGG_RS08120, reducing expression levels by 2.6-fold and 4.3-fold respectively in AMC143. Alkaline challenge of AMC010 resulted in reduced expression of the septation inhibitor protein similar to divIC (Bennett et al., 2007) (LGG_RS12050) by 3.7-fold, and minC/ftsL, genes associated with inhibition of cell division (LGG_RS06095 by 2-fold and LGG_RS06140 by −3.5-fold) ABC transport ATP binding protein (LGG_RS09570) by 4.3-fold and a YbjQ_1 domain-containing hypothetical protein homologous to LGG_RS02955 by 42.2-fold. The galactose-6-phosphate isomerase operon (LGG_RS03135, LGG_RS03140, LGG_RS03145, and LGG_RS03150) was specifically up regulated in AMC143 exposed to pH 8, while the Rpml/translation initiation factor gene cluster (LGG_RS08400 and LGG_RS08395) was down regulated in the same strain by approximately 3-fold.
Exposure to bile resulted in the differential expression of 52 genes in AMC143 and 111 genes in AMC010 with considerable overlap with alkaline treatment. Although the bsh gene has been shown to be induced by exposure to bile in Lactobacillus (Koskenniemi et al., 2011), in our study, bsh (similar to LGG_RS02395) showed a non-statistically significant (p = 0.2) 1.2-fold induction in AMC010 while AMC143 showed a marginal repression. The genes more impacted by bile exposure in AMC143 were a HU-family transcription regulator homologous to LGG_RS06660 (induced 2-fold), a septation inhibitor protein (LGG_RS12050), which was down regulated 9-fold, the glutamine ABC transport system operon composed of genes homologus to LGG_RS13875, LGG_RS13880, and LGG_RS13885, which were down regulated 2-fold, and hypothetical proteins homologous to LGG_RS01070, LGG_RS09500, and LGG_RS00415, down regulated 2.8-, 4-, and 13.9-fold respectively.
Exposure of AMC010 to Oxgall resulted in a 2.8-fold induction of a 2CRS response regulator (RR), dcuR (homologous to LGG_RS13770), similar to the RR involved in regulation of the anaerobic fumarate respiratory system in E. coli (Janausch et al., 2004), as well as a 5.3-fold induction of a major facilitator superfamily (MFS) transporter permease (LGG_RS00320) and 4.3-fold induction of the Hsp20-like molecular chaperone (LGG_RS03170). Conversely, a septation inhibitor protein (LGG_RS12050), prsA similar to LGG_RS10900, a molecular chaperone potentially involved in a late stage of protein export (Jakob et al., 2015; Jousselin et al., 2015), the NrdH-redoxin (LGG_RS07055), uracil transport (LGG_RS06995, LGG_RS07000, and LGG_RS07005) and peptide ABC transport system operon (LGG_RS07940, LGG_RS07945, LGG_RS07950, LGG_RS07955, and LGG_RS07960 were down regulated in AMC010 exposed to bile.
In Vitro Growth and Survival of L. rhamnosus Exposed to Simulated Gastrointestinal Conditions
Under normal growth conditions (MRS broth, pH 6.6 at 37°C) LGG showed the highest growth rate among strains (0.63 h−1), while AMC143 had the lowest growth rate (0.48 h−1) with no significant differences between AMC010, Lc705, and HN001 (Figure 3A). While growth rates between strains were similar under normal conditions, there were distinct differences in growth phenotypes between strains when exposed to simulated gastrointestinal conditions. Growth rates of all strains were reduced approximately 5-fold at pH 4 (adjusted with HCl); however the strains isolated from dairy products, HN001 and Lc705, showed further reduced rates than the human-derived strains, AMC010, AMC143, and LGG (Figure 3A). As expected, survival of early-log cells exposed to simulated gastric juice (pepsin, NaCl, HCl pH3) was minimally impacted. No decrease in survival was observed for LGG, Lc705, or AMC143 after 1 h of exposure, while an approximately 15% decrease was observed for HN001 and AMC010 strains. After 2 h, Lc705 showed close to 100% survival, LGG and HN001 were reduced by approximately 20%, and AMC010 and AMC143 populations were reduced by 40% (Figure 3B).
Figure 3. Growth and survival of L. rhamnosus under pH-stress conditions. (A) Growth rates were calculated for cells grown at pH 6.6 (MRS) pH4 (adjusted with HCl), or pH8 (adjusted with NaOH). (*p < 0.01). (B,C) Bacterial survival was assayed by plating and enumerating cells exposed to simulated gastric juice (NaCl, HCl, Pepsin pH3) or simulated intestinal juice (NaCl, HCl, Pancreatin pH8) for 1 or 2 h.
Exposure of L. rhamnosus strains to MRS adjusted to pH 8 with NaOH showed nearly a 50% decline in growth rate only in AMC143. HN001 and Lc705 showed increased growth rates when grown in media at higher pH 8, while AMC010 and LGG showed no significant differences in growth rates compared to pH 6.6 (Figure 3A). When early-log cells were exposed to simulated intestinal juice (pancreatin, NaCl, NaOH pH8) (Charteris et al., 1998) cell numbers increased after 1 h to then decrease in all strains, except HN001 where survival was stable (Figure 3C).
Finally, we sought to evaluate growth and survival of L. rhamnosus strains in the presence of bile. At sub-physiological concentrations of bile (0.1% w/v), the only strain that showed growth inhibition was Lc705 (60% growth rate reduction). As the concentration of bile in the growth medium increased to physiological levels (0.3% w/v), growth rates decreased by approximately 50% for AMC010, AMC143, and LGG, while HN001 and Lc705 were further inhibited (58 and 80% respectively) (Figure 4A). Lc705 showed significant decreased growth rates in all bile concentrations tested. As bile concentrations exceeded physiological levels (0.5% bile) growth rates decreased by over 60%. The growth inhibition observed at 1.0% w/v bile was greater for HN001 and Lc705 (83–92%) than for the strains of human origin (71–78%).
Figure 4. Growth and survival of L. rhamnosus under bile-stress conditions. (A) Growth rates relative to 0% bile were calculated for cells grown in MRS containing 0.1, 0.3 0.5, and 1.0% w/v oxgall. (*p > 0.05, **p > 0.01). (B) Bacterial survival was determined by plating and enumerating colonies post exposure to 0.5% w/v oxgall in MRS devoid of any carbohydrate source for 2 h. *Indicates significant difference between AMC143 and each other strain tested. (*p > 0.01, ND = No detected growth).
Exposure of L. rhamnosus to 0.5% bile resulted in a significant decrease in survival in all strains except AMC143. AMC143 exhibited nearly 20-fold higher survival rates compared to all other strains after 1 h of exposure to 0.5% bile. After 2 h, no HN001 colonies were detected, and the LGG and AMC010 populations were reduced by >95%. AMC143 however remained resistant to the antimicrobial effects of bile, with over 15% of cells surviving after 2 h of exposure (Figure 4B).
The Bile Salt Hydrolase (bsh) Gene from L. rhamnosus AMC010 and AMC143
Despite a highly conserved chromosomal architecture and gene sequence, the region upstream of the bsh gene in AMC143 contained a 208 bp deletion, eliminating a putative promoter region and ribosome binding site that is highly conserved in the other strains in the study (Figure 5B). Additionally, this deletion resulted in a truncation of the bsh gene. AMC143 showed significantly increased resistance to bile. To investigate if the deletion was responsible for the increased resistance of the strain to bile, we first assessed bsh expression levels in L. rhamnosus AMC143, AMC010, and LGG by RT-qPCR. Data showed variable expression of the gene in absence of bile. Upon bile exposure the bsh gene in AMC010 was induced approximately 4-fold, while expression of bsh homologs in AMC143 and LGG did not vary. Moreover, bsh expression was significantly lower in the presence or absence of bile in AMC143 (Figure 5A).
Figure 5. Bile salt hydrolase (bsh) expression and sequence analysis. (A) RTqPCR targeting bsh was performed on AMC010, AMC143, and LGG with and without acute bile treatment, (*p > 0.05). (B) Sequence alignment of bsh performed in CLUSTAL including annotations of coding regions and ribosome binding site.
To determine if inactivation of the bsh gene in AMC010 would replicate the bile resistant phenotype observed in AMC143, bsh was disrupted by insertion of the pFAJ-BSHi vector. Growth of the AMC010::bsh mutant under normal growth conditions was indistinguishable from the wild type strain. Survival and growth rates of the mutant were assayed under each stress condition previously performed in this study and compared to the wild type strain. When exposed to acid stress conditions, growth and survival of AMC010::bsh was identical to AMC010. Similarly, growth and survival at high pH was identical between the two strains (Supplementary Figure 1A). Growth rates of AMC010::bsh and AMC010 in varying concentrations of bile (0–1.0% w/v Oxgall) were similar in the presence of glucose (Supplementary Figure 1B); however, when exposed to high concentrations of bile (0.5% w/v Oxgall) without a carbohydrate source, survival of AMC010::bsh was enhanced 4–6-fold compared to the wild type strain (Figure 6), further suggesting a link between bsh expression and bile resistance in L. rhamnosus.
Figure 6. Survival of AMC010::bsh under bile-stress. Survival of AMC010::bsh was compared to wild type AMC010 by plating and enumerating colonies post exposure to 0.5% w/v oxgall (*p < 0.01).
The future of personalized medicine will employ defined bacterial consortia to treat or prevent diseases according to host genotype, diet, and life stage, and the essential first step to develop such consortia is the extensive characterization of gastrointestinal organisms (Faith et al., 2010).
High-throughput in silico and in vitro analyses of gene-level variation across a large array of species in the human gut are useful as screening tools; however, detailed characterization of strains is essential to determine functional implications of strain variation in the complex ecosystem of the gut (Greenblum et al., 2015). Our study confirms that strain-level variability contribute significantly to microbial diversity within complex microbial communities (Zhang and Zhao, 2016). We compared genomic and physiological data of two novel strains of L. rhamnosus (AMC010 and AMC143) (Thompson et al., 2015; Arnold et al., 2017) to three characterized strains (LGG, Lc705, and HN001) to determine their response to simulated gastrointestinal stress.
One of the most striking physiological differences between the strains in our study was the growth inhibition of AMC143 exposed to alkaline conditions. Comparative genomic analysis revealed no major differences between genes previously identified as alkaline-stress response genes, with nucleotide identities ranging from 98.5 to 99.6% between strains. Despite this growth defect, the strain had comparable survival rates when exposed to simulated intestinal juice, suggesting that the growth defect is caused by arrest of cell division as opposed to cell death. Upon exposure to alkaline stress, expression levels of genes associated with growth inhibition (divIC, ftsL, and minC) were repressed in AMC010, while expression of the same genes in AMC143 was unaffected. Conversely, expression of the translation initiation factor IF-3 was downregulated in AMC143 when exposed to alkaline conditions. We could speculate that this differential expression plays a role in the growth inhibition observed for AMC143. A previous study in Lactobacillus plantarum showed that expression of two endopeptidases (Accession numbers Q52071 and Q048X8) was reduced under alkaline conditions (Lee et al., 2011); however we failed to identify homologus genes significantly downregulated in AMC143. Further study of cellular response to alkaline pH is required to better understand the reasons for a growth defect in AMC143.
Analysis of genome sequencing data revealed a 208 bp deletion upstream of the bsh gene in AMC143 (Figure 1D). Sequence comparison between strains identified a repeat of 10 base pairs flanking the deleted region, suggesting that this deletion may have occurred in AMC143 through homologous recombination. Physiological assays showed that AMC143 was more resistant to bile-induced toxicity than the other strains tested, which led us to investigate the impact of bsh expression on bile toxicity. RTqPCR expression data showed a significantly lower bsh transcript count in AMC143 compared to LGG or AMC010, both in absence or presence of bile. A non significant induction of the bsh gene in AMC010 was also observed by mRNA sequencing data. Bile salt hydrolases have been associated to bile stress resistance in Lactobacillus and Bifidobacterium (Grill et al., 2000; McAuliffe et al., 2005; Lin et al., 2014) and implicated as a mechanism for host-microbe signaling (Zhou and Hylemon, 2014; Song et al., 2015; McMillin et al., 2016). These enzymes hydrolyze conjugated bile salts generating unconjugated bile acids, which often function as a signaling molecules to host cells and have strong anti-microbial effects (Sytnik et al., 1978; Grill et al., 2000; Schmidt et al., 2001; Kong et al., 2011), and an amino acid, which can be utilized by both host and microorganisms for protein synthesis (Ridlon et al., 2006, 2014; Patel et al., 2010). Studies have shown that BSHs do not always provide a survival advantage to bacteria exposed to bile (Fang et al., 2009), suggesting that BSH activity differs between microorganisms. In fact, lactobacilli encode variable bsh genes, each with unique substrates and activities (McAuliffe et al., 2005; Ren et al., 2011). To determine if bsh expression was correlated to lower survival upon bile exposure in our strains, we generated a mutant strain of AMC010 containing a disrupted bsh gene by site directed insertion (Lebeer et al., 2012). AMC010::bsh recapitulated the physiological phenotype observed in AMC143. These findings suggest that the unconjugated products of BSH activity were cytotoxic to the strains used in this study. Our survival assays were performed in MRS broth devoided of a carbohydrate source. In accordance to a previous report (Ziar et al., 2014), when survival experiments were performed in the presence of a carbohydrate source, bile cytotoxity was almost completely abolished for all strains (data not shown).
Strains of the same species derived from different environments are likely to have evolved different abilities to tolerate environmental stress (Douillard et al., 2013). Our data shows distinct physiological differences between strains isolated from human hosts (LGG, AMC010, and AMC143) and strains isolated from fermented dairy products (Lc705, HN001). When grown in alkaline conditions, dairy-derived strains exhibited accelerated growth, which correlates with studies done with other dairy-derived lactobacillus strains (Mojgani et al., 2015), while intestinal strains were either inhibited (AMC143) or unaffected (LGG, AMC010). Additionally, we found that while strains isolated from human hosts were able to grow with limited inhibition at sub-physiological levels of bile, the dairy isolate Lc705 was significantly inhibited even at very low bile concentrations (0.1%w/v). Moreover, survival of dairy-derived strains exposed to high bile concentrations was lower than intestinal strains, suggesting that the strain environmental origin may have driven evolution of their stress response pathways. Strains evolving in environments devoid of bile have no selective pressure to retain bile resistance genes, and these molecular mechanisms are lost over time. This phenomenon has been observed in dairy derived samples as well as in the oral microbiota (Douillard et al., 2013; de Barros et al., 2016).
Analysis of 16S rRNA sequencing data is unable to accurately resolve taxonomy at the species/sub-species level (Janda and Abbott, 2007) as single gene sequencing is insufficient to differentiate between similar species, especially in highly diverse and complex communities. Even with the best of what next generation sequencing has to offer, classical microbiology approaches are absolutely critical in understanding how microbial behavior and physiology correlate to genomic, proteomic, and transcriptomic data. As novel microbiota-derived interventions like new probiotics are being developed, it is important to keep in mind the magnitude to which seemingly minor genomic variability can result in changes to cellular physiology and thus probiotics efficacy.
JA Designed and performed experiments, and wrote the manuscript; JS and JK Performed experiments; JR Performed bioinformatics analysis; MA Designed the experiments, advised JA, contributed to bioinformatic and statistical analyses of data and wrote the manuscript; All authors read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Dr. Sarah Lebeer from the Department of Bioscience Engineering at University of Antwerp for generously supplying us with the pFAJ-5301 suicide vector for use in this study. We would also like to thank Dr. Rodolphe Barrangou from the Department of Food, Bioprocessing and Nutrition Sciences at North Carolina State University for his help with optimization of electroporation protocols. The Microbiome Core is supported in part by the NIH/National Institute of Diabetes and Digestive and Kidney Diseases grant P30 DK34987.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00242/full#supplementary-material
Supplementary Figure 1. (A) Growth rates of AMC010::bsh in MRS at pH 4, 6.6, and 8. (B) Relative growth rate of AMC010::bsh in bile (0, 0.1, 0.3, 0.5, and 1.0% w/v oxgall).
Supplementary Table 1. Expression fold change of Lactobacillus rhamnosus genes under stress conditions.
Alcántara, C., Revilla-Guarinos, A., and Zúñiga, M. (2011). Influence of two-component signal transduction systems of Lactobacillus casei BL23 on tolerance to stress conditions. Appl. Environ. Microbiol. 77, 1516–1519. doi: 10.1128/AEM.02176-10
Anderson, K. L., Roux, C. M., Olson, M. W., Luong, T. T., Lee, C. Y., Olson, R., et al. (2010). Characterizing the effects of inorganic acid and Alkaline Shock on the Staphylococcus aureus transcriptome and messenger RNA turnover. FEMS Immunol. Med. Microbiol. 60, 208–250. doi: 10.1111/j.1574-695X.2010.00736.x
Arnold, J. W., Monteagudo-Mera, A., Altermann, E., Cadenas, M. B., Thompson, A. L., and Azcarate-Peril, M. A. (2017). Genome sequences of potential probiotic Lactobacillus rhamnosus isolates from human infants. Genome Announc. 5:e00107-17. doi: 10.1128/genomeA.00107-17
Azcarate-Peril, M. A., Altermann, E., Hoover-Fitzula, R. L., Cano, R. J., and Klaenhammer, T. R. (2004). Identification and inactivation of genetic loci involved with Lactobacillus acidophilus acid tolerance. Appl. Environ. Microbiol. 70, 5315–5322. doi: 10.1128/AEM.70.9.5315-5322.2004
Azcarate-Peril, M. A., Ritter, A. J., Savaiano, D., Monteagudo-Mera, A., Anderson, C., Magness, S. T., et al. (2017). Impact of short-chain galactooligosaccharides on the gut microbiome of lactose-intolerant individuals. Proc. Natl. Acad. Sci. U.S.A. 114, E367–e375. doi: 10.1073/pnas.1606722113
Bennett, J. A., Aimino, R. M., and McCormick, J. R. (2007). Streptomyces coelicolor genes ftsL and divIC play a role in cell division but are dispensable for colony formation. J. Bacteriol. 189, 8982–8992. doi: 10.1128/JB.01303-07
Borges-Canha, M., Portela-Cidade, J. P., Dinis-Ribeiro, M., Leite-Moreira, A., and Pimentel-Nunes, F. (2015). Role of colonic microbiota in colorectal carcinogenesis: a systematic review. Rev. Esp. Enferm. Dig. 107, 659-671. doi: 10.17235/reed.2015.3830/2015
Bove, P., Capozzi, V., Garofalo, C., Rieu, A., Spano, G., and Fiocco, D. (2012). Inactivation of the ftsH gene of Lactobacillus plantarum WCFS1: effects on growth, stress tolerance, cell surface properties and biofilm formation. Microbiol. Res. 167, 187–193. doi: 10.1016/j.micres.2011.07.001
Bove, P., Russo, P., Capozzi, V., Gallone, A., Spano, G., and Fiocco, D. (2013). Lactobacillus plantarum passage through an oro-gastro-intestinal tract simulator: carrier matrix effect and transcriptional analysis of genes associated to stress and probiosis. Microbiol. Res. 168, 351–359. doi: 10.1016/j.micres.2013.01.004
Bruno-Bárcena, J. M., Andrus, J. M., Libby, S. L., Klaenhammer, T. R., and Hassan, H. M. (2004). Expression of a heterologous manganese superoxide dismutase gene in intestinal lactobacilli provides protection against hydrogen peroxide toxicity. Appl. Environ. Microbiol. 70, 4702–4710. doi: 10.1128/AEM.70.8.4702-4710.2004
Calderini, E., Celebioglu, H. U., Villarroel, J., Jacobsen, S., Svensson, B., and Pessione, E. (2017). Comparative proteomics of oxidative stress response of Lactobacillus acidophilus NCFM reveals effects on DNA repair and cysteine de novo synthesis. Proteomics 17:1600178. doi: 10.1002/pmic.201600178
Ceapa, C., Lambert, J., van Limpt, K., Wels, M., Smokvina, T., Knol, J., et al. (2015). Phenotype profiling links genotype to carbohydrate utilization signatures in the species Lactobacillus rhamnosus. Appl. Environ. Microbiol. 81, 5458–5470. doi: 10.1128/AEM.00851-15
Chapman, C. M., Gibson, G. R., and Rowland, I. (2012). In vitro evaluation of single- and multi-strain probiotics: inter-species inhibition between probiotic strains, and inhibition of pathogens. Anaerobe 18, 405–413. doi: 10.1016/j.anaerobe.2012.05.004
Charteris, W. P., Kelly, P. M., Morelli, L., and Collins, J. K. (1998). Development and application of an in vitro methodology to determine the transit tolerance of potentially probiotic Lactobacillus and Bifidobacterium species in the upper human gastrointestinal tract. J. Appl. Microbiol. 84, 759–768. doi: 10.1046/j.1365-2672.1998.00407.x
Chien, C. C., Huang, C. H., and Lin, Y. W. (2013). Characterization of a heavy metal translocating P-type ATPase gene from an environmental heavy metal resistance Enterobacter sp. isolate. Appl. Biochem. Biotechnol. 169, 1837–1846. doi: 10.1007/s12010-012-0047-4
Cusumano, Z. T., and Caparon, M. G. (2015). Citrulline protects Streptococcus pyogenes from acid stress using the arginine deiminase pathway and the F1Fo-ATPase. J. Bacteriol. 197, 1288–1296. doi: 10.1128/JB.02517-14
de Barros, J. M., Costabile, A., Charalampopoulos, D., Khutoryanskiy, V. V., and Edwards, A. D. (2016). Evaluating and optimizing oral formulations of live bacterial vaccines using a gastro-small intestine model. Eur. J. Pharm. Biopharm. 102, 115–122. doi: 10.1016/j.ejpb.2016.03.010
De Keersmaecker, S. C., Verhoeven, T. L., Desair, J., Marchal, K., Vanderleyden, J., and Nagy, I. (2006). Strong antimicrobial activity of Lactobacillus rhamnosus GG against Salmonella typhimurium is due to accumulation of lactic acid. FEMS Microbiol. Lett. 259, 89–96. doi: 10.1111/j.1574-6968.2006.00250.x
D'Haens, G. R., Sartor, R. B., Silverberg, M. S., Petersson, J., and Rutgeerts, P. (2014). Future directions in inflammatory bowel disease management. J. Crohns Colitis 8, 726–734. doi: 10.1016/j.crohns.2014.02.025
Douillard, F. P., Ribbera, A., Kant, R., Pietilä, T. E., Järvinen, H. M., Messing, M., et al. (2013). Comparative genomic and functional analysis of 100 Lactobacillus rhamnosus strains and their comparison with strain GG. PLoS Genet. 9:e1003683. doi: 10.1371/journal.pgen.1003683
Faith, J. J., Rey, F. E., O'Donnell, D., Karlsson, M., McNulty, N. P., Kallstrom, G., et al. (2010). Creating and characterizing communities of human gut microbes in gnotobiotic mice. ISME J. 4, 1094–1098. doi: 10.1038/ismej.2010.110
Fang, F., Li, Y., Bumann, M., Raftis, E. J., Casey, P. G., Cooney, J. C., et al. (2009). Allelic variation of bile salt hydrolase genes in Lactobacillus salivarius does not determine bile resistance levels. J. Bacteriol. 191, 5743–5757. doi: 10.1128/JB.00506-09
Fulde, M., Willenborg, J., Huber, C., Hitzmann, A., Willms, D., Seitz, M., et al. (2014). The arginine-ornithine antiporter ArcD contributes to biological fitness of Streptococcus suis. Front. Cell. Infect. Microbiol. 4:107. doi: 10.3389/fcimb.2014.00107
Gill, S. R., Pop, M., Deboy, R. T., Eckburg, P. B., Turnbaugh, P. J., Samuel, B. S., et al. (2006). Metagenomic analysis of the human distal gut microbiome. Science 312, 1355–1359. doi: 10.1126/science.1124234
Grill, J. P., Perrin, S., and Schneider, F. (2000). Bile salt toxicity to some bifidobacteria strains: role of conjugated bile salt hydrolase and pH. Can. J. Microbiol. 46, 878–884. doi: 10.1139/w00-066
Jakob, R. P., Koch, J. R., Burmann, B. M., Schmidpeter, P. A., Hunkeler, M., Hiller, S., et al. (2015). Dimeric structure of the bacterial extracellular foldase PrsA. J. Biol. Chem. 290, 3278–3292. doi: 10.1074/jbc.M114.622910
Janausch, I. G., Garcia-Moreno, I., Lehnen, D., Zeuner, Y., and Unden, G. (2004). Phosphorylation and DNA binding of the regulator DcuR of the fumarate-responsive two-component system DcuSR of Escherichia coli. Microbiology 150 (Pt 4), 877–883. doi: 10.1099/mic.0.26900-0
Janda, J. M., and Abbott, S. L. (2007). 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: pluses, perils, and pitfalls. J. Clin. Microbiol. 45, 2761–2764. doi: 10.1128/JCM.01228-07
Jarocki, P., Podleśny, M., Glibowski, P., and Targonski, Z. (2014). A new insight into the physiological role of bile salt hydrolase among intestinal bacteria from the genus Bifidobacterium. PLoS ONE 9:e114379. doi: 10.1371/journal.pone.0114379
Jousselin, A., Manzano, C., Biette, A., Reed, P., Pinho, M. G., Rosato, A. E., et al. (2015). The Staphylococcus aureus chaperone PrsA is a new auxiliary factor of oxacillin resistance affecting penicillin-binding protein 2A. Antimicrob. Agents Chemother. 60, 1656–1666. doi: 10.1128/AAC.02333-15
Joyce, S. A., Shanahan, F., Hill, C., and Gahan, C. G. (2014). Bacterial bile salt hydrolase in host metabolism: potential for influencing gastrointestinal microbe-host crosstalk. Gut Microbes 5, 669–674. doi: 10.4161/19490976.2014.969986
Kankainen, M., Paulin, L., Tynkkynen, S., von Ossowski, I., Reunanen, J., Partanen, P., et al. (2009). Comparative genomic analysis of Lactobacillus rhamnosus GG reveals pili containing a human- mucus binding protein. Proc. Natl. Acad. Sci. U.S.A. 106, 17193–17198. doi: 10.1073/pnas.0908876106
Kim, Y. H., Han, K. S., Oh, S., You, S., and Kim, S. H. (2005). Optimization of technical conditions for the transformation of Lactobacillus acidophilus strains by electroporation. J. Appl. Microbiol. 99, 167–174. doi: 10.1111/j.1365-2672.2005.02563.x
Kong, W., Wang, J., Xing, X., Xiao, X., Zhao, Y., Zang, Q., et al. (2011). Antifungal evaluation of cholic acid and its derivatives on Candida albicans by microcalorimetry and chemometrics. Anal. Chim. Acta 689, 250–256. doi: 10.1016/j.aca.2011.01.050
Koskenniemi, K., Laakso, K., Koponen, J., Kankainen, M., Greco, D., Auvinen, P., et al. (2011). Proteomics and transcriptomics characterization of bile stress response in probiotic Lactobacillus rhamnosus GG. Mol. Cell. Proteomics 10:M110.002741. doi: 10.1074/mcp.M110.002741
Kuroda, M., Ohta, T., and Hayashi, H. (1995). Isolation and the gene cloning of an alkaline shock protein in methicillin resistant Staphylococcus aureus. Biochem. Biophys. Res. Commun. 207, 978–984. doi: 10.1006/bbrc.1995.1281
Lebeer, S., Claes, I., Tytgat, H. L., Verhoeven, T. L., Marien, E., von Ossowski, I., et al. (2012). Functional analysis of Lactobacillus rhamnosus GG pili in relation to adhesion and immunomodulatory interactions with intestinal epithelial cells. Appl. Environ. Microbiol. 78, 185–193. doi: 10.1128/AEM.06192-11
Lee, K., Rho, B. S., Pi, K., Kim, H. J., and Choi, Y. J. (2011). Proteomic analysis of protein expression in Lactobacillus plantarum in response to alkaline stress. J. Biotechnol. 153, 1–7. doi: 10.1016/j.jbiotec.2011.02.008
Lin, J., Negga, R., Zeng, X., and Smith, K. (2014). Effect of bile salt hydrolase inhibitors on a bile salt hydrolase from Lactobacillus acidophilus. Pathogens 3, 947–956. doi: 10.3390/pathogens3040947
Makras, L., Triantafyllou, V., Fayol-Messaoudi, D., Adriany, T., Zoumpopoulou, G., Tsakalidou, E., et al. (2006). Kinetic analysis of the antibacterial activity of probiotic lactobacilli towards Salmonella enterica serovar Typhimurium reveals a role for lactic acid and other inhibitory compounds. Res. Microbiol. 157, 241–247. doi: 10.1016/j.resmic.2005.09.002
McAuliffe, O., Cano, R. J., and Klaenhammer, T. R. (2005). Genetic analysis of two bile salt hydrolase activities in Lactobacillus acidophilus NCFM. Appl. Environ. Microbiol. 71, 4925–4929. doi: 10.1128/AEM.71.8.4925-4929.2005
McMillin, M., Frampton, G., Quinn, M., Ashfaq, S., de los Santos, M. III., Grant, S., et al. (2016). Bile acid signaling is involved in the neurological decline in a murine model of acute liver failure. Am. J. Pathol. 186, 312–323. doi: 10.1016/j.ajpath.2015.10.005
Monedero, V., Revilla-Guarinos, A., and Zúñiga, M. (2017). Physiological role of two-component signal transduction systems in food-associated lactic acid bacteria. Adv. Appl. Microbiol. 99, 1–51. doi: 10.1016/bs.aambs.2016.12.002
Morel-Deville, F., Fauvel, F., and Morel, P. (1998). Two-component signal-transducing systems involved in stress responses and vancomycin susceptibility in Lactobacillus sakei. Microbiology 144 (Pt 10), 2873–2883. doi: 10.1099/00221287-144-10-2873
Nyanga-Koumou, A. P., Ouoba, L. I., Kobawila, S. C., and Louembe, D. (2012). Response mechanisms of lactic acid bacteria to alkaline environments: a review. Crit. Rev. Microbiol. 38, 185–190. doi: 10.3109/1040841X.2011.640978
Overbeck, T. J., Welker, D. L., Hughes, J. E., Steele, J. L., and Broadbent, J. R. (2017). Transient MutS-Based hypermutation system for adaptive evolution of Lactobacillus casei to low pH. Appl. Environ. Microbiol. 83:e01120-17. doi: 10.1128/AEM.01120-17
Papadimitriou, K., Alegría Á, A., Bron, P. A., de Angelis, M., Gobbetti, M., Kleerebezem, M., et al. (2016). Stress physiology of lactic acid bacteria. Microbiol. Mol. Biol. Rev. 80, 837–890. doi: 10.1128/MMBR.00076-15
Papadimitriou, K., Zoumpopoulou, G., Foligné, B., Alexandraki, V., Kazou, M., Pot, B., et al. (2015). Discovering probiotic microorganisms: in vitro, in vivo, genetic and omics approaches. Front. Microbiol. 6:58. doi: 10.3389/fmicb.2015.00058
Patel, A. K., Singhania, R. R., Pandey, A., and Chincholkar, S. B. (2010). Probiotic bile salt hydrolase: current developments and perspectives. Appl. Biochem. Biotechnol. 162, 166–180. doi: 10.1007/s12010-009-8738-1
Petrof, E. O., Gloor, G. B., Vanner, S. J., Weese, S. J., Carter, D., Daigneault, M. C., et al. (2013). Stool substitute transplant therapy for the eradication of Clostridium difficile infection: ‘RePOOPulating’ the gut. Microbiome 1:3. doi: 10.1186/2049-2618-1-3
Pithva, S., Shekh, S., Dave, J., and Vyas, B. R. (2014). Probiotic attributes of autochthonous Lactobacillus rhamnosus strains of human origin. Appl. Biochem. Biotechnol. 173, 259–277. doi: 10.1007/s12010-014-0839-9
Ren, J., Sun, K., Wu, Z., Yao, J., and Guo, B. (2011). All 4 bile salt hydrolase proteins are responsible for the hydrolysis activity in Lactobacillus plantarum ST-III. J. Food Sci. 76, M622–628. doi: 10.1111/j.1750-3841.2011.02431.x
Rossi, F., Zotta, T., Iacumin, L., and Reale, A. (2016). Theoretical insight into the heat shock response (HSR) regulation in Lactobacillus casei and L. rhamnosus. J. Theor. Biol. 402, 21–37. doi: 10.1016/j.jtbi.2016.04.029
Schmidt, E. J., Boswell, J. S., Walsh, J. P., Schellenberg, M. M., Winter, T. W., Li, C., et al. (2001). Activities of cholic acid-derived antimicrobial agents against multidrug-resistant bacteria. J. Antimicrob. Chemother. 47, 671–674. doi: 10.1093/jac/47.5.671
Sobhani, I., Tap, J., Roudot-Thoraval, F., Roperch, J. P., Letulle, S., Langella, P., et al. (2011). Microbial dysbiosis in colorectal cancer (CRC) patients. PLoS ONE 6:e16393. doi: 10.1371/journal.pone.0016393
Song, P., Rockwell, C. E., Cui, J. Y., and Klaassen, C. D. (2015). Individual bile acids have differential effects on bile acid signaling in mice. Toxicol. Appl. Pharmacol. 283, 57–64. doi: 10.1016/j.taap.2014.12.005
Tannock, G. W., Taylor, C., Lawley, B., Loach, D., Gould, M., Dunn, A. C., et al. (2014). Altered transcription of murine genes induced in the small bowel by administration of probiotic strain Lactobacillus rhamnosus HN001. Appl. Environ. Microbiol. 80, 2851–2859. doi: 10.1128/AEM.00336-14
Thompson, A. L., Monteagudo-Mera, A., Cadenas, M. B., Lampl, M. L., and Azcarate-Peril, M. A. (2015). Milk-and solid-feeding practices and daycare attendance are associated with differences in bacterial diversity, predominant communities, and metabolic and immune function of the infant gut microbiome. Front. Cell Infect. Microbiol. 5:3. doi: 10.3389/fcimb.2015.00003
Wang, X., Ren, H., Liu, D., Wang, B., Zhu, W., and Wang, W. (2013). H(+) -ATPase-defective variants of Lactobacillus delbrueckii subsp. bulgaricus contribute to inhibition of postacidification of yogurt during chilled storage. J. Food Sci. 78, M297–302. doi: 10.1111/1750-3841.12038
Yan, F., Liu, L., Dempsey, P. J., Tsai, Y. H., Raines, E. W., Wilson, C. L., et al. (2013). A Lactobacillus rhamnosus GG-derived soluble protein, p40, stimulates ligand release from intestinal epithelial cells to transactivate epidermal growth factor receptor. J. Biol. Chem. 288, 30742–30751. doi: 10.1074/jbc.M113.492397
Yu, S., Peng, Y., Chen, W., Deng, Y., and Guo, Y. (2014). Comparative genomic analysis of two-component signal transduction systems in Probiotic Lactobacillus casei. Indian J. Microbiol. 54, 293–301. doi: 10.1007/s12088-014-0456-x
Keywords: probiotics, intra-Species variability, bacterial stress, Lactobacillus rhamnosus, bile salt hydrolase (BSH)
Citation: Arnold JW, Simpson JB, Roach J, Kwintkiewicz J and Azcarate-Peril MA (2018) Intra-species Genomic and Physiological Variability Impact Stress Resistance in Strains of Probiotic Potential. Front. Microbiol. 9:242. doi: 10.3389/fmicb.2018.00242
Received: 14 November 2017; Accepted: 31 January 2018;
Published: 20 February 2018.
Edited by:Rebeca Martín, INRA Centre Jouy-en-Josas, France
Reviewed by:Giuseppe Spano, University of Foggia, Italy
Beatriz Martínez, Consejo Superior de Investigaciones Científicas (CSIC), Spain
Copyright © 2018 Arnold, Simpson, Roach, Kwintkiewicz and Azcarate-Peril. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: M. Andrea Azcarate-Peril, firstname.lastname@example.org