Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2003, p. 6750-6757, Vol. 69, No. 11
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.11.6750-6757.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada T6G 2P5,1 Department of Microbiology, University of Otago, Dunedin, New Zealand2
Received 2 June 2003/ Accepted 19 August 2003
|
|
|---|
|
|
|---|
The crop, ileum, cecum, and colon of poultry are known to harbor bacterial populations (16, 27). Recent reports have investigated the composition of the ileal (13) and cecal (35) microflora using bacteriological culture and culture-independent methods. Lactobacilli are numerous in the ileum of broilers, whereas the cecal microflora is dominated by obligately anaerobic bacteria and bacteria yet to be cultivated. From the results of culture-based studies, it has been determined that the microflora of the crop has a simple composition and is dominated by lactobacilli (16, 27). Colonization of the surface of the stratified, squamous epithelium of the crop by lactobacilli has been reported by Fuller (6) and Morishita et al. (18). Lactobacillus salivarius, Lactobacillus fermentum or Lactobacillus reuteri, and Lactobacillus acidophilus were the species most commonly detected (16, 27). These studies were conducted prior to the reclassification of L. acidophilus, which has been divided into two DNA homology groups containing six related species (5, 11, 15). DNA homology group A consists of L. acidophilus (A1), Lactobacillus crispatus (A2), Lactobacillus amylovorus (A3), and Lactobacillus gallinarum (A4); DNA homology group B consists of Lactobacillus gasseri (B1) and Lactobacillus johnsonii (B2). Even with molecular methods, it is difficult to distinguish between members of this group. Methods that have been used successfully include the following: DNA-DNA hybridization and various biochemical properties (5, 11, 15); analysis of whole-cell protein profiles by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (3, 8, 12, 21); randomly amplified polymorphic DNA (3, 8, 12, 24); sequencing of 16S and 23S ribosomal DNA (rDNA) and elongation factor Tu (2, 14, 29, 30); oligonucleotide probes and primers for species-specific hybridizations and PCR, respectively (21, 23, 24, 32); ribotyping (25); amplified fragment length polymorphism (8); and amplified ribosomal DNA restriction analysis (ARDRA) (17, 23, 31). There is a clear need, however, for simpler methods to differentiate the members of the L. acidophilus group to ensure that large-scale microbial ecological studies are logistically possible.
A detailed analysis of the crop microflora of broiler chickens using nucleic acid-based methods has not previously been reported. Our study therefore had two main aims. First, we used PCR combined with denaturing gradient gel electrophoresis (PCR-DGGE) to compare the crop microflora of birds of different ages. Second, we cultured lactobacilli from the crops of the birds and devised an ARDRA technique by which the members of the L. acidophilus group could be rapidly identified.
|
|
|---|
During the experiment, 86 birds were sampled as follows: 6 birds on day 0 (hatching day) and 10 each on 1, 3, 7, 14, 21, 28, 35, and 42 days of age. At each sampling time, one or two birds from each pen were selected randomly and in such a way that the stocking density was maintained. In a subsequent experiment, 10 birds were sampled on day 0. At each sample time, birds were euthanized and transported to the research laboratory where the crops were aseptically removed. For each bird, the crop was placed in a sterile petri plate and weighed. A section of crop tissue and contents weighing approximately 1 g or the entire crop (if it weighed less than 1 g) was transferred into a sterile 15-ml Pyrex tissue grinder and homogenized with 9 ml of 0.85% saline. For each bird, 1 ml of crop homogenate was stored at -80°C for nucleic acid-based analysis of the bacterial communities. The remaining crop homogenate was used for selective enumeration and collection of lactobacilli.
Propagation and enumeration of lactobacilli.
The crop homogenate was used to make a series of 10-fold dilutions (10-2 to 10-7) in sterile 0.85% NaCl. For each dilution, 100 µl was spread plated on Lactobacillus selective (LBS) agar (BBL) and incubated at 37°C for 48 h under anaerobic conditions (5% CO2, 10% H2, 85% N2). In the subsequent experiment using day 0 birds, 1 ml of the crop homogenate was used to make LBS pour plates. The number of CFU of presumptive lactobacilli per gram of crop for each bird was determined from the number of colonies on LBS plates. The colony morphologies on the counted LBS plates were also noted. For each crop, a total of 10 colonies, representing each colony type, were selected, restreaked on MRS agar (Difco) plates, and incubated at 37°C for 48 h under anaerobic conditions. The bacterial colonies were removed from the MRS agar plates, suspended in MRS broth containing 50% glycerol, and stored at -80°C. For routine propagation of Lactobacillus reference strains and crop isolates (Table 1), MRS broth glycerol stocks were streaked onto MRS agar and incubated at 37°C under anaerobic conditions as outlined above. Reference strains of Enterococcus and Pediococcus (Table 1) were grown on APT agar (Difco) at 30°C under aerobic conditions.
|
View this table: [in a new window] |
TABLE 1. Reference strains used in this study
|
DNA was extracted from Lactobacillus, Enterococcus, and Pediococcus type strains and crop isolates by the method of Walter et al. (34). The DNA pellet was dissolved in 20 µl of TE buffer and diluted 20-fold for PCR.
PCR-DGGE analysis of crop DNA with universal bacterial primers and primers specific for lactic acid bacteria.
PCR was conducted using either individual DNA or pooled crop DNA as the template. The pooled samples were prepared by combining the crop DNA (1 µl of each) from all 10 crops collected at the same sampling time. The V3 region of the 16S rRNA gene from the crop DNA was amplified in a GeneAmp PCR System 9700 (Applied Biosystems, Foster City, Calif.) using universal bacterial primers HDA1-GC and HDA2 (HDA PCR), following the protocol of Walter et al. (34) and subjected to DGGE (HDA PCR-DGGE). The DGGE analysis was performed with the Bio-Rad DCode Universal Mutation Detection System (Hercules, Calif.) by the method of Walter et al. (34). The V3 region of the 16S rDNA was amplified from the total crop DNA using group-specific bacterial primers Lac1 and Lac2-GC (Lac PCR) and subjected to DGGE (Lac PCR-DGGE) by the method of Walter et al. (33). Identification ladders for DGGE were prepared by combining the HDA or Lac PCR products prepared from DNA extracted from the type and reference Lactobacillus strains (Table 1).
DGGs were stained with ethidium bromide and viewed by UV transillumination. DGGE profiles were compared using Dice's similarity coefficient (Dsc) with the Bionumerics software package (Applied Maths, Austin, Tex.). When Dsc analysis was performed, only profiles within a gel were compared, not between gels. The average Dsc was calculated by adding the values of single profile comparisons for the age range stated and dividing by the total number of Dscs.
Identification of bacteria by sequencing DNA fragments.
DNA fragments generated by Lac primers were extracted from DGGs by the method of Knarreborg et al. (13). Following purification, the DNA was reamplified with the Lac1 and Lac2 (without GC clamp) primers using the PCR protocol described above. The resulting PCR products were purified with the QIAquick purification kit (Qiagen, Mississauga, Ontario, Canada), ligated into pGEM-T (Promega, Madison, Wis.), and used to transform Escherichia coli JM109. Transformants were plated on Luria-Bertani (LB) agar (26) containing ampicillin, 5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside (X-Gal), and isopropyl-ß-D-thiogalactopyranoside (IPTG) at concentrations of 100 µg/ml, 80 µg/ml, and 0.5 mM, respectively. Several colonies were subcultured, and the plasmid DNA was isolated by the alkaline lysis procedure of Sambrook et al. (26). Plasmid DNA was used as the template in a PCR with Lac1 and Lac2-GC primers, and the PCR products were analyzed by DGGE in order to compare the migration of the cloned DNA with the migration of the desired band from the original PCR-DGGE crop profile. The pGEM-T insert DNA was amplified using T7 and SP6 primers, and sequencing was conducted by the Agricultural, Food and Nutritional Science Biotech Core, University of Alberta. The sequences were compared with those in the GenBank database using the BLAST algorithm (1).
Identification of Lactobacillus isolates with ARDRA.
The total DNA extracted from each reference strain or crop isolate was used as the template for PCR amplification of either the 16S rRNA gene (16S rDNA) or the 16S rRNA gene plus the entire 16S-23S rRNA intergenic region (16-23S rDNA). The total 16S rRNA gene (1.5 kb) was amplified using SacI-POmod (5'-CCGAGCTCAACAGAGTTTGATCCTGGCTCAG-3') and SalI-T7-PC5 (5'-GGTCGACCGTTAATACGACTCACTATAGGGATACCTTGTTACGACTT-3') primers (22). The following primers were used to amplify the 16-23S rDNA (2 kb): (i) Lb16a (5'-GTGCCTAATACATGCAAGTCG-3'), which corresponds to nucleotides (nt) 17 to 36 of the 16S rDNA of L. crispatus ATCC 33820 (GenBank accession no. AF257097) (this study), and (ii) 23-1B (5'-GGGTTCCCCCATTCGGA-3'), which corresponds to nt 123 to 113 of Lactobacillus 23S rDNA and which was developed by Tannock et al. (29). PCR was performed as follows: (i) 5 min at 94°C; (ii) 25 cycles, with 1 cycle consisting of 45 s at 94°C, 30 s at 53°C, and 1.5 min at 72°C; and (iii) a final extension step of 7 min at 72°C. The PCR products were digested with HaeIII (16S rDNA and 16-23S rDNA) or MseI (16-23S rDNA) following the manufacturer's directions (Invitrogen [Burlington, Ontario, Canada] and New England Biolabs [Pickering, Ontario, Canada]). The resulting banding patterns were analyzed on a 2% agarose gel.
To confirm the species designation of the crop isolates, the V2-V3 region of the 16S rDNA was sequenced using primers Lb16a, HDA2, Lac1, and Lac2 (without GC clamp). Sequencing and analysis were conducted as outlined above.
|
|
|---|
![]() ![]() View larger version (186K): [in a new window] |
FIG. 1. PCR-DGGE profiles generated from Lactobacillus type strains used in this study. (A) HDA PCR-DGGE profiles on a 22 to 55% DGG. (B) Lac PCR-DGGE profiles on a 30 to 45% DGG. The species are indicated above the lanes. L. johnsonii* is the non-type strain L. johnsonii ATCC 11506 listed in Table 1.
|
![]() View larger version (116K): [in a new window] |
FIG. 2. PCR-DGGE profiles generated from pooled crop DNA, using primer pair HDA1-GC and HDA2 (22 to 55% DGG). The sample time is indicated above the lanes. Lane M contains the identification ladder composed of PCR products from the following reference strains of Lactobacillus: L. aviarius (Lav) ATCC 43234, L. johnsonii (Lj) ATCC 33200, L. fermentum (Lf) ATCC 14931, L. acidophilus (Lac) ATCC 4356, L. salivarius (Ls) ATCC 11741, and L. reuteri (Lr) ATCC 23272.
|
![]() View larger version (120K): [in a new window] |
FIG. 3. PCR-DGGE profiles generated from pooled crop DNA, using primer pair Lac1 and Lac2-GC (30 to 45% DGG). The sample time is indicated above the lanes. Lane M contains the identification ladder composed of PCR products from the following reference strains of Lactobacillus: L. aviarius (Lav) ATCC 43234, L. johnsonii (Lj) ATCC 33200, L. acidophilus (Lac) ATCC 4356, L. crispatus (Lc) ATCC 33820, L. salivarius (Ls) ATCC 11741, and L. reuteri (Lr) ATCC 23272. Numbered fragments were extracted and sequenced as outlined in the text.
|
Enumeration of lactobacilli in chicken crops.
For each sampling time, the number of CFU of presumptive lactobacilli per gram of crop for each bird was determined and used to calculate the mean log10 CFU per gram and standard deviation. For day 1 birds, lactobacilli were below the detection range in five birds, and a large range in log10 CFU/gram was observed for the other five birds (5.84 ± 1.23 log10 CFU/g). From days 1 to 7, the number of lactobacilli increased 1,000-fold, reaching a maximum average population of 9.00 ± 0.41 log10 CFU/g on day 7. From day 14 onwards, the average number of lactobacilli stabilized between 8.14 and 8.51 log10 CFU/g, corresponding to the period of compositional stability revealed by the PCR-DGGE results. In the initial sampling, the numbers of Lactobacillus in the day 0 crops were below the detection limit of spread plating. A subsequent experiment was conducted using an additional 10 day 0 crops using pour plating instead of spread plating. Colonies were observed for bacteria from two crops. From the bacteria from one crop, two colonies of different morphology grew; from the bacteria from the other crop, 36 colonies with the same morphology grew. Two and four colonies, respectively, were selected for further investigation.
Identification of lactobacilli and differentiation between the L. acidophilus group isolates using ARDRA.
The Lac PCR-DGGE analysis indicated that most crops contained L. reuteri, L. salivarius, and representatives of different species of the L. acidophilus group. In order to rapidly and accurately identify the Lactobacillus isolates obtained from each bird, an ARDRA technique was developed. Initially, only the 16S rRNA gene was amplified from the type strains of the L. acidophilus group and other Lactobacillus reference strains representing different species (Table 1). The 1.5-kb PCR products were digested with various restriction enzymes. HaeIII digestion generated species-specific banding patterns for L. reuteri, L. salivarius, L. fermentum, and L. aviarius as well as group-specific banding patterns for L. acidophilus group A and B species (data not shown).
Tannock et al. (29) used the sequence of the 16S-23S rRNA intergenic region to identify various members of the L. acidophilus group to the species level, and in silico restriction analysis of the intergenic region supported inclusion of this region in order to generate species-specific patterns. The 16-23S rDNA (2 kb) was amplified with primers Lb16a and 23-1B and PCR products generated from type strains were digested with HaeIII. Identical banding patterns were observed for L. amylovorus ATCC 33620 and L. gallinarum ATCC 33199. The banding pattern for L. acidophilus ATCC 4356 was similar to that obtained for ATCC 33620 and ATCC 33199 but differed slightly in the number and size of fragments below 100 bp (Fig. 4A). The banding pattern for L. crispatus ATCC 33820 differed from that of ATCC 33620, ATCC 33199, and ATCC 4356 in that the second largest fragment was
350 bp instead of
400 bp (Fig. 4A), and this difference was used to discriminate L. crispatus from the other group A species. The HaeIII banding pattern for L. johnsonii ATCC 33200 and ATCC 11506 and L. gasseri ATCC 33323 were the same but were different from those for L. crispatus and the other group A species. Unique fragmentation patterns were also evident for L. reuteri ATCC 43272, L. salivarius subsp. salivarius ATCC 11741, L. fermentum ATCC 14931, and L. aviarius ATCC 43234 (Fig. 4A). Subsequently, the HaeIII-ARDRA of the 16-23S rDNA was used to make an initial identification of the crop isolates (Fig. 4A, lanes D25 to D30). Slight variations in the HaeIII-ARDRA patterns were noted for L. reuteri (Fig. 4A, compare banding patterns of the type strain, D25, and D28) and L. crispatus (data not shown) and were characterized by the absence and/or different intensities of bands between 250 to 300 bp for L. reuteri and between 300 bp to 400 bp for L. crispatus. The identity of the crop isolates with these L. reuteri-like and L. crispatus-like patterns was confirmed by sequencing (discussed below) (Table 2). The variability observed with these patterns is likely due to strain-specific differences in the 16-23S spacer regions (10, 20).
![]() ![]() View larger version (188K): [in a new window] |
FIG. 4. ARDRA profiles of type and reference strains and crop isolates. (A) HaeIII fragmentation patterns of the 16-23S rRNA gene amplified from the type and reference strains and crop isolates as indicated. L. johnsonii* is the non-type strain L. johnsonii ATCC 11506 listed in Table 1. (B) MseI fragmentation patterns of the 16-23S rRNA gene amplified from type strains of the L. acidophilus group as indicated. In panels A and B, the 1-kb ladder (Invitrogen) was used as the molecular weight marker, and the corresponding fragment sizes are indicated at the sides of the gels.
|
|
View this table: [in a new window] |
TABLE 2. Identification of chicken crop isolates by HaeIII- and MseI-ARDRA and 16S rRNA gene sequence analysis
|
Six isolates from day 0 and 14 isolates from day 1 generated HaeIII- and MseI-ARDRA banding patterns that did not match those of the Lactobacillus reference strains. In an attempt to identify these isolates, other species of LAB were analyzed. HaeIII-ARDRA generated species-specific patterns for Enterococcus faecalis, Enterococcus faecium, Pediococcus acidilactici, and Pediococcus pentosaceus (Fig. 5). Unique patterns were also observed for Lactococcus lactis, Leuconostoc gelidum, and Carnobacterium piscicola (data not shown), suggesting this technique may be applicable for identifying species within these genera as well. Sixteen of the 20 isolates generated HaeIII banding patterns that corresponded to that of P. acidilactici (Fig. 5, strains 11-1 and 16-1); the remaining four isolates generated HaeIII banding patterns that corresponded to E. faecium (Fig. 5, strain 96-1) (Table 2).
![]() View larger version (72K): [in a new window] |
FIG. 5. HaeIII-ARDRA profiles of 16-23S rRNA gene amplified from type strains, other LAB reference strains, and crop isolates. The 1-kb ladder (Invitrogen) was used as the molecular weight marker, and the corresponding fragment sizes are indicated to the left of the gel.
|
|
|
|---|
A rapid and simple ARDRA method was developed to distinguish between the members of the L. acidophilus group. MseI-ARDRA generated unique patterns for all species of the L. acidophilus group, identifying L. crispatus, L. johnsonii, and L. gallinarum among crop isolates. Ventura et al. (31) also described an ARDRA technique in which only the 16S rDNA was amplified, which was then digested with three enzymes, Sau3AI, HinfI, and DraI. Sau3AI digestion yielded group A- and group B-specific patterns and differentiated among species of lactobacilli that were less closely related, such as L. paracasei, L. salivarius, L. reuteri, and L. fermentum. Following group designations provided by Sau3AI digestion, HinfI produced a unique banding pattern for L. acidophilus and L. gallinarum, and DraI could distinguish L. crispatus and L. amylovorus species and L. johnsonii and L. gallinarum species. All three enzymes, however, were used to distinguish between the group A species, and two enzymes were required to identify the group B species. Roy et al. (23) proposed a combined group- or species-specific PCR followed by ARDRA. Group-specific PCR was used to differentiate group A from group B lactobacilli. A second PCR was used to amplify the 16S rDNA, and L. acidophilus and L. amylovorus were identified using HinfI. Using an isoschizomer of MseI, Tru9I, Roy et al. were able to distinguish between L. crispatus and L. gallinarum species and L. gasseri and L. johnsonii species. Thus, at least two PCRs and two restriction digestions were required to identify the L. acidophilus group species. The ARDRA used in our study simplified identification in that a single MseI restriction digestion of one PCR product easily distinguished both closely related and less related lactobacilli. HaeIII-ARDRA distinguished group A L. acidophilus from group B, differentiated L. crispatus from other group A L. acidophilus, and produced unique patterns for the non-Lactobacillus LAB studied.
Our study has provided detailed knowledge of the acquisition of the Lactobacillus microflora in the broiler crop. Of particular importance was the observation of the dynamics of the crop microflora during the life of the birds, demonstrating both rapid changes during days 1 to 7 and the establishment of a stable microflora after day 14. Further, since it has been proposed that the crop microflora acts as a bacterial inoculum for the remainder of the gut (7), knowledge of the composition of this bacterial collection is critical in understanding the contribution of the microflora members to the well-being of the avian host and for selection of species for probiotics. Given the crop microflora dynamics observed in our study, it is doubtful that efficacious and scientifically valid probiotics can be derived without the use of this information because it impinges on the types of bacteria that will inoculate the digesta in the crop, suppress the multiplication of contaminating bacteria, and influence the biochemistry of the broiler gut (4, 7). The baseline information generated by this study will be essential in planning husbandry methods that utilize feed supplements other than antimicrobial drugs for the efficient production of broilers.
We thank Todd Klaenhammer (North Carolina State University) and Lynn McMullen (University of Alberta) for providing reference cultures. We also thank Matt Rawluk and Dorthe Nielsen for assistance with sampling and plating. The assistance of the animal research technicians at the Alberta Poultry Research Centre, University of Alberta, was greatly appreciated.
|
|
|---|
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»