In Silico Prediction of Horizontal Gene Transfer Events in Lactobacillus bulgaricus and Streptococcus thermophilus Reveals Protocooperation in Yogurt Manufacturing

Lactobacillus bulgaricus and Streptococcus thermophilus, used in yogurt starter cultures, are well known for their stability and protocooperation during their coexistence in milk. In this study, we show that a close interaction between the two species also takes place at the genetic level. We performed an in silico analysis, combining gene composition and gene transfer mechanism-associated features, and predicted horizontally transferred genes in both L. bulgaricus and S. thermophilus. Putative horizontal gene transfer (HGT) events that have occurred between the two bacterial species include the transfer of exopolysaccharide (EPS) biosynthesis genes, transferred from S. thermophilus to L. bulgaricus, and the gene cluster cbs-cblB(cglB)-cysE for the metabolism of sulfur-containing amino acids, transferred from L. bulgaricus or Lactobacillus helveticus to S. thermophilus. The HGT event for the cbs-cblB(cglB)-cysE gene cluster was analyzed in detail, with respect to both evolutionary and functional aspects. It can be concluded that during the coexistence of both yogurt starter species in a milk environment, agonistic coevolution at the genetic level has probably been involved in the optimization of their combined growth and interactions.Lactobacillus delbrueckii subsp. bulgaricus (Lactobacillus bulgaricus) and Streptococcus thermophilus have been used in starter cultures for yogurt manufacturing for thousands of years. Both species are known to stably coexist in a milk environment and interact beneficially. This so-called protocooperation, previously defined as biochemical mutualism, involves the exchange of metabolites and/or stimulatory factors (38). Examples of biochemical protocooperation between L. bulgaricus and S. thermophilus include the action of cell wall-bound proteases, produced by L. bulgaricus strains, and formate, required for growth of L. bulgaricus and supplied by S. thermophilus (6, 7). An overview of the interactions between the two yogurt bacteria, including the exchange of CO₂, pyruvate, folate, etc., can be found in a recently published review by Sieuwerts et al. (43). Putative genetic mechanisms underlying protocooperation, however, so far have not been studied in detail.The genomes of two L. bulgaricus strains and three S. thermophilus strains, all used in yogurt manufacturing, have been fully sequenced (3, 32, 33, 34, 39, 44, 46). The available genomic information could provide new insights into the genetic aspects of protocooperation between L. bulgaricus and S. thermophilus through the identification of putative horizontal gene transfer (HGT) events at the genome scale. HGT can be defined as the exchange of genetic material between phylogenetically unrelated organisms (23). It is considered to be a major factor in the process of environmental adaptation, for both individual species and entire microbial populations. Especially HGT events between two species existing in the same niche can reflect their interrelated activities and dependencies (13, 17). Nicolas et al. (36) predicted HGT events between Lactobacillus acidophilus and Lactobacillus johnsonii by analyzing 401 phylogenetic trees, also including the genes of L. bulgaricus. Several HGT events have been predicted in the S. thermophilus strains CNRZ1066 and LMG 18311 (3, 10, 18) as well as in L. bulgaricus ATCC 11842 (46). Moreover, a core genome of S. thermophilus and possibly acquired genes were identified by a comparative genome hybridization study of 47 strains (40).In this study, we describe an in-depth bioinformatics analysis in which we combined gene composition (GC content and dinucleotide composition) and gene transfer mechanism-associated features. Thus, we predicted horizontally transferred genes and gene clusters in the five sequenced L. bulgaricus and S. thermophilus genomes, with a focus on niche-specific genes and genes required for bacterial growth. Identification of HGT events led to a list of putative transferred genes, some of which could be important for bacterial protocooperation and the adaptation to their environment. The evolution and function of the transferred gene cluster cbs-cblB(cglB)-cysE (originally called cysM2-metB2-cysE2 in S. thermophilus), involved in the metabolism of sulfur-containing amino acids, were analyzed in detail. On the basis of our analysis, it can be concluded that both species probably agonistically coevolved at the genetic level to optimize their combined growth in a milk environment and that protocooperation thus includes both biochemical and genetic aspects.