Whole genome SNP analysis of bovine B. anthracis strains from Switzerland reflects strict regional separation of Simmental and Swiss Brown breeds in the past.

Bacillus anthracis is an evolutionarily young species that presents an extremely low genetic diversity due to its slow mode of propagation, determined by short replication phases and long sporulation periods. In our ongoing efforts to elucidate phylogenetic relationships between European B. anthracis isolates, the genomes of five strains from Switzerland belonging to lineages B.Br.CNEVA and A.Br.Aust94 were sequenced. Comparative analysis with additional, available genomes from both lineages, were used to reconstruct the substructure of these populations. Genome-wide single-nucleotide polymorphism analysis revealed two phylogeographical different groups among the Swiss B.Br.CNEVA strains (central and eastern Switzerland), that define the oldest most recent common ancestor of the B.Br.CNEVA lineage currently known. Age-old practices of livestock selection, breeding and preservation of unique traits of local breeds in Alpine valleys have likely favored differentiation of regional B. anthracis populations over centuries and the emergence of genetically distinct strains in an otherwise similar environment.

Comparative genomics of Bacillus anthracis from the wool industry highlights polymorphisms of lineage A.Br.Vollum.

BACKGROUND: With the advent of affordable next-generation sequencing (NGS) technologies, major progress has been made in the understanding of the population structure and evolution of the B. anthracis species. Here we report the use of whole genome sequencing and computer-based comparative analyses to characterize six strains belonging to the A.Br.Vollum lineage. These strains were isolated in Switzerland, in 1981, during iterative cases of anthrax involving workers in a textile plant processing cashmere wool from the Indian subcontinent. RESULTS: We took advantage of the hundreds of currently available B. anthracis genomes in public databases, to investigate the genetic diversity existing within the A.Br.Vollum lineage and to position the six Swiss isolates into the worldwide B. anthracis phylogeny. Thirty additional genomes related to the A.Br.Vollum group were identified by whole-genome single nucleotide polymorphism (SNP) analysis, including two strains forming a new evolutionary branch at the basis of the A.Br.Vollum lineage. This new phylogenetic lineage (termed A.Br.H9401) splits off the branch leading to the A.Br.Vollum group soon after its divergence to the other lineages of the major A clade (i.e. 6 SNPs). The available dataset of A.Br.Vollum genomes were resolved into 2 distinct groups. Isolates from the Swiss wool processing facility clustered together with two strains from Pakistan and one strain of unknown origin isolated from yarn. They were clearly differentiated (69 SNPs) from the twenty-five other A.Br.Vollum strains located on the branch leading to the terminal reference strain A0488 of the lineage. Novel analytic assays specific to these new subgroups were developed for the purpose of rapid molecular epidemiology. CONCLUSIONS: Whole genome SNP surveys greatly expand upon our knowledge on the sub-structure of the A.Br.Vollum lineage. Possible origin and route of spread of this lineage worldwide are discussed.

Whole Genome-Sequencing and Phylogenetic Analysis of a Historical Collection of Bacillus anthracis Strains from Danish Cattle.

Bacillus anthracis, the causative agent of anthrax, is known as one of the most genetically monomorphic species. Canonical single-nucleotide polymorphism (SNP) typing and whole-genome sequencing were used to investigate the molecular diversity of eleven B. anthracis strains isolated from cattle in Denmark between 1935 and 1988. Danish strains were assigned into five canSNP groups or lineages, i.e. A.Br.001/002 (n = 4), A.Br.Ames (n = 2), A.Br.008/011 (n = 2), A.Br.005/006 (n = 2) and A.Br.Aust94 (n = 1). The match with the A.Br.Ames lineage is of particular interest as the occurrence of such lineage in Europe is demonstrated for the first time, filling an historical gap within the phylogeography of the lineage. Comparative genome analyses of these strains with 41 isolates from other parts of the world revealed that the two Danish A.Br.008/011 strains were related to the heroin-associated strains responsible for outbreaks of injection anthrax in drug users in Europe. Eight novel diagnostic SNPs that specifically discriminate the different sub-groups of Danish strains were identified and developed into PCR-based genotyping assays.

Molecular diversity of Bacillus anthracis in the Netherlands: investigating the relationship to the worldwide population using whole-genome SNP discovery.

Bacillus anthracis, the causative agent of anthrax, has been widely described as a clonal species. Here we report the use of both canonical SNP analysis and whole-genome sequencing to characterize the phylogenetic lineages of B. anthracis from the Netherlands. Eleven strains isolated over a 25-years period (1968-1993) were paired-end sequenced using parallel sequencing technology. Five canSNP groups or lineages, i.e. A.Br.001/002 (n=6), A.Br.Aust94 (n=2), A.Br.008/011 (n=1), A.Br.011/009 (n=1) and A.Br.Vollum (n=1) were identified. Comparative analyses, with a focus on SNPs discovery, were carried out using a total of 52 B. anthracis genomes. A phylogeographic "Dutch" cluster within the dominant A.Br.001/002 group was discovered, involving isolates from a single outbreak. Diagnostic SNPs specific to the newly identified sub-groups were developed into high-resolution melting SNP discriminative assays for the purpose of rapid molecular epidemiology. Phylogenetic relationships with strains from other parts of the world are discussed.

Single-nucleotide polymorphism discrimination using high-resolution melting analysis for the genotyping of Bacillus anthracis.

High-resolution melting (HRM) is a post-PCR technique that determines with high precision the melt profile of PCR products using a new generation of double-stranded DNA-binding dyes and accurate fluorescence data acquisition over small temperature increments. The method can be used to interrogate small sets of single-nucleotide polymorphisms (SNPs). Here, we describe a simple and cost-effective HRM-based method for the screening of 14 phylogenetically informative SNPs within the genome of Bacillus anthracis that subtype the species into 13 major sublineages or subgroups. Fourteen monoplex and seven duplex SNP-discrimination assays have been designed. We detail the parameters most important for the successful application of HRM for B. anthracis genotyping.

Multiplexed genotyping of Bacillus anthracis by Luminex xMap suspension array.

The Luminex® xTAG technology is a medium to high throughput, open methodology able to test many single nucleotide polymorphisms (SNPs) in a single reaction and a minimum time. Multiplex SNPs interrogation are conducted on the Luminex xMAP system, which uses lasers to read universal tag, color-coded microspheres that attach to specific nucleic acid sequences. The present method describes a Multiplex Oligonucleotide Ligation-PCR procedure (MOL-PCR) for the simultaneous interrogation of 13 phylogenetically informative SNPs within the genome of Bacillus anthracis. The reported 13-plex assay enables efficient B. anthracis genotyping into major sublineages and groups. While cost-effective compared to other monoplex methods, the present MOL-PCR method also offers a high degree of flexibility and scalability. It can easily accommodate newly identified SNPs to increase resolving power to the canSNP typing of B. anthracis.

Draft Genomes of Three Strains Representative of the Bacillus anthracis Diversity Found in France.

We report here the draft genomes of three Bacillus anthracis strains isolated in France: 08-8_20 (A.Br.001/002), 99-100 (A.Br.011/009), and 00-82 (B.Br CNEVA). The total lengths of assemblies are 5,440,708 bp, 5,446,472 bp, and 5,436,014 bp for 08-8_20, 99-100, and 00-82, respectively.

Genotyping of French Bacillus anthracis strains based on 31-loci multi locus VNTR analysis: epidemiology, marker evaluation, and update of the internet genotype database.

BACKGROUND: Bacillus anthracis is known to have low genetic variability. In spite of this lack of diversity, multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) and single nucleotide polymorphisms (SNPs) including the canonical SNPs assay (canSNP) have proved to be highly effective to differentiate strains. Five different MLVA schemes based on a collection of 31 VNTR loci (MLVA8, MLVA15, MLVA20, MLVA25 and MLVA31) with increased resolving power have been described. RESULTS: MLVA31 was applied to characterize the French National Reference Laboratory collection. The total collection of 130 strains is resolved in 35 genotypes. The 119 veterinary and environmental strains collection in France were resolved into 26 genotypes belonging to three canSNP lineages and four MLVA clonal complexes (CCs) with particular geographical clustering. A subset of seven loci (MLVA7) is proposed to constitute a first line assay. The loci are compatible with moderate resolution equipment such as agarose gel electrophoresis and show a good congruence value with MLVA31. The associated MLVA and SNP data was imported together with published genotyping data by taking advantage of major enhancements to the MLVAbank software and web site. CONCLUSIONS: The present report provides a wide coverage of the genetic diversity of naturally occurring B. anthracis strains in France as can be revealed by MLVA. The data obtained suggests that once such coverage is achieved, it becomes possible to devise optimized first-line MLVA assays comprising a sufficiently low number of loci to be typed either in one multiplex PCR or on agarose gels. Such a selection of seven loci is proposed here, and future similar investigations in additional countries will indicate to which extend the same selection can be used worldwide as a common minimum set. It is hoped that this approach will contribute to an efficient and low-cost routine surveillance of important pathogens for biosecurity such as B. anthracis.

High-throughput sequencing of Bacillus anthracis in France: investigating genome diversity and population structure using whole-genome SNP discovery.

BACKGROUND: Single nucleotide polymorphisms (SNPs) are ideal signatures for subtyping monomorphic pathogens such as Bacillus anthracis. Here we report the use of next-generation sequencing technology to investigate the historical, geographic and genetic diversity of Bacillus anthracis in France. 122 strains isolated over a 60-years period throughout the country were whole-genome sequenced and comparative analyses were carried out with a focus on SNPs discovery to discriminate regional sub-groups of strains. RESULTS: A total of 1581 chromosomal SNPs precisely establish the phylogenetic relationships existing between the French strains. Phylogeography patterns within the three canSNP sub-lineages present in France (i.e. B.Br.CNEVA, A.Br.011/009 and A.Br.001/002) were observed. One of the more remarkable findings was the identification of a variety of genotypes within the A.Br.011/009 sub-group that are persisting in the different regions of France. The 560 SNPs defining the A.Br.011/009- affiliated French strains split the Trans-Eurasian sub-group into six distinct branches without any intermediate nodes. Distinct sub-branches, with some geographic clustering, were resolved. The 345 SNPs defining the major B.Br CNEVA sub-lineage clustered three main phylogeographic clades, the Alps, the Pyrenees, and the Massif Central, with a small Saône-et-Loire sub-cluster nested within the latter group. The French strains affiliated to the minor A.Br.001/002 group were characterized by 226 SNPs. All recent isolates collected from the Doubs department were closely related. Identification of SNPs from whole-genome sequences facilitates high-resolution strain tracking and provides the level of discrimination required for outbreak investigations. Eight diagnostic SNPs, representative of the main French-specific phylogeographic clusters, were therefore selected and developed into high-resolution melting SNP discriminative assays. CONCLUSIONS: This work has established one of the most accurate phylogenetic reconstruction of B. anthracis population structure in a country. An extensive next-generation sequencing (NGS) dataset of 122 French strains have been created that allowed the identification of novel diagnostic SNPs useful to rapidly determine the geographic origin of any strain found in France.

In silico and in vitro evaluation of PCR-based assays for the detection of Bacillus anthracis chromosomal signature sequences.

Bacillus anthracis, the causative agent of anthrax, is a zoonotic pathogen that is relatively common throughout the world and may cause life threatening diseases in animals and humans. There are many PCR-based assays in use for the detection of B. anthracis. While most of the developed assays rely on unique markers present on virulence plasmids pXO1 and pXO2, relatively few assays incorporate chromosomal DNA markers due to the close relatedness of B. anthracis to the B. cereus group strains. For the detection of chromosomal DNA, different genes have been used, such as BA813, rpoB, gyrA, plcR, S-layer, and prophage-lambda. Following a review of the literature, an in silico analysis of all signature sequences reported for identification of B. anthracis was conducted. Published primer and probe sequences were compared for specificity against 134 available Bacillus spp. genomes. Although many of the chromosomal targets evaluated are claimed to be specific to B. anthracis, cross-reactions with closely related B. cereus and B. thuringiensis strains were often observed. Of the 35 investigated PCR assays, only 4 were 100% specific for the B. anthracis chromosome. An interlaboratory ring trial among five European laboratories was then performed to evaluate six assays, including the WHO recommended procedures, using a collection of 90 Bacillus strains. Three assays performed adequately, yielding no false positive or negative results. All three assays target chromosomal markers located within the lambdaBa03 prophage region (PL3, BA5345, and BA5357). Detection limit was further assessed for one of these highly specific assays.

Genetic diversity of Bacillus anthracis in Europe: genotyping methods in forensic and epidemiologic investigations.

Bacillus anthracis, the etiological agent of anthrax, a zoonosis relatively common throughout the world, can be used as an agent of bioterrorism. In naturally occurring outbreaks and in criminal release of this pathogen, a fast and accurate diagnosis is crucial to an effective response. Microbiological forensics and epidemiologic investigations increasingly rely on molecular markers, such as polymorphisms in DNA sequence, to obtain reliable information regarding the identification or source of a suspicious strain. Over the past decade, significant research efforts have been undertaken to develop genotyping methods with increased power to differentiate B. anthracis strains. A growing number of DNA signatures have been identified and used to survey B. anthracis diversity in nature, leading to rapid advances in our understanding of the global population of this pathogen. This article provides an overview of the different phylogenetic subgroups distributed across the world, with a particular focus on Europe. Updated information on the anthrax situation in Europe is reported. A brief description of some of the work in progress in the work package 5.1 of the AniBioThreat project is also presented, including (1) the development of a robust typing tool based on a suspension array technology and multiplexed single nucleotide polymorphisms scoring and (2) the typing of a collection of DNA from European isolates exchanged between the partners of the project. The know-how acquired will contribute to improving the EU's ability to react rapidly when the identity and real origin of a strain need to be established.

Use of high-resolution melting and melting temperature-shift assays for specific detection and identification of Bacillus anthracis based on single nucleotide discrimination.

Single nucleotide polymorphisms (SNPs) are important diagnostic markers for the detection and differentiation of Bacillus anthracis. High-Resolution Melting (HRM) and Melting Temperature (Tm)-shift methods are two approaches that enable SNP detection without the need for expensive labeled probes. We evaluated the potential diagnostic capability of those methods to discriminate B. anthracis from the other members of the B. cereus group. Two assays targeting B. anthracis-specific SNPs in the plcR and gyrA genes were designed for each method and used to genotype a panel of 155 Bacilli strains. All B. anthracis isolates (n=65) were correctly and unambiguously identified. Assays also proved to be appropriate for the direct genotyping of biological samples. They could reliably detect B. anthracis in contaminated organs containing as little as 10(3)CFU/ml, corresponding to a few genome equivalents per reaction. The HRM and Tm-shift applications described here represent valuable tools for specific identification of B. anthracis at reduced cost.

Tool for quantification of staphylococcal enterotoxin gene expression in cheese.

Cheese is a complex and dynamic microbial ecosystem characterized by the presence of a large variety of bacteria, yeasts, and molds. Some microorganisms, including species of lactobacilli or lactococci, are known to contribute to the organoleptic quality of cheeses, whereas the presence of other microorganisms may lead to spoilage or constitute a health risk. Staphylococcus aureus is recognized worldwide as an important food-borne pathogen, owing to the production of enterotoxins in food matrices. In order to study enterotoxin gene expression during cheese manufacture, we developed an efficient procedure to recover total RNA from cheese and applied a robust strategy to study gene expression by reverse transcription-quantitative PCR (RT-qPCR). This method yielded pure preparations of undegraded RNA suitable for RT-qPCR. To normalize RT-qPCR data, expression of 10 potential reference genes was investigated during S. aureus growth in milk and in cheese. The three most stably expressed reference genes during cheese manufacture were ftsZ, pta, and gyrB, and these were used as internal controls for RT-qPCR of the genes sea and sed, encoding staphylococcal enterotoxins A and D, respectively. Expression of these staphylococcal enterotoxin genes was monitored during the first 72 h of the cheese-making process, and mRNA data were correlated with enterotoxin production.

Differential temporal expression of the staphylococcal enterotoxins genes during cell growth.

Staphylococcal enterotoxins (SEs) are a large family of structurally related superantigens produced by Staphylococcus aureus and responsible for staphylococcal food-poisoning (SFP). To better understand how the se genes are expressed, and especially the newly described ones (seg to ser, seu), a quantitative reverse transcription-polymerase chain reaction was developed and used to monitor their expression in a panel of 28 enterotoxigenic staphylococci including seven reference strains and 21 isolates collected from foods and SFP outbreaks in France. Kinetic mRNA studies revealed four distinct patterns of expression according to the enterotoxins genes analysed corresponding either to: (i) unchanged mRNAs abundance during bacterial growth (sea, see, sej, sek, seq and sep); (ii) slight decrease in transcript levels (seg, sei, sem, sen, seo, seu); (iii) drastic induction of expression at the end of the exponential growth phase (seb, sec, seh), or (iv) modest post-exponential increase in mRNAs level (<10-fold) (sed, ser, sel). The study demonstrates that all se containing strains are potentially able to produce SEs and that most of the newly described se genes are likely not controlled by the agr system. The rapid assessment of se transcripts levels by RT-qPCR might provide valuable clue to evaluate the poisoning risk linked to a strain.

Regulatory role of UvrY in adaptation of Photorhabdus luminescens growth inside the insect.

We report global expression profiling of a uvrY-deficient mutant of Photorhabdus luminescens. We found that the regulator moiety of the two-component regulatory system BarA/UvrY regulated more than 500 target genes coding for functions involved in the synthesis of major compartments and metabolic pathways of the cell. This regulation appeared to be in part indirect as UvrY affected the expression of several regulators. Indeed, the flagellum biosynthesis transcription activator FlhC and the flagella regulon were induced in the absence of UvrY, leading to a hyperflagellated phenotype and an increase in motility and biofilm formation. Two major regulatory systems were also altered: the type 2 quorum-sensing inducer AI-2 was activated by UvrY, and the CsrA regulator function appeared to be repressed by the increase of the small-untranslated RNA csrB, the CsrA activity inhibitor TldD and the chaperonin GroESL. Both through and independently of these systems, UvrY regulated oxidative stress resistance; bioluminescence; iron, sugar and peptide transport; proteases; polyketide synthesis enzymes and nucleobases recycling, related to insect degradation and assimilation by bacteria. As a consequence, the uvrY-deficient strain exhibited a decreased killing of insect cells and a reduced growth on insect cells culture, suggesting a UvrY role in the adaptation of P. luminescens inside the insect.

Comparison of three chromogenic media and evaluation of two molecular-based identification systems for the detection of Enterobacter sakazakii from environmental samples from infant formulae factories.

Enterobacter sakazakii is an occasional contaminant of powdered infant formula that can cause rare but severe foodborne infections in infants. To determine optimal methods for the detection and identification of E. sakazakii, 38 naturally contaminated samples from infant formulae factories were analyzed by two PCR-based methods and by a method (TS 22964/RM 210) developed by the International Organization for Standardization and the International Dairy Federation (ISO-IDF) using three different commercial chromogenic agars. The ISO-IDF method includes two enrichment steps, plating of the second enrichment broth on E. sakazakii isolation agar (a chromogenic selective agar), picking of five typical colonies for transfer onto tryptone soy agar, and subsequent confirmation of yellow-pigmented colonies by biochemical characterization. Twenty-two of the 38 samples were positive by the culture method. E. sakazakii isolation agar (ESIA; AES Laboratoires), COMPASS agar (Biokar Diagnostics), and Druggan-Forsythe-Iversen agar (Oxoid) compared favorably with violet red bile glucose agar (VRBG, a selective medium for Enterobacteriaceae), with positive predictive values of 86.96, 88, and 74.07%, respectively, in contrast to 47.83% for VRBG. One additional positive sample was detected using the nonpatented real-time PCR method evaluated, and those results were in 97.3% concordance with the ISO-IDF results. Some discrepancies between the results of the DuPont Qualicon BAX system and those of the ISO-IDF method could be explained by heterogeneity of contamination and sampling. Thus, both PCR-based systems were suitable for detecting and specifically identifying E. sakazakii within 1 to 2 days, and COMPASS agar and ESIA could be used interchangeably as a first-step medium to isolate presumptive E. sakazakii colonies.

A robotic DNA purification protocol and real-time PCR for the detection of Enterobacter sakazakii in powdered infant formulae.

BACKGROUND: Enterobacter sakazakii is the causative agent of rare but severe food-borne infections associated with meningitis, necrotizing enterocolitis and sepsis in infants. Rehydrated powdered infant formulae have been implicated as the source of infection in several outbreaks and sporadic cases. In this work, a real time fluorescence resonance energy transfer PCR assay incorporating an internal amplification control (IAC) was developed for the specific detection of E. sakazakii in foods. Performance of the assay, coupled to an automated DNA extraction system and the E. sakazakii ISO-IDF (TS 22964/RM 210) enrichment procedure, was evaluated on infant formulae and samples from production environment. RESULTS: The real-time PCR assay had 100% specificity as assessed using 35 E. sakazakii and 184 non-E. sakazakii strains. According to the E. sakazakii strains tested, the detection limits ranged from 5 to 25 genomic copies. Assays on pure cultures (including real-time PCR and DNA extraction) gave a sensitivity of about 102 to 103 CFU/ml. Out of 41 naturally contaminated infant formulae and environmental samples analysed for the presence of E. sakazakii, 23 were positive by real-time PCR and 22 by the conventional culture method, giving 97.5% concordance with the ISO-IDF reference method. CONCLUSION: This method, combining specific real-time PCR, automated DNA extraction and ISO-IDF standard enrichments, provides a useful tool for rapid screening of E. sakazakii in food and environmental matrices.

Pleiotropic role of quorum-sensing autoinducer 2 in Photorhabdus luminescens.

Bacterial virulence is an integrative process that may involve quorum sensing. In this work, we compared by global expression profiling the wild-type entomopathogenic Photorhabdus luminescens subsp. laumondii TT01 to a luxS-deficient mutant unable to synthesize the type 2 quorum-sensing inducer AI-2. AI-2 was shown to regulate more than 300 targets involved in most compartments and metabolic pathways of the cell. AI-2 is located high in the hierarchy, as it controls the expression of several transcriptional regulators. The regulatory effect of AI-2 appeared to be dose dependent. The luxS-deficient strain exhibited decreased biofilm formation and increased type IV/V pilus-dependent twitching motility. AI-2 activated its own synthesis and transport. It also modulated bioluminescence by regulating the synthesis of spermidine. AI-2 was further shown to increase oxidative stress resistance, which is necessary to overcome part of the innate immune response of the host insect involving reactive oxygen species. Finally, we showed that the luxS-deficient strain had attenuated virulence against the lepidopteran Spodoptera littoralis. We concluded that AI-2 is involved mainly in early steps of insect invasion in P. luminescens.

Proteome analysis of the phenotypic variation process in Photorhabdus luminescens.

Photorhabdus luminescens is an insect pathogen associated with specific soil nematodes. The bacterium has a complex life cycle with a symbiotic stage in which bacteria colonize the intestinal tract of the nematodes, and a pathogenic stage against susceptible larval-stage insect. Symbiosis-"deficient" phenotypic variants (known as secondary forms) arise during prolonged incubation. Correspondence analysis of the in silico proteome translated from the genome sequence of strain TT01 identified two major biases in the amino acid composition of the proteins. We analyzed the proteome, separating three classes of extracts: cellular, extracellular, and membrane-associated proteins, resolved by 2-DE. Approximately 450 spots matching the translation products of 231 different coding DNA sequences were identified by PMF. A comparative analysis was performed to characterize the protein content of both variants. Differences were evident during stationary growth phase. Very few proteins were found in variant II supernatants, and numerous proteins were lacking in the membrane-associated fraction. Proteins up-regulated by the phenotypic variation phenomenon were involved in oxidative stress, energy metabolism, and translation. The transport and binding of iron, sugars and amino acids were also affected and molecular chaperones were strongly down-regulated. A potential role for H-NS in phenotypic variation control is discussed.

Proteome analysis of Streptococcus thermophilus grown in milk reveals pyruvate formate-lyase as the major upregulated protein.

We investigated the adaptation to milk of Streptococcus thermophilus LMG18311 using a proteomic approach. Two-dimensional electrophoresis of cytosolic proteins were performed after growth in M17 medium or in milk. A major modification of the proteome concerned proteins involved in the supply of amino acids, like the peptidase PepX, and several enzymes involved in amino acid biosynthesis. In parallel, we observed the upregulation of the synthesis of seven enzymes directly involved in the synthesis of purines, as well as formyl-tetrahydrofolate (THF) synthetase and serine hydroxy-methyl transferase, two enzymes responsible for the synthesis of compounds (THF and glycine, respectively) feeding the purine biosynthetic pathway. The analysis also revealed a massive increase in the synthesis of pyruvate formate-lyase (PFL), the enzyme which converts pyruvate into acetyl coenzyme A and formate. PFL has been essentially studied for its role in mixed-acid product formation in lactic acid bacteria during anaerobic fermentation. However, formate is an important methyl group donor for anabolic pathway through the formation of folate derivates. We hypothesized that PFL was involved in purine biosynthesis during growth in milk. We showed that PFL expression was regulated at the transcriptional level and that pfl transcription occurred during the exponential growth phase in milk. The complementation of milk with formate or purine bases was shown to reduce pfl expression, to suppress PFL synthesis, and to stimulate growth of S. thermophilus. These results show a novel regulatory mechanism controlling the synthesis of PFL and suggest an unrecognized physiological role for PFL as a formate supplier for anabolic purposes.