Skip to main content

Genomic analysis of group B Streptococcus from milk demonstrates the need for improved biosecurity: a cross-sectional study of pastoralist camels in Kenya



Streptococcus agalactiae (Group B Streptococcus, (GBS)) is the leading cause of mastitis (inflammation of the mammary gland) among dairy camels in Sub-Saharan Africa, with negative implications for milk production and quality and animal welfare. Camel milk is often consumed raw and presence of GBS in milk may pose a public health threat. Little is known about the population structure or virulence factors of camel GBS. We investigated the molecular epidemiology of camel GBS and its implications for mastitis control and public health.


Using whole genome sequencing, we analysed 65 camel milk GBS isolates from 19 herds in Isiolo, Kenya. Six sequence types (STs) were identified, mostly belonging to previously described camel-specific STs. One isolate belonged to ST1, a predominantly human-associated lineage, possibly as a result of interspecies transmission. Most (54/65) isolates belonged to ST616, indicative of contagious transmission. Phylogenetic analysis of GBS core genomes showed similar levels of heterogeneity within- and between herds, suggesting ongoing between-herd transmission. The lactose operon, a marker of GBS adaptation to the mammary niche, was found in 75 % of the isolates, and tetracycline resistance gene tet(M) in all but two isolates. Only the ST1 isolate harboured virulence genes scpB and lmb, which are associated with human host adaptation.


GBS in milk from Kenyan camel herds largely belongs to ST616 and shows signatures of adaptation to the udder. The finding of similar levels of within- and between herd heterogeneity of GBS in camel herds, as well as potential human-camel transmission highlights the need for improved internal as well as external biosecurity to curb disease transmission and increase milk production.

Peer Review reports


In the arid and semi-arid lands of the Horn of Africa, camels are valuable assets for household income and food security, and commonly kept for milk production. In these areas, where desertification and prolonged droughts are constraints for food production, camel pastoralism is an integral part of the sustenance of the inhabitants, and in some areas camel milk and meat contribute to more than 50 % of the diet [1, 2]. Camels are well-adapted to harsh conditions and can continue to produce milk despite limited access to feed and water, which sets them apart from other types of livestock [3]. Kenya has the third largest camel population globally [4]. The majority of these camels are kept by pastoralists, who adhere to nomadic husbandry traditions, browsing camels over large areas and milking by hand [5]. These pastoralists live in close contact with their animals, and milk is consumed without prior pasteurization, both practices that contribute to zoonotic disease transmission risks [2].

Mastitis, inflammation of the mammary gland, is a frequently-occurring problem among dairy camels with negative implications for production and animal welfare [6, 7]. Streptococcus agalactiae, or group B Streptococcus (GBS), has been identified as a common cause of both clinical mastitis (CM) and subclinical mastitis (SCM) in camels [8, 9], often resulting in chronic infections [10], reduced milk yield [11] and high bacterial counts in milk [12]. In Kenya, a lack of involvement of veterinary services and poor compliance with recommended dosing regimens in camels [13], in combination with the common use of substandard antimicrobial products [14, 15], may serve as drivers in the development of antimicrobial resistance. Resistance has been reported in bacteria from camel milk within the region, including GBS [8, 16, 17].

In the bovine dairy industry, the molecular epidemiology of GBS largely depends on the local context. Studies in dairy herds in northern Europe (Denmark, Finland, Norway) [18,19,20] and in Australia [21] show limited strain diversity of GBS at farm level, with one genotype predominating in each herd. Differently, in India and Colombia, extensive within-herd heterogeneity has been observed, with an overlap in sequence types (STs) isolated from cattle and people [22, 23]. Variation in access to veterinary services [24], milking hygiene [8], and other biosecurity practices, [25] affecting pathogen transmission and mastitis control may contribute to such differences.

Investigations of the global GBS population have revealed the existence of clonal complexes (CC) associated with specific host species or niches (e.g. the mammary gland) and generalist lineages [22, 26, 27]. Niche adaptation, such as lactose fermentation, has been described in mastitis-causing isolates, and the lactose operon (Lac.2) has been identified as responsible for lactose metabolism in the vast majority of bovine isolates [19, 28, 29]. In addition to genes encoding metabolic properties, virulence genes, including adhesins and invasins, can affect pathogens’ ability to cause infections. For example, scpB (C5a peptidase) and lmb (laminin-binding protein) have been strongly associated with disease in humans, but not animals [29].

In camels, little is known about the population structure and epidemiology of GBS. Fischer et al. [17] found that isolates from camel mastitis primarily belonged to a common genotype, ST616, but with limited description of spatial or social relations between animals or herds under investigation. Combined molecular and epidemiological investigation is needed to increase the understanding of sources and transmission routes for GBS in pastoralist camel herds, to expand the knowledge base underpinning mastitis control strategies, and to explore the potential threat to public health. Here, we investigate the molecular epidemiology of GBS in dairy camel herds in Kenya using genomic and phylogenetic analysis, as well as spatial mapping.


Study area and selection of herds

This cross-sectional study was undertaken around Isiolo town, located in Isiolo County in central Kenya. Isiolo County is classified as arid or semi-arid with an annual rainfall of approximately 150–250 mm [30]. The camel population in the region was estimated to comprise of 45,309 individuals [31] distributed over 2,050 camel milk producers [32] out of a total human population of 268,002 individuals [33]. Camel keeping households produce milk for household consumption, but there is also an expanding camel milk market with milk being sold commercially along an informal milk value chain.

Herds were selected based on the following inclusion criteria: Pastoralist herd practicing extensive browsing and selling milk, willingness among camel owners to participate, and accessibility of the herds [8]. All selected herds were visited once during the period from February to April 2017, which corresponded to the end of dry season or early wet season. Herd data regarding management was collected digitally on a tablet using free open-source software [34]. All herds belonged to the Mlango-Ngarendare-Burat camel milk cluster [35] and sold milk along the informal milk value chain with end markets in Isiolo town and central Nairobi.

Milk sampling, bacteriological culture and isolate selection

Milk was sampled and cultured as described elsewhere [8]. In brief, 20 pastoralist camel herds were visited and a subset of 7 to 13 camels per herd were screened for evidence of mastitis (udder inflammation). Initially, udders were examined clinically and checked for the presence of blind (non-milk producing) quarters and signs of inflammation, such as swelling, increased temperature, pain and redness or palpatory findings of induration of udder tissue. Quarter milk was subjected to the California Mastitis Test (CMT), which provides semi-quantitative measurement of milk leucocyte content as indicator of the degree of inflammation [36]. After discarding the first streaks of milk, approximately 10 mL of milk from each quarter was milked in to the CMT-paddle. Milk was assessed visually for colour, viscosity, presence of blood and clots and then mixed with an equal amount of test reagent. The liquids were mixed by gently rotating the paddle and the viscosity of the solution was scored according to the Scandinavian scoring system (scale from 1 to 5) where 1 represents no change in viscosity, and 5 represents gel formation with a distinct peak [37]. A CMT-score of ≥ 3 was considered indicative of mastitis. For herds with fewer than 20 lactating camels, all animals were sampled. In herds with 20 or more camels every second camel was sampled until a target of 10 camels per herd was reached. Milk samples were frozen at -18 °C to -20 °C for 1 to 7 days prior to culturing. All samples were cultured on blood agar (CM0271, Oxoid, Thermo Fisher Scientific, Waltham, MA) containing 5 % defibrinated sheep blood, and Edwards agar (Oxoid, CM0027) and incubated aerobically at 37℃. After 24 to 48 h, primary species identification was based on colony morphology and catalase testing. Species confirmation was conducted using MALDI-ToF mass spectrometry (MS). Out of 804 quarter milk samples collected across all herds, 154 samples from 65 animals in 19 herds were GBS-positive. One isolate per animal was arbitrarily selected for further analysis. Most (n = 53 of 65) isolates were collected from CMT-positive udder quarters, with the remainder (n = 12) collected from CMT-negative quarters. Forty quarters had SCM (CMT-positive but no palpatory or visual changes of the udder or milk) and 13 had signs of CM (visible or palpable abnormalities in the milk or udder). Five CM cases were classified as acute CM (ACM) based on swelling, pain, redness and/or abnormal milk), with the remainder classified as chronic CM (CCM) (based on induration of the udder). Distribution of mastitis categories is shown in Fig. 1. One isolate per camel (1 to 8 isolates per herd, additional file 1) was arbitrarily selected for sequencing and further analysis (additional file 2).

Fig. 1
figure 1

Flowchart of sample and mastitis classification: normal milk (CMT < 3), mastitis (defined as a CMT-score ≥ 3), subclinical mastitis (defined as a CMT- score ≥ 3 and no palpatory or visual changes of the udder or milk), clinical mastitis (defined as CMT-score ≥ 3 and abnormalities in the milk or udder), acute clinical mastitis (ACM) (defined as a CMT-score ≥ 3 and swelling, pain, redness or abnormal milk), and chronic clinical mastitis (CCM) (defined as a CMT-score ≥ 3 and induration of the udder)

Potential associations between categorical variables were investigated using Pearson’s chi2 test and Fisher’s exact test. A p-value of < 0.05 was considered significant. Statistical analyses were performed in Stata (Stata; Stata Statistical Software, release 13.1; StataCorp LP, College Station, TX).

DNA extraction and sequencing

For DNA extraction, GBS colony material was collected with a calibrated loop (1 µl) and suspended in 600 µl nuclease free water (Sigma-Aldrich, St Louis, MO, USA). The suspension was mixed with 0.1 mm silica beads (BioSpec Products Inc., Bartlesville, USA) and added to the FastPrep24 (MP Biomedicals LLC, Irvine, CA, USA) and then run at 6.5 m/s for three 2-minute cycles. DNA was extracted from 200 µl samples using the IndiMag Pathogen kit (Indical Bioscience GmbH, Leipzig, Germany) and eluted in nuclease free water. DNA concentration was measured using the Invitrogen Qubit 3.0 Fluorometer (ThermoFisher Scientific Inc., Waltham, MA, USA), and adjusted to 7.5 ng/µL. Library preparation and whole genome sequencing were performed by Clinical Genomics, Science for Life Laboratory (Clinical Genomics, Solna, Sweden) on the Illumina NovaSeq (Illumina, Inc. CA, US) resulting in paired end reads of 150 bp in length.

Antimicrobial susceptibility and virulence testing

Phenotypic antimicrobial susceptibility testing was carried out as previously described [8]. Minimum inhibitory concentrations were determined using broth microdilution. Testing was performed according to the recommendations of the Clinical and Laboratory Standards Institute [38] using VetMIC panels (SVA, Uppsala, Sweden). SRST2 v0.2.0 [39] was used to detect antimicrobial resistance (AMR) genes from raw sequence reads with the ARG-ANNOT v3 database [40].

Lactose fermentation was assessed phenotypically by inoculating each of the selected isolates onto bromocresol purple lactose agar (SVA, Uppsala, Sweden). A yellow colour change of the colonies indicated lactose fermentation whereas purple colonies were classified as negative for lactose fermentation. Escherichia coli ATCC35218 and Proteus mirabilis CCUG26767 were used as positive and negative controls. Plates were incubated aerobically at 37℃ and checked for colour change at 24, 48 and 72 h. Assembled genomes (see below) were scanned for the presence of the lactose operon (Lac.2) [28] with a BLASTn v2.9.0 search [41] based on a database of four known Lac.2 genotypic variants [27, 42]. Minimum thresholds for identity (ID) and query coverage (QC) were both set at 90 %. Lac.2-negative isolates based on BLAST searches were further scanned for annotations related to Lac.2 in files obtained using Prokka v1.14.6 [43]. To confirm presence/absence of the Lac.2 operon in the chromosome, a PCR targeting a ≈ 2.5-kbp region straddling lacEG was used. Positive and negative controls were selected from the study isolates based on genomic detection of Lac.2.

Assembled genomes were scanned for the presence of human-associated virulence genes scpB and lmb [44] using tBLASTn. Capsular serotyping was conducted in silico using a standard method [45].

Phylogenetic and statistical analysis

Reads were filtered for quality and trimmed with ConDeTri suite v2.3 [46]. De novo assembly was performed using SPAdes v3.13.1 [47]. Assembly quality was checked with QUAST v5.0.2 [48] and species identity was confirmed with KmerFinder v3.2 [49]. All assembled genomes passed quality control. Multi locus sequence typing (MLST) was carried out with SRST2 and new alleles were submitted for allele number and ST assignment through pubMLST [50].

A core genome alignment was obtained with Snippy v4.6.0 [51] using ILRI112, an ST617 isolate from a Kenyan camel, as reference genome (accession HF952106). A maximum likelihood tree was inferred with RAxML-NG v0.9.0 [52] under a GTR + G model. A map of herd coordinates was created with ggplot in RStudio, with R (v4.0.). All figures were edited using Inkscape [53]).

Pairwise single nucleotide polymorphism (SNP)-distances between ST616 genomes were calculated using pairsnp v0.2.0 [54] Within-herd pairwise distances and between-herd pairwise distances were plotted using matplotlib v3.3.2 [55].


Sequence types and serotypes associated with camel mastitis

All 65 sequenced isolates were confirmed as GBS. Among them, only six STs were identified. These included three novel STs that were single locus variants (SLVs) of STs already described (ST1652, SLV of ST617; ST1653 and ST1654, SLVs of ST616). The vast majority of isolates (n = 54) belonged to ST616, which was found in all herds but one, followed by ST1652 (n = 6) found in five herds, ST612 (n = 2) found in two herds, ST1653 (n = 1), and ST1654 (n = 1), all belonging to previously described camel-associated clonal complexes. One isolate, however, belonged to a completely unrelated clonal complex (CC1) and was identified as ST1. In most herds (n = 12), isolates belonged to a single ST, but multiple STs were detected in the remaining seven herds. In five herds (herds A, B, D, J and N), two STs were present and in two herds (P and R), isolates belonged to three STs.

Four capsular serotypes were detected in silico (III to VI). The most common serotype was serotype III (n = 56), followed by serotype VI (n = 6), serotype IV (n = 2) and serotype V (n = 1). There was perfect concordance between STs and serotypes, with ST616 and its SLVs ST1653 and ST1654 all belonging to serotype III, ST612 isolates to serotype IV and ST1652 to serotype VI. Serotype III was significantly associated with SCM (36 of 56 serotype III isolates) (Fisher’s exact test; p = 0.04).

Antimicrobial susceptibility and virulence testing

Phenotypic susceptibility testing revealed tetracycline resistance (minimum inhibitory concentration (MIC) above 1 µg/mL) in 63 of 65 GBS isolates (MIC-values for tetracycline shown in Additional file 2, for details on the susceptibility testing, please see [8]). The tet(M) gene was found in all phenotypically tetracycline-resistant isolates and was absent from the two isolates that were phenotypically susceptible to tetracycline. Apart from the tet(M) gene, no more genes coding for antibiotic resistance were found within the GBS genomes.

Lactose fermentation was detected phenotypically in 75 % (n = 49) of the isolates. PCR-results were in 100 % agreement with the findings of phenotypic lactose fermentation and the presence of a lactose operon in the genomic analysis. There was an association between ST and lactose fermentation (Fisher’s exact test; p = 0.002) with the majority of ST616 (n = 45) and ST1653 (n = 1) being lactose fermenters. Three variants of the lactose operon (Lac.2) were detected, including two known genotypic variants (Lac.2b, n = 8; Lac.2d, n = 26) and a new variant (n = 15) that was named Lac.2e. Lac.2e (length = 9,535 bp) has the same gene arrangement as Lac.2a, with the exception of an additional gene in Lac.2e, a glucokinase (glk, length = 951 bp), upstream lacA. The Lac.2-variants were located at several insertion sites (Lac.2b, n = 2; Lac.2d, n = 1; Lac.2e n = 2), which largely agreed with their integrase types (additional file 3). Multiple STs were associated with each Lac.2 type and multiple Lac.2 types were detected within STs (Table 1) but there was no association between serotype and lactose operon. For mastitis isolates, i.e. those associated with the presence of an inflammatory response, Lac.2d was overrepresented (25 of 41 mastitis-derived isolates), (Pearson’s chi2-test; p = 0.001) and the same association was found for isolates from quarters with SCM (19 of 31 SCM isolates), (Pearson’s chi2-test; p = 0.001).

Table 1 Distribution of lactose operons across sequence types (STs) in 65 GBS isolates from camel milk

The human-associated virulence genes scpB and lmb were detected in a single genome assembly, which belonged to ST1, as part of a known composite transposon [56]. Two scpB variants were found within this mobile element, with scpB1 (1,071 bp) upstream scpB2 (2,214 bp), followed by lmb (additional file 4).

Phylogenetic analysis

After initial inspection, the isolate belonging to ST1 was removed from the phylogenetic tree to better visualise the relatedness among the other isolates. Three main lineages were observed in the core genome phylogenetic tree (Fig. 2). Two lineages corresponded to a single ST (ST612, ST1652), whereas the third and largest lineage included ST616 and its SLVs, ST1653 and ST1654. 

Fig. 2
figure 2

Maximum likelihood phylogenetic tree based on a core genome alignment of 65 camel group B Streptococcus (GBS) isolates. Leaf colours correspond to herd of origin (n = 19), whose locations are indicated on the map. Sequence types (ST) are shown on the branches, with the ST616 nomenclature also comprising its two single locus variants (ST1653 and ST1654). Grey bars indicate Lac.2 genotypic variants. A single genome assembly belonging to ST1 was removed from the tree to facilitate visualisation. Tree has been rooted at midpoint. ILRI112 (ST617; accession HF952106) was used as reference genome. The map was created with ggplot in RStudio, with R (v4.0.)

Within the major lineage, pairwise genetic distances (i.e. number of single nucleotide polymorphisms) between isolates did not differ for within-herd (mean = 51.33, standard deviation = 24.01) versus between-herd (mean = 57.86, standard deviation = 20.29) comparisons (Fig. 3).

Fig. 3
figure 3

Frequency distribution of pairwise genetic distances calculated between genomes from the largest group B Streptococcus lineage from camel milk (sequence types 616, 1653 and 1654) belonging to the same herd (within-herd diversity) or to different herds (between-herd diversity)


Genomic analyses show the existence of a predominant niche-adapted strain in camel milk

Here, we describe the genomic diversity of GBS from camel milk in Kenyan dairy herds and demonstrate that the main strain responsible for mastitis shows genetic signatures of adaptation to the mammary gland. ST616 serotype III was the predominant lineage among milk isolates, similar to previous findings [17]. The presence of Lac.2, which always corresponded to phenotypic lactose fermentation in GBS, suggests niche-adaptation to the mammary gland and has been described for GBS [19] as well as for Gram-negative mastitis pathogens [57]. The Lac.2 operon in cattle GBS shows signatures of genetic mobility (e.g. integrase) [28] and several possible insertion sites [42], and in camel GBS, Lac.2 showed similar mobility features. Isolates belonging to STs other than ST616 were less likely to ferment lactose which could indicate infection by strains from extramammary sources, a route of exposure that has also been suggested for dairy cattle [20]. Some non-ST616 isolates encoded Lac.2, which may have been acquired through horizontal gene transfer (HGT) as a means of adaptation to the mammary gland, again as seen in cattle GBS isolates [28]. To further explore infection dynamics within camel herds, within host-diversity needs to be investigated as well as the prevalence and molecular characteristics of extramammary GBS isolates in camels.

Tetracycline resistance was common among the investigated isolates and the tet(M) gene was found in all phenotypically resistant isolates. The presence of tet(M) has also been described in GBS populations in humans, fish, and, to a much lesser extent, cattle [26, 58, 59]. It has been considered a marker or even driver of expansion of GBS in the human population, starting in the 1940s. Its presence in piscine and bovine isolates has been interpreted as evidence of a human origin of animal GBS [42], but only in isolates from lineages that are shared between host species. Camel GBS is unique and, with the exception of the ST1 isolate reported here, there is no evidence of strain sharing between camels and humans. Considering that tetracycline is commonly used in pastoralist camel herds [13], its presence is likely due to acquisition of mobile genetic elements under selective pressure [17, 60].

Camel GBS isolates from different herds are genetically related

This is the first report on GBS genomic variation within camel herds. Within-herd heterogeneity of GBS was detected both at ST level (based on ca. 3500 nucleotides from 7 partial housekeeping genes) and within STs (based on single nucleotide polymorphisms detected through whole genome sequencing), with similar heterogeneity within and between herds. Heterogeneity may result from within-herd evolution, as described for GBS in dairy cattle [18]. Considering the finding of GBS in camels with evidence of chronic mastitis such as induration, and the high within-herd prevalence [8], within-camel and within-herd evolution are plausible. In addition, heterogeneity may reflect the occurrence of multiple transmission events between herds. In dairy cattle farms in high income countries, GBS is commonly introduced through the purchase of infected animals [61]. In Kenya, lactating female camels are rarely sold due to their high economic value as milk producers [62]. However, camel dairy herds are not closed, and females are frequently moved between groups depending on reproductive status. Additionally, herds could come into contact with each other during the daytime while browsing or at common watering points. Moreover, herdsmen in commercial herds may move between herds or share duties when camping together, increasing the likelihood of transmission of GBS from one herd to another, especially if not following proper hygiene practices (e.g. washing hands) (personal communication, Yussuf Maalim, Kenya Camel Association). In addition to between-herd biosecurity, within-herd biosecurity is affected by poor hygiene, leading to transmission of GBS between camels. Lack of water and inadequate milking hygiene are common issues in camel pastoralist herds and significant risk factors for mastitis [63, 64]. Fundamental milking hygiene practices, such as cleaning of hands before milking or using gloves, washing of the udder prior to milking or using post-milking teat disinfectant, were lacking in all of the herds included in this study [8]. Some herds would use a milking order, which could be a mastitis-reducing intervention [65], however, this alone is insufficient to halt contagious transmission of GBS without other hygiene measures, which would also require awareness of the existence of SCM, which, by definition, is not visible with the naked eye. Suckling calves may cause contamination of the camel udder [62] but they are unlikely to contribute to camel-to-camel transmission because cross-suckling is thought to be uncommon among camels [66].

Potential interspecies transmission

This is the first reported case of GBS ST1 isolated from a camel. ST1 is commonly isolated from humans, including healthy carriers, although there are no reports on distribution of human carriage strains in Kenya [67, 68]. This observation may represent interspecies transmission as also observed for cattle [19, 22, 27] and fish [26]. Contamination of the milk sample with a human GBS carriage isolate during the sampling procedure cannot be ruled out, although the isolate was collected from an udder quarter with signs of mastitis (CMT4) and using aseptic sampling technique, suggesting that the bacteria was present intramammarily. The fact that the scpB1, scpB2 and lmb-genes were found only in the ST1 isolate further points towards this isolate being of human origin, as those genes are specifically associated with the human host [29]. The possibility of interspecies transmission could interfere with camel mastitis control strategies as well as pose a risk to the consumers of camel milk. The camel milk industry in Kenya is currently undergoing rapid changes [69] and urban herds kept under more intensive management are becoming more common. In the aquaculture sector, intensification has been associated with the emergence of zoonotic GBS [26] and this risk may also exist in other production animal industries. The direct exposure to GBS of people consuming unpasteurized camel milk warrants further investigation into the potential human health hazard of camel GBS. To investigate potential interspecies transmission, isolates should be collected and compared from camels and from humans living in close contact with camels.


We found that GBS isolated from camel milk collected in and around Isiolo town (Isiolo County, Kenya) belonged to one predominant sequence type, ST616 or its single locus variants, of which a large proportion showed signs of niche adaptation. In light of this, it is likely that mastitis-causing GBS strains in camels are largely udder-bound and this demonstrates the need for improved milking hygiene. The similarity in heterogeneity of isolates within and across herds suggests that internal and external biosecurity measures would be needed to reduce within- and between herd-transmission, respectively. Our finding of one likely human-derived isolate (ST1) in an infected udder strongly suggests that human to camel transmission is possible, and this potential risk should be further explored. Biosecurity is a cornerstone of disease control; in order for it to be feasible, sustainable and thus efficiently implemented in pastoralist settings, locally-appropriate interventions should be devised and tested.

Availability of data and materials

All data analysed during this study is included in this published article [and its supplementary information files 1–4]. Sequence data generated for this study is available at the European Nucleotide Archive (ENA) under study PRJEB43245 (accession numbers from ERS5810822 to ERS5810886).



Acute clinical mastitis


Antimicrobial resistance


Clonal complex


Chronic clinical mastitis


Clinical mastitis


California mastitis test


Group B Streptococcus


Minimum inhibitory concentration


Lactose operon

lmb :

Laminin-binding protein


Matrix assisted laser desorption ionization-time of flight mass spectrometry

scpB :

C5a peptidase


Subclinical mastitis


Sequence type


National Veterinary Institute


  1. Amwata DA, Nyariki DM, Musimba NRK. Factors influencing pastoral and agropastoral household vulnerability to food insecurity in the drylands of Kenya: a case study of Kajiado and Makueni Counties. J Int Dev. 2016;28(5):771–87.

    Article  Google Scholar 

  2. Elhadi YA, Nyariki DM, Wasonga OV. Role of camel milk in pastoral livelihoods in Kenya: contribution to household diet and income. Pastoralism. 2015 Apr 20;5(1):8.

  3. Bekele T, Lundeheim N, Dahlborn K. Milk production and feeding behavior in the camel (Camelus dromedarius) during 4 watering regimens. J Dairy Sci. 2011 Mar 1;94(3):1310–7.

  4. FAOSTAT [Internet]. [cited 2021 Mar 2]. Available from:

  5. Nyariki DM, Amwata DA. The value of pastoralism in Kenya: Application of total economic value approach. Pastoralism. 2019 Jul 17;9(1):9.

  6. Farah Z, Mollet M, Younan M, Dahir R. Camel dairy in Somalia: Limiting factors and development potential. Livest Sci. 2007 Jun 1;110(1):187–91.

  7. Obied AI, Bagadi HO, Mukhtar MM. Mastitis in Camelus dromedarius and the somatic cell content of camels’ milk. Res Vet Sci. 1996 Jul;61(1):55–8.

  8. Seligsohn D, Nyman A-K, Younan M, Sake W, Persson Y, Bornstein S, et al. Subclinical mastitis in pastoralist dairy camel herds in Isiolo, Kenya: Prevalence, risk factors, and antimicrobial susceptibility. J Dairy Sci. 2020 May;103(5):4717–31.

  9. Toroitich KC, Gitau GK, Kitala PM, Gitao GC. The prevalence and causes of mastitis in lactating traditionally managed one-humped camels (Camelus dromedarius) in West Pokot County, Kenya. Livest Res Rural Dev. 2017;29(4):62.

  10. Younan M. Parenteral treatment of Streptococcus agalactiae mastitis in Kenyan camels (Camelus dromedarius). Rev Elev Médecine Vét Pays Trop. 2002;55(3):177–81.

    Article  Google Scholar 

  11. Wahinya MM, Karuku KJ, Richard K, Muli MV, Younan M. Effect of subclinical mastitis caused by beta-haemolytic streptococci on milk yields in Kenyan camels (Camelus dromedarius). Bull Anim Health Prod Afr. 2014;62(3):267–74.

    Google Scholar 

  12. Farah Z, Fischer A. Milk and meat from the camel, Handbook on products and processing. Zürich, Switzerland: vdf Hochschulverlag AG and der ETH Zürich; 2004.

    Google Scholar 

  13. Lamuka PO, Njeruh FM, Gitao GC, Abey KA. Camel health management and pastoralists’ knowledge and information on zoonoses and food safety risks in Isiolo County, Kenya. Pastoralism. 2017 Aug 2;7(1):20.

  14. Koech LC, Irungu BN, Ng’ang’a MM, Ondicho JM, Keter LK. Quality and brands of amoxicillin formulations in Nairobi, Kenya. BioMed Res Int. 2020 Jun 13;2020:E7091278.

  15. Muthiani M. Factors influencing the influx of counterfeit medicines in Kenya: A survey of pharmaceutical importing small and medium enterprises within Nairobi. Int J Bus Soc Sci. 2012;3(11):87–96.

    Google Scholar 

  16. Omwenga I, Aboge GO, Mitema ES, Obiero G, Ngaywa C, Ngwili N, et al. Antimicrobial usage and detection of multidrug-resistant Staphylococcus aureus, including methicillin-resistant strains in raw milk of livestock from northern Kenya. Microb Drug Resist [Internet]. 2020 [cited 2020 Dec 10]; Available from:

  17. Fischer A, Liljander A, Kaspar H, Muriuki C, Fuxelius H-H, Bongcam-Rudloff E, et al. Camel Streptococcus agalactiae populations are associated with specific disease complexes and acquired the tetracycline resistance gene tetM via a Tn916-like element. Vet Res. 2013 Oct 1;44(1):86.

  18. Mahmmod YS, Klaas IC, Katholm J, Lutton M, Zadoks RN. Molecular epidemiology and strain-specific characteristics of Streptococcus agalactiae at the herd and cow level. J Dairy Sci. 2015 Oct;98(10):6913–24.

  19. Lyhs U, Kulkas L, Katholm J, Waller KP, Saha K, Tomusk RJ, et al. Streptococcus agalactiae serotype IV in humans and cattle, northern Europe. Emerg Infect Dis. 2016 Dec;22(12):2097–103.

  20. Jørgensen HJ, Nordstoga AB, Sviland S, Zadoks RN, Sølverød L, Kvitle B, et al. Streptococcus agalactiae in the environment of bovine dairy herds – rewriting the textbooks? Vet Microbiol. 2016 Feb 29;184:64–72.

  21. Baseggio N, Mansell PD, Browning JW, Browning GF. Strain differentiation of isolates of streptococci from bovine mastitis by pulsed-field gel electrophoresis. Mol Cell Probes. 1997 Oct 1;11(5):349–54.

  22. Cobo-Angel CG, Jaramillo-Jaramillo AS, Palacio-Aguilera M, Jurado-Vargas L, Calvo-Villegas EA, Ospina-Loaiza DA, et al. Potential group B Streptococcus interspecies transmission between cattle and people in Colombian dairy farms. Sci Rep. 2019 Dec;9(1):14025.

  23. Shome BR, Bhuvana M, Mitra SD, Krithiga N, Shome R, Velu D, et al. Molecular characterization of Streptococcus agalactiae and Streptococcus uberis isolates from bovine milk. Trop Anim Health Prod. 2012 Dec;44(8):1981–92.

  24. Heffernan C, Misturelli F. The delivery of veterinary services to the poor: findings from Kenya [Internet]. Reading, UK: Veterinary and Economics Research Unit, Department of Agriculture, University of Reading.; 2000 [cited 2020 Dec 18]. Available from:

  25. Wolff C, Abigaba S, Sternberg Lewerin S. Ugandan cattle farmers’ perceived needs of disease prevention and strategies to improve biosecurity. BMC Vet Res. 2019 Jun 21;15(1):208.

  26. Barkham T, Zadoks RN, Azmai MNA, Baker S, Bich VTN, Chalker V, et al. One hypervirulent clone, sequence type 283, accounts for a large proportion of invasive Streptococcus agalactiae isolated from humans and diseased tilapia in Southeast Asia. PLoS Negl Trop Dis. 2019 Jun 27;13(6):e0007421.

  27. Sørensen UBS, Klaas IC, Boes J, Farre M. The distribution of clones of Streptococcus agalactiae (group B streptococci) among herdspersons and dairy cows demonstrates lack of host specificity for some lineages. Vet Microbiol. 2019 Aug 1;235:71–9.

  28. Richards VP, Lang P, Pavinski Bitar PD, Lefébure T, Schukken YH, Zadoks RN, et al. Comparative genomics and the role of lateral gene transfer in the evolution of bovine adapted Streptococcus agalactiae. Infect Genet Evol. 2011 Aug 1;11(6):1263–75.

  29. Sørensen UBS, Poulsen K, Ghezzo C, Margarit I, Kilian M. Emergence and global dissemination of host-specific Streptococcus agalactiae clones. mBio [Internet]. 2010 Aug 31 [cited 2020 Dec 18];1(3). Available from:

  30. Mati BM, Muchiri JM, Njenga K, Penning de Vries F, Merrey DJ. Assessing water availability under pastoral livestock systems in drought-prone Isiolo District, Kenya. Colombo, Sri Lanka: International Water Management Institute; 2005.

    Google Scholar 

  31. Isiolo County government, Isiolo County Sectoral Plan 2018–2027 [Internet]. Nairobi, Kenya: Department of Agriculture, Livestock & Fisheries Development, Republic of Kenya; 2019 [cited 2021 Mar 10]. Available from:

  32. Musinga M, Kimenye D, Kivolonzi P. The camel milk industry in Kenya. Resource Mobilization Centre, Nairobi, Kenya: Netherlands Development Organisation (SNV); 2008.

  33. Kenya National Bureau of Statistic. Kenya Population and Housing Census Volume I: Population by County and Sub-County 2019 [Internet]. Nairobi, Kenya: Ministry of Finance and Planning; 2019 [cited 2021 Mar 10]. Available from:

  34. KoBoToolbox | Data Collection Tools for Challenging Environments [Internet]. KoBoToolbox. [cited 2021 Apr 29]. Available from:

  35. Lamuka PO, Njeruh FM, Gitao GC, Matofari J, Bowen R, Abey KA. Prevalence of bovine and avian tuberculosis in camel herds and associated public health risk factors in Isiolo County, Kenya. Trop Anim Health Prod. 2018 Jun 1;50(5):937–45.

  36. Schalm OW, Noorlander DO. Experiments and observations leading to development of the California mastitis test. J Am Vet Med Assoc. 1957 Mar 1;130(5):199–204.

  37. Klastrup O, Madsen PS. Nordiske rekommendationer vedrorende mastitisundersogelser af kirtelprover (Nordic recommendations concerning mastitis control of quarter samples). Nord Veterinaermed. 1974;26:197–204.

    CAS  Google Scholar 

  38. Clinical and Laboratory Standards Institute (CLSI). Performance standards for antimicrobial susceptibility testing: M100, S27. 28th ed. Wayne, PA; 2017.

  39. Inouye M, Dashnow H, Raven L-A, Schultz MB, Pope BJ, Tomita T, et al. SRST2: Rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 2014 Nov 20;6(11):90.

  40. Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother. 2014 Jan 1;58(1):212–20.

  41. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009 Dec 15;10(1):421.

  42. Crestani C, Forde TL, Lycett SJ, Holmes MA, Fasth C, Persson-Waller K, et al. The fall and rise of group B Streptococcus in dairy cattle: reintroduction due to human-to-cattle host jumps? Preprint at (2021).

  43. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinforma Oxf Engl. 2014 Jul 15;30(14):2068–9.

  44. Morach M, Stephan R, Schmitt S, Ewers C, Zschöck M, Reyes-Velez J, et al. Population structure and virulence gene profiles of Streptococcus agalactiae collected from different hosts worldwide. Eur J Clin Microbiol Infect Dis. 2018 Mar 1;37(3):527–36.

  45. Metcalf BJ, Chochua S, Gertz RE, Hawkins PA, Ricaldi J, Li Z, et al. Short-read whole genome sequencing for determination of antimicrobial resistance mechanisms and capsular serotypes of current invasive Streptococcus agalactiae recovered in the USA. Clin Microbiol Infect. 2017 Aug 1;23(8):574.e7-574.e14.

  46. Smeds L, Künstner A. ConDeTri - A content dependent read trimmer for Illumina data. PLOS ONE. 2011 Oct 19;6(10):e26314.

  47. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012 Apr 16;19(5):455–77.

  48. Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013 Apr 15;29(8):1072–5.

  49. Cineros JLB, Lund O. KmerFinderJS: A client-server method for fast species typing of bacteria over slow Internet connections. bioRxiv. 2017 Jun 2;145284.

  50. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the website and their applications. Wellcome Open Res. 2018 Sep 24;3:124.

  51. Seemann T. tseemann/snippy [Internet]. 2021 [cited 2021 Apr 29]. Available from:

  52. Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019 Nov 1;35(21):4453–5.

  53. Draw Freely | Inkscape [Internet]. [cited 2021 Apr 29]. Available from:

  54. Tonkin-Hill G. gtonkinhill/pairsnp [Internet]. 2021 [cited 2021 Apr 29]. Available from:

  55. Hunter JD. Matplotlib: A 2D graphics environment. Comput Sci Eng. 2007 May;9(3):90–5.

  56. Franken C, Haase G, Brandt C, Weber-Heynemann J, Martin S, Lämmler C, et al. Horizontal gene transfer and host specificity of beta-haemolytic streptococci: the role of a putative composite transposon containing scpB and lmb. Mol Microbiol. 2001;41(4):925–35.

  57. Holt KE, Wertheim H, Zadoks RN, Baker S, Whitehouse CA, Dance D, et al. Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc Natl Acad Sci. 2015 Jul 7;112(27):E3574–81.

  58. Da Cunha V, Davies MR, Douarre P-E, Rosinski-Chupin I, Margarit I, Spinali S, et al. Streptococcus agalactiae clones infecting humans were selected and fixed through the extensive use of tetracycline. Nat Commun. 2014 Aug 4;5(1):4544.

  59. Duarte RS, Bellei BC, Miranda OP, Brito MAVP, Teixeira LM. Distribution of antimicrobial resistance and virulence-related genes among Brazilian group B streptococci recovered from bovine and human sources. Antimicrob Agents Chemother. 2005 Jan 1;49(1):97–103.

  60. León-Sampedro R, Novais C, Peixe L, Baquero F, Coque TM. Diversity and evolution of the Tn5801-tet(M)-like integrative and conjugative elements among Enterococcus, Streptococcus, and Staphylococcus. Antimicrob Agents Chemother. 2016 Mar 1;60(3):1736–46.

  61. Agger J, Priou C, Huda A, Aagaard K. Risk factors for transmission of Streptococcus agalactiae infection between Danish dairy herds: a case control study. Vet Res. 1994;25(2–3):227–34.

    CAS  PubMed  Google Scholar 

  62. Noor IM, Guliye AY, Tariq M, Bebe BO. Assessment of camel and camel milk marketing practices in an emerging peri-urban production system in Isiolo County, Kenya. Pastor Res Policy Pract. 2013 Dec 2;3(1):28.

  63. Kashongwe OB, Bebe BO, Matofari JW, Huelsebusch CG. Associations between milking practices, somatic cell counts and milk postharvest losses in smallholder dairy and pastoral camel herds in Kenya. Int J Vet Sci Med. 2017 Jun 1;5(1):57–64.

  64. Odongo N, Matofari J, Abong G, Lamuka P, Abey K. Knowledge and practices of food hygiene and safety among camel milk handlers in the pastoral camel value chain in Kenya. Afr J Food Agric Nutr Dev. 2017 Mar 23;17:11803–21.

  65. Hutton CT, Fox LK, Hancock DD. Mastitis control practices: differences between herds with high and low milk somatic cell counts. J Dairy Sci. 1990 Apr 1;73(4):1135–43.

  66. Packer C, Lewis S, Pusey A. A comparative analysis of non-offspring nursing. Anim Behav. 1992 Feb 1;43(2):265–81.

  67. Huber CA, McOdimba F, Pflueger V, Daubenberger CA, Revathi G. Characterization of invasive and colonizing isolates of Streptococcus agalactiae in East African Adults. J Clin Microbiol. 2011 Oct 1;49(10):3652–5.

  68. Seale AC, Koech AC, Sheppard AE, Barsosio HC, Langat J, Anyango E, et al. Maternal colonization with Streptococcus agalactiae and associated stillbirth and neonatal disease in coastal Kenya. Nat Microbiol. 2016 May 23;1(7):1–10.

  69. Anderson DM, Elliott H, Kochore HH, Lochery E. Camel herders, middlewomen, and urban milk bars: the commodification of camel milk in Kenya. J East Afr Stud. 2012 Aug 1;6(3):383–404.

Download references


The authors would like to extend their gratitude to the pastoralists that participated in the study. We thank Mattias Myrenås (Department of Animal Health and Antimicrobial Strategies, National Veterinary Institute, Sweden) for assistance with PCR-analysis and invaluable help. We also thank Ann Nyman (Växa Sweden, Stockholm, Sweden) for statistical support.


This study was funded by the Swedish Research Council [grant number VR2015_03583]. The funding body was not involved in the study design, data collection and analysis, interpretation of data or in writing the manuscript. Open Access funding provided by Swedish University of Agricultural Sciences.

Author information

Authors and Affiliations



DS, CC, EC, TF and RZ designed the study. DS performed the field work, the laboratory analyses and the phenotypic characterization. CC carried out WGS analysis. DS drafted the manuscript with input from CC on genomic and phylogenetic analysis. All authors read, critically reviewed and approved the final manuscript.

Authors’ information

Not applicable.

Corresponding author

Correspondence to Dinah Seligsohn.

Ethics declarations

Ethics approval and consent to participate

The isolates used in this study were collected as a part of a prevalence study [8] which was approved by the National Commission for Science, Technology and Innovation, Nairobi, Kenya (Permit number: NACOSTI/P/19/84995/13088).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Number of sampled camels and udder quarters and the number of Group B Streptococcus (GBS)-positive camels and quarters in 20 dairy camel herds in Isiolo county, Kenya. One isolate per camel was selected for whole genome sequencing.

Additional file 2.

Authors’ original data for 65 GBS-isolates collected from camel milk in Isiolo county, Kenya, February-May, 2017, including accession numbers available at European Nucleotide Archive (ENA).

Additional file 3: Figure S1.

Maximum likelihood phylogenetic tree of 49 Lac.2 integrase amino-acid sequences extracted from group B Streptococcus genomes from Kenyan camels. Sequences were aligned using MAFFT v7.475 and phylogeny was estimated with PhyML v3.3.20190909. Isolate names are shown. Leaf colours correspond to the insertion site where Lac.2 is integrated (light green: deoD, red: yxdL, blue: hypothetical, purple: ClbS/DfsB family four-helix bundle protein). For 30 isolates, the Lac.2 insertion site could not be determined because their integrases were found at the edge of a contig (dark green leaves).

Additional file 4: Figure S2.

Organisation of genes of the scpB-lmb mobile transposon as found in group B Streptococcus isolate P4 (sequence type 1) in milk from a camel (Camelus dromedarius) in Kenya. In this transposon, two variants of C5a peptidase gene scpB are present (scpB1 and scpB2) upstream the laminin binding protein gene lmb.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Seligsohn, D., Crestani, C., Forde, T.L. et al. Genomic analysis of group B Streptococcus from milk demonstrates the need for improved biosecurity: a cross-sectional study of pastoralist camels in Kenya. BMC Microbiol 21, 217 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: