Molecular Microbial Analysis of Lactobacillus Strains Isolated from the Gut of Calves for Potential Probiotic Use
The intestinal microbiota has an influence on the growth and health status of the hosts. This is of particular interest in animals reared using intensive farming practices. Hence, it is necessary to know more about complexity of the beneficial intestinal microbiota. The use of molecular methods has revolutionized microbial identification by improving its quality and effectiveness. The specific aim of the study was to analyze predominant species of Lactobacillus in intestinal microbial ecosystem of young calves. Forty-two lactic acid bacteria (LAB) isolated from intestinal tract of young calves were characterized by: Amplified Ribosomal DNA Restriction Analysis (ARDRA), by using Hae III, Msp I, and Hinf I restriction enzymes, and 16S rDNA gene sequencing. ARDRA screening revealed nine unique patterns among 42 isolates, with the same pattern for 29 of the isolates. Gene fragments of 16S rDNA of 19 strains representing different patterns were sequenced to confirm the identification of these species. These results confirmed that ARDRA is a good tool for identification and discrimination of bacterial species isolated from complex ecosystem and between closely related groups. This paper provides information about the LAB species predominant in intestinal tract of young calves that could provide beneficial effects when administered as probiotic.
The natural microbiota of the gastrointestinal tract has an influence on the biochemistry, immunology, physiology, and nonspecific host's resistance against infectious diseases . Therefore, the role of the intestinal microbiota is of vital importance in the nutritional status of the host, and particularly in farm animals that are reared in intensive systems . Because of this it is necessary to determine the complexity of the intestinal flora and recognize the different microorganisms that compose it. This is particularly relevant in the probiotic therapy field where it is necessary to distinguish between probiotics and autochthonous microbiota .
Lactobacilli are part of the normal human gastrointestinal microbiota and may also be found in other mammalian species [4–7] and birds . It has been reported that some Lactobacillus species have probiotic properties and that they are “live micro-organisms which when administered in adequate amounts confer a health benefit on the host” .
The first step in the probiotic production is the isolation and identification of the normal components of the gut microbiota, because one of the desirable characteristics of strains used as probiotics is that they should be autochthonous to the ecosystem of which they will be part once ingested . Then, we must assess the probiotic and technological properties of the strains  in order to select the best examples that will form the probiotic inoculum. The inocula can be either monostrain or multistrain . The latter is more effective because it can use the complementary and synergistic effects of each microorganism .
To analyse and rapidly identify bacteria from microbial communities, classical physiological and biochemical tests are not adequate because the bacterial populations involved often have similar nutritional requirements and grow under similar environmental conditions. Currently, there is a wide variety of molecular strategies, such as PCR with specific primers, DGGE, RAPD, PFGE, FISH, RFLP, and PCR-ARDRA, among others , which are available to determine the species diversity of Lactobacillus .
The comparison of sequences of the 16S rDNA gene is a very reliable method for sorting and identifying bacterial species . Because these genes are highly conserved and are present in large numbers of copies within each bacterial cell, their use as a molecular target has increased in the recent years .
The ARDRA technique is a highly discriminatory method, simple and quick to identify Gram positive nonspore bacteria. Many authors have shown that this method is suitable for the discrimination of different species of Lactobacillus [8, 16, 17]. In addition, many LAB used as starters or probiotics have been identified with the ARDRA methodology .
The aim of this study was to analyse the predominant species of Lactobacillus that constitute the intestinal microbial ecosystem of young calves, by means of isolating and identifying strains through the application of the ARDRA technique and 16S rDNA gene sequencing, as a prior step to the design of a probiotic inoculum for cattle.
2. Materials and Methods
2.1. Bacterial Isolation
Isolates were taken from the mucosa of cecum and jejunum of six young calves reared in intensive conditions. For this, a selective Lactobacillus Anaerobic MRS broth with Vancomycin and Bromocresol green (LAMVAB, 7) was used. Forty-one colonies were multiplied in MRS broth for 24 hours at 37°C. For preservation, the cultures were frozen at 80°C with the addition of glycerol 25% v/v.
2.2. DNA Isolation
An aliquot of 2 mL of each 24 hours culture was centrifuged at 14000 g (for 5 minutes). The sediment was frozen at 20°C for 24 hours to facilitate the breaking of the cells. The DNA was extracted according to Marmur  modified by Kurzak et al.  and then resuspended in 50 L of TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8). An aliquot of 5 L of this template DNA was added directly to the PCR tube. The amount of DNA obtained was quantified by measuring it in an UV spectrum (260 nm) and its integrity was visualised by agarose gel electrophoresis to 0.7% w/v, by staining with ethidium bromide and visualising under UV light.
2.3. 16S rDNA Amplification
The 16S rDNA gene was amplified by PCR with a thermal cycler (MJ Research). DNA fragments of approximately 1.5 kpb were amplified using the primers 27F (-AGAGTTTGATCCTGGCTCAG-) and 1492R (-GGYTACCTTGTTACGACTT-). Each PCR tube (50 L) contained a reaction mix of 10 L 5X PCR buffer for Taq polymerase (Promega), 1.5 mM MgCl2, 200 M of each deoxynucleotide triphosphate (Promega), 0.4 M of each primer and 2 U of Taq Polymerase (Promega) and 5 L of template DNA. The termocycle programme was as follows: 94°C for 5 minutes; 30 cycles of 94°C for 1 minute, 55°C for 1 minute and 72°C for 1 minute; and a final extension step at 72°C for 7 minute. After cycling, the PCR products were visualised by electrophoresis on a 1% w/v agarose gel (40 minute, 75 V), by staining with ethidium bromide (0.5 g/mL) and visualising under UV light (DyNA Light UV Transilluminator, LabNet, UV light source wavelength 302 nm).
In order to achieve complete digestion, restriction mixes (20 L of final volume) were carried out for 4 hours at 37°C. Each reaction tube contained 2 L of 10X incubation buffer, 0.2 L of bovine serum albumin, 6 U of the respective restriction enzyme, 2.5 L of bidistilled water and 15 L of PCR product. Three restriction enzymes were used: Hae III, Msp I and Hinf I (Promega). The resulting digestion products were visualised under UV-light (LabNet Transilluminator, UV light source wavelength 302 nm), after agarose gel electrophoresis 3% w/v (90 minutes, 75 V) by staining with ethidium bromide (0.5 g/mL). Restriction patterns identical to the sequenced strains led to the identification of the corresponding species .
2.5. In Silico Study
For this study, Nebcutter software testing protocols (http://tools.neb.com/NEBcutter2/index.php) were used. The theoretical restriction profiles of the 16S rDNA sequence of each species, which had a high percentage of identity in the alignment of the BLAST algorithm, were compared with profiles of the isolates in this study. Besides, theoretical restriction profiles of the 16S rDNA gene sequences were obtained from other species of Lactobacillus and Enterococcus to determine the power of the ARDRA technique to discriminate from other species.
The PCR products of 19 representative strains of each restriction group were purified with the Wizard PCR SV Gel & PCR Clean-Up System kit (Promega) and sequenced. The sequences were compared with the sequences deposited in the GenBank database using the BLAST algorithm (http://www.ncbi.nlm.nih.gov/BLAST/; 1).
2.7. Nucleotide Sequence Accession Numbers
The sequences were deposited in the GenBank database using the web-based data submission tool, BankIt (http://www.ncbi.nlm.nih.gov/BankIt, 1).
3.1. Identification of LAB Isolates by ARDRA
Lactic acid bacteria isolated from calves’ intestinal tract samples yielded nine unique ARDRA patterns among the 42 isolates tested (Figure 1). One ARDRA pattern clearly dominated the samples, accounting for 29 of the 42 colonies tested. The other ARDRA patterns from the isolated bacteria were present at a low frequency (Table 2). Most of the ARDRA patterns derived from lactobacilli. Although the isolation medium was specific for Lactobacillus spp., two of the patterns found belonged to Enterococcus spp.
The restriction of the amplified fragment of the 16S rDNA gene with Hae III generated six different profiles. Lactobacillus plantarum, Weissella paramesenteroides, L. salivarius, L. ruminis and L. mucosae presented specific profiles for each of these species. Instead, Pediococcus acidilactici, L. farciminis, L. curvatus and E. hirae showed restriction fragments different from the species listed above but not distinguishable among them (Figure 1(a)).
The enzyme Msp I also showed six different restriction profiles. Species that showed characteristic profiles were: L. salivarius, L. curvatus, L. mucosae and E. hirae. It was not possible to differentiate between P. acidilactici and W. paramesenteroides and between L. plantarum, L. ruminis and L. farciminis (Figure 1(b)).
Hinf I produced seven restriction profiles, five of which were typical of L. ruminis, L. curvatus, L. farciminis, L. mucosae and E. hirae. The restriction profiles produced by the species L. plantarum and P. acidilactici, and by W. paramesenteroides and L. salivarius were not able to distinguish between them (Figure 1(c)).
The restriction profile of each isolate and its association with the concerned species are detailed in Table 1.
3.2. In Silico Study
The size of the fragments obtained by the theoretical restriction of the sequences obtained from GenBank that had a higher percentage of identity with the isolations coincided with the restriction fragments obtained in the in vitro study.
On the other hand, some nonisolated species that belonged to the same genus or phylogenetic group as that of the isolates were distinguishable with the in silico study, in most cases by restriction with the Hinf I enzyme (Table 3).
3.3. Identification by Sequencing of the 16S rRNA Gene
Nineteen representative clones of the ARDRA profiles observed were selected for sequencing.
The sequences of the gene fragments obtained from the 16S rDNA were aligned with those from GenBank using the BLAST algorithm. Table 1 shows the percentage of identity of the isolated strains in relation to those found in the database and the access number to GenBank for each of the sequences obtained.
We found that the isolates DSPV 322T, 324T, 325T, 327T, 329T, 333T, 340T, 344T, and 355T represented the ARDRA patterns that were observed most frequently (29 times) among the 42 isolates tested, and that their 16S sequence was most closely related to Lactobacillus salivarius.
This species was found in all the calves studied. The other species were represented by one, two or three isolates and were found in one or two calves depending on the case. On the other hand, with the exception of P. acidilactici, which was found in the jejunum and in the cecum, the species isolated in the large intestine were different from those isolated in the small intestine (Table 2).
The identification of microbial species through the use of phenotypic methods can sometimes be uncertain, complicated and time-consuming. The use of molecular methods has revolutionised their identification, by improving the quality and effectiveness of this identification. Some of these methodologies use either the rDNA spacer region or its target. These techniques are useful for both the identification and reliable detection of different bacterial species as well as the monitoring of the species . In this way, members of a probiotic multistrain inoculum can be identified and distinguished from strains that share the same environment such as starters in foods (yogurt, cheese, etc.).
The use of species-specific primers or probes is not applicable in environments where there are several Lactobacillus species because prior knowledge of them is required. In these cases, more general molecular tools should be applied . The techniques used to identify Lactobacillus species in different environments are the comparison of total or partial sequences of 16S rDNA, ARDRA patterns of 16S rDNA or the intergenic region of the 16S-23S rDNA [5, 8, 22].
While the use of 16S DNA sequencing methods gives a high resolution of the diversity of microbial species in an environment, it is very time-consuming and too costly to be used for routine screening of samples. Methods for the initial analysis of faecal samples should be rapid and able to give a broad view of the microbial ecology. ARDRA has been used to compare bacterial isolates within a wide range of microbial communities. The advantages of ARDRA are that it is rapid, reproducible, relates to microbial diversity, and will be invaluable in analysing a greater number of samples together with experimental objectives such as dietary interventions .
In the present work, ARDRA allowed us to differentiate Enterococcus hirae from the rest of the Lactobacillus spp. isolates. This differentiation was observed by restricting with any of the three enzymes used.
The Lactobacillus isolated belonged to two groups: the L. casei-Pediococcus group and the Leuconostoc group. The latter includes the species Weissella and Lactobacillus paramesenteroides, which can be differentiated from the L. casei-Pediococcus group by the typical profile obtained with the restriction enzymes Hae III and Hinf I. This methodology also allowed the distinction between phylogenetically related species belonging to the L. casei-Pediococcus group. These species were L. ruminis, L. salivarius, L. curvatus, P. acidilactici, L. farciminis, L. plantarum, and L. mucosae, which have a 16S rDNA homology of 90.3 to 99% .
The similarity between the profiles obtained by the in silico study of the sequences of the GenBank and the isolates revealed that the strains of the same species had similar profiles. This result proved to be another tool for the identification of the species. The possibility to obtain these profiles, distinguishable between the isolates, together with the differentiation of these isolates from other related species (Table 3), shows that this technique allows the distinction of species with high homology. Such is the case of E. faecium and E. faecalis, which were also found in the intestines of calves  and can be distinguished from E. hirae by the restriction enzyme Hinf I. There are species within the same phylogenetic groups, such as L. fermentum and L. reuteri, which have higher homology than others and are most closely related to L. mucosae . Despite these similarities, in the present work the in silico study showed that the latter could be distinguished from the first two by the ARDRA methodology (Table 3). These results show that the ARDRA technique is a tool that highly discriminates between LAB species and seems to group the isolates by species and then sequence some exponents of each group. This may save both time and money when it is necessary to analyse large numbers of isolates.
L. salivarius was the predominant species in the gastrointestinal tract of calves. It was found in the cecum of all individuals (Table 2) and in some animals it was the only species isolated. This species was also detected by Schneider et al.  in calves reared in the same geographical area. L. salivarius is an inhabitant of the gastrointestinal tract of other species such as chickens , pigs , and humans [4, 5]. Many strains that correspond to this species have been studied to evaluate their probiotic properties. Some strains isolated from infants have shown antimicrobial capacity against pathogens , and, in particular, L. salivarius CTC2197 was able to prevent the colonization of Salmonella enteritidis in chickens .
The probiotic properties of microorganisms are characteristic of each strain. Therefore, belonging to a species is not sufficient to guarantee the possession of such properties. For a strain to be used as a probiotic, it should be considered GRAS, that is, possessing probiotic effects and technological capabilities suitable for its propagation and preservation over time. Therefore, in order to select the best specimens, in future works we aim at evaluating the probiotic properties of each isolate obtained in this study (in vitro: aggregation, coaggregation with pathogens, production of inhibitory substances, bile and pH resistance; in vivo: effect on calves performance, challenge with pathogen). The knowledge of such properties will allow the development of an inoculum for young calves to improve their performance in intensive farming systems.
This study is part of the CAI+D Project financed by Universidad Nacional del Litoral. L. S. Frizzo and L. P. Soto are fellows at Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET, Argentina).
S. Salminem and A. von Wright, Lactic Acid Bacteria: Microbiology and Functional Aspects, Marcel Dekker, Nueva York, NY, USA, 2nd edition, 1998.
M. Rosmini, G. Sequeira, I. Guerrero-Legarreta et al., “Producción de probióticos para animales de abasto: importancia del uso de la microbiota intestinal indígena,” Revista Mexicana de Ingeniería Quimica, vol. 3, pp. 187–197, 2004.View at: Google Scholar
A. L. McCartney, “Application of molecular biological methods for studing probiotics and the gut flora,” British Journal of Nutrition, vol. 88, supplement 1, pp. S29–S37, 2002.View at: Google Scholar
S. Silvi, M. C. Verdenelli, C. Orpianesi, and A. Cresci, “EU project Crownalife: functional foods, gut microflora and healthy ageing: isolation and identification of Lactobacillus and Bifidobacterium strains from faecal samples of elderly subjects for a possible probiotic use in functional foods,” Journal of Food Engineering, vol. 56, no. 2-3, pp. 195–200, 2003.View at: Publisher Site | Google Scholar
H. G. H. J. Heilig, E. G. Zoetendal, E. E. Vaughan, P. Marteau, A. D. L. Akkermans, and W. M. de Vos, “Molecular diversity of Lactobacillus spp. and other lactic acid bacteria in the human intestine as determined by specific amplification of 16S ribosomal DNA,” Applied and Environmental Microbiology, vol. 68, no. 1, pp. 114–123, 2002.View at: Publisher Site | Google Scholar
R. Schneider, M. Rosmini, M. Ehrmann, and R. Vogel, “Identificación de bacterias lácticas componentes de la microbiota típica de los terneros criados en condiciones artificiales,” Revista FAVE—Ciencias Veterinarias, vol. 3, no. 1-2, pp. 7–15, 2004, (Argentina).View at: Google Scholar
L. L. Guan, K. E. Hagen, G. W. Tannock, D. R. Korver, G. M. Fasenko, and G. E. Allison, “Detection and identification of Lactobacillus species in crops of broilers of different ages by using PCR-denaturing gradient gel electrophoresis and amplified ribosomal DNA restriction analysis,” Applied and Environmental Microbiology, vol. 69, no. 11, pp. 6750–6757, 2003.View at: Publisher Site | Google Scholar
FAO/WHO, “Evaluation of health and nutritional properties of probiotics in food including powder milk with live lactic acid bacteria. Expert consultation report,” Food and Agriculture Organization of the United Nations and World Health Organization, Córdoba, Argentina, October 2001.View at: Google Scholar
L. Soto, L. Frizzo, E. Bertozzi et al., “Milk evaluation as growth and cold preservation medium of a probiotic inoculum for young calves,” Journal of Animal and Veterinary Advances, vol. 8, no. 7, pp. 1353–1360, 2009.View at: Google Scholar
C. M. Lee, C. C. Sieo, C. M. V. L. Wong, N. Abdullah, and Y. W. Ho, “Sequence analysis of 16S rRNA gene and 16S-23S rRNA gene intergenic spacer region for differentiation of probiotics Lactobacillus strains isolated from the gastrointestinal tract of chicken,” Annals of Microbiology, vol. 58, no. 1, pp. 133–140, 2008.View at: Google Scholar
G. Tannock, “Molecular assessment of intestinal microflora,” American Journal of Clinical Nutrition, vol. 73, no. 2, supplement, pp. 410S–414S, 2001.View at: Google Scholar
M. Ventura, M. Elli, R. Reniero, and R. Zink, “Molecular microbial analysis of Bifidobacterium isolates from different environments by the species-specific amplified ribosomal DNA restriction analysis (ARDRA),” FEMS Microbiology Ecology, vol. 36, no. 2-3, pp. 113–121, 2001.View at: Publisher Site | Google Scholar
J. Marmur, “A procedure for the isolation of deoxyribonucleic acid from microorganisms,” Journal of Molecular Biology, vol. 13, pp. 208–218, 1961.View at: Google Scholar
P. Kurzak, M. A. Ehrmann, and R. F. Vogel, “Diversity of lactic acid bacteria associated with ducks,” Systematic and Applied Microbiology, vol. 21, no. 4, pp. 588–592, 1998.View at: Google Scholar
G. Blaiotta, V. Fusco, D. Ercolini, M. Aponte, O. Pepe, and F. Villani, “Lactobacillus strain diversity based on partial hsp60 gene sequences and design of PCR-restriction fragment length polymorphism assays for species identification and differentiation,” Applied and Environmental Microbiology, vol. 74, no. 1, pp. 208–215, 2008.View at: Publisher Site | Google Scholar
P. Vandamme, B. Pot, M. Gillis, P. De Vos, K. Kersters, and J. Swings, “Polyphasic taxonomy, a consensus approach to bacterial systematics,” Microbiological Reviews, vol. 60, no. 2, pp. 407–438, 1996.View at: Google Scholar
S. Roos, F. Karner, L. Axelsson, and H. Jonsson, “Lactobacillus mucosae sp. nov., a new species with in vitro mucus-binding activity isolated from pig intestine,” International Journal of Systematic and Evolutionary Microbiology, vol. 50, no. 1, pp. 251–258, 2000.View at: Google Scholar
Y.-J. Lee, W.-K. Yu, and T.-R. Heo, “Identification and screening for antimicrobial activity against Clostridium difficile of Bifidobacterium and Lactobacillus species isolated from healthy infant faeces,” International Journal of Antimicrobial Agents, vol. 21, no. 4, pp. 340–346, 2003.View at: Publisher Site | Google Scholar
M. Pascual, M. Hugas, J. I. Badiola, J. M. Monfort, and M. Garriga, “Lactobacillus salivarius CTC2197 prevents Salmonella enteritidis colonization in chickens,” Applied and Environmental Microbiology, vol. 65, no. 11, pp. 4981–4986, 1999.View at: Google Scholar