BACKGROUND
the streptococcus bovis/streptococcus equinus complex  is a highly diverse group of bacteria that includes human and animal commensals, opportunistic pathogens and organisms that contribute to traditional food fermentations  <cit> .
the complex currently comprises the species streptococcus gallolyticus subsp.
gallolyticus , streptococcus gallolyticus subsp.
pasteurianus , streptococcus gallolyticus subsp.
macedonicus , streptococcus infantarius subsp.
infantarius , streptococcus lutetiensis , streptococcus alactolyticus and a remaining group of strains originally described as s. equinus and s. bovis which are allocated to s. equinus  <cit> .
sbsec are commensal colonizers of the human and animal gastrointestinal tract .
the carriage rate in humans is estimated at  <dig>  % in neonates in the uk  <cit>  and around 5 % in adults in france and the uk, which corresponds to the decreasing relative abundance of streptococci in the gut microbiota during aging .
sbsec are also highly prevalent among most domesticated and wild animals, including ruminants such as cattle, goats, sheep, deer and camels.
their habitat also extends to bears, piglets, rodents, dogs, sea otters and birds  <cit> .
the sbsec has been associated with a variety of diseases such as infective endocarditis , bacteremia, biliary tract and prosthetic joint infections as well as meningitis and diarrhea in humans or ruminal acidosis, bloat and laminitis in animals  <cit> .
these diseases are associated with specific species within the sbsec such as sgg for ie and bacteremia or sgp and s. lutetiensis for infant meningitis, biliary and urinary tract infections  <cit> .
in addition, s. bovis has been linked to cancer; specifically s. bovis biotype i  has been linked to colorectal cancer  in humans  <cit>  and possibly s. lutetiensis to non-colonic cancer  <cit> .
in contrast, sgm and sii play a predominant role in traditional fermented food products of animal and plant origin in southern europe, africa, asia and north america, indicating the large spectrum of roles and potential public health risks by members of the sbsec to cause diseases  <cit> .
sbsec-related factors to cause disease, their associated virulence mechanisms, infection routes, population structure and epidemiology are however not yet elucidated.
pilus proteins encoded in a pil1-operon of sgg are an important virulence factor responsible for adhesion to damaged tissues of heart valves and adenocarcinomas in the colon .
furthermore, potential pro-inflammatory proteins  <cit>  were identified in s. infantarius nctc <dig> , which suggests involvement in crc development  <cit> .
shared strain lineages between food products, animals and humans suggest a zoonotic potential and possible infection route via food and fecal-oral transmission  <cit> , which agrees with increased sbsec-related incidences in rural areas  <cit> .
however, a lack of reliable molecular epidemiological tools and marker genes hinder accurate differentiation of sbsec subspecies and their individual risk assessments, disease association and the differentiation between commensal and pathogenic subspecies or strain lineages.
this difficulty is mainly due to the high degree of gene conservation among members of the sbsec in combination with evidence of horizontal gene transfer between sbsec and other streptococci  <cit> , the ability of sbsec to be naturally competent  <cit> , rapidly changing taxonomy and a broad range of strains from potential food-grade, to commensal and pathogenic organisms fitting the emerging pathobiont concept  <cit> .
the sbsec therefore requires a holistic approach to elucidate the phylogeny, epidemiology and pathogenicity of its members in relation to their animal, human and food habitats.
multi locus sequence typing  has been used to assess global epidemiology and strain lineages based on the combined analysis of short dna sequences of housekeeping genes  <cit> .
recently, two such mlst schemes were developed for sgg on different sets of seven housekeeping genes  <cit> .
however, these mlst schemes were not designed to capture the entire sbsec.
in this study, we aimed to develop an overarching sbsec mlst scheme based on  <dig> housekeeping genes.
the novel mlst scheme was clustered with phenotypic data of strains, including their survival in simulated gastric conditions, adhesion to extracellular matrix proteins of the human intestinal epithelium and nf-kb activation capacity, all of which are important parameters to evaluate intestinal bacteria and characterize potential pathogens.
methods
bacterial strains and growth conditions
sbsec strains  were used in this study.
of these,  <dig> originated from the american type culture collection , the university of gothenburg , deutsche sammlung für mikroorganismen und zellkulturen  and donated by other researchers while  <dig> strains were analyzed in silico using genome sequences available on genbank  .
african strains were selected from our own culture collection .
these strains were isolated from fermented dairy products of cow, camel and goat origin as predominant organisms at  <dig> colony forming units  ml− <dig>  .
all enterococci and streptococci strains were cultured aerobically overnight in brain heart infusion  or m <dig> broth  at 37 °c, except streptococcus thermophilus, which was cultured in m <dig> broth  at 42–43 °c.
lactobacilli were grown anaerobically overnight in mrs tween  <dig> broth  supplemented with  <dig>  % w/w l-cysteine hydrochloride .
mrs-c agar was incubated at 37 °c for 2 days.
the purity of strains was checked by microscopy, streak plating and rep-pcr fingerprinting as previously described  <cit> .
stocks were kept in 33 %  glycerol solution at −80 °c.
all chemicals, proteins and enzymes were obtained from sigma-aldrich  unless noted otherwise.
dna isolation and pcr assays for bacterial identification
dna was isolated from single colonies on agar media via a short lysis in triton x-100-based buffer at 95 °c  <cit> .
all primers for pcr amplification and sanger sequencing were obtained from microsynth .
pcr reactions, visualization via gel electrophoresis and analyses of rep-pcr fingerprints were performed as previously described  <cit> .
sanger sequencing was performed at gatc  and microsynth .
general 16s rrna gene amplification was performed using bak11w and bak <dig> primers  <cit> .
all strains originating from our culture collection were previously identified by rep-pcr fingerprinting, 16s rrna gene assay and by a pcr assay targeting genes groes/groel for some representative strains  <cit> .
all strains from external sources were identified using a novel groel-sequencing approach, with primers groel-fw and groel-rev for amplification and sequencing that were designed based on clustalw-aligned groel sequences derived from all available sbsec genomes .
these novel primers target conserved binding sites within groel and amplify a 1167-bp product.
the groel pcr assay consisted of 2 min at 95 °c, followed by 35 cycles of 30 s at 95 °c, 30 s at 60 °c, 60 s at 72 °c and final replication for 7 min at 72 °c.
presence of lacs/lacz as marker genes for adapted lactose metabolism
lacs and lacz with high sequence identity to s. thermophilus/s.
salivarius were shown to be marker genes for dairy adaptation for african variants of sii  <cit> .
the presence of these two genes in sbsec strains was assessed and assigned to mlst clusters.
primer pairs lacz <dig> /lacz-17rev and lacs-8/lacs- <dig> rev  were used as previously described  <cit> , resulting in the amplification of approx.
1 kb of the lacz and 501 bp of the lacs gene  <cit> .
multi locus sequence typing
selection of target genes/loci
the mlst scheme candidate target genes were first selected and evaluated based on genes used in the mlst-schemes for streptococcus pneumoniae, streptococcus suis, streptococcus agalactiae, streptococcus uberis, streptococcus zooepidemicus and the salivarius group in addition to several novel candidate loci  .
thirty potential loci for initial in silico analysis from sbsec genomes were identified using clc genomic workbench .
sequences for each locus were aligned using bioedit, and compared to establish a phylogenetic tree and a sequence identity matrix using mega  <dig>   <cit> .
ten target loci were chosen on the basis of maximal nucleotide divergence  within a highly conserved sequence  to enable intra- and inter-species differentiation .
potential forward and reverse primer binding sites were selected within conserved regions and a distance of approximately 400–500 bp.
all primers were designed with a length of 18-22-bp .
all primers were also used for sanger sequencing of the amplified dna fragments.
primers for glucokinase gene amplification were optimized for s. alactolyticus using glucokinase sequences of related streptococci to develop degenerate primers.
a uniform pcr amplification protocol was used, with an initial denaturation at 95 °c for 3 min, followed by 35 cycles of 30 s each at 95 °c, 59 °c and 72 °c.
final polymerization was performed at 72 °c for 7 min.
establishment of allelic profiles and sequence types for sbsec isolates
sequencing reactions for all ten loci were performed on both strands using the same primers as for the initial pcr amplification .
sequencing chromatograms were proofread, corrected and assembled using clc genomic workbench.
whole genome sequences of sbsec isolates obtained from genbank  were processed in silico using the corresponding sbsec mlst primers to locate and extract the desired gene sequences.
sequences for each gene were then aligned in mega  <dig>  using the clustalw multiple alignment algorithm and trimmed to equal lengths.
for each locus, the sequences of all tested isolates were compared and allele numbers were assigned to each unique sequence using clc genomic workbench.
every isolate was hereby defined by a series of ten integers, constituting its allelic profile.
each unique allelic profile was subsequently assigned as the multilocus sequence type  of a strain.
the diversity within each allele was visualized using the splitstree  <dig> software package and aligned sequence files processed by uncorrectedp and neighbor-net algorithms  <cit> .
determination of genetic relationships among isolates by computational analysis
the relatedness of isolates was analyzed using the start <dig> software  <cit> .
start <dig> was further used to determine the number of polymorphic nucleotide sites, calculate dn/ds ratios and construct profile-based dendrograms using the neighbor-joining  algorithm.
sequence-based trees were constructed in mega  <dig>  using the nj algorithm and  <dig> bootstrap replications.
the classical and standardized index of association  were calculated in start <dig> using only one strain per st to avoid bias  <cit> .
additionally, the test of sawyer   <cit>  was applied to the synonymous polymorphic sites  within the alleles at each locus.
related sts were grouped using the start <dig> and eburst v <dig>  <cit>  to form clonal complexes .
ccs were defined as groups of two or more independent isolates where each isolate had identical alleles at seven or more loci with at least one other member of the group  <cit> .
this definition was also used to number phylogenetic tree clades and branches.
the allele profiles per locus were used to calculate the simpson’s index of diversity   <cit> , where a value close to  <dig> indicates high diversity and low diversity for values close to  <dig>  twice the standard deviation of the sid value was used to obtain an approximation for the 95 % confidence intervals as previously described  <cit> .
assay for simulated gastric conditions in simulated gastric juice
for simulated gastric juice  assays, strains were grown in m <dig> broth , containing lactose as carbon source to mimic the carbon source of milk.
s. alactolyticus dsm20728t, s. bovis dsm20480t and s. equinus dsm20558t were incubated in bhi broth  containing glucose instead of lactose to promote growth.
bacteria cultures were standardized to an od <dig> of  <dig>  using 5 ml of phosphate-buffered saline , which is near the optimal ph for growth of sbsec  <cit> .
the standardized bacteria suspensions were then centrifuged at 6000 × g for 6 min at 20 °c.
the pellet was resuspended in 5 ml pbs  and stored on ice for 30 min.
sgj was prepared using 1 g of nacl ,  <dig>  g of pepsin from porcine source , and  <dig>  g peptone from casein  dissolved in 200 ml distilled h2o  <cit> .
ph was adjusted to ph  <dig>  or ph  <dig>  using 5 m hcl.
these ph values were selected within the range of human gastric conditions   <cit>  allowing for adequate time-dependent survival assays over the 15 min experiment duration.
all sgj were filter-sterilized  and stored at 4 °c.
phosphate buffered saline  was prepared containing 40 g l− <dig> nacl , 1 g l− <dig> kcl ,  <dig>  g l− <dig> na2hpo <dig> ,  <dig>  g l− <dig> kh2po <dig>  in distilled h2o, and the ph was adjusted to  <dig>  using 5 m hcl for low ph-stress conditions.
pbs was autoclaved and stored at 4 °c.
the assay was performed in sterile 96-well plates  using 270 μl of sgj or pbs  per well.
the plates were pre-heated to 37 °c and inoculated with 30 μl of the standardized test culture in parallel for all conditions.
plates were incubated in air-tight plastic containers to reduce evaporation.
samples were analyzed after  <dig>   <dig>   <dig> and 15 min of incubation at 37 °c.
enumeration of surviving bacteria was done after serial dilution from 10− <dig> to 10− <dig> and microspotting of 15 μl onto agar medium  with a detection limit of  <dig>  log <dig> cfu ml− <dig>  the assays were performed with two biological replications each comprising three technical replications.
sii cj <dig> and sii ccug43820t were used as control strains in each replication.
relative survival rates per strain were calculated where the input is set at  <dig> log <dig> cfu ml− <dig> and reduction is expressed in negative log-values  <cit> .
for ease of comparison between strains at a given condition, relative survival rates of each strain were normalized using the formula /sd, where xi is the relative survival rate of a strain, xmean the mean relative survival rate of all strains for one condition and sd the standard deviation.
assay for bacterial adhesion to extracellular matrix proteins
bhi broth overnight cultures were standardized to an optical density at 600 nm  of  <dig>  using pbs  and centrifuged for 10 min at 3000 × g. the supernatant was discarded and the pellet resuspended in the original volume of pbs .
the bacteria sample was then divided into two equal aliquots and centrifugation was repeated.
the supernatant was discarded and one bacteria aliquot per culture was resuspended in pbs  and acidified pbs , respectively.
aliquots were kept on ice until use.
extracellular matrix proteins  used in this study were collagen type i  and type iv , fibrinogen and fibronectin  and the glyco-protein mucin type ii .
bovine serum albumin  was used as a control protein .
stock solutions of ecms were prepared at a concentration of 10 μg ml− <dig> in pbs  except bsa and mucin type ii at  <dig>  mg ml− <dig> in  <dig>  m tris hcl .
maxisorptm 96-well plates  were treated with 100 μl of an ecm stock solution per well  <cit> .
one 96-well plate was filled with 100 μl of pbs  as a negative adhesion control of the plastic surface.
the plates were incubated overnight at 4 °c.
subsequently, the liquid was poured off by plate inversion and dried for 10 min at 65 °c.
the plate surface was blocked with 100 μl pbs  containing 1 % tween  <dig>  incubated for 1 h at 37 °c and subsequently washed  <dig> times with 100 μl pbs  containing  <dig>  % tween  <dig>  coated plates were stored at 4 °c with 100 μl pbs  added per well.
bacteria cells  resuspended in pbs ph  <dig>  or ph  <dig>  were then added to the coated 96-well plates.
the inoculated plates were gently centrifuged for 10 min at 400 × g at room temperature to enable physical contact of the bacteria with the tested proteins, and the plates were then incubated at 37 °c for 1 h. the wells were then washed three times with 100 μl pbs  containing  <dig>  % tween  <dig>  and adhering bacteria were fixed for 20 min at 65 °c.
staining was performed using crystal violet  for 45 min at room temperature, followed by washing three times with 100 μl pbs .
bacteria-associated crystal violet was then solubilized by adding 100 μl citrate buffer  to each well and incubating for 1 h at 37 °c under constant agitation at 400–500 rpm.
the absorbance of crystal violet was measured at 595 nm .
the normalized cell adhesion of a strain  was calculated from the mean of three independent biological repetitions, using the formula /sd  <cit> , where xmedian and sd were obtained from the adhesion measurements of all strains under one condition.
values larger than  <dig> from the normalized  <dig> were considered to be adherent.
assay for inflammatory potential via nf-kb activation
the activation of nuclear factor k-light-chain-enhancer of activated b cells  was assessed using thp1-blue™ cells .
these cells express all known toll-like receptors and enable quantification of the nf-kb activation via the seap reporter gene.
thp1-blue™ cells were grown in rpmi  <dig> supplemented with 2 mm l-glutamine, 50 iu ml− <dig> penicillin,  <dig> μg ml− <dig> streptomycin and 10 % heat-inactivated fetal calf serum  in a humidified 10 % co <dig> atmosphere at 37 °c.
all culture media were supplied by life technologies .
the assay was performed in 96-well plates containing 90 μl of thp1-blue™ cells at a concentration of 2 ×  <dig> cells per well.
bacteria cells were prepared by centrifugation for 5 min at  <dig>  × g of 1-ml aliquots of bacteria cultures standardized to an od <dig> of  <dig> .
the supernatant was transferred into a new tube and both supernatant and pellet were stored at −20 °c until use.
thp1-blue™ cells were activated using 10 μl of the supernatant or resuspended pellet in 1 ml pbs.
the plates were incubated for  <dig> to 24 h. nf-kb activation was quantified by adding 10 μl of supernatant of the stimulated thp1-blue™ cells to 100 μl quanti-blue™ .
after  <dig> to 8 h of incubation at 37 °c, seap activity was quantified at 655 nm as described by the supplier.
normalization was performed as described for the adhesion assay.
statistics
statistical analysis was performed in jmp <dig>  .
data was first analyzed for normal distribution using the shapiro-wilk-w-test.
non-normal distribution of log-transformed values required the utilization of a kruskal-wallis-test, which was upon rejection of h <dig> further compared in a post-hoc pair-wise wilcoxon-test without correction for multiple comparisons  <cit>  in order to reduce type ii errors.
this was complemented by graphical analysis and detection of outliers by calculating the interquartile range iq and the defined limitations of lower fence of q1- <dig> *iq and upper fence q3 +  <dig> *iq.
strains displaying values outside of the upper or lower fence were considered outliers.
RESULTS
identification of sbsec species by groel and 16s rrna gene sequencing
partial groel sequencing was performed on all sbsec strains to confirm species and subspecies status .
all groel nucleotide sequences were trimmed to an equal length of 772 bp for comparison and neighbor-joining tree construction .
sii and s. lutetiensis had a maximum identity of  <dig> – <dig>  %, whereas none of the s. gallolyticus branch members displayed more than 91 % identity to sii.
the s. gallolyticus branch was more conserved, yielding  <dig> – <dig>  % or  <dig> – <dig>  % between sgm and sgg or sgp, respectively.
intraspecies variations were generally low, yielding  <dig> – <dig>  % sequence identity among strains of a single subspecies.
s. equinus jb <dig> showed highest identity to sii  by groel and  <dig>  % to the s. lutetiensis by 16s rrna gene sequence.
in addition, ccug <dig> , which was initially classified as s. bovis and subsequently reclassified as sii  <cit> , was again reclassified in this study as s. equinus.
similarly, s. bovis atcc <dig> was reclassified as sgp based on groel and 16s rrna gene sequences .
strain j <dig> 40– <dig> isolated from fermented milk from bangladesh was identified as sgm by groel sequencing, whereas 16s rrna gene sequence and rep-pcr fingerprinting indicated s. infantarius/s.
equinus branch .
mlst loci characterization and key performance identifiers
the mlst assay was designed as an overarching scheme to target the entire sbsec and validate species assignment .
aligned and trimmed sequences were between 393 bp and 573 bp in length .
in total,  <dig> sequence types  were defined for the  <dig> strains analyzed.
the number of alleles ranged from  <dig> for tpia to  <dig> for pyre and muts <dig>  while the number of polymorphic sites  ranged from  <dig> for tpia to  <dig> for muts .
the number of alleles and pms were lower for individual species of the sbsec .
tpia showed the highest degree of conservation whereas muts <dig> and pyre displayed the highest variability, which was also reflected in the total number of alleles defined.table  <dig> key characteristics of the sbsec-mlst scheme for the overall sbsec, sii and s. lutetiensis
index of association1
ia
s
ddla
gki
glna
muts
muts2
phes
pros
pyre
thrs
tpia
1index of association based on one strain per st: ia: calculated using classical maynard smith approach, ia
s: standardized according to haubold
2sbsec: includes also 3 s. gallolyticus subsp.
pasteurianus and 1 s. alactolyticus strains not listed separately
3
a allele
4
pms polymorphic sites
5sid : simpson’s index of diversity with 95 % confidence interval  calculated from sid+/− 2*standard deviation
6no dn/ds calculation due to limited allele numbers and differences
values were calculated for the entire sbsec and individual species, including all  <dig> target genes and simpson’s index of diversity  as an indication for the differentiation power among strains
index of association1
ia
s
ddla
gki
glna
muts
muts2
phes
pros
pyre
thrs
tpia
1index of association based on one strain per st: ia: calculated using classical maynard smith approach, ia
s: standardized according to haubold
2sbsec: includes also  <dig> s. gallolyticus subsp.
pasteurianus and  <dig> s. alactolyticus strains not listed separately
3
a allele
4
pms polymorphic sites
5sid : simpson’s index of diversity with 95 % confidence interval  calculated from sid+/− 2*standard deviation
6no dn/ds calculation due to limited allele numbers and differences
7nan: maximum diversity, no ci calculcated.
values were calculated for the entire sbsec and individual species, including all  <dig> target genes and simpson’s index of diversity  as an indication for the differentiation power among strains
splitstree analysis of each individual locus provided network-like structures for all ten loci .
least-square tree fitting resulted in values > <dig>  for all trees.
within each splitstree network, separate clusters were distinguishable for the s. gallolyticus and s. infantarius branch, with s. alactolyticus as the most distant species.
the index of association  and its standardized version  for the overall sbsec were  <dig>  and  <dig> , respectively .
ias values calculated for individual groups were  <dig>  ,  <dig>  ,  <dig>  ,  <dig>  ,  <dig>  , and  <dig>  , while sgp comprised too few strains for calculation.
except for sgm, all other groups including sbsec displayed significant linkage disequilibrium, suggesting recombinatory evolution.
simpson’s index of diversity , which indicates the discrimination power of the method, was  <dig>   for the overall sbsec approach and  <dig>   for sii .
definition of clonal complexes
clonal complexes  are an important classification tool for mlst schemes.
ccs were defined for strains that share  <dig> out of  <dig> loci using eburst v <dig>  and these were depicted in a profile-based neighbor-joining tree .
a total of  <dig> ccs were determined, which grouped  <dig> strains and left  <dig> singletons.
cc <dig> and cc <dig> were both composed of west african dairy strains, whereas kenyan and somali strains were both grouped in cc <dig> and cc <dig> .
cc <dig> represented the main human pathogen strain cluster that delineated in close relationship with ivorian dairy cluster cc <dig>  small ccs comprising two strains each were found for sgp  and s. lutetiensis .
cc <dig> comprised most of the dairy sgm strains from africa and greece, suggesting close relationships among strains of this subspecies.
s. bovis  and s. equinus  cannot yet be considered as real ccs, as each cc contained only two versions of the corresponding type strain obtained from different sources.
the number of strains did not allow the definition of ccs for sgg or the prediction of a founding strain per cc.
the defined ccs 1– <dig> suggest clear regional relationships among sii dairy strains  and even across continents , but also pathogenic relationships concentrated in cc <dig>
mlst-based phylogeny of sbsec and individual species groups
the sbsec mlst scheme provides the first detailed overview of all species within the sbsec .
all branches 1– <dig> were numbered sequentially for ease of comparison using the cc definition of  <dig> or more identical alleles out of  <dig> fig.
<dig> sequence-based phylogenetic mlst tree of the sbsec.
the mlst tree comprised commensal animal and human strains, human pathogenic strains and food-derived strains of all sbsec species.
the tree was calculated from the concatenated sequences of the  <dig> mlst loci of each strain using the neighbor joining algorithm and  <dig> bootstrap replications.
the tree was rooted to s. alactolyticus dsm20728t.
corresponding sts, alleles and prevalence of dairy adaptation marker genes , origin and isolation source are indicated.
branches were defined and numbered according to the clonal complex specifications
primary differentiation established an s. gallolyticus clade  and an s. equinus/s.
lutetiensis/sii clade .
within each clade, species-specific clades were delineated.
the s. equinus clade  harbored both the s. equinus and s. bovis type strains that originally defined sbsec .
s. equinus/s.bovis strains were distributed between a ruminal/human clade  and an equine branch .
s. equinus jb <dig>  a ruminant isolate, delineated a separate branch  with closer association with the s. lutetiensis clade than the major s. equinus branch, suggesting reclassification of this strain based on mlst analysis.
the s. lutetiensis clade reflected a major division between human clade  and a potential side branch  comprising an animal strain .
the human clade was further divided into human blood and fecal isolates.
the s. gallolyticus clade showed subdivision into a clade comprising sgg and sgp, and a separate clade of sgm .
within the sgg/sgp clade, separate clades were delineated for human pathogenic strains of sgg  and sgp .
animal commensal strains seemed more heterogeneous, forming two separate branches, one for the sgg type strain dsm16831t  and one even further related for the calf isolate lmg <dig> , indicating different lineages of commensal and pathogenic origins.
among the sgm strains, the divisions between greek , east african  and west african  strains were fully reflected, suggesting a main african lineage that was further differentiated into east and west african lineages.
the sii clade was analyzed in more detail .
several dairy and human clades were delineated, suggesting a general trend of different lineages.
strains of human clinical cases were concentrated in only two specific clades comprising branches  <dig> – <dig>  and branch  <dig> regardless of geographic origin.
the majority of dairy strains were found in two dairy clades comprising  east african dairy isolates  and  west african and asian dairy isolates , suggesting different lineages within the main dairy lineage of sii.
east african dairy strains formed a single clade subdivided into two main clades of dairy isolates only  and mixed human and dairy isolates comprising the sii type strain .
the majority of west african and the asian dairy isolates were grouped in a single dairy clade comprising clades  <dig> and  <dig>  the delineation of clades  <dig> and  <dig> was directly linked to thrs alleles  <dig>   <dig> and  <dig> featuring highest sequence identities of  <dig> – <dig>  % to s. thermophilus instead of other sii or sbsec members as for the other  <dig> alleles.
however, west african dairy strains were more heterogeneous than those of east africa , and also featured in addition to clade  <dig> several shared or related clades with human blood isolates  and east african dairy isolates .
prevalence of lacs/lacz genes as marker genes for dairy adaptation
the prevalence of s. thermophilus-like lacs/lacz genes was assessed as marker genes for dairy adaptations among all sbsec, including sii strains, using specific pcr assays and previously obtained data  <cit> .
only sii strains were found to harbor lacs/lacz, with highest prevalence among the major dairy branches  <dig> ,  <dig> ,  <dig>  and  <dig>  .
in contrast, dairy strains more closely related to human strains, such as sii 150a  and the four west african dairy strains found in branches  <dig> and  <dig>  did not harbor the lacs/lacz genes.
resistance to simulated gastric juice
the resistance of sbsec strains  to sgj was tested in sgj at ph  <dig>  or ph  <dig>  and pbs at ph  <dig>  as control .
all differences in this section were tested for significance at p <  <dig>  unless mentioned otherwise.
generally, relative cell survival was significantly different for the three conditions ph  <dig>  , ph  <dig>   and pbs ph  <dig>   over 15 min.
significantly different survival was also observed among strains at ph  <dig>  or ph  <dig>   whereas no significant difference was observed among the strains tested under control conditions of pbs ph <dig>   .
significantly different survival rates were determined between dairy sii clades  <dig>   <dig> and  <dig> featuring higher survival compared to sii commensal clade  <dig> and predominantly blood isolate clade  <dig>  s. lutetiensis clade  <dig>  s. gallolyticus clades  <dig> and  <dig> as well as s. alactolyticus branch  <dig> .
survival rates were significantly dependent on time-ph combinations.
conditions of ph  <dig>  yielded low survival rates, often below detection limit after 15 min .
this strong inhibition resulted in a few significant differences between strains and mlst branches with the exception of outliers in clade  <dig> and branch  <dig> .
significantly highest relative survival in comparison to all other strains was determined for the four west african dairy sii strains 6by-11b, 13af, 12ay- <dig> and p2vc1- <dig>  clustered in mlst branches  <dig>  and  <dig>  as well as the human blood isolate s. lutetiensis  <dig>  and all s. thermophilus strains .
significantly higher relative survival was therefore determined for west african  compared to east african sii strains in branches  <dig> and  <dig>  no significant correlation was observed between lacs/lacz prevalence, habitat and relative survival under the conditions tested .
adhesion to extracellular matrix proteins
adhesion of sbsec strains to various ecms was investigated to estimate their interaction with cell surfaces of the gi tract and blood vessels.
due to heavy tailed distributed values, all adhesion values were normalized and analyzed using a median-based calculation .fig.
<dig> analysis of the adhesion ability of sbsec members to extracellular matrix proteins and activation of the nf-kb signaling pathway.
adhesion and nf-kb activation values of individual strains displayed in a heat map.
the phylogenetic tree is based on the mlst-tree but not drawn to scale for better visualization.
adhesion measurements were normalized for each strain based on the adhesion median of all strains.
values larger than  <dig> sd from the normalized  <dig> adhesion were considered to be adherent.
outliers were marked using an asterix*
lb.
plantarum wcfs <dig> and lb.
reuteri dsm20016t were used as positive adhesion controls showing significant adhesion towards fibronectin, fibrinogen, collagen i and mucin type ii that are known factors contributing to their persistence in the human gi tract.
as expected, lb.
plantarum nz <dig>  the sortase knockout mutant of wcsf <dig> depleted of its adhering surface proteins, showed no adhesion and served as negative adhesion control.
adhesion values obtained at experimental conditions of ph  <dig>  showed better reproducibility in contrast to those values obtained at ph  <dig>  and were therefore chosen for analysis.
adhesion characteristics were separated by mlst clustering and isolate origin.
generally, adhesion was most prevalent to the two ecms collagen i and mucin type ii, with  <dig>  and  <dig>  % of sbsec isolates showing adhesion , respectively, followed by adhesion to fibronectin, fibrinogen and collagen iv by  <dig>  % of sbsec isolates.
among human sii, adhesion to mucin type ii and unspecific adhesion was observed for  <dig>  % and adhesion to all remaining ecms for 33 % of human sii strains .
among these human strains, comparison of blood isolates  and potentially commensal isolates  showed that significant adhesion was exclusively present among sii blood isolates comprised mainly in branches  <dig>  and  <dig>  human fecal isolate sii ccug43820t featured no significant adhesion abilities to fibronectin, fibrinogen or collagen i and only weak adhesion to collagen iv and mucin type ii.
similarly, s. lutetiensis blood isolate  <dig>  displayed significant wide range adhesion abilities whereas human fecal isolates of s. lutetiensis  displayed no significant adhesion abilities, suggesting low prerequisite or other mechanisms contributing to the persistence of these possibly commensal sii and s. lutetiensis strains in the human gi tract .
similarly, dairy sii  generally had a low prevalence of adhesion abilities to the ecms tested and those dairy sii isolates showing adhesion, featured mainly adhesion to collagen iv and mucin type ii detected in  <dig>  and  <dig>  % of strains, respectively, with the exception of cji <dig>  cj <dig>  cj <dig> and ab2vb2- <dig> showing also adhesion to collagen i or mucin type ii.
clustered according to mlst and geography, west african dairy sii of mlst branches  <dig> ,  <dig>  and  <dig> displayed low adhesion abilities to the ecms tested.
adhesion was more prevalent among east african dairy sii strains with significant adhesion to at least one ecm for strains cji <dig>  cj <dig>  cj <dig> and cj <dig> as well as west african dairy sii strain ab2vb2- <dig>  which all share the same mlst clade 1– <dig> .
the significantly adherent west african dairy strain sii av2a <dig> was the exception to cluster with the otherwise low-adherent strains of the west african main dairy clade  <dig> .
among the s. gallolyticus clade 14– <dig>  high adhesion abilities towards collagen i, iv and fibrinogen was observed.
dairy sgm branches of east africa  and greece  showed adhesion towards mucin type ii and fibronectin.
in contrast, sgp  did not display any significant adhesion abilities allowing clear separation by mlst clustering from sgg .
determination of immune stimulatory potential by nf-kb activation
activation of the transcription factor nf-kb was used as marker for the immune stimulation of the sbsec strains studied  using the thp1-blue reporter cell line that expresses all tlrs that detect bacterial motifs in the intestine.
nf-kb activation was generally low for the strain panel tested with significantly higher values obtained for cell pellets  in contrast to cell-free supernatants .
cell pellets, significant inflammation was detected only for individual strains such as sii ccug <dig> , s. equinus ccug <dig>  and sgg dsm 16831t .
the majority of dairy isolates and also human commensal and blood isolates did not significantly activate nf-kb, particularly in comparison to control strain e. faecalis jh2- <dig>  among cell-free supernatants, highest nf-kb activation was measured for dairy isolates sii ab2vb2- <dig> , sii cj <dig>  and sgm cji <dig> , as well as human blood isolates sii  <dig>  and sii ccug <dig> , displayed significant nf-kb activation.
again, no significant nf-kb activation was detected among supernatants of other human commensal strains of s. lutetiensis, sii and even pathogenic s. lutetiensis, sii, sgg and sgp.
discussion
the sbsec is a highly diverse group of bacteria colonizing multiple habitats such as food, human and animal gi tracts.
although the 16s rrna gene sequence provides sufficient power to differentiate the main branches of s. gallolyticus, s. infantarius and s. alactolyticus, differentiation at subspecies level requires sequencing of additional genes such as groel  <cit> .
we showed that groel provided stringent sequence data and discrimination power for sbsec subspecies identification.
interestingly, sii j <dig> 40- <dig> seems to be an exception, 16s rrna gene sequence suggested an s. infantarius branch member while groel indicates highest similarity with sgm; however, the identity in the sbsec mlst scheme was confirmed to be sii.
the ongoing horizontal gene transfer  among streptococci and sbsec  <cit>  in combination with natural competence  <cit>  may have changed even genes considered as housekeeping genes.
further indications for hgt between dairy sii and s. thermophilus were observed among thrs alleles of strains comprised in mlst clades  <dig> and  <dig> representing asian dairy strain sii j <dig> 40- <dig> and the major west african clade.
most strains in clade  <dig> harbor lacs and lacz which serve as marker genes for s. thermophilus-like dairy adaptation.
in combination with the thrs alleles present, this suggests a common ancestor for clades  <dig> and  <dig> that underwent multiple hgt with s. thermophilus.
ultimately, these observations of hgt and natural competence support the need for a multi gene approach such as mlst for sbsec species to ensure accurate identification as well as inter- and intraspecies phylogenetics.
the sbsec mlst scheme developed in this study provides the first multi-gene-based comparison across the entire sbsec.
in contrast to previous sgg mlst schemes  <cit> , this scheme provides a tool for clustering all currently described sbsec species and subspecies from human, animal and food sources.
the sid of  <dig> – <dig>  is within the range and confidence interval of that of the sgg mlst   <cit> .
therefore, our novel sbsec mlst assay provides equally high differentiation power as the sgg-specific assays without the restriction of being limited to a single species but instead being applicable to all species of the sbsec.
furthermore, extraction and phylogenetic comparison of mlst alleles from whole genome sequences as performed for several strains in this study ensure the application of this sbsec mlst scheme also in combination with next generation sequencing of bacterial genomes as data source instead of traditional sanger sequencing.
the sbsec mlst scheme provided highly comparable phylogenetic tree delineation for sgg as those obtained via the sgg-specific mlst scheme.
dn/ds ratios of all loci in the sbsec mlst scheme were < <dig> and comparable to those of the two sgg-specific schemes  <cit> .
number and percentage of pms of most loci were in the range of dumke et al.
at  <dig> – <dig>  % and shibata et al.
at  <dig> – <dig>  %.
the ia of  <dig> – <dig>  for individual species was slightly lower than the ia of  <dig>  determined by dumke et al.
whereas shibata et al.
did not provide a value for ia.
the ia values for the overall sbsec, s. equinus, sgg, sii and s. lutetiensis indicated significant linkage disequilibrium between the  <dig> selected loci as previously observed for the loci set of sgg  <cit> .
linkage disequilibrium is the effect of non-random association of the selected alleles in an mlst scheme  <cit> , which indicates that recombination was a factor in the evolution of the sbsec as seen also by network-like structures in splitstree analyses   <cit> .
sgm was the only subspecies where linkage disequilibrium was not detected.
sgm also displayed the lowest genetic diversity despite the analysis of greek, somali, ivorian and mali strains, suggesting that sgm strains from these different origins are highly related.
this could support the theory that sgm is actually a subspecies of dairy origin within the sbsec  <cit> , enabled through clear delineation from other sbsec members in this mlst scheme.
genome sequencing of sgm strains will be required to further assess their evolution in relation to the presence of putative sgm dairy adaptation markers such as duplicated lacg <dig> described for sgm aca-dc  <dig> and decayed or active virulence factors  <cit> .
sgg and s. lutetiensis showed a clear differentiation between presumptive human commensal lineages and those of infectious nature isolated from feces, blood or heart samples of humans and animals.
in particular, commensal strains such as koala bear isolate sgg dsm16831t and calf isolate sgg lmg <dig> were clearly distinct from human pathogenic lineages, suggesting the existence of potentially commensal and pathogenic lineages among many sbsec species.
this requires evaluation in a one health approach, encompassing animal and human health, food, the zoonotic potential of sbsec members and their status as pathobionts  <cit> .
the subspecies sii that originated from human and animal-derived food sources had the highest representation in this mlst assay.
clear regional differentiation was observed for the majority of east and west african dairy strains.
differentiation by host and source was observed for most human and dairy strains, suggesting different lineages depending on fecal or blood origin.
the mlst phylogeny suggests that human fecal and especially blood isolates cluster by isolation source , and not by geographical origin in contrast to most dairy strains.
this possibly reflects different adaptation and risk levels.
dairy strains clustered in distinct lineages featuring a prominent west african and asian main branch and an east african branch.
the main lineage of east african dairy strains displayed remarkable conservation in contrast to the higher diversity among west african strains and especially those collected in côte d’ivoire.
the regions of mandera and garissa in kenya share a common cultural background with somalia.
large somali populations are represented in both kenyan regions, with a long tradition in livestock trading, mobile pastoral lifestyle and migration across borders  <cit> , possibly explaining the spreading of similar strain lineages.
certain lineages among the west african dairy sii without the marker genes lacs/lacz clustered separately from the main dairy lineages and closer to human fecal and pathogenic isolates, suggesting the existence of overlapping lineages of sii or contamination of foods by fecal material.
this raises concerns about potential health risks of these strain lineages and the origin of these lineages in food.
exposure of strains to sgj was used to assess the survival of bacteria in stomach conditions and delivery potential to the colon.
generally, strains of mlst clusters showed similar relative survival being affected by sgj.
however, at ph  <dig>  and ph  <dig>  several dairy sii clades and mainly four dairy strains showed significantly higher relative survival than human blood sii isolates/clades, respectively, suggesting different adaptation to acidic stress in relation to strain origin and mlst clustering.
mlst-based clustering was further used to assess the ability of individual strains to adhere to ecm and activate the nf-kb pathway in correlation with mlst clades and phylogeny.
adherence is an important first step in colonization of the gi tract or damaged body tissues such as heart valves.
mlst clusters provided clear differentiation of highly adhesive sii/s.
lutetiensis blood isolates, the s. gallolyticus clade, less adhesive east african dairy strains and nearly non-adhesive west african strains, suggesting a generally lower tendency among west african dairy clades and commensal s. lutetiensis to adhere to the ecms tested.
the genetic factors contributing to adhesion were so far only related to a specific pilus loci for collagen binding termed pil <dig>  pil <dig> and pil <dig> present in the sgg genomes analyzed  <cit> .
only pil <dig> was found in sii, sgm and sgp, whereas the crucial pil <dig> factor of sgg was absent  <cit> .
comparative genomics and further functional analyses are needed to determine the responsible factors in highly adhering strains of sii, s. lutetiensis, sgm and sgg, and possibly identify new virulence factors to enable an in-depth risk assessment for the different sbsec species in relation to human and animal diseases.
the nf-kb activation by bacteria attaching to the intestinal mucosa plays a central role in immune signaling of the host, particularly in the intestinal tract.
for most strains, nf-kb activation was positively correlated with significant adhesion abilities, suggesting health implications and the need for further functional analyses for such strains.
several blood isolates as well as pathogenic sgp dsm 15351t did not yield significant nf-kb activation, which could enable them to infect the host without inducing a strong immune response  <cit> .
CONCLUSIONS
the overarching mlst scheme for the entire sbsec provides an epidemiologic and diagnostic tool enabling in-depth inter- and intraspecies comparisons of all species within this emerging human and animal pathogenic complex.
the mlst scheme based on  <dig> housekeeping genes allowed delineating all species within the sbsec.
it enabled division by host and habitat and furthermore possibly pathogenic lineages from commensal and food lineages.
this approach will be useful in enhancing taxonomic classifications of human, animal and food isolates as well as their relationships.
our data suggest that the parallel evolution of east and west african dairy lineages is currently in the process of differentiating from the human lineage including genetic exchange with other streptococci.
however, enhanced adhesion abilities of certain lineages similar to those of sgg and the ability to cause inflammation or evade immune response are certainly of concern and require further investigations into pathogenicity mechanisms shared among the sbsec and those unique to only specific species.
this will greatly assist in designing appropriate treatment strategies for patients infected with this emerging pathogenic group depending on mlst clustering.
furthermore, accurately collecting clinically relevant epidemiologic data in relation to transmission, zoonosis, host specificity and identification of virulence mechanisms will greatly contribute to an enhanced risk assessment of pathogenic, commensal and food lineages or the potential role as pathobionts in an overarching one health approach.
abbreviations
cc, clonal complex; crc, colorectal cancer; ecm, extracellular matrix protein; gi tract, gastrointestinal tract; hgt, horizontal gene transfer; ia, index of association; ias, standardized index of association; ie, infective endocarditis; mlst, multi locus sequence typing; pbs, phosphate-buffered saline; pms, polymorphic sites; sbsec, streptococcus bovis/streptococcus equinus complex; sgg, streptococcus gallolyticus subsp.
gallolyticus; sgj, simulated gastric juice; sgm, streptococcus gallolyticus subsp.
macedonicus; sgp, streptococcus gallolyticus subsp.
pasteurianus; sii, streptococcus infantarius subsp.
infantarius; sid, simpson’s index of diversity; st, sequence type
additional files
additional file 1: title of data: strains used in this study.
description of data: strains used in this study.
additional file 2: title of data: primers developed and used in this study.
description of data: primers developed and used in this study for the sii/sbsec mlst scheme of housekeeping genes and identification based on groel, lacs and lacz genes.
additional file 3: title of data: mlst candidate loci used in this study.
description of data: all  <dig> candidate loci which were used in this study were selected from literature and mlst.net .
conserved housekeeping genes were selected on the basis of maximal nucleotide divergence for intraspecies differentiation.
all loci which are marked with an asterisk were used for the final sbsec mlst scheme.
additional file 4: title of data: sequence identity matrix of 772-bp groel fragments of sbsec species.
description of data: sequence identity matrix of 772-bp groel fragments of sbsec species
additional file 5: title of data: groel-based phylogenetic tree of sbsec strains  of this study.
description of data: phylogenetic tree of all  <dig> sbsec strains investigated this study based on neighbor-joining  <dig> calculation of a 772-bp sequence fragment of the groel gene.
the scale bar below the tree  <dig> indicates the evolutionary distance using the number of base substitutions per site as units.
<dig> species abbreviations: s. infantarius subsp.
infantarius , s. gallolyticus subsp.
gallolyticus  <dig> , s, gallolyticus subsp.
macedonicus , s. gallolyticus subsp.
pasteurianus .
additional file 6: title of data: 16s rrna gene-based phylogenetic tree of selected sbsec members.
description of data: phylogenetic tree of selected sbsec strains  of all sbsec species based on a 930-bp fragment of the 16s rrna gene.
calculations were performed using the neighbor-joining algorithm.
the scale bar below the tree indicates the evolutionary distance using the number of base substitutions per site as units.
species abbreviations: s. infantarius subsp.
infantarius , s. gallolyticus subsp.
gallolyticus , s, gallolyticus subsp.
macedonicus , s. gallolyticus subsp.
pasteurianus .
additional file 7: title of data: splitstree visualization of the  <dig> mlst loci.
description of data: the trees were obtained using the defined allele sequences and calculation by neighbor-joining algorithm.
corresponding alleles of selected reference strains are highlighted: sii ccug43820t = bright green, sgg dsm16831t = red, sgm dsm15879t = ice blue, s. lutetiensis ccug43822 = yellow, sgp dsm15351t = brown, s. bovis dsm20480t = pink, s. equinus dsm20558t=dark blue, s. alactolyticus dsm20728t = orange.
additional file 8: title of data: profile-based mlst tree of all species within the sbsec.
description of data: the strain pool comprises commensal animal and human strains, human pathogenic strains and food-derived strains.
the tree was calculated by neighbor joining algorithm using the allele profiles of  <dig> housekeeping genes.
corresponding sequence types , alleles and prevalence of dairy adaptation marker genes  including origin and isolation source are indicated.
species abbreviations: s. infantarius subsp.
infantarius , s. gallolyticus subsp.
gallolyticus , s, gallolyticus subsp.
macedonicus , s. gallolyticus subsp.
pasteurianus .
additional file 9: title of data: sequence-based phylogenetic mlst tree of s. infantarius subsp.
infantarius .
description of data: the mlst tree comprised commensal animal and human strains, human pathogenic strains and food-derived strains of all sbsec species.
the tree was calculated by neighbor joining algorithm using the concatenated partial sequences of  <dig> housekeeping genes.
trees were rooted to s. equinus jb <dig>  corresponding sequence types , alleles and prevalence of dairy adaptation marker genes  including origin and isolation source are indicated.
branches were defined and numbered according to the clonal complex specifications.
additional file 10: title of data: relative survival of sbsec members under simulated gastric conditions at ph  <dig> .
description of data: relative survival of sbsec strains after  <dig>   <dig> and 15 min incubation in sgj at ph  <dig>  performed in two biological replications.
relative survival values in log <dig> cfu ml− <dig> were normalized according to the mean of all measurements of one condition to allow comparison between strains.
positive normalized values indicate relative survival higher than the mean of all strains.
statistical distribution indicators: median: large black dotted line; q <dig> and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line .
the phylogenetic tree is based on the mlst-tree but not drawn to scale.
additional file 11: title of data: relative survival of sbsec members under simulated gastric conditions in pbs ph  <dig> .
description of data: relative survival of sbsec strains at  <dig>   <dig> and 15 min in pbs ph  <dig>  performed in  <dig> biological replications.
input of  <dig> log <dig> cfu ml− <dig> was used as basis from which relative survival is expressed.
relative survival values were normalized to allow comparison between strains where positive values indicate relative survival higher than the mean of all strains.
statistical distribution indicators: median: large black dotted line; q <dig> and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line .
the phylogenetic tree is based on the mlst-tree but not drawn to scale.
additional file 12: title of data: relative survival of sbsec members under simulated gastric conditions at ph  <dig> .
description of data: relative survival of sbsec strains at  <dig>   <dig> and 15 min in in simulated gastric conditions at ph  <dig>  performed in  <dig> biological replications.
input of  <dig> log <dig> cfu ml− <dig> was used as basis from which relative survival is expressed.
relative survival values were normalized to allow comparison between strains where positive values indicate relative survival higher than the mean of all strains.
statistical distribution indicators: median: large black dotted line; q <dig> and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line .
the phylogenetic tree is based on the mlst-tree but not drawn to scale.
the authors greatly acknowledge the sbsec strains received from ingolf nes , tsilia lazarovitch , taku miyamoto , rosa del campo/beatriz romero hernández , marisa haenni/jean-yves madec .
the authors thank peter bron and michiel kleerebezem  for the gift of lactobacillus plantarum wcfs <dig> and lactobacillus plantarum nz <dig>
this study and all authors were funded by the ubs optimus foundation, switzerland and the erafrica program theme “new ideas”.
funding
this project was financed by the ubs optimus foundation, switzerland and erafrica “new ideas”.
the funders had no role in study design, data collection, analysis, interpretation, decision to publish, or preparation of the manuscript.
all authors have read and approved the final manuscript.
availability of data and material
the data set supporting the results of this article are available in the genbank and pubmlst.org repositories under accession number ranges  and  and sbsec mlst database   <cit> , respectively.
authors’ contributions
cj, dwmk, ja, db, sv, ts and ji performed the experiments to type bacteria isolates, develop and validate the mlst assay and contributed to data analysis and interpretation.
ja and tw conducted adhesion and inflammation assays.
fk conducted gastric survival assays.
ja, tw and fk contributed to data analysis and interpretation of the respective assays.
cj, tw, cl and lm designed and supervised the experiments.
cj, tw, bb, cl, dwmk, wkm and lm performed analyses and interpretation, and wrote and revised the manuscript.
all authors read and approved the final manuscript.
competing interests
the authors declare that they have no competing interests.
consent for publication
not applicable.
ethical approval and consent to participate
not applicable.

 the streptococcus bovis/streptococcus equinus complex is a highly diverse group of bacteria that includes human and animal commensals, opportunistic pathogens and organisms that contribute to traditional food fermentations. the complex currently comprises the species streptococcus gallolyticus subsp. gallolyticus, streptococcus gallolyticus subsp. pasteurianus, streptococcus gallolyticus subsp. macedonicus, streptococcus infantarius subsp. infantarius, streptococcus lutetiensis, streptococcus alactolyticus and a remaining group of strains originally described as s. equinus and s. bovis which are allocated to s. equinus. sbsec are commensal colonizers of the human and animal gastrointestinal tract. the carriage rate in humans is estimated at  % in neonates in the uk and around 5 % in adults in france and the uk, which corresponds to the decreasing relative abundance of streptococci in the gut microbiota during aging. sbsec are also highly prevalent among most domesticated and wild animals, including ruminants such as cattle, goats, sheep, deer and camels. their habitat also extends to bears, piglets, rodents, dogs, sea otters and birds. the sbsec has been associated with a variety of diseases such as infective endocarditis, bacteremia, biliary tract and prosthetic joint infections as well as meningitis and diarrhea in humans or ruminal acidosis, bloat and laminitis in animals. these diseases are associated with specific species within the sbsec such as sgg for ie and bacteremia or sgp and s. lutetiensis for infant meningitis, biliary and urinary tract infections. in addition, s. bovis has been linked to cancer; specifically s. bovis biotype i has been linked to colorectal cancer in humans and possibly s. lutetiensis to non-colonic cancer. in contrast, sgm and sii play a predominant role in traditional fermented food products of animal and plant origin in southern europe, africa, asia and north america, indicating the large spectrum of roles and potential public health risks by members of the sbsec to cause diseases. sbsec-related factors to cause disease, their associated virulence mechanisms, infection routes, population structure and epidemiology are however not yet elucidated. pilus proteins encoded in a pil1-operon of sgg are an important virulence factor responsible for adhesion to damaged tissues of heart valves and adenocarcinomas in the colon. furthermore, potential pro-inflammatory proteins were identified in s. infantarius nctc, which suggests involvement in crc development. shared strain lineages between food products, animals and humans suggest a zoonotic potential and possible infection route via food and fecal-oral transmission, which agrees with increased sbsec-related incidences in rural areas. however, a lack of reliable molecular epidemiological tools and marker genes hinder accurate differentiation of sbsec subspecies and their individual risk assessments, disease association and the differentiation between commensal and pathogenic subspecies or strain lineages. this difficulty is mainly due to the high degree of gene conservation among members of the sbsec in combination with evidence of horizontal gene transfer between sbsec and other streptococci, the ability of sbsec to be naturally competent, rapidly changing taxonomy and a broad range of strains from potential food-grade, to commensal and pathogenic organisms fitting the emerging pathobiont concept. the sbsec therefore requires a holistic approach to elucidate the phylogeny, epidemiology and pathogenicity of its members in relation to their animal, human and food habitats. multi locus sequence typing has been used to assess global epidemiology and strain lineages based on the combined analysis of short dna sequences of housekeeping genes. recently, two such mlst schemes were developed for sgg on different sets of seven housekeeping genes. however, these mlst schemes were not designed to capture the entire sbsec. in this study, we aimed to develop an overarching sbsec mlst scheme based on housekeeping genes. the novel mlst scheme was clustered with phenotypic data of strains, including their survival in simulated gastric conditions, adhesion to extracellular matrix proteins of the human intestinal epithelium and nf-kb activation capacity, all of which are important parameters to evaluate intestinal bacteria and characterize potential pathogens. methods bacterial strains and growth conditions sbsec strains were used in this study. of these, originated from the american type culture collection, the university of gothenburg, deutsche sammlung für mikroorganismen und zellkulturen and donated by other researchers while strains were analyzed in silico using genome sequences available on genbank. african strains were selected from our own culture collection. these strains were isolated from fermented dairy products of cow, camel and goat origin as predominant organisms at colony forming units ml−. all enterococci and streptococci strains were cultured aerobically overnight in brain heart infusion or m broth at 37 °c, except streptococcus thermophilus, which was cultured in m broth at 42–43 °c. lactobacilli were grown anaerobically overnight in mrs tween broth supplemented with  % w/w l-cysteine hydrochloride. mrs-c agar was incubated at 37 °c for 2 days. the purity of strains was checked by microscopy, streak plating and rep-pcr fingerprinting as previously described. stocks were kept in 33 % glycerol solution at −80 °c. all chemicals, proteins and enzymes were obtained from sigma-aldrich unless noted otherwise. dna isolation and pcr assays for bacterial identification dna was isolated from single colonies on agar media via a short lysis in triton x-100-based buffer at 95 °c. all primers for pcr amplification and sanger sequencing were obtained from microsynth. pcr reactions, visualization via gel electrophoresis and analyses of rep-pcr fingerprints were performed as previously described. sanger sequencing was performed at gatc and microsynth. general 16s rrna gene amplification was performed using bak11w and bak primers. all strains originating from our culture collection were previously identified by rep-pcr fingerprinting, 16s rrna gene assay and by a pcr assay targeting genes groes/groel for some representative strains. all strains from external sources were identified using a novel groel-sequencing approach, with primers groel-fw and groel-rev for amplification and sequencing that were designed based on clustalw-aligned groel sequences derived from all available sbsec genomes. these novel primers target conserved binding sites within groel and amplify a 1167-bp product. the groel pcr assay consisted of 2 min at 95 °c, followed by 35 cycles of 30 s at 95 °c, 30 s at 60 °c, 60 s at 72 °c and final replication for 7 min at 72 °c. presence of lacs/lacz as marker genes for adapted lactose metabolism lacs and lacz with high sequence identity to s. thermophilus/s. salivarius were shown to be marker genes for dairy adaptation for african variants of sii. the presence of these two genes in sbsec strains was assessed and assigned to mlst clusters. primer pairs lacz /lacz-17rev and lacs-8/lacs- rev were used as previously described, resulting in the amplification of approx. 1 kb of the lacz and 501 bp of the lacs gene. multi locus sequence typing selection of target genes/loci the mlst scheme candidate target genes were first selected and evaluated based on genes used in the mlst-schemes for streptococcus pneumoniae, streptococcus suis, streptococcus agalactiae, streptococcus uberis, streptococcus zooepidemicus and the salivarius group in addition to several novel candidate loci. thirty potential loci for initial in silico analysis from sbsec genomes were identified using clc genomic workbench. sequences for each locus were aligned using bioedit, and compared to establish a phylogenetic tree and a sequence identity matrix using mega. ten target loci were chosen on the basis of maximal nucleotide divergence within a highly conserved sequence to enable intra- and inter-species differentiation. potential forward and reverse primer binding sites were selected within conserved regions and a distance of approximately 400–500 bp. all primers were designed with a length of 18-22-bp. all primers were also used for sanger sequencing of the amplified dna fragments. primers for glucokinase gene amplification were optimized for s. alactolyticus using glucokinase sequences of related streptococci to develop degenerate primers. a uniform pcr amplification protocol was used, with an initial denaturation at 95 °c for 3 min, followed by 35 cycles of 30 s each at 95 °c, 59 °c and 72 °c. final polymerization was performed at 72 °c for 7 min. establishment of allelic profiles and sequence types for sbsec isolates sequencing reactions for all ten loci were performed on both strands using the same primers as for the initial pcr amplification. sequencing chromatograms were proofread, corrected and assembled using clc genomic workbench. whole genome sequences of sbsec isolates obtained from genbank were processed in silico using the corresponding sbsec mlst primers to locate and extract the desired gene sequences. sequences for each gene were then aligned in mega using the clustalw multiple alignment algorithm and trimmed to equal lengths. for each locus, the sequences of all tested isolates were compared and allele numbers were assigned to each unique sequence using clc genomic workbench. every isolate was hereby defined by a series of ten integers, constituting its allelic profile. each unique allelic profile was subsequently assigned as the multilocus sequence type of a strain. the diversity within each allele was visualized using the splitstree software package and aligned sequence files processed by uncorrectedp and neighbor-net algorithms. determination of genetic relationships among isolates by computational analysis the relatedness of isolates was analyzed using the start software. start was further used to determine the number of polymorphic nucleotide sites, calculate dn/ds ratios and construct profile-based dendrograms using the neighbor-joining algorithm. sequence-based trees were constructed in mega using the nj algorithm and bootstrap replications. the classical and standardized index of association were calculated in start using only one strain per st to avoid bias. additionally, the test of sawyer was applied to the synonymous polymorphic sites within the alleles at each locus. related sts were grouped using the start and eburst v to form clonal complexes. ccs were defined as groups of two or more independent isolates where each isolate had identical alleles at seven or more loci with at least one other member of the group. this definition was also used to number phylogenetic tree clades and branches. the allele profiles per locus were used to calculate the simpson’s index of diversity, where a value close to indicates high diversity and low diversity for values close to twice the standard deviation of the sid value was used to obtain an approximation for the 95 % confidence intervals as previously described. assay for simulated gastric conditions in simulated gastric juice for simulated gastric juice assays, strains were grown in m broth, containing lactose as carbon source to mimic the carbon source of milk. s. alactolyticus dsm20728t, s. bovis dsm20480t and s. equinus dsm20558t were incubated in bhi broth containing glucose instead of lactose to promote growth. bacteria cultures were standardized to an od of using 5 ml of phosphate-buffered saline, which is near the optimal ph for growth of sbsec. the standardized bacteria suspensions were then centrifuged at 6000 × g for 6 min at 20 °c. the pellet was resuspended in 5 ml pbs and stored on ice for 30 min. sgj was prepared using 1 g of nacl,  g of pepsin from porcine source, and  g peptone from casein dissolved in 200 ml distilled h2o. ph was adjusted to ph  or ph  using 5 m hcl. these ph values were selected within the range of human gastric conditions allowing for adequate time-dependent survival assays over the 15 min experiment duration. all sgj were filter-sterilized and stored at 4 °c. phosphate buffered saline was prepared containing 40 g l− nacl, 1 g l− kcl,  g l− na2hpo,  g l− kh2po in distilled h2o, and the ph was adjusted to using 5 m hcl for low ph-stress conditions. pbs was autoclaved and stored at 4 °c. the assay was performed in sterile 96-well plates using 270 μl of sgj or pbs per well. the plates were pre-heated to 37 °c and inoculated with 30 μl of the standardized test culture in parallel for all conditions. plates were incubated in air-tight plastic containers to reduce evaporation. samples were analyzed after and 15 min of incubation at 37 °c. enumeration of surviving bacteria was done after serial dilution from 10− to 10− and microspotting of 15 μl onto agar medium with a detection limit of log cfu ml− the assays were performed with two biological replications each comprising three technical replications. sii cj and sii ccug43820t were used as control strains in each replication. relative survival rates per strain were calculated where the input is set at log cfu ml− and reduction is expressed in negative log-values. for ease of comparison between strains at a given condition, relative survival rates of each strain were normalized using the formula /sd, where xi is the relative survival rate of a strain, xmean the mean relative survival rate of all strains for one condition and sd the standard deviation. assay for bacterial adhesion to extracellular matrix proteins bhi broth overnight cultures were standardized to an optical density at 600 nm of using pbs and centrifuged for 10 min at 3000 × g. the supernatant was discarded and the pellet resuspended in the original volume of pbs. the bacteria sample was then divided into two equal aliquots and centrifugation was repeated. the supernatant was discarded and one bacteria aliquot per culture was resuspended in pbs and acidified pbs, respectively. aliquots were kept on ice until use. extracellular matrix proteins used in this study were collagen type i and type iv, fibrinogen and fibronectin and the glyco-protein mucin type ii. bovine serum albumin was used as a control protein. stock solutions of ecms were prepared at a concentration of 10 μg ml− in pbs except bsa and mucin type ii at  mg ml− in  m tris hcl. maxisorptm 96-well plates were treated with 100 μl of an ecm stock solution per well. one 96-well plate was filled with 100 μl of pbs as a negative adhesion control of the plastic surface. the plates were incubated overnight at 4 °c. subsequently, the liquid was poured off by plate inversion and dried for 10 min at 65 °c. the plate surface was blocked with 100 μl pbs containing 1 % tween incubated for 1 h at 37 °c and subsequently washed times with 100 μl pbs containing  % tween coated plates were stored at 4 °c with 100 μl pbs added per well. bacteria cells resuspended in pbs ph  or ph  were then added to the coated 96-well plates. the inoculated plates were gently centrifuged for 10 min at 400 × g at room temperature to enable physical contact of the bacteria with the tested proteins, and the plates were then incubated at 37 °c for 1 h. the wells were then washed three times with 100 μl pbs containing  % tween and adhering bacteria were fixed for 20 min at 65 °c. staining was performed using crystal violet for 45 min at room temperature, followed by washing three times with 100 μl pbs. bacteria-associated crystal violet was then solubilized by adding 100 μl citrate buffer to each well and incubating for 1 h at 37 °c under constant agitation at 400–500 rpm. the absorbance of crystal violet was measured at 595 nm. the normalized cell adhesion of a strain was calculated from the mean of three independent biological repetitions, using the formula /sd, where xmedian and sd were obtained from the adhesion measurements of all strains under one condition. values larger than from the normalized were considered to be adherent. assay for inflammatory potential via nf-kb activation the activation of nuclear factor k-light-chain-enhancer of activated b cells was assessed using thp1-blue™ cells. these cells express all known toll-like receptors and enable quantification of the nf-kb activation via the seap reporter gene. thp1-blue™ cells were grown in rpmi supplemented with 2 mm l-glutamine, 50 iu ml− penicillin, μg ml− streptomycin and 10 % heat-inactivated fetal calf serum in a humidified 10 % co atmosphere at 37 °c. all culture media were supplied by life technologies. the assay was performed in 96-well plates containing 90 μl of thp1-blue™ cells at a concentration of 2 ×  cells per well. bacteria cells were prepared by centrifugation for 5 min at  × g of 1-ml aliquots of bacteria cultures standardized to an od of. the supernatant was transferred into a new tube and both supernatant and pellet were stored at −20 °c until use. thp1-blue™ cells were activated using 10 μl of the supernatant or resuspended pellet in 1 ml pbs. the plates were incubated for to 24 h. nf-kb activation was quantified by adding 10 μl of supernatant of the stimulated thp1-blue™ cells to 100 μl quanti-blue™. after to 8 h of incubation at 37 °c, seap activity was quantified at 655 nm as described by the supplier. normalization was performed as described for the adhesion assay. statistics statistical analysis was performed in jmp. data was first analyzed for normal distribution using the shapiro-wilk-w-test. non-normal distribution of log-transformed values required the utilization of a kruskal-wallis-test, which was upon rejection of h further compared in a post-hoc pair-wise wilcoxon-test without correction for multiple comparisons in order to reduce type ii errors. this was complemented by graphical analysis and detection of outliers by calculating the interquartile range iq and the defined limitations of lower fence of q1- *iq and upper fence q3 +  *iq. strains displaying values outside of the upper or lower fence were considered outliers.  identification of sbsec species by groel and 16s rrna gene sequencing partial groel sequencing was performed on all sbsec strains to confirm species and subspecies status. all groel nucleotide sequences were trimmed to an equal length of 772 bp for comparison and neighbor-joining tree construction. sii and s. lutetiensis had a maximum identity of –  %, whereas none of the s. gallolyticus branch members displayed more than 91 % identity to sii. the s. gallolyticus branch was more conserved, yielding –  % or –  % between sgm and sgg or sgp, respectively. intraspecies variations were generally low, yielding –  % sequence identity among strains of a single subspecies. s. equinus jb showed highest identity to sii by groel and  % to the s. lutetiensis by 16s rrna gene sequence. in addition, ccug, which was initially classified as s. bovis and subsequently reclassified as sii, was again reclassified in this study as s. equinus. similarly, s. bovis atcc was reclassified as sgp based on groel and 16s rrna gene sequences. strain j 40– isolated from fermented milk from bangladesh was identified as sgm by groel sequencing, whereas 16s rrna gene sequence and rep-pcr fingerprinting indicated s. infantarius/s. equinus branch. mlst loci characterization and key performance identifiers the mlst assay was designed as an overarching scheme to target the entire sbsec and validate species assignment. aligned and trimmed sequences were between 393 bp and 573 bp in length. in total, sequence types were defined for the strains analyzed. the number of alleles ranged from for tpia to for pyre and muts while the number of polymorphic sites ranged from for tpia to for muts. the number of alleles and pms were lower for individual species of the sbsec. tpia showed the highest degree of conservation whereas muts and pyre displayed the highest variability, which was also reflected in the total number of alleles defined.table key characteristics of the sbsec-mlst scheme for the overall sbsec, sii and s. lutetiensis index of association1 ia s ddla gki glna muts muts2 phes pros pyre thrs tpia 1index of association based on one strain per st: ia: calculated using classical maynard smith approach, ia s: standardized according to haubold 2sbsec: includes also 3 s. gallolyticus subsp. pasteurianus and 1 s. alactolyticus strains not listed separately 3 a allele 4 pms polymorphic sites 5sid: simpson’s index of diversity with 95 % confidence interval calculated from sid+/− 2*standard deviation 6no dn/ds calculation due to limited allele numbers and differences values were calculated for the entire sbsec and individual species, including all target genes and simpson’s index of diversity as an indication for the differentiation power among strains index of association1 ia s ddla gki glna muts muts2 phes pros pyre thrs tpia 1index of association based on one strain per st: ia: calculated using classical maynard smith approach, ia s: standardized according to haubold 2sbsec: includes also s. gallolyticus subsp. pasteurianus and s. alactolyticus strains not listed separately 3 a allele 4 pms polymorphic sites 5sid: simpson’s index of diversity with 95 % confidence interval calculated from sid+/− 2*standard deviation 6no dn/ds calculation due to limited allele numbers and differences 7nan: maximum diversity, no ci calculcated. values were calculated for the entire sbsec and individual species, including all target genes and simpson’s index of diversity as an indication for the differentiation power among strains splitstree analysis of each individual locus provided network-like structures for all ten loci. least-square tree fitting resulted in values > for all trees. within each splitstree network, separate clusters were distinguishable for the s. gallolyticus and s. infantarius branch, with s. alactolyticus as the most distant species. the index of association and its standardized version for the overall sbsec were and, respectively. ias values calculated for individual groups were,,,,, and, while sgp comprised too few strains for calculation. except for sgm, all other groups including sbsec displayed significant linkage disequilibrium, suggesting recombinatory evolution. simpson’s index of diversity, which indicates the discrimination power of the method, was for the overall sbsec approach and for sii. definition of clonal complexes clonal complexes are an important classification tool for mlst schemes. ccs were defined for strains that share out of loci using eburst v and these were depicted in a profile-based neighbor-joining tree. a total of ccs were determined, which grouped strains and left singletons. cc and cc were both composed of west african dairy strains, whereas kenyan and somali strains were both grouped in cc and cc. cc represented the main human pathogen strain cluster that delineated in close relationship with ivorian dairy cluster cc small ccs comprising two strains each were found for sgp and s. lutetiensis. cc comprised most of the dairy sgm strains from africa and greece, suggesting close relationships among strains of this subspecies. s. bovis and s. equinus cannot yet be considered as real ccs, as each cc contained only two versions of the corresponding type strain obtained from different sources. the number of strains did not allow the definition of ccs for sgg or the prediction of a founding strain per cc. the defined ccs 1– suggest clear regional relationships among sii dairy strains and even across continents, but also pathogenic relationships concentrated in cc mlst-based phylogeny of sbsec and individual species groups the sbsec mlst scheme provides the first detailed overview of all species within the sbsec. all branches 1– were numbered sequentially for ease of comparison using the cc definition of or more identical alleles out of fig. sequence-based phylogenetic mlst tree of the sbsec. the mlst tree comprised commensal animal and human strains, human pathogenic strains and food-derived strains of all sbsec species. the tree was calculated from the concatenated sequences of the mlst loci of each strain using the neighbor joining algorithm and bootstrap replications. the tree was rooted to s. alactolyticus dsm20728t. corresponding sts, alleles and prevalence of dairy adaptation marker genes, origin and isolation source are indicated. branches were defined and numbered according to the clonal complex specifications primary differentiation established an s. gallolyticus clade and an s. equinus/s. lutetiensis/sii clade. within each clade, species-specific clades were delineated. the s. equinus clade harbored both the s. equinus and s. bovis type strains that originally defined sbsec. s. equinus/s.bovis strains were distributed between a ruminal/human clade and an equine branch. s. equinus jb a ruminant isolate, delineated a separate branch with closer association with the s. lutetiensis clade than the major s. equinus branch, suggesting reclassification of this strain based on mlst analysis. the s. lutetiensis clade reflected a major division between human clade and a potential side branch comprising an animal strain. the human clade was further divided into human blood and fecal isolates. the s. gallolyticus clade showed subdivision into a clade comprising sgg and sgp, and a separate clade of sgm. within the sgg/sgp clade, separate clades were delineated for human pathogenic strains of sgg and sgp. animal commensal strains seemed more heterogeneous, forming two separate branches, one for the sgg type strain dsm16831t and one even further related for the calf isolate lmg, indicating different lineages of commensal and pathogenic origins. among the sgm strains, the divisions between greek, east african and west african strains were fully reflected, suggesting a main african lineage that was further differentiated into east and west african lineages. the sii clade was analyzed in more detail. several dairy and human clades were delineated, suggesting a general trend of different lineages. strains of human clinical cases were concentrated in only two specific clades comprising branches – and branch regardless of geographic origin. the majority of dairy strains were found in two dairy clades comprising east african dairy isolates and west african and asian dairy isolates, suggesting different lineages within the main dairy lineage of sii. east african dairy strains formed a single clade subdivided into two main clades of dairy isolates only and mixed human and dairy isolates comprising the sii type strain. the majority of west african and the asian dairy isolates were grouped in a single dairy clade comprising clades and the delineation of clades and was directly linked to thrs alleles and featuring highest sequence identities of –  % to s. thermophilus instead of other sii or sbsec members as for the other alleles. however, west african dairy strains were more heterogeneous than those of east africa, and also featured in addition to clade several shared or related clades with human blood isolates and east african dairy isolates. prevalence of lacs/lacz genes as marker genes for dairy adaptation the prevalence of s. thermophilus-like lacs/lacz genes was assessed as marker genes for dairy adaptations among all sbsec, including sii strains, using specific pcr assays and previously obtained data. only sii strains were found to harbor lacs/lacz, with highest prevalence among the major dairy branches,, and. in contrast, dairy strains more closely related to human strains, such as sii 150a and the four west african dairy strains found in branches and did not harbor the lacs/lacz genes. resistance to simulated gastric juice the resistance of sbsec strains to sgj was tested in sgj at ph  or ph  and pbs at ph  as control. all differences in this section were tested for significance at p <  unless mentioned otherwise. generally, relative cell survival was significantly different for the three conditions ph , ph  and pbs ph  over 15 min. significantly different survival was also observed among strains at ph  or ph  whereas no significant difference was observed among the strains tested under control conditions of pbs ph. significantly different survival rates were determined between dairy sii clades and featuring higher survival compared to sii commensal clade and predominantly blood isolate clade s. lutetiensis clade s. gallolyticus clades and as well as s. alactolyticus branch. survival rates were significantly dependent on time-ph combinations. conditions of ph  yielded low survival rates, often below detection limit after 15 min. this strong inhibition resulted in a few significant differences between strains and mlst branches with the exception of outliers in clade and branch. significantly highest relative survival in comparison to all other strains was determined for the four west african dairy sii strains 6by-11b, 13af, 12ay- and p2vc1- clustered in mlst branches and as well as the human blood isolate s. lutetiensis and all s. thermophilus strains. significantly higher relative survival was therefore determined for west african compared to east african sii strains in branches and no significant correlation was observed between lacs/lacz prevalence, habitat and relative survival under the conditions tested. adhesion to extracellular matrix proteins adhesion of sbsec strains to various ecms was investigated to estimate their interaction with cell surfaces of the gi tract and blood vessels. due to heavy tailed distributed values, all adhesion values were normalized and analyzed using a median-based calculation.fig. analysis of the adhesion ability of sbsec members to extracellular matrix proteins and activation of the nf-kb signaling pathway. adhesion and nf-kb activation values of individual strains displayed in a heat map. the phylogenetic tree is based on the mlst-tree but not drawn to scale for better visualization. adhesion measurements were normalized for each strain based on the adhesion median of all strains. values larger than sd from the normalized adhesion were considered to be adherent. outliers were marked using an asterix* lb. plantarum wcfs and lb. reuteri dsm20016t were used as positive adhesion controls showing significant adhesion towards fibronectin, fibrinogen, collagen i and mucin type ii that are known factors contributing to their persistence in the human gi tract. as expected, lb. plantarum nz the sortase knockout mutant of wcsf depleted of its adhering surface proteins, showed no adhesion and served as negative adhesion control. adhesion values obtained at experimental conditions of ph  showed better reproducibility in contrast to those values obtained at ph  and were therefore chosen for analysis. adhesion characteristics were separated by mlst clustering and isolate origin. generally, adhesion was most prevalent to the two ecms collagen i and mucin type ii, with and  % of sbsec isolates showing adhesion, respectively, followed by adhesion to fibronectin, fibrinogen and collagen iv by  % of sbsec isolates. among human sii, adhesion to mucin type ii and unspecific adhesion was observed for  % and adhesion to all remaining ecms for 33 % of human sii strains. among these human strains, comparison of blood isolates and potentially commensal isolates showed that significant adhesion was exclusively present among sii blood isolates comprised mainly in branches and human fecal isolate sii ccug43820t featured no significant adhesion abilities to fibronectin, fibrinogen or collagen i and only weak adhesion to collagen iv and mucin type ii. similarly, s. lutetiensis blood isolate displayed significant wide range adhesion abilities whereas human fecal isolates of s. lutetiensis displayed no significant adhesion abilities, suggesting low prerequisite or other mechanisms contributing to the persistence of these possibly commensal sii and s. lutetiensis strains in the human gi tract. similarly, dairy sii generally had a low prevalence of adhesion abilities to the ecms tested and those dairy sii isolates showing adhesion, featured mainly adhesion to collagen iv and mucin type ii detected in and  % of strains, respectively, with the exception of cji cj cj and ab2vb2- showing also adhesion to collagen i or mucin type ii. clustered according to mlst and geography, west african dairy sii of mlst branches, and displayed low adhesion abilities to the ecms tested. adhesion was more prevalent among east african dairy sii strains with significant adhesion to at least one ecm for strains cji cj cj and cj as well as west african dairy sii strain ab2vb2- which all share the same mlst clade 1–. the significantly adherent west african dairy strain sii av2a was the exception to cluster with the otherwise low-adherent strains of the west african main dairy clade. among the s. gallolyticus clade 14– high adhesion abilities towards collagen i, iv and fibrinogen was observed. dairy sgm branches of east africa and greece showed adhesion towards mucin type ii and fibronectin. in contrast, sgp did not display any significant adhesion abilities allowing clear separation by mlst clustering from sgg. determination of immune stimulatory potential by nf-kb activation activation of the transcription factor nf-kb was used as marker for the immune stimulation of the sbsec strains studied using the thp1-blue reporter cell line that expresses all tlrs that detect bacterial motifs in the intestine. nf-kb activation was generally low for the strain panel tested with significantly higher values obtained for cell pellets in contrast to cell-free supernatants. cell pellets, significant inflammation was detected only for individual strains such as sii ccug, s. equinus ccug and sgg dsm 16831t. the majority of dairy isolates and also human commensal and blood isolates did not significantly activate nf-kb, particularly in comparison to control strain e. faecalis jh2- among cell-free supernatants, highest nf-kb activation was measured for dairy isolates sii ab2vb2-, sii cj and sgm cji, as well as human blood isolates sii and sii ccug, displayed significant nf-kb activation. again, no significant nf-kb activation was detected among supernatants of other human commensal strains of s. lutetiensis, sii and even pathogenic s. lutetiensis, sii, sgg and sgp. discussion the sbsec is a highly diverse group of bacteria colonizing multiple habitats such as food, human and animal gi tracts. although the 16s rrna gene sequence provides sufficient power to differentiate the main branches of s. gallolyticus, s. infantarius and s. alactolyticus, differentiation at subspecies level requires sequencing of additional genes such as groel. we showed that groel provided stringent sequence data and discrimination power for sbsec subspecies identification. interestingly, sii j 40- seems to be an exception, 16s rrna gene sequence suggested an s. infantarius branch member while groel indicates highest similarity with sgm; however, the identity in the sbsec mlst scheme was confirmed to be sii. the ongoing horizontal gene transfer among streptococci and sbsec in combination with natural competence may have changed even genes considered as housekeeping genes. further indications for hgt between dairy sii and s. thermophilus were observed among thrs alleles of strains comprised in mlst clades and representing asian dairy strain sii j 40- and the major west african clade. most strains in clade harbor lacs and lacz which serve as marker genes for s. thermophilus-like dairy adaptation. in combination with the thrs alleles present, this suggests a common ancestor for clades and that underwent multiple hgt with s. thermophilus. ultimately, these observations of hgt and natural competence support the need for a multi gene approach such as mlst for sbsec species to ensure accurate identification as well as inter- and intraspecies phylogenetics. the sbsec mlst scheme developed in this study provides the first multi-gene-based comparison across the entire sbsec. in contrast to previous sgg mlst schemes, this scheme provides a tool for clustering all currently described sbsec species and subspecies from human, animal and food sources. the sid of – is within the range and confidence interval of that of the sgg mlst. therefore, our novel sbsec mlst assay provides equally high differentiation power as the sgg-specific assays without the restriction of being limited to a single species but instead being applicable to all species of the sbsec. furthermore, extraction and phylogenetic comparison of mlst alleles from whole genome sequences as performed for several strains in this study ensure the application of this sbsec mlst scheme also in combination with next generation sequencing of bacterial genomes as data source instead of traditional sanger sequencing. the sbsec mlst scheme provided highly comparable phylogenetic tree delineation for sgg as those obtained via the sgg-specific mlst scheme. dn/ds ratios of all loci in the sbsec mlst scheme were < and comparable to those of the two sgg-specific schemes. number and percentage of pms of most loci were in the range of dumke et al. at –  % and shibata et al. at –  %. the ia of – for individual species was slightly lower than the ia of determined by dumke et al. whereas shibata et al. did not provide a value for ia. the ia values for the overall sbsec, s. equinus, sgg, sii and s. lutetiensis indicated significant linkage disequilibrium between the selected loci as previously observed for the loci set of sgg. linkage disequilibrium is the effect of non-random association of the selected alleles in an mlst scheme, which indicates that recombination was a factor in the evolution of the sbsec as seen also by network-like structures in splitstree analyses. sgm was the only subspecies where linkage disequilibrium was not detected. sgm also displayed the lowest genetic diversity despite the analysis of greek, somali, ivorian and mali strains, suggesting that sgm strains from these different origins are highly related. this could support the theory that sgm is actually a subspecies of dairy origin within the sbsec, enabled through clear delineation from other sbsec members in this mlst scheme. genome sequencing of sgm strains will be required to further assess their evolution in relation to the presence of putative sgm dairy adaptation markers such as duplicated lacg described for sgm aca-dc and decayed or active virulence factors. sgg and s. lutetiensis showed a clear differentiation between presumptive human commensal lineages and those of infectious nature isolated from feces, blood or heart samples of humans and animals. in particular, commensal strains such as koala bear isolate sgg dsm16831t and calf isolate sgg lmg were clearly distinct from human pathogenic lineages, suggesting the existence of potentially commensal and pathogenic lineages among many sbsec species. this requires evaluation in a one health approach, encompassing animal and human health, food, the zoonotic potential of sbsec members and their status as pathobionts. the subspecies sii that originated from human and animal-derived food sources had the highest representation in this mlst assay. clear regional differentiation was observed for the majority of east and west african dairy strains. differentiation by host and source was observed for most human and dairy strains, suggesting different lineages depending on fecal or blood origin. the mlst phylogeny suggests that human fecal and especially blood isolates cluster by isolation source, and not by geographical origin in contrast to most dairy strains. this possibly reflects different adaptation and risk levels. dairy strains clustered in distinct lineages featuring a prominent west african and asian main branch and an east african branch. the main lineage of east african dairy strains displayed remarkable conservation in contrast to the higher diversity among west african strains and especially those collected in côte d’ivoire. the regions of mandera and garissa in kenya share a common cultural background with somalia. large somali populations are represented in both kenyan regions, with a long tradition in livestock trading, mobile pastoral lifestyle and migration across borders, possibly explaining the spreading of similar strain lineages. certain lineages among the west african dairy sii without the marker genes lacs/lacz clustered separately from the main dairy lineages and closer to human fecal and pathogenic isolates, suggesting the existence of overlapping lineages of sii or contamination of foods by fecal material. this raises concerns about potential health risks of these strain lineages and the origin of these lineages in food. exposure of strains to sgj was used to assess the survival of bacteria in stomach conditions and delivery potential to the colon. generally, strains of mlst clusters showed similar relative survival being affected by sgj. however, at ph  and ph  several dairy sii clades and mainly four dairy strains showed significantly higher relative survival than human blood sii isolates/clades, respectively, suggesting different adaptation to acidic stress in relation to strain origin and mlst clustering. mlst-based clustering was further used to assess the ability of individual strains to adhere to ecm and activate the nf-kb pathway in correlation with mlst clades and phylogeny. adherence is an important first step in colonization of the gi tract or damaged body tissues such as heart valves. mlst clusters provided clear differentiation of highly adhesive sii/s. lutetiensis blood isolates, the s. gallolyticus clade, less adhesive east african dairy strains and nearly non-adhesive west african strains, suggesting a generally lower tendency among west african dairy clades and commensal s. lutetiensis to adhere to the ecms tested. the genetic factors contributing to adhesion were so far only related to a specific pilus loci for collagen binding termed pil pil and pil present in the sgg genomes analyzed. only pil was found in sii, sgm and sgp, whereas the crucial pil factor of sgg was absent. comparative genomics and further functional analyses are needed to determine the responsible factors in highly adhering strains of sii, s. lutetiensis, sgm and sgg, and possibly identify new virulence factors to enable an in-depth risk assessment for the different sbsec species in relation to human and animal diseases. the nf-kb activation by bacteria attaching to the intestinal mucosa plays a central role in immune signaling of the host, particularly in the intestinal tract. for most strains, nf-kb activation was positively correlated with significant adhesion abilities, suggesting health implications and the need for further functional analyses for such strains. several blood isolates as well as pathogenic sgp dsm 15351t did not yield significant nf-kb activation, which could enable them to infect the host without inducing a strong immune response.  the overarching mlst scheme for the entire sbsec provides an epidemiologic and diagnostic tool enabling in-depth inter- and intraspecies comparisons of all species within this emerging human and animal pathogenic complex. the mlst scheme based on housekeeping genes allowed delineating all species within the sbsec. it enabled division by host and habitat and furthermore possibly pathogenic lineages from commensal and food lineages. this approach will be useful in enhancing taxonomic classifications of human, animal and food isolates as well as their relationships. our data suggest that the parallel evolution of east and west african dairy lineages is currently in the process of differentiating from the human lineage including genetic exchange with other streptococci. however, enhanced adhesion abilities of certain lineages similar to those of sgg and the ability to cause inflammation or evade immune response are certainly of concern and require further investigations into pathogenicity mechanisms shared among the sbsec and those unique to only specific species. this will greatly assist in designing appropriate treatment strategies for patients infected with this emerging pathogenic group depending on mlst clustering. furthermore, accurately collecting clinically relevant epidemiologic data in relation to transmission, zoonosis, host specificity and identification of virulence mechanisms will greatly contribute to an enhanced risk assessment of pathogenic, commensal and food lineages or the potential role as pathobionts in an overarching one health approach. abbreviations cc, clonal complex; crc, colorectal cancer; ecm, extracellular matrix protein; gi tract, gastrointestinal tract; hgt, horizontal gene transfer; ia, index of association; ias, standardized index of association; ie, infective endocarditis; mlst, multi locus sequence typing; pbs, phosphate-buffered saline; pms, polymorphic sites; sbsec, streptococcus bovis/streptococcus equinus complex; sgg, streptococcus gallolyticus subsp. gallolyticus; sgj, simulated gastric juice; sgm, streptococcus gallolyticus subsp. macedonicus; sgp, streptococcus gallolyticus subsp. pasteurianus; sii, streptococcus infantarius subsp. infantarius; sid, simpson’s index of diversity; st, sequence type additional files additional file 1: title of data: strains used in this study. description of data: strains used in this study. additional file 2: title of data: primers developed and used in this study. description of data: primers developed and used in this study for the sii/sbsec mlst scheme of housekeeping genes and identification based on groel, lacs and lacz genes. additional file 3: title of data: mlst candidate loci used in this study. description of data: all candidate loci which were used in this study were selected from literature and mlst.net. conserved housekeeping genes were selected on the basis of maximal nucleotide divergence for intraspecies differentiation. all loci which are marked with an asterisk were used for the final sbsec mlst scheme. additional file 4: title of data: sequence identity matrix of 772-bp groel fragments of sbsec species. description of data: sequence identity matrix of 772-bp groel fragments of sbsec species additional file 5: title of data: groel-based phylogenetic tree of sbsec strains of this study. description of data: phylogenetic tree of all sbsec strains investigated this study based on neighbor-joining calculation of a 772-bp sequence fragment of the groel gene. the scale bar below the tree indicates the evolutionary distance using the number of base substitutions per site as units. species abbreviations: s. infantarius subsp. infantarius, s. gallolyticus subsp. gallolyticus, s, gallolyticus subsp. macedonicus, s. gallolyticus subsp. pasteurianus. additional file 6: title of data: 16s rrna gene-based phylogenetic tree of selected sbsec members. description of data: phylogenetic tree of selected sbsec strains of all sbsec species based on a 930-bp fragment of the 16s rrna gene. calculations were performed using the neighbor-joining algorithm. the scale bar below the tree indicates the evolutionary distance using the number of base substitutions per site as units. species abbreviations: s. infantarius subsp. infantarius, s. gallolyticus subsp. gallolyticus, s, gallolyticus subsp. macedonicus, s. gallolyticus subsp. pasteurianus. additional file 7: title of data: splitstree visualization of the mlst loci. description of data: the trees were obtained using the defined allele sequences and calculation by neighbor-joining algorithm. corresponding alleles of selected reference strains are highlighted: sii ccug43820t = bright green, sgg dsm16831t = red, sgm dsm15879t = ice blue, s. lutetiensis ccug43822 = yellow, sgp dsm15351t = brown, s. bovis dsm20480t = pink, s. equinus dsm20558t=dark blue, s. alactolyticus dsm20728t = orange. additional file 8: title of data: profile-based mlst tree of all species within the sbsec. description of data: the strain pool comprises commensal animal and human strains, human pathogenic strains and food-derived strains. the tree was calculated by neighbor joining algorithm using the allele profiles of housekeeping genes. corresponding sequence types, alleles and prevalence of dairy adaptation marker genes including origin and isolation source are indicated. species abbreviations: s. infantarius subsp. infantarius, s. gallolyticus subsp. gallolyticus, s, gallolyticus subsp. macedonicus, s. gallolyticus subsp. pasteurianus. additional file 9: title of data: sequence-based phylogenetic mlst tree of s. infantarius subsp. infantarius. description of data: the mlst tree comprised commensal animal and human strains, human pathogenic strains and food-derived strains of all sbsec species. the tree was calculated by neighbor joining algorithm using the concatenated partial sequences of housekeeping genes. trees were rooted to s. equinus jb corresponding sequence types, alleles and prevalence of dairy adaptation marker genes including origin and isolation source are indicated. branches were defined and numbered according to the clonal complex specifications. additional file 10: title of data: relative survival of sbsec members under simulated gastric conditions at ph . description of data: relative survival of sbsec strains after and 15 min incubation in sgj at ph  performed in two biological replications. relative survival values in log cfu ml− were normalized according to the mean of all measurements of one condition to allow comparison between strains. positive normalized values indicate relative survival higher than the mean of all strains. statistical distribution indicators: median: large black dotted line; q and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line. the phylogenetic tree is based on the mlst-tree but not drawn to scale. additional file 11: title of data: relative survival of sbsec members under simulated gastric conditions in pbs ph . description of data: relative survival of sbsec strains at and 15 min in pbs ph  performed in biological replications. input of log cfu ml− was used as basis from which relative survival is expressed. relative survival values were normalized to allow comparison between strains where positive values indicate relative survival higher than the mean of all strains. statistical distribution indicators: median: large black dotted line; q and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line. the phylogenetic tree is based on the mlst-tree but not drawn to scale. additional file 12: title of data: relative survival of sbsec members under simulated gastric conditions at ph . description of data: relative survival of sbsec strains at and 15 min in in simulated gastric conditions at ph  performed in biological replications. input of log cfu ml− was used as basis from which relative survival is expressed. relative survival values were normalized to allow comparison between strains where positive values indicate relative survival higher than the mean of all strains. statistical distribution indicators: median: large black dotted line; q and q3: small black dotted line; lower and upper outlier fence: large yellow dotted line. the phylogenetic tree is based on the mlst-tree but not drawn to scale. the authors greatly acknowledge the sbsec strains received from ingolf nes, tsilia lazarovitch, taku miyamoto, rosa del campo/beatriz romero hernández, marisa haenni/jean-yves madec. the authors thank peter bron and michiel kleerebezem for the gift of lactobacillus plantarum wcfs and lactobacillus plantarum nz this study and all authors were funded by the ubs optimus foundation, switzerland and the erafrica program theme “new ideas”. funding this project was financed by the ubs optimus foundation, switzerland and erafrica “new ideas”. the funders had no role in study design, data collection, analysis, interpretation, decision to publish, or preparation of the manuscript. all authors have read and approved the final manuscript. availability of data and material the data set supporting the results of this article are available in the genbank and pubmlst.org repositories under accession number ranges and and sbsec mlst database, respectively. authors’ contributions cj, dwmk, ja, db, sv, ts and ji performed the experiments to type bacteria isolates, develop and validate the mlst assay and contributed to data analysis and interpretation. ja and tw conducted adhesion and inflammation assays. fk conducted gastric survival assays. ja, tw and fk contributed to data analysis and interpretation of the respective assays. cj, tw, cl and lm designed and supervised the experiments. cj, tw, bb, cl, dwmk, wkm and lm performed analyses and interpretation, and wrote and revised the manuscript. all authors read and approved the final manuscript. competing interests the authors declare that they have no competing interests. consent for publication not applicable. ethical approval and consent to participate not applicable.
changing taxonomy, shared habitats, natural competence and evidence for horizontal gene transfer pose difficulties for determining their phylogeny, epidemiology and virulence mechanisms. commensal and pathogenic sbsec strains of human, animal and food origin from europe, asia, america and africa were used in the mlst scheme yielding sequence types and clonal complexes differentiated into distinct habitat-associated and mixed lineages. the streptococcus bovis/streptococcus equinus complex comprises seven species classified as human and animal commensals, emerging opportunistic pathogens and food fermentative organisms. an sbsec overarching multi locus sequence type scheme targeting housekeeping genes was developed, validated and combined with host-related properties of adhesion to extracellular matrix proteins, activation of the immune responses via nf-kb and survival in simulated gastric juice. they highlighted the health risk among shared lineages and establish the basis to elucidate transmission, host specificity, virulence mechanisms and enhanced risk assessment as pathobionts in an overarching one health approach. nf-kb activation was observed for a limited number of dairy and blood isolates suggesting the potential of some pathogenic strains for reduced immune activation. adhesion to ecms collagen i and mucin type ii was a common characteristic followed by adhesion to fibronectin and fibrinogen. high adhesion values shared between certain dairy and blood strains as well as the behavior of nf-kb activation are concerning for specific lineages. combining phylogenetic and functional analyses via sbsec mlst enabled the clear delineation of strain clades to unravel the complexity of this bacterial group. strains from dairy mlst clades displayed the highest relative survival to sgj independently of dairy adaptation markers lacs/lacz.
 the streptococcus bovis/streptococcus equinus complex comprises seven species classified as human and animal commensals, emerging opportunistic pathogens and food fermentative organisms. changing taxonomy, shared habitats, natural competence and evidence for horizontal gene transfer pose difficulties for determining their phylogeny, epidemiology and virulence mechanisms. thus, novel phylogenetic and functional classifications are required. an sbsec overarching multi locus sequence type scheme targeting housekeeping genes was developed, validated and combined with host-related properties of adhesion to extracellular matrix proteins, activation of the immune responses via nf-kb and survival in simulated gastric juice.  commensal and pathogenic sbsec strains of human, animal and food origin from europe, asia, america and africa were used in the mlst scheme yielding sequence types and clonal complexes differentiated into distinct habitat-associated and mixed lineages. adhesion to ecms collagen i and mucin type ii was a common characteristic followed by adhesion to fibronectin and fibrinogen. high adhesion abilities were found for east african dairy and human blood isolate branches whereas commensal fecal sbsec displayed low adhesion. nf-kb activation was observed for a limited number of dairy and blood isolates suggesting the potential of some pathogenic strains for reduced immune activation. strains from dairy mlst clades displayed the highest relative survival to sgj independently of dairy adaptation markers lacs/lacz.  combining phylogenetic and functional analyses via sbsec mlst enabled the clear delineation of strain clades to unravel the complexity of this bacterial group. high adhesion values shared between certain dairy and blood strains as well as the behavior of nf-kb activation are concerning for specific lineages. they highlighted the health risk among shared lineages and establish the basis to elucidate transmission, host specificity, virulence mechanisms and enhanced risk assessment as pathobionts in an overarching one health approach. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords streptococcus infantariusstreptococcus macedonicusstreptococcus gallolyticusone healthadhesioninflammationcolorectal cancerfoodborne diseaseinfective endocarditispathobiontubs optimus foundationerafricaissue-copyright-statement© the author 2016 