Isolation and purification
About forty-seven LAB were isolated from fifteen raw cow milk samples collected from different areas in Babylon City, Iraq. These bacterial genera were initially purified based on the morphological and microscopically characteristics (cocci or bacilli depending on gram stain).
Results showed that all the isolated bacteria were gram-positive, microscope examination revealed that most of the LAB including 52 isolates were found cocci, whereas the 22 LAB isolates were found to be bacilli. As a result, the LAB isolates were classified into five grouped genera that were: (29.85%)
Lactobacillus, (27.99%)
Lactococcus, (11.36%)
Pediococcus and the rest distributed among
Leuconostocs and Enterococcus.
In some previous studies conducted by
El-Ghany et al. (2016) and
Allaith et al. (2022) found that LABs have cocci or bacilli shapes, are homo-fermentative or hetero-fermentative and have positive or negative gram stain.
In the current study,
Lactobacillus was the prevailing genus among the other genera isolated from the collected raw cow milk samples with a ratio (of 29.85%). These results agreed with
Badis et al., (2004) who indicated that the
Lactobacillus species were the most dominant species isolated from raw goat’s milk involving
L. plantarum, L. helveticus, L. brevis and
L. delbrueckii subsp. Bulgaricus.
Lactobacillus species were isolated and identified from traditional fermented milk by
El-Ghany et al. (2016) and
Allaith et al. (2022). Species of Lactobacillus, Leuconostoc and Lactococcus were identified as prevailing isolates in Dhan traditional butter by
Guetouache et al. (2014).
Holzapfel (1998) and
Osmanagaoglu (2010) reported that the LAB species within the
Pediococcus,
Lactobacillus and
Lactococcus are commonly utilized as probiotics for animals and humans and to produce and preserve various foods.
In Egypt,
El-Ghany et al. (2016) and
Allaith et al., (2022) lactobacillus species were isolated and identified from traditional fermented milk. As shown by
Guetouache et al., (2014) species of
Lactobacillus, Leuconostoc and
Lactococcus were identified as prevailing isolates in Dhan traditional Butter.
Holzapfel, (1998) and
Osmanagaoglu, (2010) reported that the LAB species within the
Pediococcus, Lactobacillus and
Lactococcus are commonly utilized as probiotics for animals and humans and to produce and preserve various foods.
Molecular identification of LAB isolates
Seventy-four PCR amplified of
16S rRNA products for the LAB isolates with a band size of approximately 1500 bp compared with the 1500 bp for the marker as displayed in (Fig 1). In this regard,
(Stackebrandt et al., 1983; El-Ghany et al., 2016;
Allaith, 2018;
Allaith et al., 2022; AL-Shammari et al., 2023) reported the significant role of the
16S rRNA gene in diagnosing the LAB species because the
16S rRNA gene is considered a universal and conserved region.
Sequencing of PCR products amplified and alignment with NCBI
PCR product amplified from each LAB isolate (Fig 1) sequenced in two directions to identify the isolated bacteria at the species level. The nucleotide sequence obtained from each bacterial isolate was analyzed using a GenBank BLASTN alignment (
El-Ghany et al., 2016;
Allaith et al., 2022). Results showed that all the identified isolates belong to LAB distributed among the species
Lactobacillus rhamnosus (12%),
Lactobacillus kefiri (8%),
Lactobacillus curvatus (8%),
Lactococcus lactis (33%),
Pediococcus pentosaceus (2%) and
Pediococcus acidilactici (8%).
The similarity % among local isolates and various strains in NCBI
The aligned nucleotide sequences were obtained from the bacterial isolates with the same strain in NCBI that reached complete similarity (100%) (
El-Ghany et al., 2016). The nucleotide sequences of the
L. rhamnosus isolates in this study were analyzed and compared with those
L. rhamnosus isolates previously registered in NCBI (Fig 2). The genetic similarity percentage of nucleotide sequences ranged between 99-100%. It was the closest to many isolates previously recorded such as those identified in China (CP053619 and CP019305), South Korea (MT322928) and Thailand (MN435970) with a nucleotide sequence similarity of 100%; whereas the lowest sequence homology (99%)was noticed with some
L. rhamnosus isolates previously identified in many countries in the world such those identified in South Korea (MN533908 and USA (HQ811818). Similarly,
El-Ghany et al., (2016) found in their study that
L. rhamnosus strain was used as a control to compare with the rest of the LAB isolates that were identified by
16S rRNA and alignment with NCBI.
L. curvatus isolates showed the highest similarity in the nucleotide sequence, with a match of 100% to the
L. curvatus LC129555 strain in China (Fig 3). Whereas the
L. curvatus isolate had the lowest similarity percentage with the L. curvatus strains identified in many countries, such as South Korea (MN372055.1), Peru (MG031212) and China (EU855223). Additionally, the highest nucleotide sequence similarity reaching 100% was found with the
L. kefiri isolates in Turkey (MH549129) and Italy (NR_042230) while the lowest similarity nucleotide sequence of the
L. kefiri isolate was identified in China (HM218366) reaching 99% (Fig 4). These results agreed with
Guetouache et al. (2014) that the species of Lactobacillus, Leuconostoc and Lactococcus were identified as prevailing isolates in fermented food by 16S rRNA and alignment with NCBI.
Additionally
, P. acidilactici isolates had the highest similarity of the nucleotide sequence with a match of 100%
P. acidilactici AF515229 strain in Japan (Fig 5). Whereas the
P. acidilactici isolate had the lowest similarity percentage with the
L. curvatus CP006854 strains identified in Denmark reaching 98%. The results displayed as shown in (Fig 6) comparison of the nucleotide sequence obtained from
P. pentosaceus, identified in this study, showed that there is a 100% similarity with most of the
P. pentosaceus isolates, deposited in NCBI, such as the
P. pentosaceus MH045191 in India, (MT510516, MT516146 and MT464432) in China and (CP028269) in Korea. In this regard,
El-Ghany et al., (2016) reported that
P. parvulus F1030 was isolated from fermented food and identified by 16SrRNA sequencing by alignment with NCBI. Similarly,
Immerstrand et al. (2010) found in their study that isolates and diagnoses the
P. parvulus 2.6 by 16S DNA sequencing and alignment with NCBI.
Morphological and molecular identification results also demonstrated that 25 out of 74 isolates obtained in this study belong to the species
L. lactis. The alignment analysis of the sequence for
L. lactis had 100% similarity with most of L. lactis strains recorded previously in NCBI such as
L. lactis (MN749817) in Korea, (MN749817 and MH549138) in Morocco, (HE805077) in Thailand, (MH549123 and MH549135) in Turkey and (NR-040955) in the USA. The other isolates available in NCBI gave a 99% similarity with
L. lactis strains in NCBI such as (NR_116443) in the USA and (AM406671) in Ireland as shown in (Fig 7). These results agreed with
Nomura et al., (2006) that the
L. lactis were isolated from milk and plant and identified by
16S rRNA and alignment with NCBI, subsequently studying the genetic and phenotypic traits for it.
The phylogenetic tree analysis
An unrooted phylogenetic tree was generated using the Neighbor-Joining method based on comparing the nucleotide sequences generated from (LAB) isolates including (
L. curvatus, L. kefiri, L. rhamnosus) and
L. plantarum (as a control) as shown in (Fig 8) (
Saitou and Nei, 1987). The most accurate tree had a total branch length with evolution distances estimated following the maximum likelihood approach
(Tamura et al., 2013) indicating the number of base changes per site. The results of the phylogenetic analysis proved the tree had two clades distributed among 22 isolates, the first clade consisted of Polyphyletic groups that had various species that lacked a common ancestor, including (1,2,3,4,5,6,7,8,9 isolates) belonging to
L. rhamnosus, (10,11,12,13,14,15 isolates) belonging to
L. kefiri and (22) isolate belonging to
L. plantarum. Besides, the second clade consisted of monophyletic groups including (16, 17, 18, 19, 20 and 21) isolates belonging to
L. curvatus.
On the other hand, an unrooted A phylogenetic tree was generated using the Neighbor-Joining method based on a comparison of the nucleotide sequences generated from the isolates
P. acidilactici and
P. pentosaceus as shown in (Fig 9) (
Saitou and Nei, 1987).
The most accurate tree had a total branch length with evolution distances estimated following the maximum likelihood approach
(Tamura et al., 2013), indicating the number of base changes per site. The results of the phylogenetic analysis proved the tree had two clades of monophyletic groups (lacked a common ancestor) distributed among 8 isolates including (23, 24, 25, 26, 27 and 28 isolates) belonging to
P. acidilactici and 29, 30 isolates belonging to
P. pentosaceus.
Moreover, an unrooted phylogenetic tree was generated using the Neighbor-Joining method based on comparing the nucleotide sequences generated from
L. lactis isolates (Fig 10) (
Saitou and Nei, 1987;
Alhissnawi et al., 2024). The most accurate tree had a total branch length with evolution distances estimated following the maximum likelihood approach
(Tamura et al., 2013), indicating the number of base changes per site. The results of the phylogenetic analysis proved the tree had two clades of paraphyletic groups (lacked a common ancestor) distributed among 24 isolates including (31-55 isolates) belonging to
L. lactis.
The advantages of studying the phylogenetic tree relationship among LAB species are in their interrelation-ships. Furthermore, environmental factors affected bacterial evolution that led to the detection of polymorphism in LAB because of the variation sequencing of nucleotides in the
16S rRNA gene (
Romanenko et al., 2008;
Jasim and Maaroof, 2017;
Mnati et al., 2021; Al-Shugeairy et al., 2021).