Comparative genomic analysis of Brevibacterium strains: insights into key genetic determinants involved in adaptation to the cheese habitat
BMC Genomics volume 18, Article number: 955 (2017)
Brevibacterium strains are widely used for the manufacturing of surface-ripened cheeses, contributing to the breakdown of lipids and proteins and producing volatile sulfur compounds and red-orange pigments. The objective of the present study was to perform comparative genomic analyses in order to better understand the mechanisms involved in their ability to grow on the cheese surface and the differences between the strains.
The genomes of 23 Brevibacterium strains, including twelve strains isolated from cheeses, were compared for their gene repertoire involved in salt tolerance, iron acquisition, bacteriocin production and the ability to use the energy compounds present in cheeses. All or almost all the genomes encode the enzymes involved in ethanol, acetate, lactate, 4-aminobutyrate and glycerol catabolism, and in the synthesis of the osmoprotectants ectoine, glycine-betaine and trehalose. Most of the genomes contain two contiguous genes encoding extracellular proteases, one of which was previously characterized for its activity on caseins. Genes encoding a secreted triacylglycerol lipase or involved in the catabolism of galactose and D-galactonate or in the synthesis of a hydroxamate-type siderophore are present in part of the genomes. Numerous Fe3+/siderophore ABC transport components are present, part of them resulting from horizontal gene transfers. Two cheese-associated strains have also acquired catecholate-type siderophore biosynthesis gene clusters by horizontal gene transfer. Predicted bacteriocin biosynthesis genes are present in most of the strains, and one of the corresponding gene clusters is located in a probable conjugative transposon that was only found in cheese-associated strains.
Brevibacterium strains show differences in their gene repertoire potentially involved in the ability to grow on the cheese surface. Part of these differences can be explained by different phylogenetic positions or by horizontal gene transfer events. Some of the distinguishing features concern biotic interactions with other strains such as the secretion of proteases and triacylglycerol lipases, and competition for iron or bacteriocin production. In the future, it would be interesting to take the properties deduced from genomic analyses into account in order to improve the screening and selection of Brevibacterium strains, and their association with other ripening culture components.
Microbial communities from rinds of surface-ripened cheeses are composed of various bacteria, yeasts and molds, which contribute to the flavor, texture and appearance of the final products. These microorganisms may come from the milk, the ripening environment or from ripening cultures that are widely used in the cheese industry. The function of the ripening cultures is to provide specific organoleptic properties, to ensure a better regularity of manufacturing, and to outcompete pathogens or spoilage microorganisms . However, strains from ripening cultures frequently do not establish themselves in cheeses . Even if they are massively inoculated, these strains are sometimes outcompeted by the resident “house flora” due to insufficient fitness in the cheese surface habitat. The ability to grow on the cheese surface depends on various properties such as efficient salt tolerance and iron acquisition systems, or on the ability to use the energy compounds present in the cheese . In addition, growth is influenced by the other microorganisms present at the cheese surface, with which they may have positive or negative interactions.
One example of a ripening culture component that may have problematic growth in cheeses is Brevibacterium [4,5,6,7,8,9]. This microorganism contributes to the breakdown of lipids and proteins, and produces volatile sulfur compounds that are key aroma impact compounds, as well as red-orange pigments [10,11,12]. For a long time, Brevibacterium linens was considered to be the major Brevibacterium species in cheeses. In 2004, it was broken down into three species: B. linens, B. antiquum and B. aurantiacum . These three species, together with B. casei and other not-yet described Brevibacterium species, have been isolated from cheeses .
In order to improve the strategies for selecting Brevibacterium strains for ripening cultures, it is important to better understand the mechanisms involved in their ability to grow on the cheese surface. This can be investigated by genomic analyses. For example, the study of the genomes of the cheese strains Glutamicibacter arilaitensis Re117 (formerly Arthrobacter arilaitensis Re117) and Corynebacterium variabile DSM 44702 revealed several metabolic capabilities that were considered to play roles in growth on cheese [15, 16]. In addition, a recent study provided evidence of extensive horizontal gene transfer (HGT) in cheese-associated bacteria, including Brevibacterium strains, for which genes involved in iron acquisition were particularly abundant in the transferred islands .
The aim of the present study was to investigate, in Brevibacterium strains, key genetic determinants known to be important for growth in cheese: the catabolism of energy compounds present in cheeses, iron acquisition, salt tolerance and bacteriocin production. For that purpose, we sequenced the genome of 13 Brevibacterium strains, including eleven strains isolated from cheeses. We performed comparative analyses of these genomes and of ten other Brevibacterium genomes from strains isolated from diverse environments and already present in the Integrated Microbial Genomes (IMG) database .
Growth conditions and DNA extraction
The Brevibacterium strains were grown under aerobic conditions (rotary shaker at 150 rpm) for three days at 25 °C in 50-ml conical flasks containing 10 ml of brain heart infusion broth (Biokar Diagnostics, Beauvais, France). Bacterial cells were recovered by centrifugation at 4500 x g for 10 min from 5 ml of culture, washed once with 5 ml of TE buffer (Tris-HCl 10 mM, EDTA 1 mM, pH 8) and resuspended in 500 μl of the same buffer. Seventy-five μl of lysis solution containing lysozyme (40 mg/ml) and lyticase (1333 U/ml) were added and the suspensions were incubated for 30 min at 37 °C. After addition of 60 μl of 0.5 M EDTA pH 8.0, 20 μl of proteinase K (20 mg/ml) and 100 μl of 20% SDS, the samples were incubated for 1 h at 55 °C and subsequently transferred to 2-ml bead-beating tubes containing 100 mg of 0.1 mm-diameter zirconium beads and 100 mg of 0.5 mm-diameter zirconium beads. After cooling on ice, 500 μl of phenol-chloroform-isoamyl alcohol (25:24:1; saturated with 10 mM Tris, pH 8.0, and 1 mM EDTA) were added and the tubes were shaken in a bead beater (FastPrep-24, MP Biomedicals, Illkirch, France) using two 45-s mixing sequences at a speed of 6.0 m/s. The tubes were cooled on ice for 5 min after each mixing. The content of the tubes was transferred to Phase Lock Gel Heavy tubes (5 PRIME, Hilden, Germany), which were then centrifuged at 18,500 × g for 30 min at 20 °C. The aqueous phases were subsequently transferred to new Phase Lock Gel tubes. After adding 500 μl of phenol-chloroform-isoamyl alcohol and gentle mixing, centrifugation was performed at 18,500 × g for 20 min at 20 °C. Five hundred μl of chloroform were then added to the aqueous phase, and the tubes were centrifuged at 18,500 × g for 20 min at 20 °C after gently mixing. The aqueous phases (approximately 200 μl) were recovered, mixed with 1 μl of RNase A (20 mg/ml), and incubated for 30 min at 37 °C. DNA was then precipitated by adding 200 μl of cold isopropanol and the tubes were incubated for 10 min at 4 °C. The DNA was recovered by centrifugation for 15 min at 18,500 x g and 4 °C, and the pellets were subsequently washed two times with 500 μl of 70% (vol/vol) ethanol. The pellets were then dried for 15 min at 42 °C and dissolved in 100 μl of water. Finally, DNA cleanup was performed using the DNeasy Blood & Tissue Kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s instructions. DNA yield and purity (absorbance ratio at 260/280 nm) were determined using a NanoDrop ND-1000 spectrophotometer (Labtech, Palaiseau, France). DNA integrity was verified by electrophoresis in a 1% agarose gel (Qbiogene, Illkirch, France) in 1X TAE buffer (40 mM Tris, 20 mM acetic acid, 1 mM EDTA, pH 8.3) stained with SYBR® Safe 1X (Invitrogen, Carlsbad, CA, USA).
Genome sequencing, assembly and annotation
Libraries were generated using the TruSeq DNA Sample Preparation Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Sequencing was carried out on an Illumina MiSeq apparatus at the INRA GeT-PlaGE platform (http://get.genotoul.fr) in order to generate paired-end reads (250 bases in length). For each strain, the paired-end reads were merged using FLASH . De novo assembly was performed using SPAdes version 3.1.1 . Only contigs with length > 1000 bp were considered for further study. Gene predictions and annotations were performed automatically using the Integrated Microbial Genomes (IMG) database and comparative analysis system , as described in the corresponding standard operating procedure . PHASTER  was used to predict prophages. Sequences involved in the production of secondary metabolites and bacteriocins were searched using antiSMASH 3.0 , BAGEL3  and BACTIBASE .
Comparative genomic analysis and phylogenetic classification
Comparative genomic analyses were performed considering the 13 genomes sequenced in the present study as well as the ten other Brevibacterium genomes present in the IMG database in September 2016 (Table 1). A phylogenetic analysis was performed for the 23 Brevibacterium genomes using the sequences of 40 marker genes, as described by Mende et al. . The genomes of Glutamicibacter arilaitensis Re117  (Project accession number PRJEA50353) and Corynebacterium casei LMG S-19264  (Project accession number PRJNA186910) were used as outgroups. An in-house database of the 40 marker genes present in 388 bacterial strains, which included 117 strains isolated from dairy products , was used to detect the 40 marker genes in the Brevibacterium genomes using tBlastN . The best hits were selected with at least 60% sequence identity and 80% coverage. Each marker gene was translated into an amino acid sequence using T-Coffee  and aligned using MUSCLE . The 40 individual alignments were then concatenated to a single one, which was used to build the tree using FastTree 2  with the following parameters: –gamma –pseudo –mlacc 3 –slownni and the default bootstrap procedure (1000 resamples). The tree was visualized and annotated using MEGA7 .
Homologous gene families were computed using the OrthoMCL procedure implemented in GET-HOMOLOGUES software . Amino acid sequences of the CDSs from the 23 Brevibacterium genomes were grouped into clusters using 75% identity and 75% coverage thresholds with a BlastP cutoff E value <1e-05. The inflation index of OrthoMCL algorithm was set to 1.5, as recommended by Li et al. . Strain ATCC 9175T was set as the first reference genome, the subsequent ones were randomly chosen by GET-HOMOLOGUES. Functional category assignment to each cluster was done according to the Clusters of Orthologous Groups (COG) database . The Average Nucleotide Identity (ANI), implemented in GET-HOMOLOGUES software, was computed on homologous genes for all possible pairs of genomes.
General genomic features
Assembly and annotation metrics for the 13 newly sequenced Brevibacterium genomes and those concerning the ten additional genomes (which included one cheese isolate) for the comparative analysis are detailed in Additional file 1. Brevibacterium genomes show considerable size heterogeneity, ranging from about 2.3 to 4.5 Mbp. There is no clear relationship between genome size and the habitat from which the strains were isolated. However, genome size is quite similar among the 12 cheese isolates, 4 Mbp on average (min: 3.7; max: 4.5), and these genomes harbor an average of 3712 genes (min: 3395; max: 4154). The two smallest genomes corresponded to human-associated strains: B. ravenspurgense 5401308T and B. mcbrellneri ATCC 49030T (2.3 and 2.6 Mbp, respectively). The G + C content varied from 58.0% (B. mcbrellneri ATCC49030T) to 70.9% (B. album DSM 18261T), and was between 62 and 65% for most of the strains. Clustered regularly interspaced short palindromic repeat (CRISPR) candidates were found in most genomes (19 out of 23). Complete CRISPR-Cas systems consist of an array of CRISPRs interspaced by spacers and an adjacent cas gene cluster. Such a complete structure was observed only in the genomes of B. casei S18 (scaffold ID: S272_Contig15.15; 17 spacers), B. album DSM 18261T (scaffold ID: K318DRAFT_scaffold00001.1; 2 + 1 + 6 + 35 + 12 spacers) and B. ravenspurgense 5401308T (scaffold ID: Y1ADRAFT_CAJD01000011_1.11; 36 spacers). They all belong to the Type I CRISPR-Cas system but their overall gene content is variable. The PHASTER tool identified only one complete prophage region, in strain B. sp. Mu109. This 11.4 kb region (scaffold ID: Ga0063700_1029; locus tag Ga0063700_02875 to Ga0063700_02886) contains attachment sites (attL and attR) and CDSs encoding putative transposase, recombinase, integrase and capsid scaffolding proteins. However, because most genomes are draft genomes, other prophage regions might be present as regions split on several contigs and, consequently, difficult to detect.
Phylogenomic analyses and orthology
The genomic data from the 23 Brevibacterium strains were used to assess their intra-genus phylogenetic relationships. The phylogenomic tree partitioned the strains into two major lineages (Fig. 1). Lineage 1 contains two Brevibacterium strains isolated from human-associated samples, i.e., B. mcbrellneri ATCC49030T and B. ravenspurgense 5401308T. Lineage 2, containing the 21 other strains, is composed of two branches corresponding to the groups 2.A and 2.B. Group 2.A contains four strains isolated from diverse habitats: B. album DSM 18261T (saline soil), B. senegalense JC43T (human stool), B. sp. Mu109 (cheese) and B. jeotgali SJ5-8T (fermented seafood). Group 2.B contains all the Brevibacterium strains isolated from cheese except B. sp. Mu109. It can be divided into three clades: 2.B.1, 2.B.2 and 2.B.3. Clade 2.B.1 contains nine strains, including eight strains isolated from cheese, and B. sandarakinum DSM 22082T, which was isolated from a house wall. Clade 2.B.2 consists of six strains, including two from cheeses and three from freshwater environments. Clade 2.B.3 contains two strains belonging to the B. casei species, i.e., CIP 102111T, isolated from cheese and S18, isolated from human skin. In the following, clades 2.B.1, 2.B.2, 2.B.3 will be referred to as group “aurantiacum/sandarakinum/antiquum”, group “linens/siliguriense/iodinum” and group “casei”, respectively. The taxonomy of the Brevibacterium genus is still in process of reclassification . According to Stackebrandt et al. , 16S rRNA (rrs) gene sequences and DNA-DNA hybridization (DDH) should be considered as molecular criteria for species delineation. However, the Average Nucleotide Identity (ANI) has recently been proposed to replace DDH values [39,40,41]. In the present study, we used two criteria for delineating species: two strains were considered to belong to the same species if (i) their rrs gene sequence identity was ≥98%  and (ii) their ANI was ≥95% . The results are presented in Additional file 2. Based on these criteria, two strains isolated from cheese products, i.e. 239c and Mu109 could not be linked to any sequenced species. In addition, B. linens AE038–8 and B. linens ATCC 9172T have an ANI value of only 91.2%, suggesting that strain AE038–8 does not belong to the B. linens species.
Pan-genome analysis of the 23 Brevibacterium strains resulted in 25,376 orthologous gene clusters, containing a total of 78,702 protein-coding genes (Additional files 3 and 4). Of these clusters, only 263 (1%) are shared by all the strains (the core genome), reflecting a high intragenus genomic variability. We further investigated the B. aurantiacum species, in which all five of the sequenced strains were isolated from cheese. The pan-genome of this group contains 5988 orthologous gene clusters, of which 2684 are core and 3304 are accessory, i.e., variable among strains (Fig. 2A). Functional prediction revealed an overrepresentation of some COG categories in the accessory genome relative to the core genome, especially the categories [V] (Defense mechanisms) and [X] (Mobilome: prophages, transposons), in which the ratio accessory/core was 4.4 and 45.3, respectively (Fig. 2B).
Catabolism of energy compounds present in cheeses
Catabolism of lactose, galactose and D-galactonate
During the manufacturing of cheeses, lactose is consumed by lactic acid bacteria, but some lactose may still be present at the beginning of ripening. When the lactic starter culture contains Streptococcus thermophilus strains, some galactose is produced from lactose, and this compound may be present in the cheese curd for several weeks . No beta-galactosidase gene was identified in the 23 Brevibacterium genomes, which is consistent with the fact that most Brevibacterium strains are not able to consume lactose . However, the genomes of twelve of the 23 Brevibacterium strains encode the four enzymes of the Leloir pathway for galactose utilization (Additional file 5). The corresponding strains belong to the phylogenetic groups aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum. A complete pathway for D-galactonate catabolism is present in strains from the phylogenetic groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei. The genes are organized in a cluster, that encodes a transcriptional regulator, a 2-dehydro-3-deoxyphosphogalactonate aldolase [EC 22.214.171.124], a 2-dehydro-3-deoxygalactonokinase [EC 126.96.36.199], a galactonate dehydratase [EC 188.8.131.52] and one or two D-galactonate importers.
Catabolism of lactate, acetate, ethanol and citrate
Lactate, which is produced from lactose by lactic acid bacteria, is an important energy substrate for most aerobic cheese ripening bacteria . At least two predicted lactate permease genes are present in the genomes of the 23 Brevibacterium strains (Additional file 5). The oxidation of lactate by lactate dehydrogenase generates pyruvate that is later catabolized through the TCA cycle. There are two types of lactate dehydrogenases in bacteria, NAD-dependant lactate dehydrogenases (nLDHs) and NAD-independent lactate dehydrogenases (iLDHs), the latter generally being considered as the enzymes responsible for lactate oxidation in bacteria . Four types of iLDHs were identified in the genomes of the Brevibacterium strains: a quinone-dependent D-iLDH, a quinone-dependent L-iLDH, a Dld II-type quinone or cytochrome-dependent D-iLDH  and a three subunit quinone or cytochrome-dependent L-iLDH complex LldEFG . Dld II-type quinone or cytochrome-dependent D-iLDHs were present only in group aurantiacum/sandarakinum/antiquum. A lactate permease gene is adjacent to the quinone-dependent L-iLDH gene in 14 strains. In cheese, acetate is produced from lactose by heterofermentative lactic acid bacteria, from citrate by citrate-utilizing lactic acid bacteria, or from lactate by Pediococcus and Propionibacterium strains. Twenty of the 23 Brevibacterium genomes encoded the monocarboxylic acid transporter MctC, which is an uptake system for acetate, propionate and pyruvate  (Additional file 5). Genes encoding acetyl-CoA synthase, which channels acetate toward the TCA cycle, are present in all the 23 strains. In cheese, ethanol is produced from lactose by yeasts such as Kluyveromyces lactis and K. marxianus, and by heterofermentative lactic acid bacteria. It is a potential energy substrate for aerobic microorganisms during cheese ripening but, to our knowledge, this has never been investigated. There are numerous candidate genes encoding enzymes involved in the catabolism of ethanol to acetate in the genomes of the Brevibacterium strains: between one and nine for the alcohol dehydrogenase and between nine and 31 genes for the acetaldehyde dehydrogenase (Additional file 5). In Corynebacterium glutamicum, the genes encoding the alcohol dehydrogenase (adhA) and the acetaldehyde dehydrogenase (ald) are responsible for ethanol catabolism . All of the 23 Brevibacterium genomes encode an ortholog of adhA (65 to 69% identity at the amino acid level) and an ortholog of ald (67 to 82% identity), and in 21 strains, these genes form a cluster. Milk contains citrate (~8 mM in cow’s milk), which may persist during cheese ripening if the lactic starter culture does not contain citrate-utilizing lactic acid bacteria. Citrate importers belonging to the Citrate-Mg2+:H+ (CitM)/Citrate-Ca2+:H+ (CitH) Symporter (CitMHS) family (TC no. 2.A.11) and to the 2-HydroxyCarboxylate Transporter (2-HCT) family (TC no. 2.A.24) were detected in 16 Brevibacterium strains (Additional file 5). The genomes of the two B. casei strains did not encode any citrate transporter, and those of the six strains from the phylogenetic group linens/siliguriense/iodinum did not encode 2-HCT family citrate transporters.
Catabolism of lipids and glycerol
Milk contains large amounts of triglycerides (~35 g/l in cow’s milk), which can be used by cheese microorganisms as an energy source. Lipid catabolism involves the release of free fatty acids and glycerol and the subsequent breakdown of these compounds. The Brevibacterium genomes encode between 8 and 37 proteins with putative lipase or esterase activity (Additional file 6). Two types of secreted lipases / esterases were identified. The first type is a triacylglycerol lipase [EC 184.108.40.206], that was identified in 15 genomes, corresponding to eight of the nine strains of the phylogenetic group aurantiacum/sandarakinum/antiquum, to all six strains of the group linens/siliguriense/iodinum and to only one of the eight other strains (the cheese-associated strain Mu109). The second type is a glycerophosphoryl diester phosphodiesterase [EC 220.127.116.11], which was identified only in the strains of the phylogenetic groups aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum. Uptake of fatty acids can be done by passive diffusion through the membrane lipid bilayer or by protein-facilitated transfer . Genomic analysis revealed the presence of short-chain fatty acid uptake (AtoE) family proteins (TC no. 2.A.73) in all of the 23 Brevibacterium strains and it is noteworthy that Brevibacterium strains have a greater number of transporters from this family than other Actinobacteria strains (1.7 vs. 0.1 genes per genome), based on assignations to the COG2031 for the 5726 Actinobacteria strains present in April 2017 in the IMG database. A complete beta-oxidation pathway for fatty acid degradation was identified in all of the investigated strains, except in strains ATCC 49030T and 5401308T, for which genes encoding L-3-hydroxyacyl-CoA dehydrogenase or 3-ketoacyl-CoA thiolase seemed to be lacking. Beta-oxidation of odd-chain-length fatty acids yields propionyl-CoA in addition to acetyl-CoA. Propionyl-CoA may be catabolized via the methylcitrate cycle, which oxidizes it to pyruvate. Three enzymes are characteristic of this cycle: the methylcitrate synthase, the methylcitrate dehydratase and the 2-methylisocitrate lyase, encoded by prpC, prpD and prpB, respectively . In all of the investigated Brevibacterium genomes, these genes form a cluster (prpDBC). It is noteworthy that in all of the strains of the phylogenetic groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei, prpDBC is located upstream of the genes involved in the catabolism of glycerol (glycerol uptake protein, glycerol kinase and glycerol-3-phosphate dehydrogenase). This glycerol pathway occurs in aerobic conditions because glycerol-3-phosphate dehydrogenase reduces quinones of the respiratory chain . The genome of the two B. casei strains also encodes a glycerol dehydrogenase and a dihydroxyacetone kinase, which constitute another glycerol utilization pathway.
Catabolism of proteins and amino acids
Cheese contains a large amount of proteins, mainly caseins, which can be degraded by various proteolytic microorganisms. The amino acids resulting from proteolysis can be used as an energy source by the cheese ripening microorganisms and they are also precursors of key flavor compounds. Genomic analysis of the 23 Brevibacterium genomes by the MEROPS peptidase BLAST search tool  and manual curation of the results revealed the presence of between four and 16 putative excreted enzymes with proteolytic activities, depending on the strain (Additional file 7). About 15% of these enzymes have a LPXTG motif (TIGR01167) at their C-terminus and are thus probably cell-wall-associated. Interestingly, in many cases the corresponding genes are contiguous. In addition, there was a complete identity between the sequence of one of the two predicted cell-wall-associated proteases of B. aurantiacum ATCC 9174 (locus tag BlinB01003410) and the first 20 N-terminal amino acid sequence of an extracellular protease purified from this strain . The presence of amino acid degradation pathways was inferred from the annotation of the 23 genomes (Additional file 7). All of them encode a bifunctional proline dehydrogenase/L-glutamate gamma-semialdehyde dehydrogenase [EC 18.104.22.168 / EC 22.214.171.124], which catalyzes oxidation of proline to glutamate using a membrane-bound quinone and NAD as the electron acceptor. They also encode a NAD-specific glutamate dehydrogenase [EC 126.96.36.199] that produces NADH and alpha-ketoglutarate, an intermediate of the TCA cycle, and the enzymes involved in the degradation of threonine and serine. The histidine and alanine catabolic pathways seem to be present in all of the strains except B. senegalense JC43T (histidine) and B. ravenspurgense 5401308T (alanine). Phenylalanine, tyrosine, methionine and arginine degradation pathways were identified only in the strains of the phylogenetic groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei. Gamma-aminobutyrate is a four-carbon non-protein amino acid produced from glutamate by lactic acid bacteria during the ripening of cheeses . All the genomes encode a 4-aminobutyrate transaminase [EC 188.8.131.52] and a succinate semialdehyde dehydrogenase [EC 184.108.40.206]. These enzymes convert Gamma-aminobutyrate into succinate, an intermediate of the TCA cycle.
Genes encoding putative Mn2+ or Fe2+ transporters are present in all of the 23 Brevibacterium genomes (Additional file 8). The EfeUOB transporter, which is a high-affinity uptake system for both Fe2+ and Fe3+ , is present in seven strains. Five genomes encode a putative ABC-type iron transport system whose closest homolog in Haemophilus influenza (FbpABC) is required for acquiring iron from transferrin . Fe3+/siderophore transport components are present in all of the 23 genomes, varying from nine components for B. senegalense JC43T, up to 31 components for B. sp. Mu109. Comparison of the abundance of these components based on the analysis of the COG1120, COG4604, COG0609, COG4605, COG4606, COG4779, COG0614, COG4592 and COG4607, showed that the average number was higher for Brevibacterium (17.6 genes per genome) than for the other Actinobacteria (8.4 genes per genome, calculated from 5726 Actinobacteria genomes). The Brevibacterium genomes also encode siderophore interacting proteins (mean number of 3.6 genes per genome), which are required for iron release from Fe3+/siderophore complexes. A putative hydroxamate-type siderophore biosynthesis cluster, which encodes a lysine N6-hydroxylase [EC 220.127.116.11], a siderophore synthetase component and, occasionally, a L-2,4-diaminobutyrate decarboxylase [EC 18.104.22.168], is present in 15 strains. All these strains belong to the phylogenetic groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei (Additional file 8). The genome of strain Mu109 contains a cluster (locus tag Ga0063700_02161 to Ga0063700_02193) encoding eight Fe3+/siderophore transport components, one siderophore interacting protein and four proteins that are probably involved in the biosynthesis of a catecholate or a mixed catecholate/hydroxamate siderophore: 4′-phosphopantetheinyl transferase EntD, MbtH protein, non-ribosomal siderophore peptide synthetase component and L-ornithine N5-oxygenase [EC 22.214.171.124]. The closest homologs of the 4′-phosphopantetheinyl transferase EntD and the non-ribosomal siderophore peptide synthetase component are found in Streptomyces species. The genome of strain ATCC 9174 contains a cluster (locus tag BlinB01002486 to BlinB01002496) encoding three Fe3+/siderophore transport components, one siderophore interacting protein, one siderophore exporter and seven genes that are probably involved in the biosynthesis of a catecholate siderophore: non-ribosomal siderophore peptide synthetase component, glycosyltransferase IroB, 2,3-dihydro-2,3-dihydroxybenzoate dehydrogenase [EC 126.96.36.199], isochorismate synthase [EC 188.8.131.52], 2,3-dihydroxybenzoate-AMP ligase [EC 184.108.40.206], isochorismatase [EC 220.127.116.11], and a putative transferase component of siderophore synthetase. These siderophore biosynthesis genes probably result from HGTs since many of their closest homologs are present, either in Gram-negative species or in Streptomyces or Paenibacillus species. In addition, comparison of the flanking regions of the siderophore biosynthesis cluster of strain ATCC 9174 to the CNRZ 920 genome (in which the cluster is absent) revealed that the cluster corresponded to an insertion that occurred in an ancestor of strain ATCC 9174 at the end of a tRNA-gly gene (locus tag BlinB_R0152 in ATCC 9174 and Ga0063691_00673 in CNRZ 920), which is followed in CNRZ 920 by a protein of unknown function (locus tag Ga0063691_00672, which is an ortholog of BlinB01002497 in ATCC 9174). For strains ATCC 9172T and Mu101, there is little evidence for siderophore biosynthesis, even if it cannot be excluded, since their genome encodes a lysine decarboxylase that is clustered with a siderophore interacting protein and a Fe3+/siderophore binding component. With the exception of isochorismate synthase, which is an enzyme that is also involved in menaquinone biosynthesis, no genes involved in siderophore biosynthesis were identified in the genomes of strains DSM 18261T, SJ5-8T, 5401308T, ATCC 49030T and JC43T. In summary, the genomic analyses indicate that 16 of the 23 Brevibacterium strains are probably able to produce siderophores. Hydroxamate-type siderophore genes are found in 14 strains, catecholate-type siderophore genes in strain Mu109, and both types in strain ATCC 9174. Siderophore biosynthesis is predicted to occur in all of the nine strains of the phylogenetic group aurantiacum/sandarakinum/antiquum, in four of the six strains of the group linens/siliguriense/iodinum, in the two strains of the group casei, but only in one of the six other Brevibacterium strains, which corresponded to the cheese isolate Mu109. Interestingly, analysis of the genomes present in the IMG database also revealed that four gene clusters involved in iron acquisition were shared between Brevibacterium strains isolated from cheeses and cheese isolates belonging to other genera (Glutamicibacter, Microbacterium and Corynebacterium). This corresponded to recent HGT events since the percentages of identity at the amino acid level between the genes in Brevibacterium and the genes in the other genera were typically ~95–100%. In most cases, transposase genes were located close to the clusters. These clusters were denoted Iron-Brev1, Iron-Brev2, Iron-Brev3 and Iron-Brev4 (Additional file 9). The cluster Iron-Brev1 corresponds to the ActinoRUSTI region that was recently described .
Osmotic stress tolerance
Cheeses are salted by applying salt to their surface or by submerging them in saturated brine. One mechanism to overcome osmotic stress that results from high salt concentration is the accumulation of osmoprotectants. Brevibacterium strains are known to produce the osmoprotectant ectoine . It is synthesized from L-aspartate-semialdehyde by the action of three enzymes: diaminobutyrate aminotransferase (EctB) [EC 18.104.22.168], diaminobutyrate acetyltransferase (EctA) [EC 22.214.171.124] and ectoine synthase (EctC) [EC 126.96.36.199]. Ectoine can be further converted to hydroxyectoine by the action of ectoine dioxygenase (EctD). Except for strains ATCC 49030T and 5401308T, all the investigated Brevibacterium genomes contained the ectABC cluster (Additional file 10). The ectD gene was present in 13 of the 15 strains from the phylogenetic groups aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum, but it was absent in the two strains of the group casei, and only present in one of the six other Brevibacterium strains. Glycine-betaine is another osmoprotectant that can be synthesized by Brevibacterium. The genomic analysis revealed the presence of two possible pathways. In the first one, choline is converted to glycine betaine by the combined action of a choline dehydrogenase [EC 188.8.131.52] and a betaine aldehyde dehydrogenase [EC 184.108.40.206] whereas in the second one, this conversion is catalyzed by a single enzyme, choline oxidase [EC 220.127.116.11]. Twelve of the strains have both pathways, nine have only the choline oxidase pathway, one has only the choline dehydrogenase pathway, and one has none of them. When present, the choline oxidase gene is located in a cluster containing a betaine/carnitine/choline transporter (BCCT) and a betaine-aldehyde dehydrogenase. Choline is present in milk and cheese , making it available for glycine-betaine biosynthesis. Interestingly the Brevibacterium genomes generally exhibit a greater number of BCCT than the other Actinobacteria genomes (6.3 vs. 1.1 genes per genome, based on assignations to COG1292) (Additional file 10). However, they do not have more ABC transport components involved in the transport of glycine betaine or related osmolytes (4.7 vs. 6.1). Trehalose is a non-reducing sugar that plays a physiological role in energy storage and also as a compatible solute [59, 60]. The genome of the 23 strains encodes a trehalose synthase [EC 18.104.22.168], which catalyzes the synthesis of trehalose from maltose (Additional file 10). Except for strain ATCC 49030T, trehalose can also be produced from UDP-glucose and glucose-6-phosphate via the trehalose-6-phosphate synthase [EC 22.214.171.124]/trehalose-6-phosphate phosphatase [EC 126.96.36.199] pathway. Two operons encoding multisubunit (Na+)(K+)/proton antiporters (Mrp systems) were identified in all the Brevibacterium genomes. These Mrp systems are composed of six subunits and the genes are organized as “group 2” mrp operons . In other Actinobacteria genomes, the mean number of Mrp systems, based on assignation to COG1006, is about four times lower (0.5 vs. 2.0). The abundance of the other (Na+)(K+)/proton antiporters was similar in Brevibacterium and in the other Actinobacteria genomes (mean values of 3.9 and 4.0 genes per genome, respectively, based on assignations to the COG0025, COG1055, COG1757, COG3004, COG3067 and COG3263). It is noteworthy that Brevibacterium strains have a greater number of transporters from the Sodium Solute Symporter (SSS) family than other Actinobacteria strains (7.3 vs. 1.7 genes per genome, based on the number of proteins matching the PF00474 Hidden Markov Model).
Bacteriocines and phenazines
In cheese, bacteriocin producers can inhibit other microbial groups that share the same ecological niche, which confers them a selective advantage. Genomic analysis predicted the production of several bacteriocins in Brevibacterium strains (Additional file 11). Linocin M18-related bacteriocins (PF04454) were identified in 15 genomes, including five of the nine strains of the phylogenetic group aurantiacum/sandarakinum/antiquum, all six strains of the group linens/siliguriense/iodinum and four of the eight other strains (B. album DSM 18261T, B. senegalense JC43T, B. jeotgali SJ5-8T and B. sp. Mu109). Five groups of ribosomally synthesized and post-translationally modified peptides (RiPPs)  were predicted in the 23 Brevibacterium strains. The first group corresponds to lanthipeptides, which are characterized by the presence of lanthionine (Lan) and/or methyllanthionine (MeLan) and/or labionin (Lab) residues [62, 63]. Many lanthipeptides have an antimicrobial activity, mainly against Gram-positive bacteria , and are referred to as lantibiotics. Lanthipeptide gene clusters were predicted in seven strains, of which six are cheese-associated (Additional file 11). Even if they have different structures, these clusters contain one or two putative lanthipeptide synthetases, one or two putative precursor peptides, one to three putative ABC transport system components and, in some cases, a prolyl oligopeptidase (PF00326) (Fig. 3A). Protein sequence analysis of the predicted lanthipeptide synthetases suggests that all of them belong to class III, corresponding to a trifunctional enzyme LanKC . The sequences of all putative precursor peptides (Fig. 3B) contain two characteristic and conserved Ser/Ser(Thr)/Cys motifs necessary for the formation of (Me)Lan and/or (Me)Lab. One of them also contains a B-A-C-Leu-Gln motif in its N-terminal part (where B is Ile, Leu or Val; A is Phe or Leu; and C is Glu or Asp), which is highly conserved in class III precursor peptides and essential for the enzymatic processing of the labyrinthopeptin A2 . To our knowledge, as of this time, the Kribbella flavida DSM 17836 protein FlaP, which is a proline-specific oligopeptidase, is the only characterized protease involved in the removal of the leader peptide of a class III lanthipeptide . Genes with a weak homology with flaP were identified in six of the seven predicted lanthipeptide gene clusters, but only one of the six predicted precursor peptides contains a Pro residue in the leader peptide, which could serve as the primary cleavage site for the predicted peptidase (Fig. 3B). Interestingly, the putative lanthipeptide gene cluster from the four cheese-associated strains B. antiquum CNRZ 918, B. antiquum P10, B. aurantiacum ATCC 9174 and B. linens ATCC 9172T is located in a ~96 kb genomic island, which is absent in the other Brevibacterium genomes investigated, but which is present in the genome of Corynebacterium casei LMG S-19264, a strain isolated from a smear-ripened cheese . This genomic island, which we denoted as BreLI (Brevibacterium Lanthipeptide Island), probably corresponds to an Integrative and Conjugative Element (ICE, or conjugative transposon). ICEs are typically composed of three core genetic modules involved in: (i) integration and excision; (ii) conjugation; and (iii) regulation . Functional prediction of genes in BreLI revealed assignations to these three functions (Fig. 3C and Additional file 11). Integration of this ICE occurred at the 3′ end of a gene encoding a class Ib ribonucleotide reductase beta subunit, resulting in a 12-bp perfect repeat sequence, present at the two borders of BreLI. It is noteworthy that, in comparison to the four other strains, a segment containing three genes of the lanthipeptide gene cluster region is lacking in the ATCC 9174 BreLI Island (compare clusters (a) and (b) in Fig. 3A). The second group of predicted RiPPs is related to lactococcin 972 (PF09683), which has been characterized in Lactococcus lactis IPLA 972 [69, 70]. Genes with a weak homology with the lactococcin 972 precursor peptide were identified in six cheese-associated strains from the phylogenetic group aurantiacum/sandarakinum/antiquum, corresponding to two different putative prepeptides, one of them being present only in strain 239c, and the other in five B. aurantiacum strains (Additional file 11). In all these six strains, the bacteriocin structural gene forms an operon with three genes encoding two transmembrane proteins and a putative ABC transport systeme ATP-binding protein. Comparison of the flanking regions of the lactococcin 972-related bacteriocin biosynthesis gene cluster in strains 6(3) and 8(6) with the CNRZ 918 genome (in which the cluster is absent) revealed that the gene cluster is located in a genomic island that is inserted at the end of a tRNA-val gene (locus tag Ga0063697_02943 in strain 6(3), Ga0063698_00389 in strain 8(6) and Ga0063689_00412 in strain CNRZ 918), which is followed in CNRZ 918 by an uracil-xanthine permease (locus tag Ga0063689_00413, which is an ortholog of Ga0063697_02919 in 6(3) and Ga0063698_00365 in 8(6)). The lactococcin 972-related bacteriocin biosynthesis gene cluster in strains ATCC 9174 and CNRZ 920 is also flanked by a tRNA-val gene. The third group of predicted RiPPs corresponds to linear azol(in)e-containing peptides (LAPs), which are characterized by the presence of multiple thiazole and (methyl)oxazole heterocycles, and sometimes by their corresponding 2-electron reduced azoline state. LAP gene clusters were predicted in four B. aurantiacum cheese-associated strains (Additional file 11). One cluster is present in strains 6(3), ATCC 9175T and ATCC 9174, and another is present only in strain 8(6). All these clusters contain putative precursor peptides, a SagB-type dehydrogenase domain-containing protein (TIGR03605), one or two YcaO cyclodehydratases (PF02624) and a cyclodehydratase-docking fusion protein (TIGR03882), representing the critical components in LAP biosynthesis [62, 71]. In addition, the gene cluster in strains ATCC 9175T, ATCC 9174 and 6(3) also contains a methyltransferase, which may be involved in the methylation of the peptide . Comparison of the flanking regions of the predicted LAP gene cluster in strain 6(3) with the CNRZ 920 genome (in which the cluster is absent) revealed that the gene cluster is located in a genomic island whose insertion is between an aspartate racemase (orthologs with locus tag Ga0063697_02599 in strain 6(3) and Ga0063691_01440 in strain CNRZ 920) and a short-chain dehydrogenase (orthologs with locus tag Ga0063697_02614 in strain 6(3) and Ga0063691_01438 in strain CNRZ 920). The LAP gene cluster in strain 8(6) contains two putative ABC transport system components. This type of component is frequently associated with LAP biosynthesis gene clusters [62, 72]. The fourth group of predicted RiPPs corresponds to the recently discovered linaridin family, which is characterized by the presence of thioether crosslinks, like for the lanthipeptide family, but whose post-translational modifications are unrelated to the lanthipeptide biosynthesis pathway . Linaridin precursor peptides were predicted in the two cheese-associated strains 239c and 8(6). These two prepeptides showed significant similarity to the characterized cypemycin (51 to 53% identity) and grisemycin (58% identity) prepeptides [73, 74]. The putative linaridin gene clusters of strains 239c and 8(6) are very similar (Additional file 11). They contain homologs of the cypemycin biosynthetic gene cluster, i.e. cypH, cypL, cypD and cypT . The fifth group of predicted RiPPs is related to the sporulation delaying protein (SDP), which has been characterized in Bacillus subtilis . SDP is produced by the sdpABC operon in which sdpC encodes the SDP precursor peptide (pro-SdpC) and sdpA and sdpB are essential for the production of the active SDP toxin . Genes with a weak homology with sdpABC were identified in three B. aurantiacum cheese-associated strains (Additional file 11). Comparison of the corresponding flanking regions in strains ATCC 9175T and 6(3) with the 8(6) genome (in which the cluster is absent) revealed that this gene cluster is located in a genomic island that is inserted at the end of a NAD(P)H-dependent flavin oxidoreductase (locus tag Ga0063690_01635 in strain ATCC 9175T, Ga0063697_00339 in strain 6(3) and Ga0063698_01092 in strain 8(6)), which is followed in 8(6) by a LuxR family transcriptional regulator (locus tag Ga0063698_01093, which is an ortholog of Ga0063690_01646 in ATCC 9175T and Ga0063697_00350 in strain 6(3)).
Phenazines are heterocyclic compounds that are substituted at different points around their rings . They have various colors and biological activities, such as antibiotic and intercellular signaling activities. Brevibacterium iodinum ATCC 49514T, isolated from milk, is known to produce purple extracellular crystals of the phenazine iodinin . Its genome contains the phzGFEDCB gene cluster, which encodes one enzyme involved in chorismate synthesis (PhzG) and the five enzymes required for the generation of the “core” phenazines, phenazine-1,6- dicarboxylic acid or phenazine-1-carboxylic acid, which are precursors for strain-specific phenazine derivates (Additional file 12). This cluster is also present in four other strains, all isolated from cheeses (ATCC 9175T, ATCC 9172T, P10, and ATCC 9174), which are thus probably able to produce phenazine derivatives.
Microbial cultures are widely used in the cheesemaking industry. There is a high degree of expertise today in terms of the design of lactic starter cultures. This expertise includes several important properties such as the ability to grow in milk, the resistance to bacteriophages, and the generation of adequate acidifying and texturizing activities. The situation is more problematic for ripening cultures for surface-ripened cheese, which contain strains that sometimes do not grow well in cheese. This may be explained by an insufficient fitness in the cheese surface habitat in comparison to the resident “house flora”. Indeed, contrary to lactic starter cultures where growth is very fast (several hours) and facilitated by a massive inoculation in milk that contains only a low concentration of microbial cells, growth at the surface of cheeses takes a longer time (several weeks) and therefore offers more opportunities for the development of well-adapted adventitious strains.
Cheese-associated strains of Brevibacterium belong to different species, showing that adaptation to the cheese habitat occurred independently in different lineages. The strains sequenced in the present study belong to the species B. aurantiacum, B. antiquum, B. linens, B. casei, as well as to two not-yet described species (B. sp. Mu109 and 239c). In some cases, the genomes of cheese-associated strains are closely related to other strains, suggesting a recent adaptation. For example, the ANI value between strains S18 (human skin) and CIP 102111T (cheese) is higher than the values between the different cheese isolates belonging to the species B. aurantiacum (98.22% vs. 97.48–98.04%). Interestingly, a complete prophage was found in B. sp. Mu109, which indicates that cheese-associated Brevibacterium strains may undergo phage attacks, even if, to our knowledge, there is no study in the scientific literature concerning the sensitivity to phages of ripening cultures containing Brevibacterium strains. It is possible that the impact of phages on these cultures has been overlooked.
In this study, genomic analyses revealed that Brevibacterium strains show differences in the ability to use energy compounds present in cheeses. Indeed, the galactose catabolism pathway is predicted to occur only in the phylogenetic groups aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum, and the D-galactonate catabolism pathway only in the phylogenetic groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei. The strains from the aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum groups have more lactate dehydrogenases than the strains from the other groups, suggesting that they have a better ability to consume lactate, which may be useful during growth in cheese. A predicted extracellular triacylglycerol lipase was found in 14 of the 15 genomes from the aurantiacum/sandarakinum/antiquum and linens/siliguriense/iodinum groups and in the cheese-associated strain Mu109, suggesting that many Brevibacterium strains are likely to contribute to cheese lipolysis. The resulting fatty acids constitute a potential energy source for the lipolitic Brevibacterium strains as well as for other microbial populations living on the cheese surface. A cluster of two genes encoding two secreted proteases with LPXTG motif is present in 16 of the 23 investigated Brevibacterium genomes. In B. aurantiacum ATCC 9174, one of these two proteases has been purified and characterized . This protease is active on casein and secreted during the growth of the strain, which indicates that it probably contributes to cheese proteolysis [79, 80]. The amino acids produced from caseins and the non-protein amino acid 4-aminobutyrate are major energy substrates for the bacteria growing at the cheese surface. The 17 strains belonging to the groups aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei are well equipped in amino acid catabolism pathways and have a similar gene content. The six other Brevibacterium strains have a lower catabolic potential since they do not encode the enzymes involved in arginine, phenylalanine, methionine and tyrosine degradation.
Microorganisms growing at the cheese surface have to be able to withstand the osmotic stress due to the presence of salt. Most Brevibacterium genomes encode the pathways for the production of the osmoprotectants ectoine or hydroxyectoine from aspartate, glycine-betaine from choline, and trehalose. In addition, in comparison to the other Actinobacteria genomes, they also encode a large number of Betaine/Carnithine/Choline family Transporters (BCCT), of multisubunit (Na+)(K+)/proton antiporters, and of Sodium Solute Symporters (SSS). The two last systems take advantage of Na+ gradients to regulate the intracellular pH and to import nutrients. These properties probably contribute to the good resistance of Brevibacterium strains to the cheese salt and to their stimulation in the presence of 4% NaCl . No major differences in the gene content concerning osmotic stress resistance was observed between the cheese-associated or the other strains, and between the investigated phylogenetic groups, except for B. ravenspurgense 5401308T and B. mcbrellneri ATCC 49030T, which have a smaller number of BCCT and SSS systems. These two strains have the smallest genome sizes of the 23 investigated strains and have probably undergone a reductive genome evolution.
The present study confirms the central role of iron metabolism in bacteria from cheese surface microbial communities. Brevibacterium strains produce and utilize siderophores , and the growth of typical aerobic bacteria at the cheese surface is limited by the availability of iron . The low iron content of milk, the presence of sequestering compounds such as lactoferrin, the presence of oxygen and the high pH of the cheese matrix during the growth of the acido-sensitive bacteria contribute to restricting the availability of iron. One strategy for improving iron acquisition in the cheese habitat is to produce siderophores, and this capacity can be acquired or improved by HGT. Hydroxamate-type siderophore biosynthesis seems to be widespread in Brevibacterium strains since gene clusters involved in their production were found in most genomes of the aurantiacum/sandarakinum/antiquum, linens/siliguriense/iodinum and casei groups. Clusters encoding catecholate-type siderophores were present only in two strains and resulted from horizontal transfers from other genera. Interestingly, both strains were isolated from cheese, and one of them (ATCC 9174) also possesses a gene cluster involved in the production of an acetohydroxamate-type siderophore. At least one other cheese-associated strain, Glutamicibacter arilaitensis Re117, has two siderophore biosynthesis clusters, including a catecholate siderophore resulting from HGT . It is likely that siderophore production by cheese microorganisms also results in biotic interactions. In a recent study, it was observed that a Staphylococcus equorum strain, which was a weak competitor against other closely related Staphylococcus species in model cheese experiments, became dominant in the presence of the fungus Scopulariopsis . This effect was attributed to fungal siderophore production, which may relieve S. equorum of the costly production of the siderophore staphyloferrin B and potentially provide an iron source through cross-feeding. This type of cross-feeding may also occur in microbial communities containing Brevibacterium strains that do not have the potential to produce siderophores, such as B. linens ATCC 9172T and B. linens Mu101. The Fe3+/siderophore complexes are imported into the cells by ABC transport systems, whether or not these siderophores are produced by the same strain. Numerous ABC-type Fe3+/siderophore components are present in Brevibacterium strains. These systems allow the cells to take advantage of different types of siderophores available in the medium. Horizontal gene transfers concerning ABC-type Fe3+/siderophore components were observed in the cheese-associated strains Glutamicibacter arilaitensis Re117 , in Corynebacterium variabile DSM 44702 , and a recent study provided evidence of extensive HGTs concerning Fe3+/siderophore acquisition in very diverse cheese-associated bacteria . One of these regions, known as ActinoRUSTI, was present in two of the 23 Brevibacterium strains investigated in the present study. Three other islands involved in iron acquisition that are present both in cheese-associated Brevibacterium strains and in cheese-associated strains belonging to other genera (Glutamicibacter arilaitensis Re117, Corynebacterium casei UCMA 3821, Corynebacterium variabile DSM 44702 and Microbacterium gubbeenense DSM 15944) were detected in the investigated genomes. Interestingly, three of these four islands were present in the cheese-associated strain B. sp. Mu109, which is also the strain with the highest number of ABC-type Fe3+/siderophore components (31 components).
Bacteriocin production is thought to play a critical role in mediating the microbial population or community interactions . It may thus be assumed that it could have a significant impact on the cheese surface where microbial density may exceed 1010 cells per cm2. In the present study, we detected several putative bacteriocin gene clusters in Brevibacterium genomes. They corresponded to six groups of bacteriocins, i.e., linocin M18-related bacteriocins and five groups of RiPPs, and were highly variable among strains. Linocin M18, which has been characterized in the red smear cheese bacterium B. linens M18, is an antilisterial and wide-spectrum bacteriocin . Its activity against Listeria spp. has also been demonstrated in a model cheese . In our study, linocin M18-related bacteriocins were detected in 15 of the 23 investigated genomes, belonging to different phylogenetic groups and isolated from different habitats. This result is consistent with the fact that the structural gene lin encoding linocin M18 is widely distributed in coryneform bacteria . To our knowledge, up until now, there has been no experimental evidence about the production of RiPPs in Brevibacterium. However, it cannot be excluded that some of the antibacterial substances characterized from Brevibacterium strains, such as the Linecin A from B. aurantiacum ATCC 9175T  and/or the Linenscin OC2 from B. linens OC2 , are in fact RiPPs. Almost all the RiPP gene clusters detected in our study were found in cheese-associated strains, except for a lanthipeptide gene cluster in B. mcbrellneri ATCC 49030T, which is a human-associated strain. Interestingly, many of them seem to result from HGTs. In this study, we identified the BreLI island, which is probably a ~96-kb integrative and conjugative element (ICE) encoding for class III lanthipeptides. This island is present in four cheese-associated Brevibacterium strains as well as in a cheese-associated strain from another genus (Corynebacterium casei LMG S-19264). Taken together, our results consolidate the hypothesis that bacteriocin production may provide an ecological advantage to cheese-associated bacteria. It would be interesting in further studies to examine whether these bacteriocin gene clusters are functional, to investigate the activity spectra of these bacteriocins and to determine the influence of environmental conditions on their biosynthesis. Such information would be useful for the design of surface-ripened cheese cultures in order to improve their competitiveness against adventitious strains or to prevent growth of pathogens and spoilage microorganisms. It would also be interesting to investigate the role of phenazines, which have a broad-spectrum antibiotic activity and whose biosynthesis is predicted to occur in four cheese-associated strains. These strains belong to three different species, but the presence of the corresponding gene clusters cannot be explained by recent HGT events in view of the divergence in gene sequences.
The present study confirms that there are HGT events between microorganisms growing at the surface of cheeses. Acquisition of genes involved in siderophore biosynthesis, iron import and bacteriocin production can probably improve the fitness of the strains in the cheese habitat. It can be hypothesized that these gene transfers exert an influence on the balance between the resident “house flora” and the strains from the ripening cultures. For the latter type of strains, it is, in most cases, the same strains that are massively inoculated in all the manufacturing runs. For the adventitious “house flora”, it would be beneficial to acquire genes that improve their competitiveness to the detriment of the inoculated strains, such as genes governing bacteriocin production or iron acquisition. It would thus be interesting to examine the extent to which gene acquisition by the adventitious strains impacts the growth and stability of the components from ripening cultures used for the production of surface-ripened cheeses.
Some properties deduced from genome analyses are similar in all the investigated strains, such as the ability to catabolize ethanol or 4-aminobutyrate. This is also observed for the ability to catabolize glycerol and for osmotolerant biosynthesis, except for strains 5401308T and ATCC 49030T, which are two strains that have a smaller genome size. Other properties are mainly correlated to the phylogenetic position of the strains, whether they were isolated or not from cheese, such as the ability to catabolize galactose, lactate, amino acids, or to secrete triacylglycerol lipases. The ability to catabolize D-galactonate is present in part of the strains, and this property does not seem to be correlated to the phylogenetic position or to the habitat of the strains. There are great differences in the number of Fe3+/siderophore ABC transport components. Some of these genes are in clusters that are also present in cheese-associated bacteria belonging to other genera, indicating that these genes are disseminated by HGTs among strains living on the cheese surface. Two Brevibacterium strains isolated from cheeses also acquired a catecholate-type siderophore biosynthesis gene cluster by HGT. Bacteriocin biosynthesis genes are present in most of the strains, and one of the corresponding gene clusters is located in a probable conjugative transposon of ~96 kb (BreLI), which is present in four cheese-associated Brevibacterium strains as well as in Corynebacterium casei LMG S-19264, a strain isolated from a smear-ripened cheese. Brevibacterium strains thus show differences in genetic determinants involved in the growth on the cheese surface. Some of them are correlated to the phylogenetic position and others are the result of gene transfers. Part of these properties contributes to biotic interactions between strains. In the future, it would be interesting to take this information into account in order to improve the screening and selection of Brevibacterium strains and their association with other ripening culture components.
Average Nucleotide Identity
Betaine/Carnithine/Choline family Transporter
Brevibacterium Lanthipeptide Island
Coding DNA sequence
Clustered Regularly Interspaced Short Palindromic Repeats
Horizontal Gene transfer
Hidden Markov Model
Integrative and Conjugative Element
NAD-independent lactate dehydrogenase
Linear Azol(in)e-containing Peptides
NAD-dependant lactate dehydrogenase
Ribosomally synthesized and Post-translationally modified Peptide
Sporulation Delaying Protein
Sodium Solute Symporter
Bockelmann W. Secondary cheese starter cultures. In: law BA, Tamime AY, editors. Technol. Cheesemaking second Ed. Oxford: Wiley-Blackwell; 2010. p. 193–230.
Irlinger F, Layec S, Hélinck S, Dugat-Bony E. Cheese rind microbial communities: diversity, composition and origin. FEMS Microbiol Lett. 2015;362:1–11.
Monnet C, Landaud S, Bonnarme P, Swennen D. Growth and adaptation of microorganisms on the cheese surface. FEMS Microbiol Lett. 2015;362:1–9.
Brennan NM, Ward AC, Beresford TP, Fox PF, Goodfellow M, Cogan TM. Biodiversity of the bacterial flora on the surface of a smear cheese. Appl Environ Microbiol. 2002;68:820–30.
Feurer C, Vallaeys T, Corrieu G, Irlinger F. Does smearing inoculum reflect the bacterial composition of the smear at the end of the ripening of a French soft, red-smear cheese? J Dairy Sci. 2004;87:3189–97.
Goerges S, Mounier J, Rea MC, Gelsomino R, Heise V, Beduhn R, et al. Commercial ripening starter microorganisms inoculated into cheese milk do not successfully establish themselves in the resident microbial ripening consortia of a south german red smear cheese. Appl Environ Microbiol. 2008;74:2210–7.
Gori K, Ryssel M, Arneborg N, Jespersen L. Isolation and identification of the microbiota of Danish farmhouse and industrially produced surface-ripened cheeses. Microb Ecol. 2013;65:602–15.
Mounier J, Gelsomino R, Goerges S, Vancanneyt M, Vandemeulebroecke K, Hoste B, et al. Surface microflora of four smear-ripened cheeses. Appl Environ Microbiol. 2005;71:6489–500.
Rea MC, Görges S, Gelsomino R, Brennan NM, Mounier J, Vancanneyt M, et al. Stability of the biodiversity of the surface consortia of Gubbeen, a red-smear cheese. J Dairy Sci. 2007;90:2200–10.
Rattray FP, Fox PF. Aspects of enzymology and biochemical properties of Brevibacterium linens relevant to cheese ripening: a review. J Dairy Sci. 1999;82:891–909.
Forquin-Gomez M-P, Weimer BC, Sorieul L, Kalinowski J, Vallaeys T. The family Brevibacteriaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. Prokaryotes Actinobacteria. 4th ed. Berlin, Heidelberg: Springer; 2014. p. 141–53.
Onraedt A, Soetaert W, Vandamme E. Industrial importance of the genus Brevibacterium. Biotechnol Lett. 2005;27:527–33.
Gavrish EI, Krauzova VI, Potekhina NV, Karasev SG, Plotnikova EG, Altyntseva OV, et al. Three new species of brevibacteria, Brevibacterium antiquum sp. nov., Brevibacterium aurantiacum sp. nov. and Brevibacterium permense sp. nov. Mikrobiologiia. 2004;73:218–25.
Montel M-C, Buchin S, Mallet A, Delbes-Paus C, Vuitton DA, Desmasures N, et al. Traditional cheeses: rich and diverse microbiota with associated benefits. Int J Food Microbiol. 2014;177:136–54.
Schröder J, Maus I, Trost E, Tauch A. Complete genome sequence of Corynebacterium Variabile DSM 44702 isolated from the surface of smear-ripened cheeses and insights into cheese ripening and flavor generation. BMC Genomics. 2011;12:545.
Monnet C, Loux V, Gibrat J-F, Spinnler E, Barbe V, Vacherie B, et al. The arthrobacter arilaitensis Re117 genome sequence reveals its genetic adaptation to the surface of cheese. PLoS One. 2010;5:e15489.
Bonham KS, Wolfe BE, Dutton RJ. Extensive horizontal gene transfer in cheese-associated bacteria. elife. 2017;6
Markowitz VM, I-MA C, Chu K, Pati A, Ivanova NN, Kyrpides NC. Ten years of maintaining and expanding a microbial genome and metagenome analysis system. Trends Microbiol. 2015;23:730–41.
Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinforma Oxf Engl. 2011;27:2957–63.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. J. Comput Mol Cell Biol. 2012;19:455–77.
Markowitz VM, I-MA C, Palaniappan K, Chu K, Szeto E, Grechkin Y, et al. IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res. 2012;40:D115–22.
Huntemann M, Ivanova NN, Mavromatis K, Tripp HJ, Paez-Espino D, Palaniappan K, et al. The standard operating procedure of the DOE-JGI microbial genome annotation pipeline (MGAP v.4). Stand Genomic Sci. 2015;10:86.
Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44:W16–21.
Weber T, Blin K, Duddela S, Krug D, Kim HU, Bruccoleri R, et al. antiSMASH 3.0-a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 2015;43:W237–43.
van Heel AJ, de Jong A, Montalbán-López M, Kok J, Kuipers OP. BAGEL3: Automated identification of genes encoding bacteriocins and (non-)bactericidal posttranslationally modified peptides. Nucleic Acids Res. 2013;41:W448–53.
Hammami R, Zouhir A, Le Lay C, Ben Hamida J, Fliss I. BACTIBASE second release: a database and tool platform for bacteriocin characterization. BMC Microbiol. 2010;10:22.
Mende DR, Sunagawa S, Zeller G, Bork P. Accurate and universal delineation of prokaryotic species. Nat Methods. 2013;10:881–4.
Walter F, Albersmeier A, Kalinowski J, Rückert C. Complete genome sequence of Corynebacterium Casei LMG S-19264T (=DSM 44701T), isolated from a smear-ripened cheese. J Biotechnol. 2014;189:76–7.
Almeida M, Hébert A, Abraham A-L, Rasmussen S, Monnet C, Pons N, et al. Construction of a dairy microbial genome catalog opens new perspectives for the metagenomic analysis of dairy fermented products. BMC Genomics. 2014;15:1101.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Notredame C, Higgins DG, Heringa J. T-coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302:205–17.
Edgar RC. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Price MN, Dehal PS, Arkin AP. FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33:1870–4.
Contreras-Moreira B, Vinuesa P. GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis. Appl Environ Microbiol. 2013;79:7696–701.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6.
Stackebrandt E, Frederiksen W, Garrity GM, Grimont PAD, Kämpfer P, Maiden MCJ, et al. Report of the ad hoc committee for the re-evaluation of the species definition in bacteriology. Int J Syst Evol Microbiol. 2002;52:1043–7.
Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A. 2005;102:2567–72.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57:81–91.
Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106:19126–31.
Stackebrandt E, Goebel BM. Taxonomic note: a place for DNA-DNA Reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. Int J Syst Evol Microbiol. 1994;44:846–9.
Michel V, Martley FG. Streptococcus Thermophilus in cheddar cheese--production and fate of galactose. J Dairy Res. 2001;68:317–25.
Mounier J, Rea MC, O’Connor PM, Fitzgerald GF, Cogan TM. Growth characteristics of Brevibacterium, Corynebacterium, microbacterium, and staphylococcus spp. isolated from surface-ripened cheese. Appl Environ Microbiol. 2007;73:7732–9.
Jiang T, Gao C, Ma C, Microbial XP. Lactate utilization: enzymes, pathogenesis, and regulation. Trends Microbiol. 2014;22:589–99.
Pinchuk GE, Rodionov DA, Yang C, Li X, Osterman AL, Dervyn E, et al. Genomic reconstruction of Shewanella oneidensis MR-1 metabolism reveals a previously uncharacterized machinery for lactate utilization. Proc Natl Acad Sci U S A. 2009;106:2874–9.
Jolkver E, Emer D, Ballan S, Krämer R, Eikmanns BJ, Marin K. Identification and characterization of a bacterial transport system for the uptake of pyruvate, propionate, and acetate in Corynebacterium glutamicum. J Bacteriol. 2009;191:940–8.
Auchter M, Arndt A, Eikmanns BJ. Dual transcriptional control of the acetaldehyde dehydrogenase gene ald of Corynebacterium glutamicum by RamA and RamB. J Biotechnol. 2009;140:84–91.
Hajri T, Abumrad NA. Fatty acid transport across membranes: relevance to nutrition and metabolic pathology. Annu Rev Nutr. 2002;22:383–415.
Muñoz-Elías EJ, Upton AM, Cherian J, McKinney JD. Role of the methylcitrate cycle in mycobacterium tuberculosis metabolism, intracellular growth, and virulence. Mol Microbiol. 2006;60:1109–22.
Bott M, Niebisch A. The respiratory chain of Corynebacterium glutamicum. J Biotechnol. 2003;104:129–53.
Rawlings ND, Barrett AJ, Finn R. Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2016;44:D343–50.
Rattray FP, Bockelmann W, Fox PF. Purification and characterization of an extracellular proteinase from Brevibacterium linens ATCC 9174. Appl Environ Microbiol. 1995;61:3454–6.
Nomura M, Kimoto H, Someya Y, Furukawa S, Suzuki I. Production of gamma-aminobutyric acid by cheese starters during cheese ripening. J Dairy Sci. 1998;81:1486–91.
Miethke M, Monteferrante CG, Marahiel MA, van Dijl JM. The Bacillus Subtilis EfeUOB transporter is essential for high-affinity acquisition of ferrous and ferric iron. Biochim Biophys Acta. 2013;1833:2267–78.
Khan AG, Shouldice SR, Kirby SD, Yu R, Tari LW, Schryvers AB. High-affinity binding by the periplasmic iron-binding protein from Haemophilus influenzae is required for acquiring iron from transferrin. Biochem J. 2007;404:217–25.
Bernard T, Jebbar M, Rassouli Y, Himdi-Kabbab S, Hamelin J, Blanco C. Ectoine accumulation and osmotic regulation in Brevibacterium linens. Microbiology. 1993;139:129–36.
Zeisel SH, Mar M-H, Howe JC, Holden JM. Concentrations of choline-containing compounds and betaine in common foods. J Nutr. 2003;133:1302–7.
Frings E, Kunte HJ, Galinski EA. Compatible solutes in representatives of the genera Brevibacterium and Corynebacterium: occurrence of tetrahydropyrimidines and glutamine. FEMS Microbiol Lett. 1993;109:25–32.
Kempf B, Bremer E. Uptake and synthesis of compatible solutes as microbial stress responses to high-osmolality environments. Arch Microbiol. 1998;170:319–30.
Swartz TH, Ikewada S, Ishikawa O, Ito M, Krulwich TA. The Mrp system: a giant among monovalent cation/proton antiporters. Extremophiles. 2005;9:345–54.
Arnison PG, Bibb MJ, Bierbaum G, Bowers AA, Bugni TS, Bulaj G, et al. Ribosomally synthesized and post-translationally modified peptide natural products: overview and recommendations for a universal nomenclature. Nat Prod Rep. 2013;30:108–60.
Meindl K, Schmiederer T, Schneider K, Reicke A, Butz D, Keller S, et al. Labyrinthopeptins: a new class of carbacyclic lantibiotics. Angew Chem Int Ed Engl. 2010;49:1151–4.
Bierbaum G, Sahl H-G. Lantibiotics: mode of action, biosynthesis and bioengineering. Curr Pharm Biotechnol. 2009;10:2–18.
Zhang Q, Yu Y, Vélasquez JE, van der Donk WA. Evolution of lanthipeptide synthetases. Proc Natl Acad Sci U S A. 2012;109:18361–6.
Müller WM, Ensle P, Krawczyk B, Süssmuth RD. Leader peptide-directed processing of labyrinthopeptin A2 precursor peptide by the modifying enzyme LabKC. Biochemistry (Mosc). 2011;50:8362–73.
Völler GH, Krawczyk B, Ensle P, Süssmuth RD. Involvement and unusual substrate specificity of a prolyl oligopeptidase in class III lanthipeptide maturation. J Am Chem Soc. 2013;135:7426–9.
Bi D, Xu Z, Harrison EM, Tai C, Wei Y, He X, et al. ICEberg: a web-based resource for integrative and conjugative elements found in bacteria. Nucleic Acids Res. 2012;40:D621–6.
Martínez B, Fernández M, Suárez JE, Rodríguez A. Synthesis of lactococcin 972, a bacteriocin produced by Lactococcus lactis IPLA 972, depends on the expression of a plasmid-encoded bicistronic operon. Microbiol Read Engl. 1999;145(Pt 11):3155–61.
Sánchez C, Hernández de Rojas A, Martínez B, Argüelles ME, Suárez JE, Rodríguez A, et al. Nucleotide sequence and analysis of pBL1, a bacteriocin-producing plasmid from Lactococcus lactis IPLA 972. Plasmid. 2000;44:239–49.
Haft DH, Basu MK, Mitchell DA. Expansion of ribosomally produced natural products: a nitrile hydratase- and Nif11-related precursor family. BMC Biol. 2010;8:70.
Lee SW, Mitchell DA, Markley AL, Hensler ME, Gonzalez D, Wohlrab A, et al. Discovery of a widely distributed toxin biosynthetic gene cluster. Proc Natl Acad Sci U S A. 2008;105:5879–84.
Claesen J, Bibb M. Genome mining and genetic analysis of cypemycin biosynthesis reveal an unusual class of posttranslationally modified peptides. Proc Natl Acad Sci U S A. 2010;107:16297–302.
Claesen J, Bibb MJ. Biosynthesis and regulation of grisemycin, a new member of the linaridin family of ribosomally synthesized peptides produced by Streptomyces Griseus IFO 13350. J Bacteriol. 2011;193:2510–6.
Liu W-T, Yang Y-L, Xu Y, Lamsa A, Haste NM, Yang JY, et al. Imaging mass spectrometry of intraspecies metabolic exchange revealed the cannibalistic factors of Bacillus Subtilis. Proc Natl Acad Sci U S A. 2010;107:16286–90.
Pérez Morales TG, Ho TD, Liu W-T, Dorrestein PC, Ellermeier CD. Production of the cannibalism toxin SDP is a multistep process that requires SdpA and SdpB. J Bacteriol. 2013;195:3244–51.
Price-Whelan A, Dietrich LEP, Newman DK. Rethinking “secondary” metabolism: physiological roles for phenazine antibiotics. Nat Chem Biol. 2006;2:71–8.
Whitman W, Goodfellow M, Kämpfer P, Busse H-J, Trujillo M, Ludwig W, et al. Bergey’s manual of systematic bacteriology (second edition) volume 5: the Actinobacteria. New York: Springer Science & Business Media; 2012.
Rattray FP, Fox PF, Healy A. Specificity of an extracellular proteinase from Brevibacterium linens ATCC 9174 on bovine alpha s1-casein. Appl Environ Microbiol. 1996;62:501–6.
Rattray FP, Fox PF, Healy A. Specificity of an extracellular proteinase from Brevibacterium linens ATCC 9174 on bovine beta-casein. Appl Environ Microbiol. 1997;63:2468–71.
Masoud W, Jakobsen M. The combined effects of pH, NaCl and temperature on growth of cheese ripening cultures of Debaryomyces Hansenii and coryneform bacteria. Int Dairy J. 2005;15:69–77.
Noordman WH, Reissbrodt R, Bongers RS, Rademaker JLW, Bockelmann W, Smit G. Growth stimulation of Brevibacterium sp. by siderophores. J Appl Microbiol. 2006;101:637–46.
Monnet C, Back A, Irlinger F. Growth of aerobic ripening bacteria at the cheese surface is limited by the availability of iron. Appl Environ Microbiol. 2012;78:3185–92.
Kastman EK, Kamelamela N, Norville JW, Cosetta CM, Dutton RJ, Wolfe BE. Biotic interactions shape the ecological distributions of staphylococcus species. MBio. 2016;7:e01157–16.
Riley MA, Wertz JE. Bacteriocins: evolution, ecology, and application. Annu Rev Microbiol. 2002;56:117–37.
Valdés-Stauber N, Scherer S. Isolation and characterization of Linocin M18, a bacteriocin produced by Brevibacterium linens. Appl Environ Microbiol. 1994;60:3809–14.
Eppert I, Valdés-Stauber N, Götz H, Busse M, Scherer S. Growth reduction of listeria spp. caused by undefined industrial red smear cheese cultures and bacteriocin-producing Brevibacterium linens as evaluated in situ on soft cheese. Appl Environ Microbiol. 1997;63:4812–7.
Valdes-Stauber N, Scherer S. Nucleotide sequence and taxonomical distribution of the bacteriocin gene lin cloned from Brevibacterium linens M18. Appl Environ Microbiol. 1996;62:1283–6.
Kato F, Eguchi Y, Nakano M, Oshima T, Murata A. Purification and characterization of Linecin-a, a Bacteriocin of Brevibacterium linens. Agric Biol Chem. 1991;55:161–6.
Maisnier-Patin S, Richard J. Activity and purification of linenscin OC2, an antibacterial substance produced by Brevibacterium linens OC2, an orange cheese coryneform bacterium. Appl Environ Microbiol. 1995;61:1847–52.
Kämpfer P, Schäfer J, Lodders N, Busse H-J. Brevibacterium sandarakinum sp. nov., isolated from a wall of an indoor environment. Int J Syst Evol Microbiol. 2010;60:909–13.
Maizel D, Utturkar SM, Brown SD, Ferrero MA, Rosen BP. Draft genome sequence of Brevibacterium linens AE038-8, an extremely arsenic-resistant bacterium. Genome Announc. 2015;3:e00316–5.
Roux V, Robert C, Gimenez G, Raoult D. Draft genome sequence of Brevibacterium massiliense strain 541308T. J Bacteriol. 2012;194:5151–2.
Kokcha S, Ramasamy D, Lagier J-C, Robert C, Raoult D, Fournier P-E. Non-contiguous finished genome sequence and description of Brevibacterium senegalense sp. nov. Stand Genomic Sci. 2012;7:233–45.
Ganesan B, Seefeldt K, Weimer BC. Fatty acid production from amino acids and alpha-keto acids by Brevibacterium linens BL2. Appl Environ Microbiol. 2004;70:6385–93.
We are grateful to the INRA MIGALE bioinformatics platform (http://migale.jouy.inra.fr) for providing computational resources, and the IMG-ER system (https://img.jgi.doe.gov/cgi-bin/mer/main.cgi) for generating annotations.
N-P.P. is the recipient of a doctoral fellowship from the AgroParisTech Foundation and the INRA MEM metaprogram. This work was supported by France Génomique National infrastructure, funded as part of “Investissement d’avenir” program managed by Agence Nationale pour la Recherche (contrat ANR-10-INBS-09).
Availability of data and materials
Draft genome sequences of the 13 Brevibacterium strains sequenced in this study were deposited at the European Nucleotide Archive (http://www.ebi.ac.uk/ena). Accession numbers are provided in Table 1.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genome statistics. General features of the Brevibacterium genomes. (XLSX 15 kb)
16S–ANI. (XLSX 15 kb)
Orthology. (XLSX 6663 kb)
Orthology (Fig). Number of orthologous gene clusters in the 23 Brevibacterium genomes. (PDF 17 kb)
Galactose-galactonate-lactate-acetate-ethanol-citrate. (XLSX 36 kb)
Lipid-glycerol. (XLSX 27 kb)
Excreted proteases-aminoacids. (XLSX 43 kb)
Iron. (XLSX 53 kb)
Iron (Fig). (PDF 98 kb)
Osmotolerance. (XLSX 24 kb)
Bacteriocines. (XLSX 140 kb)
Phenazines. (XLSX 17 kb)
About this article
Cite this article
Pham, NP., Layec, S., Dugat-Bony, E. et al. Comparative genomic analysis of Brevibacterium strains: insights into key genetic determinants involved in adaptation to the cheese habitat. BMC Genomics 18, 955 (2017). https://doi.org/10.1186/s12864-017-4322-1
- Horizontal gene transfer
- Comparative genomics
- Cheese rind
- Cheese ripening
- Iron acquisition