Previous Article | Next Article ![]()
Journal of Bacteriology, July 2007, p. 4624-4634, Vol. 189, No. 13
0021-9193/07/$08.00+0 doi:10.1128/JB.00337-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Genomic Sciences Graduate Program, North Carolina State University, Raleigh, North Carolina 27695,1 Department of Food Science, North Carolina State University, Raleigh, North Carolina 27695,2 Southeast Dairy Foods Research Center, North Carolina State University, Raleigh, North Carolina 276953
Received 7 March 2007/ Accepted 16 April 2007
|
|
|---|
|
|
|---|
Bile salts are one of the major components of bile and consist of a cholesterol-derived ring structure that is amide linked to an amino acid, either glycine or taurine. The amphipathic nature of conjugated bile salts allows them to act as emulsifiers that facilitate the dietary absorption of lipids, but this property also gives them the ability to emulsify the lipid membranes of bacterial cells (7). Some microbes in the gastrointestinal tract have the ability to hydrolyze this amide linkage, creating deconjugated bile salts that can damage bacterial membranes and lead to cell death (18). Some bile salts have the ability to cross lipid membranes into the cellular cytoplasm, where they can damage DNA and proteins (20, 25). Bile has been shown to induce transcription of molecular chaperones, such as GroESL, and other proteins involved in DNA and protein repair (20). While bile-specific defense mechanisms are not wholly understood, a diverse group of them have been elucidated, including those that export bile from the cell and chemically modify bile constituents (7, 27, 32, 34).
Mechanisms to sense the presence of bile and alter transcription are not well characterized but may involve a two-component regulatory system (2CRS). 2CRSs allow bacteria to sense and respond to changes in their environment after receiving an environmental signal through transmembrane sensing domains of the histidine protein kinase (HPK). Once it receives a signal input, the HPK autophosphorylates through ATP hydrolysis. The phosphoryl group is then transferred to the regulatory domain of the response regulator (RR), which in turn promotes a transcriptional response through its DNA binding domain (28). Previous studies have shown the ability of bacteria to alter the transcription of necessary genes in the presence of bile or bile constituents, including potential induction through the RR of a 2CRS (26), and have indicated the increased expression of genes in a 2CRS in response to bile (20, 30). However, to date, no system has been functionally characterized in any gram-positive organism.
The complete genome sequence of L. acidophilus NCFM allows examination of the genomic elements that confer probiotic properties on the strain (1). cDNA microarrays have proven a key tool for the examination of bacterial transcriptomes (2, 6). This study employed a whole-genome microarray representing 97.4% of the open reading frames (ORFs) to investigate how L. acidophilus NCFM reacts to the presence of bile in its environment and to identify genes that contribute to bile tolerance.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Strains and plasmids
|
|
View this table: [in a new window] |
TABLE 2. Primers used to generate gene deletions
|
|
View this table: [in a new window] |
TABLE 3. Primers used in RT-QPCR
|
cDNA probe preparation and microarray hybridization. NCFM cells were grown to an OD600 between 0.2 and 0.3 (representing early log phase). The cells were then pelleted at room temperature and resuspended in either MRS or MRS with 0.5% (wt/vol) oxgall (Difco) and incubated for 30 min at 37°C. The cells were then cold shocked in an ethanol-dry ice bath, pelleted at 4°C, flash frozen, and stored at 80°C. Total RNA was isolated using Trizol (Invitrogen)-chloroform extraction, and residual genomic DNA was removed using DNase I, amplification grade (Invitrogen), according to the manufacturer's directions. The RNA concentration and quality were determined by spectrophotometry and agarose gel electrophoresis. Twenty-five micrograms of RNA were aminoallyl labeled by reverse transcription with random hexamers in the presence of aminoallyl dUTP (Sigma Chemical Co.) using Superscript II reverse transcriptase (Life Technologies), followed by fluorescent labeling with N-hydroxysuccinimide-activated Cy3 or Cy5 ester (Amersham Biotechnology). The labeled cDNA probes were purified using a PCR purification kit (QIAGEN). Coupling of the Cy3 and Cy5 dyes to the aminoallyl-dUTP-labeled cDNA and hybridization of samples to microarrays were performed according to the protocols outlined by The Institute for Genome Research (http://www.tigr.org/tdb/microarray/protocolsTIGR.shtml). Cy5- and Cy3-labeled cDNA probes were hybridized to the arrays for 16 h at 42°C. After hybridization, the slides were washed twice in high-stringency buffer (1x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate] with 0.2% sodium dodecyl sulfate[SDS]) for 5 min each. The first wash was performed at 42°C and the second at room temperature. The slides were then washed in a low-stringency buffer (0.1x SSC containing 0.2% SDS) for 5 min at room temperature and finally in 0.1x SSC for two 2.5-min washes at room temperature.
Data normalization and gene expression analysis. Immediately after the slides were washed, fluorescence intensities were acquired at 10-µm resolution using a ScanArray 4000 microarray scanner (Packard Biochip BioScience; Biochip Technologies LLC) and stored as TIFF images. Signal intensities were quantified, background was subtracted, and the data were normalized using the QuantArray 3.0 software package (Perkin-Elmer). Two slides per experiment, each containing triplicate arrays, were hybridized reciprocally to Cy3- and Cy5-labeled probes (dye swapping). Spots were analyzed by adaptive quantitation. The data were median normalized. When the local background intensity was higher than the spot signal, no data were recorded for those spots. The mean of the six ratios per gene was recorded, and the ratio between the average absolute pixel values for the replicated spots of each gene, with or without treatment, represented the change in gene expression. Genes in potential operons were considered for analysis if at least one gene of the operon showed significant expression changes and the remaining genes showed trends toward that expression. Confidence intervals and P values for the change were also calculated with the use of a two-sample t test. P values of < 0.05 were considered significant.
Site-specific integration and deletion in L. acidophilus NCFM. Two noncontiguous fragments flanking an internal region of the target ORF were amplified using PCR primers listed in Table 2. The fragments were cloned into pORI28 and then transformed into L. acidophilus NCK1392 containing the temperature-sensitive helper plasmid pTRK669. Selection of integrants was then carried out as described previously (29). Upon successful integration of the plasmid into the genome, a single integrant colony was propagated in the absence of antibiotic selection and replica plated onto MRS agar and MRS agar containing Em. Ems cells were screened for a deletion mutation using PCR with primers flanking the targeted region, and deletions were confirmed by Southern blotting.
Bile survival assays. To test deletion mutants for survival in the presence of bile, cells were grown to early log phase (OD600 = 0.2 to 0.3), serially diluted, and plated on MRS agar containing 0%, 0.5%, 1.0%, and 2.0% (wt/vol) oxgall (Difco) or bile salts (Sigma-Aldrich, St. Louis, MO) using a Whitley Automatic Spiral Plater. The concentrations of oxgall were chosen to emulate the concentrations found in the human intestinal tract, which range from 0.2 to 2% (16).
MICs of oxgall, glycocholic acid (GCA), taurocholic acid (TCA), glycodeoxycholic acid (GDCA), taurodeoxycholic acid (TDCA), and sodium chloride were determined by spotting approximately 104 early-log-phase cells (OD600 = 0.2 to 0.3) onto MRS agar plates containing increasing concentrations of these compounds. Inhibitory concentrations of Triton X-100 and SDS (Fisher Scientific, Fair Lawn, NJ) were determined by inoculating 200-µl MRS broth cultures containing increasing amounts of these compounds. Inhibition was reported as the lowest concentration tested that prevented the growth of the strain. The increments of concentrations tested were as follows: bile salts, 1 µg/ml; salt, 5 µg/ml (1%); TX-100, 0.5%; SDS, 0.02%.
RT-PCR and RT-QPCR. RNA samples were collected and treated as described for microarray experiments. Reverse transcriptase-quantitative PCR (RT-QPCR) primers are shown in Table 3. PCR was performed using the Quanti-Tect Reverse Transcriptase PCR kit (QIAGEN) according to the manufacturer's directions, but the total reaction size was reduced to 20 µl. PCR was performed on a Bio-Rad My-IQ single-color detection system (Bio-Rad Laboratories, Hercules, CA). Absolute mRNA copy numbers were determined using a standard curve of known concentrations, and transcripts were enumerated using the MyIQ software. Copy numbers per cell were determined by normalizing the transcript number of the queried ORF to the transcript number of LBA383 (DNA polymerase III, delta subunit).
RT-PCR was carried out to determine intergenic transcription using the Superscript II kit (Invitrogen) according to the manufacturer's directions. cDNAs were generated using random primers. PCR was performed using primers specific to the intergenic regions of the hypothetical cotranscript.
Northern blotting.
Northern blotting of total RNA was carried out using the NorthernMax Gly kit (Ambion). Briefly, RNA was electrophoresed on a 1% agarose gel and transferred to a Magnacharge nylon transfer membrane (MSI). The membrane was probed with a 692-bp fragment of LBA1430 labeled with [
-32P]CTP. Transfer, hybridization, washing, and exposure to films (Eastman Kodak Company, Rochester, NY) were carried out according to standard methods and the manufacturer's instructions.
|
|
|---|
2-fold upregulated or downregulated) and met the criteria for statistical significance (P < 0.05). Overall, 78 genes (3.9% of the genome) were shown to be induced while 168 (8.4%) were repressed. |
View this table: [in a new window] |
TABLE 4. L. acidophilus NCFM ORFs differentially expressed upon exposure to 0.5% oxgalla
|
![]() View larger version (14K): [in a new window] |
FIG. 1. Distribution of ORFs within COG classifications. The black bars represent the numbers of upregulated ORFs (total, 78), and the gray bars represent downregulated ORFs (total, 168).
|
Among genes not differentially expressed upon exposure to bile were the two bile salt hydrolases characterized previously. In some species, including Listeria monocytogenes and Lactobacillus plantarum, these genes have been shown to play a role in bile tolerance. In L. acidophilus NCFM, knockout mutations of the two bsh genes showed no differences in survival in the presence of bile, despite encoding active enzymes that are able to deconjugate the specific bile salts GCA, GDCA, TCA, and TDCA, as well as glycochenodeoxycholic acid and taurochenodeoxycholic acid (24).
Operon configuration and annotations. A cluster of genes seemingly unrelated by function (LBA1425 to LBA1432) was significantly induced in the presence of bile, as indicated by the microarray experiments (ranging from 2.1- to 3.0-fold upregulation). To date, this group of genes has not been shown to be differentially expressed under any test conditions examined (2, 3, 6). This cluster was particularly interesting because it contained genes for a type IIIA HPK (LBA1430) and an RR of the OmpR family (LBA1431), both with uncharacterized functions (1). Its induction in the presence of bile indicated that this 2CRS may potentially regulate transcription under these conditions. LBA1425 and LBA1426 were annotated as a hypothetical protein and an unknown protein, respectively. Neither protein contained a COG that could give a clue to its function, although BLAST analysis showed that LBA1425 had 46% identity to a putative cell surface hydrolase in L. plantarum WCFS1. LBA1427 was annotated as a putative oxidoreductase but showed less than 30% similarity to oxidoreductases from other species. It had similarity to COG0656, aldo/keto reductases (E value, 4e11). LBA1428 showed homology (58% identity) to a predicted hypothetical protein in Lactobacillus gasseri, and it contained similarity to a COG of predicted redox proteins (E value, 3e5). LBA1429 was annotated as a transporter based on its similarity to the major facilitator superfamily COG. This protein showed homology to a multidrug transport protein (L. plantarum WCFS1; 46% identity). LBA1432, a hypothetical protein, showed similarity to the RelA/SpoT COG (E value, 5e21); however, this protein showed stronger similarity to COG2357, an uncharacterized protein domain of unknown function conserved in bacteria.
Because of their upregulation in the presence of oxgall and their physical proximity in the genome, ORFs LBA1425 to LBA1432 were analyzed as a potential cotranscript. Clone Manager Software (Scientific and Educational Software, Cary, NC) revealed a region of dyad symmetry downstream of LBA1425 with a Gibbs free-energy value of 14.9 kcal/mol, followed closely by a T-rich region (Fig. 2). According to de Hoon et al. (13), this free-energy value falls below the average predicted free-energy value for terminators in NCFM (13.8 kcal/mol), indicating good potential for the area to form a rho-independent terminator in an RNA transcript. RT-PCR and Northern analysis were employed to functionally determine the length of transcription in this region. These experiments indicated that the cotranscription of genes in the region extended from LBA1432 through LBA1425, for an overall transcript length of 7 kb.
![]() View larger version (15K): [in a new window] |
FIG. 2. Organization of a bile-induced operon in L. acidophilus NCFM. The nucleotide sequence of the upstream promoter is shown, with the ribosomal binding site and 10 and 35 boxes in boldface type. Predicted rho-independent terminators are shown with the predicted free energy of folding. ORFs upregulated in the presence of bile are shown in gray. Putative annotations for each ORF are shown. Hypo., hypothetical.
|
![]() View larger version (10K): [in a new window] |
FIG. 3. Survival of early-log-phase L. acidophilus NCFM (), NCK 1871 ( LBA1427; oxidoreductase) ( ), NCK 1873 ( LBA1428; hypothetical) ( ), NCK 1875 ( LBA1429; transporter) ( ), NCK 1877 ( LBA1430; HPK) ( ), NCK1879 ( LBA1431; RR) ( ), and NCK 1881 ( LBA1432; hypothetical) ( ) on MRS agar plates with increasing concentrations of oxgall. The error bars represent the standard errors of the mean for triplicate counts.
|
![]() View larger version (13K): [in a new window] |
FIG. 4. Percent recovery of early-log-phase cells on MRS agar supplemented with 0.2% GDCA and TDCA compared to MRS agar. Bars (left to right): L. acidophilus NCFM, NCK 1871 ( LBA1427; oxidoreductase), NCK 1873 ( LBA1428; hypothetical), NCK 1875 ( LBA1429; transporter), NCK 1877 ( LBA1430; HPK), NCK1879 ( LBA1431; RR), and NCK 1881 ( LBA1432; hypothetical). The error bars represent the standard errors of the mean for triplicate counts.
|
|
View this table: [in a new window] |
TABLE 5. Inhibitory concentrations of selected compounds on three independent cultures of L. acidophilus NCFM strains
|
![]() View larger version (10K): [in a new window] |
FIG. 5. Effects of 2CRS gene mutations on induction of LBA1429 (black bars) and LBA1432 (gray bars) in the presence of 0.5% oxgall. The error bars represent standard errors of the mean for four biological replicates.
|
|
|
|---|
Analysis of a previously uncharacterized operon upregulated in the presence of bile suggested that the 2CRS within it was concerned with induction of genes involved in response to bile. HPKs and RRs that are in close proximity typically act together in sensing and regulatory capacities. Analysis of these mutant strains on oxgall plates initially indicated that this 2CRS was no exception, with both the HPK and RR deletion mutant strains showing a marked decrease in recovery in the presence of increasing levels of oxgall. The increased survival and increased growth in the presence of the specific bile salt TDCA seen in the HPK mutant strain were surprising.
Some 2CRSs control their own expression, so the transcriptional effects of mutating the HPK and the RR were examined on genes in the operon in the presence and absence of oxgall. Mutation in the RR led to an enhanced level of induction, indicating that the RR acts in the capacity of a repressor of this operon. Since the majority of the differentially regulated genes in the genome (Fig. 1) were repressed in the presence of bile, the role of this RR in the bile transcriptome could be substantial. Mutation in the HPK had no effect on induction of the operon in the presence of bile. This observation, along with the differences in phenotype between the HPK and RR mutant strains, raises two possibilities regarding this 2CRS. The first is that these two proteins do not form a cognate 2CRS but that they interact with other 2CRS components in the genome to promote a transcriptional response. Through comparative analysis, Grebe and Stock (15) demonstrated a correlation between the sequence subfamilies of HPKs and RRs that were known to act cognately. Previous analysis of the 2CRS of NCFM (1) showed that LBA1430 is a type IIIA HPK. LBA1431 belongs to the OmpR subfamily of RRs, which type IIIA HPKs have been shown to interact with exclusively. The second possibility regarding this 2CRS is that the RR is activated independently of phosphorylation by its cognate HPK. There are several reports of 2CRSs in which acetyl phosphate acts as a phospho-donor to the RR in the absence of the HPK or where the HPK and acetyl phosphate act synergistically to promote activation (4, 12, 22, 23, 38). The biological purpose for this alternative activation is uncertain, but it has been shown to play a role in transcriptional control during acid stress (4). Since there are no orphan HPK genes present in the genome that might indicate RR cross talk, and sequence analysis indicates the high probability that LBA1430 and LBA1431 form a cognate 2CRS, it is likely that this RR was activated by another phosphate-donating member, such as acetyl phosphate, in the absence of its HPK. Since deletion of the HPK did result in a decrease in survival in the presence of oxgall, it is most likely that this protein does play a role in signal transduction. Also, since the 2CRS within the operon does not act to induce it in the presence of bile, this single system is most likely a player in a more complex regulatory network that controls bile-influenced transcription. Experiments to examine the global transcriptional effects of both the HPK and the RR mutations in the presence of bile are in progress.
In addition to the 2CRS, other genes were functionally shown to play a role in bile tolerance. Decreased survival of the LBA1429 mutant in bile and growth inhibition by the cholic acid-based bile salts, but not salt or detergents, supports the hypothesis that this transporter interacts with bile salts. Transport proteins, including those belonging to the major facilitator superfamily, have been shown to act in bile salt transport in other species, including E. coli, L. monocytogenes, L. plantarum, and Eubacterium sp. (14, 21, 32, 34). Such transporters have been shown to act in various capacities in species, including influx, efflux, and antiport, and are believed to both import salts for modification and export them to prevent their deleterious effects. The differences in the roles of the transporters can lead to varying effects of mutations on bile tolerance. The exact role that LBA1429 plays is unclear from these experiments, although the LBA1429 mutant strain's lowered growth inhibition on GCA and TCA suggest interaction with cholic acid. These experiments introduce the first transporter gene in NCFM that has been implicated in bile tolerance. Sequence analysis of hypothetical protein LBA1432 does not suggest a role that it plays in bile tolerance, although weak similarity to RelA/SpoT domains might suggest that the protein is associated with gene regulation in times of stress or starvation (10). Increased recovery of LBA1427 (oxidoreductase) and LBA1428 (hypothetical) mutant strains in the presence of bile, compared to the wild type, was surprising considering that these genes are members of an operon induced by bile. Since both proteins contain homology to oxidoreductases, they could play a role in bile salt modification, as it is known that some species of intestinal bacteria possess the ability to modify bile salts through oxidation and reduction (7, 9). These proteins may act as previously unrecognized bile salt modifiers, although a role for modifiers of this type in bile tolerance or sensitivity has not been established. Many bacterial interactions with bile, such as those that lead to modification of bile salts, are poorly understood. It is possible that these proteins play an important role with bile in the gastrointestinal system but that any effects cannot be observed in vitro. Induction of the operon, along with the repression by LBA1431 (RR), could indicate a need for fine control of these genes and suggest that these proteins provide a necessary benefit to the cell.
Other findings from the microarray experiments showed upregulation of genes involved in carbohydrate metabolism but downregulation of genes involved in many other aspects of growth, indicating that bile slows growth in the cells. The induction of the lactose and galactose operons was an intriguing finding, as this strain is widely utilized in yogurt and other fermented dairy products. It might be speculated that as some lactobacilli evolved in the small intestines of mammals, a frequently encountered energy source was lactose from milk. In the gastrointestinal tract, bile may act as an environmental signal for the cells to increase production of proteins needed to metabolize this energy source. This finding could suggest that consumption of this Lactobacillus strain in dairy foods might promote survival in the gastrointestinal tract, as it is better prepared to metabolize those available carbohydrates.
Bile induction of potential adhesion genes was also identified by microarray experiments and confirmed using RT-QPCR. Adhesion of bacterial cells to intestinal cells is considered an important attribute for lactobacilli in the gastrointestinal tract (8, 35). Induction of such adhesion determinants suggests that bile may act to signal arrival in the intestinal tract. While bile's effect on bacterial cells has been indicated as largely detrimental, induction of these genes could indicate a new role for bile as a location indicator in the cellular environment.
Examination of the global transcriptional effects of bile led to the discovery of an inducible operon that encodes proteins involved in both bile tolerance and, surprisingly, sensitivity. This duality illustrates the complexity of transcriptional regulation and reflects how little bile tolerance is currently understood. A 2CRS that aids bile tolerance is a new finding. The fact that L. acidophilus NCFM utilizes one of its nine 2CRSs in bile tolerance shows the importance of this characteristic in species that travel through the intestinal tract.
We thank E. Durmaz, M. Miller, T. Duong, and P. Ruse for their helpful discussions and technical advice.
Published ahead of print on 20 April 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»