All relevant data are within the manuscript and its Supporting Information files.

Introduction {#sec004}
============

In 2015, there was an unprecedented outbreak of invasive disease due to the group B *Streptococcus* (GBS; *Streptococcus agalactiae*) in Singapore. Disease, including septic arthritis and meningitis, was associated with a GBS belonging to serotype III, subtype 4 (serotype III-4), and multilocus sequence type (MLST) 283 (ST283) \[[@pntd.0007421.ref001]\]. Human disease was associated with consumption of raw, farmed, freshwater fish: an official public advisory was issued, and bacteraemia rates promptly fell \[[@pntd.0007421.ref002], [@pntd.0007421.ref003]\]. Whereas GBS sepsis in neonates, post-partum adults, and older adults with co-morbidities is well known, the ST283 outbreak was different as it affected non-pregnant, younger adults with fewer co-morbidities \[[@pntd.0007421.ref001]\], suggesting greater virulence. Although GBS colonisation has previously been associated with fish consumption \[[@pntd.0007421.ref004]\], and the origin of GBS in late onset neonatal disease, while uncertain, may be enteric \[[@pntd.0007421.ref005], [@pntd.0007421.ref006]\], this was the first report of an invasive GBS outbreak associated with foodborne transmission.

Previous reports of ST283 were limited to 20 cases in adults in Hong Kong between 1993 and 2003 \[[@pntd.0007421.ref007]\], two cases in a survey of 119 osteoarticular GBS infections in France between 2002 and 2007 \[[@pntd.0007421.ref008]\], and an infected tilapia (*Oreochromis* sp.) amongst samples collected in Thailand between 2000 and 2010 \[[@pntd.0007421.ref009]\]. Southeast (SE) Asia was, however, under-represented in the literature, with MLST data reported for only a handful of human GBS from SE Asia amongst thousands in global studies.

We hypothesised that fish-borne ST283 might be a regional problem, causing disease in multiple SE Asian countries. As a prelude to studying transmission, we searched for ST283 and its variants, collectively called clonal complex 283 (CC283), in GBS collections from humans and aquaculture. Our aim was achieved, as we verified the wide prevalence of ST283 in Southeast Asia, though whether human disease is fish-borne outside Singapore remains an open question.

Methods {#sec005}
=======

The study was approved by the Institutional Review Board of Tan Tock Seng Hospital (TTSH), Singapore, NHG DSRB 2016/00202.

Literature search {#sec006}
-----------------

We identified articles in English with searches of Medline, PubMed, and references from articles, with the terms "*Streptococcus agalactiae*" OR "ST283" AND "meningitis" OR "invasive", up to December 2017.

GBS collections {#sec007}
---------------

Institutions, selected through personal contacts and recommendations, were invited to contribute GBS and/or datasets. There were no specified criteria other than the availability of isolates or DNA from invasive GBS, or GBS from high vaginal swabs, over time, with metadata if possible. Invasive disease was defined as isolation of GBS from normally sterile sites. Indications for testing patients and accompanying characterisation data were not standardised or complete; rates of meningitis, osteoarthritis, and endocarditis may therefore be greater than recorded. Country-specific regulations resulted in some contributors sending GBS isolates to Singapore for processing while others contributed extracted DNA, whole genome sequencing (WGS) data, or other data sets. GBS were screened with a ST283-specific PCR (Lezhava A., Sarma S., Chen S. and Barkham T.M.S. A method for the detection of Group B Streptococcus. (2017) PCT Patent Application PCT/SG2017/050579). This PCR has been evaluated against WGS on over 660 invasive GBS representing 27 MLST types collected over 18 years in Singapore and found to be 100% sensitive and 99.8% specific. GBS identification, determination of serotype and WGS were performed as previously described \[[@pntd.0007421.ref001]\]. All new sequencing data were deposited in GenBank under BioProject PRJNA293392. We expressed ST283 prevalence as simple proportions of all GBS in each collection ([Fig 1](#pntd.0007421.g001){ref-type="fig"}).

![Prevalence of clonal complex 283 in human and animal collections of group B *Streptococcus* showing their host, geographic origin, and period of collection.\
This figure represents our new data, as well as the literature on *Streptococcus agalactiae* (GBS) that includes multi locus sequencing typing data up to December 2017. The vertical bars on the left indicate human or animal origin and the geographical region, Southeast Asia (SEA) or outside SEA (Ex- SE Asia), where collections of GBS originated. The horizontal bars delineate the time period of each collection of GBS, with reference to the central time bar; associated text shows the host, the country of origin, and number of ST283/all GBS in each collection, except where other STs are shown. The figure shows the lack of reports, from humans and animals, in SEA compared with outside SEA, both in terms of time periods and absolute numbers of GBS studied. It also shows that GBS CC283 is rare outside SEA, with only four human and no animal CC283 reported, despite the large number of GBS studied. In contrast, CC283 is prevalent in all human and animal GBS collections from SEA. \*One of these 34 from a tilapia in Vietnam is ST1311, a double locus variant of ST283. \*\* Australia, Ghana, South America, North America, Israel and Kuwait. Abbreviation: ND = not determined.](pntd.0007421.g001){#pntd.0007421.g001}

Human GBS {#sec008}
---------

### Singapore {#sec009}

Invasive bacterial isolates are routinely saved: 331 invasive GBS isolates from 2001 to 2010 were retrieved from freezers in TTSH and Changi General Hospital. Data and sequences for 408 invasive GBS ST283 isolated from 2011--2015 were publicly available \[[@pntd.0007421.ref001]\]. Sequences of five GBS, isolated in Singapore General Hospital in 1998 from cases of meningitis, are now available at GenBank accessions SRR6282417-SRR6282421 (under BioProject PRJNA417692).

### Thailand {#sec010}

139 invasive GBS isolates were available from population-based surveillance for blood stream infections conducted in all 20 hospitals in Nakhon Phanom and Sa Kaeo provinces from 2007--2015; these include paediatric patients. GBS isolates from the isolate repository were screened with the ST283 PCR described above. A selection of positive isolates were subjected to WGS on the MiSeq platform using the Nextera XT DNA Library Prep Kits and the MiSeq Reagent Kit v3. In addition, sequences of seven GBS ST283 from invasive human cases in Bangkok in 2015 were publicly available \[[@pntd.0007421.ref001]\].

### Lao PDR {#sec011}

This set comprised 38 GBS, available as isolates or as DNA extracts, obtained from blood cultures and CSF samples at the Microbiology Laboratory of Mahosot Hospital, Vientiane, Lao PDR between 2000--2017 during studies of the aetiology of fever and central nervous system infections.

### Vietnam {#sec012}

13 invasive GBS isolated between 2015 and 2017 were retrieved from the Hospital for Tropical Diseases (HTD), Ho Chi Minh City (HCMC). HTD is a 550-bed hospital that serves as a main primary and secondary facility for the surrounding local population in HCMC and a tertiary referral center for infectious diseases for the southern provinces of Vietnam. Nearly 70% of HTD admissions live in HCMC, with the remainder residing in the surrounding provinces. Neonates, patients without infectious diseases, including those with surgical requirements, tuberculosis, cancer, primary hematological disorders or immunosuppression (other than HIV) are referred to other hospitals within HCMC. HIV-infected children are often referred to local paediatric hospitals. Blood cultures were performed for patients in whom an infection was suspected on the basis of a fever (\>38°C) or who had evidence of sepsis on the basis of the presence of two or more of the following features: fever (\>38°C) or low temperature (\<36°C); tachycardia (exact level according to age); tachypnea (exact level according to age); an elevated white cell count (\>12,000 cells/mm3) or depressed white cell count (\<4,000 cells/mm3). There was no systematic change in the application of these criteria during the time course of the study. All data originating from consecutive patients admitted to the hospital who had a blood culture performed for suspected bloodstream infection between 2015 and 2017 were included. In addition, 38 GBS representing the most recent high vaginal swab GBS isolates from women with colpitis, three of whom were pregnant, were collected from outpatients at the National Hospital for Obstetrics & Gynaecology (Phu San Hospital), Hanoi, between September 2016 and May 2017.

### Cambodia {#sec013}

DNA from eleven GBS isolated from skin and umbilical swabs collected from children presenting to Angkor Hospital for Children, a non-governmental paediatric referral hospital in Siem Reap, between 2012 and 2016 were included.

### Britain {#sec014}

WGS data were available from 1,017 invasive human GBS submitted to the Respiratory and Vaccine Preventable Bacteria Reference Unit, Public Health England (PHE), as part of the British Society for Antimicrobial Chemotherapy Resistance Surveillance Project over a one year period bridging 2014 and 2015.

Piscine GBS {#sec015}
-----------

### Singapore {#sec016}

Sequences of GBS ST283 isolated from fish in 2015 were publicly available \[[@pntd.0007421.ref001], [@pntd.0007421.ref010]\].

### Malaysia {#sec017}

We sequenced 28 GBS, isolated between 2007 and 2008, from brain, eye or kidney samples from 28 tilapia (*Oreochromis* sp.) from farms suffering streptococcosis outbreaks. Ten of these 28 tilapia did not show any external or internal signs of disease. The fish were obtained from nine farms separated by 20--250 km in Kedah and Terengganu states, Peninsular Malaysia.

### Vietnam {#sec018}

We sequenced 34 GBS isolated in 2016 from the brains of sick fish from five farms in the Mekong River running through An Giang Province and Can Tho District. GBS was isolated directly from the brains of red or black tilapia (*Oreochromis* sp. and *O*. *niloticus*) that showed abnormal behaviour (erratic swimming pattern) or "pop eye" (exophthalmus), both of which are recognized clinical signs of streptococcosis.

### Other datasets {#sec019}

Other GBS MLST datasets were contributed by collaborators. All whole genome data sets annotated as *S*. *agalactiae* in GenBank deposited as assemblies (as of June 07, 2017) and in the SRA database (as of November 10, 2017) were downloaded.

### Genomic methods {#sec020}

We performed Illumina sequencing on 145 CC283 GBS strains. All new sequencing data was deposited in GenBank under BioProject PRJNA293392. All assembled genomes (as of June 7, 2017) and short read data sets (as of November 10, 2017) annotated as *S*. *agalactiae* in GenBank were also downloaded; this provided another 82 CC283 isolates, for a total of 227 CC283 isolates amongst 7,468 GBS. All primary sequence analysis was performed by the Genome Institute of Singapore Efficient Rapid Microbial Sequencing (GERMS) platform (<https://www.a-star.edu.sg/gis/Our-Science/Technology-Platforms/GERMS>). Reference-based analyses were performed using the SG-M1 genome \[[@pntd.0007421.ref011]\] as the reference. FASTQ files were mapped using bwa (version 0.7.10) \[[@pntd.0007421.ref012]\]; indel realignment and single nucleotide polymorphism (SNP) calling were performed using Lofreq\* (version 2.1.2) with default parameters \[[@pntd.0007421.ref013]\]. SNP positions for assembled genomes were inferred by using nucmer and show-snps from the MUMmer package (version 3.23) \[[@pntd.0007421.ref014]\]. MLST and resistance gene predictions were made using SRST2 0.2.0 \[[@pntd.0007421.ref015]\] for Illumina sequenced strains or using a custom BLASTN \[[@pntd.0007421.ref016]\] based script for fully assembled reference sequences, using the recommended MLST database (<http://pubmlst.org/sagalactiae/>) \[[@pntd.0007421.ref017]\] and the ARG-ANNOT resistance gene database \[[@pntd.0007421.ref018]\] included with SRST2.

### Recombination analysis {#sec021}

From an initial neighbour-joining tree made from all SNPs called relative to the SG-M1 reference, a set of strains containing all CC283 strains and strains from the nearest non-CC283 clade (consisting of 92 strains within CC10) was taken. We excluded those strains which had no metadata for isolation date, leaving a total of 273 strains (215 from the ST283 clade, 58 from the non-ST283 clade). We reconstructed an aligned genome sequence for each of these 273 strains by introducing that strain's SNPs into the SG-M1 genome; gaps were also inserted where mapping coverage was below 10 as per Lofreq default parameters. The BRATNextGen \[[@pntd.0007421.ref019]\] (using a clustering cut off of 0.2 and a significance cut off of 0.05) and ClonalFrameML \[[@pntd.0007421.ref020]\] (using default parameters) programs were used to call recombination on this set of 273 aligned genomes. We took the union of all recombined segments called by either program in any of the 273 strains as the maximal recombination set (in total, 608 kb or 28.7% of the total chromosome); these were removed for the subsequent Bayesian analysis. We also took the union of all recombined segments called by either program only in the 215 CC283 strains as a CC283-specific recombination set (totalling 62 kb or 3.0% of the total chromosome); removal of these segments from the reconstructed genomes of the original 227 CC283 strains (including those that had no available year or country metadata) resulted in a recombination-free CC283 alignment which was used to generate the maximum-likelihood tree shown in [Fig 2B](#pntd.0007421.g002){ref-type="fig"}.

![Phylogenetic analysis of group B *Streptococcus* (GBS) with emphasis on clonal complex (CC) 283 genomes.\
(A) Approximately maximum-likelihood phylogenetic tree of 1,236 GBS strains. The ST283 isolate SG-M1 was used as a reference sequence. The scale bar is shown on the x-axis, in mutations/nucleotide. The bootstrap values for selected branches (supporting the difference between different CCs) are indicated by black circles; all are 1.000 except where indicated. Major clonal complexes are indicated with coloured branches and a matching coloured label. The clade containing all CC283 isolates is highlighted in blue. All ST283 variants discussed in the text (ST491, ST739, ST751, and ST1311) are highlighted with black branches and an adjacent label in black text. (B) Approximately maximum-likelihood phylogenetic tree of 227 GBS CC283 genomes. Reconstructed genome sequences (based on the SG-M1 reference sequence) of the isolates indicated as CC283 in (A) were used, after excluding redundant isolates, defined as identical sequences from the same site, based on SNP calls. Bootstrap support is indicated for selected branches by black circles. Arrows show predicted events resulting in loss of tetracycline resistance. (C) For each CC283 isolate, the metadata are indicated at the same horizontal position (i.e. directly to the right of the phylogenetic tree tip) as in panel B. Host and country are represented by coloured rectangles, as indicated in the legend. Different values for Host and Country are further offset horizontally for clarity. Asterisks to the right of the Year box indicate isolates for which individual loss of tetracycline resistance appears to have occurred.](pntd.0007421.g002){#pntd.0007421.g002}

### Phylogenetic trees {#sec022}

Approximately maximum-likelihood SNP trees were created using FastTree 2.1.8 with the--gtr and--nt command line options \[[@pntd.0007421.ref021]\] on reconstructed genome sequences (relative to the SG-M1 reference genome). [Fig 2A](#pntd.0007421.g002){ref-type="fig"} was constructed from 1,236 strains (all CC10 and CC283 strains; all strains described in Da Cunha, et al 2014 \[[@pntd.0007421.ref022]\] that were available in GenBank; and 23 complete genomes in the GenBank Refseq database as of August 4, 2016) using all called SNPs (i.e., no removal of recombination). [Fig 2B](#pntd.0007421.g002){ref-type="fig"} was created from all CC283 strains after removal of CC283-specific recombination as described above (Recombination analysis). All phylogenetic trees were visualised with GGTREE 3.2 \[[@pntd.0007421.ref023]\] in R (3.2.2) ([https://www.R-project.org](https://imsva91-ctp.trendmicro.com/wis/clicktime/v1/query?url=https:%2f%2fwww.R-project.org&umid=21DE4B13-6DEE-4005-A39F-479A096547E6&auth=6e3fe59570831a389716849e93b5d483c90c3fe4-63fb27d76ec781c1338b476c0cb562dc5de69afd)). Origins of the GBS CC283 WGS sequences used for phylogenetic analysis are given in [Table 1](#pntd.0007421.t001){ref-type="table"}.

10.1371/journal.pntd.0007421.t001

###### New and Existing group B *Streptococcus* (GBS) clonal complex (CC) 283 sequences used for phylogenetic analysis.

This gives an overview of the origins and numbers of CC283. Further details are in Supporting Information.

![](pntd.0007421.t001){#pntd.0007421.t001g}

  Host                     Country       Site or reference                                    Collection Year   Sample type   GBS isolates screened for ST283 (No.)   GBS isolates with WGS data (No.)            CC283 isolates with WGS data (No.)
  ------------------------ ------------- ---------------------------------------------------- ----------------- ------------- --------------------------------------- ------------------------------------------- ------------------------------------------
  **New sequences**                                                                                                                                                                                               
  Human                    Singapore     Singapore General Hospital                           1998              B             5                                       5                                           5
  Human                    Singapore     Tan Tock Seng Hospital & Changi General Hospital     2001--2010        B, CSF        331                                     331                                         21
  Human                    Thailand      Nakhon Phanom and Sa Kaeo provinces                  2007--2015        B             139                                     22 [^a^](#t001fn001){ref-type="table-fn"}   22
  Human                    Lao PDR       Mahosot Hospital, Vientiane                          2000--2017        B, CSF        38                                      38                                          30[^b^](#t001fn002){ref-type="table-fn"}
  Human                    Vietnam       Hospital for Tropical Diseases, Ho Chi Minh City     2015--2017        B             13                                      13                                          4
  Human                    Vietnam       Phu San Hospital, Hanoi                              2016--2017        HVS           38                                      0 [^a^](#t001fn001){ref-type="table-fn"}    0
  Human                    Cambodia      Angkor Hospital for Children, Siem Reap              2012--2016        S             11                                      0 [^a^](#t001fn001){ref-type="table-fn"}    0
  Human                    Britain       Multiple hospitals                                   2015              B             ND                                      1,017                                       1
  Tilapia                  Malaysia      Kedah and Terengganu states                          2007--2008        Br, E, K      28                                      28                                          28
  Tilapia                  Vietnam       An Giang and Can Tho                                 2016              Br            34                                      34                                          34
  **Existing sequences**                                                                                                                                                                                          
  Human                    Singapore     \[[@pntd.0007421.ref001]\]                           2011--2015        B, CSF        ND                                      145                                         40
  Human                    Thailand      \[[@pntd.0007421.ref001]\]                           2015              B             ND                                      7                                           6
  Human                    Hong Kong     \[[@pntd.0007421.ref024]\]                           1993--2012        SS            ND                                      11                                          11
  Human                    Netherlands   GenBank                                              ND                ND            ND                                      1                                           1
  Fish                     Singapore     \[[@pntd.0007421.ref001], [@pntd.0007421.ref010]\]   2015              S, O, M       ND                                      20                                          20
  Fish                     Thailand      \[[@pntd.0007421.ref009]\]                           2000--2010        ND            ND                                      1                                           1
  Fish                     Vietnam       \[[@pntd.0007421.ref009]\]                           2006              ND            ND                                      1                                           1
  Fish                     Thailand      \[[@pntd.0007421.ref025]\]                           2012--2014        ND            ND                                      1                                           1
  Frog                     China         GenBank                                              2014              L             ND                                      1                                           1

^a^ = WGS was only performed on GBS positive by an ST283 specific PCR.

^b^ = Included one technical replicate. Where sequences from the same site were identical (based on SNP calls), only one was used. Abbreviations: No. = number, B = blood, CSF = cerebro-spinal fluid, HVS = high vaginal swab, S = superficial swabs, ND = not determined, SS = sterile site, Br = brain, E = eye, K = kidney, O = organs, M = muscle, L = liver.

### BEAST analysis {#sec023}

This was performed on all CC10 and CC283 strains described above in the recombination section. Analysis of divergence times was performed using BEAST 2.4.8 \[[@pntd.0007421.ref026]\] using a combination of site substitution models (HKY, GTR4), clock models (strict, relaxed exponential, relaxed lognormal), and population change models (coalescent constant, exponential, Bayesian skyline, and extended Bayesian skyline). MCMC runs were done in triplicate with 200,000,000 steps, using the first 20,000,000 as burn-in, sampling every 10,000 steps. Logcombiner was used to combine the three runs. Models were compared using AICM analysis in Tracer V1.6 \[[@pntd.0007421.ref027]\], from which the GTR4 substitution model with a strict clock and exponential model was chosen as the best set of parameters. Treeannotator was used to generate a final tree. Visualization was done in FigTree 1.4.3 (<http://tree.bio.ed.ac.uk/software/figtree/>). The divergence time reported was calculated after removing the maximal recombination set; since this was a significant fraction of the genome, we repeated the analysis using a data set from which only the ST283-specific recombination regions were removed. The 95% highest posterior density for the maximal recombination data set fell entirely within that for the ST283-specific recombination set, suggesting that the removal of additional recombination regions did not drastically bias this result. However, despite this technical agreement, the opportunistic nature of the source data for these predictions must be kept in mind when interpreting these results.

Results {#sec024}
=======

Human GBS CC283 in Asia {#sec025}
-----------------------

GBS ST283 accounted for 11% to 76% of invasive GBS per country, with the earliest known example collected in 1995, in Hong Kong. Five GBS ST11, single locus variants (SLV) of ST283, isolated from meningitis cases in Singapore in 1998 \[[@pntd.0007421.ref028], [@pntd.0007421.ref029]\] were recently corrected to ST283 \[[@pntd.0007421.ref030]\], leaving no known examples of ST11. In total, 29% of GBS were identified as ST283 ([Fig 1](#pntd.0007421.g001){ref-type="fig"}, [Table 2](#pntd.0007421.t002){ref-type="table"}). The majority of patients with ST283 (345/357 (97%) were adults, and 36% to 80% did not have comorbidities. Meningitis, endocarditis, and septic arthritis were noted in 10% to 35%, 4.5% to 10%, and 23% to 39% of adult patients with ST283, respectively ([S1 Table](#pntd.0007421.s001){ref-type="supplementary-material"}).

10.1371/journal.pntd.0007421.t002

###### Asian[^a^](#t002fn001){ref-type="table-fn"} Human group B *Streptococcus* (GBS) clonal complex (CC) 283 reported up to December 2017 as a proportion of invasive GBS by location and year.

All CC283 in this table are sequence type (ST) 283: other examples of CC283 were not found. This table shows that ST283 was found in the first year of all newly described GBS collections in Southeast Asia, so it may have predated these collections, and that both GBS numbers available and ST283 proportions vary from year to year.

![](pntd.0007421.t002){#pntd.0007421.t002g}

          Vientiane, Lao PDR   Nakhon Phanom and Sa Kaeo provinces, Thailand   Ho Chi Minh City, Vietnam   Singapore [^b^](#t002fn002){ref-type="table-fn"}   Total                                                                       
  ------- -------------------- ----------------------------------------------- --------------------------- -------------------------------------------------- ------- -------- ----- ---------------------------------------------- ----- ----------
  2000    1                    1 (100)                                         ND                          ND                                                 ND      ND       ND    ND                                             1     1
  2001    0                    0                                               ND                          ND                                                 ND      ND       16    3 (19)                                         16    3 (19)
  2002    0                    0                                               ND                          ND                                                 ND      ND       13    2 (15)                                         13    2 (15)
  2003    1                    1 (100)                                         ND                          ND                                                 ND      ND       28    2 (7)                                          29    3 (10)
  2004    0                    0                                               ND                          ND                                                 ND      ND       32    7 (22)                                         32    7 (22)
  2005    0                    0                                               ND                          ND                                                 ND      ND       23    0 (0)                                          23    0 (0)
  2006    3                    2 (67)                                          ND                          ND                                                 ND      ND       36    2 (6)                                          38    4 (11)
  2007    5                    5 (100)                                         7                           3 (43)                                             ND      ND       41    2 (5)                                          52    10 (19)
  2008    0                    0                                               10                          6 (60)                                             ND      ND       46    1 (2)                                          54    7 (13)
  2009    2                    1 (50)                                          1                           0                                                  ND      ND       53    1 (2)                                          55    2 (4)
  2010    7                    5 (71)                                          4                           3 (75)                                             ND      ND       43    1 (2)                                          52    9 (17)
  2011    1                    0 (0)                                           5                           4 (80)                                             ND      ND       38    0 (0)                                          44    4 (9)
  2012    1                    1 (100)                                         13                          8 (62)                                             ND      ND       43    3 (7)                                          57    12 (21)
  2013    3                    3 (100)                                         31                          27 (87)                                            ND      ND       58    9 (16)                                         92    39 (42)
  2014    2                    1 (50)                                          66                          50 (76)                                            ND      ND       77    15 (19)                                        145   66 (46)
  2015    4                    3 (75)                                          2                           1 (50)                                             4       1 (25)   110   61 (55)[\*](#t002fn003){ref-type="table-fn"}   120   66 (55)
  2016    7                    5 (71)                                          ND                          ND                                                 8       3 (38)   50    15 (30)                                        15    8 (53)
  2017    1                    1 (100)                                         ND                          ND                                                 1       0 (0)    ND    ND                                             2     1 (50)
  Total   38                   29 (76)                                         139                         102 (73)                                           13      4 (31)   707   124 (18)                                       840   244 (29)

^a^ Similarly detailed data for Hong Kong were not found.

^b^ This Singapore data is from the only two institutions that had systematically saved GBS; Tan Tock Seng Hospital and Changi General Hospital.

\*The outbreak year. Abbreviations: ND = not determined.

Piscine GBS CC283 in Asia {#sec026}
-------------------------

All newly sequenced GBS isolated from tilapia in Peninsular Malaysia and Vietnam, during outbreaks of streptococcosis, were CC283 ([Fig 1](#pntd.0007421.g001){ref-type="fig"}, [Table 3](#pntd.0007421.t003){ref-type="table"}). Notably, ten of the 28 fish sampled in Malaysia did not have signs of disease. ST283 was previously reported from 17/57 (30%) GBS isolated from seven freshwater fish species, including tilapia, on sale for human consumption in Singapore in 2015: these 17 ST283 were found in only 6/586 (1%) fish from ports, but in 11/39 (28%) fish from markets \[[@pntd.0007421.ref010]\]. Two previously reported ST283 were reported from farmed tilapia in Thailand between 2000 and 2014 \[[@pntd.0007421.ref009], [@pntd.0007421.ref025]\]. One ST491 and one ST1311, (a new SLV of ST491, and also a double locus variant of ST283) were isolated from tilapia in Vietnam in 2006 \[[@pntd.0007421.ref009]\] and 2016 respectively. An ST739, another SLV of ST283, was isolated in 2014 from a diseased tiger frog (*Hoplobatrachus chinensis*) farmed for human consumption in Guangdong, China, about 100Km from Hong Kong; although frogs are not piscine, they are often farmed in close proximity to fish and all known GBS infections in frogs are caused by clades that are also found in fishes \[[@pntd.0007421.ref009]\].

10.1371/journal.pntd.0007421.t003

###### Asian Piscine group B *Streptococcus* (GBS) clonal complex (CC) 283 reported up to December 2017.

All these CC283 were ST283, except for two variants, as indicated.

![](pntd.0007421.t003){#pntd.0007421.t003g}

  Country/fish                                                                                                                                                      Data source                  Date range   Locations   Sample type   GBS No.   CC283 No. (%)
  ----------------------------------------------------------------------------------------------------------------------------------------------------------------- ---------------------------- ------------ ----------- ------------- --------- -----------------------------------------------------
  Singapore; Asian bighead carp, red tilapia, black tilapia, giant snakehead, common snakehead, grass carp and silver carp [^a^](#t003fn001){ref-type="table-fn"}   \[[@pntd.0007421.ref010]\]   2015         Ports       I, S          27        6 (22)
  Markets                                                                                                                                                           I, S                         30           11 (36)                             
  Peninsular Malaysia; red hybrid tilapia                                                                                                                           New data                     2007--2008   9 farms     I             28        28 (100)
  Vietnam, tilapia                                                                                                                                                  \[[@pntd.0007421.ref009]\]   2006         ND          ND            ND        1 [^b^](#t003fn002){ref-type="table-fn"} (ND)
  Vietnam, Can Tho and An Giang; tilapia                                                                                                                            New data                     2016         5 farms     I             34        34 [^**c**^](#t003fn003){ref-type="table-fn"} (100)
  Thailand; tilapia                                                                                                                                                 \[[@pntd.0007421.ref009]\]   2000--2010   ND          I             8         1 (12·5)
  Thailand; tilapia                                                                                                                                                 \[[@pntd.0007421.ref025]\]   2012--2014   ND          ND            ND        1 (ND)

^a^ A glossary of fish scientific names is in [S2 Table](#pntd.0007421.s002){ref-type="supplementary-material"}.

^b^ This is ST491, a single locus variant of ST283.

^C^ One of these 34 is ST1311, a double locus variant of ST283 (See [S3 Table](#pntd.0007421.s003){ref-type="supplementary-material"}). Abbreviations: I = invasive (brain, organs, muscle), S = superficial, ND = not determined.

GBS CC283 outside Asia {#sec027}
----------------------

Only four CC283, of human origin, were found: two from France, one from Britain, and one from the Netherlands ([Fig 1](#pntd.0007421.g001){ref-type="fig"}, [Table 4](#pntd.0007421.t004){ref-type="table"}). All four were ST283 and none had associated epidemiological data. CC283 was otherwise absent from studies of over 4,000 human GBS from nine countries, of over 1,200 bovine GBS from five continents, and from the aquaculture literature from 6 continents, as of December 2017 ([Fig 1](#pntd.0007421.g001){ref-type="fig"} and [S4 Table](#pntd.0007421.s004){ref-type="supplementary-material"}) \[[@pntd.0007421.ref009], [@pntd.0007421.ref022], [@pntd.0007421.ref031]--[@pntd.0007421.ref054]\]. Three human ST751, although an SLV of ST283 and ST10, were distinguished from ST283 by over 2,000 SNPs, based on phylogenetic analysis, and clustered outside of CC283, in CC10 ([Fig 2](#pntd.0007421.g002){ref-type="fig"}), so ST751 is excluded from tables and figures referring to CC283. Other SLVs of ST283 listed in the MLST database are ST690, a human isolate, and ST160, of unknown origin, but we could not find WGS data for them. Three GBS serotype III-4 were not included as ST data were not available, although ST283 is the only published ST within serotype III-4: they were reported, without epidemiological data, amongst a collection from Australia and New Zealand in the original description of this subtype \[[@pntd.0007421.ref055]\].

10.1371/journal.pntd.0007421.t004

###### Human group B *Streptococcus* (GBS) clonal complex (CC) 283 outside Asia.

![](pntd.0007421.t004){#pntd.0007421.t004g}

  ST, place/country, isolate No.   Data source                                          Date range   Sample type   GBS No.   ST283 No.   ST283%
  -------------------------------- ---------------------------------------------------- ------------ ------------- --------- ----------- --------
  ST283, Britain; ERR1742070       GenBank [^**a**^](#t004fn001){ref-type="table-fn"}   2014         Invasive      1017      1           \<0·1
  ST283, France; S80 & S81         \[[@pntd.0007421.ref008]\]                           2002--2007   Invasive      119       2           2
  ST283, Netherlands; ERR1659855   GenBank                                              ND           ND            ND        1           ND

^**a**^ Originally Public Health England (V. Chalker, personal communication). Abbreviations: ST = sequence type, ND = not determined.

Genomic analysis {#sec028}
----------------

A previous report, using Bayesian analysis of a smaller subset of mostly Singaporean isolates, estimated the time of emergence of CC283 as 1994 (95% highest posterior density (HPD) 1991--1997) \[[@pntd.0007421.ref001]\]. Analysis of our expanded set confirms that CC283 forms a monophyletic clade that appears to have arisen from within CC10 ([Fig 2A](#pntd.0007421.g002){ref-type="fig"}). Bayesian analysis of our current set of CC283 and closely related CC10 strains suggests a slightly earlier predicted emergence date of 1985 (95% HPD 1980--1990). Given that the original Bayesian analysis was done on a data set that consisted largely of strains isolated from Singapore during 2015, the shift in the predicted emergence date is perhaps not surprising. However, our current data set was also not systematically collected, either by time, geography, or host, and these results should be interpreted in this context. Additional systematic sampling may provide further insights into emergence dates as well as the likely geographical and host origin of CC283. We found very little recombination in CC283 isolates, 62 kb (3.0%) of the total chromosome, and little variation in genome content, as indicated in the phylogenetic tree ([Fig 2](#pntd.0007421.g002){ref-type="fig"}) and single nucleotide polymorphism (SNP) distances between CC283 pairs ([S5 Table](#pntd.0007421.s005){ref-type="supplementary-material"}).

Discussion {#sec029}
==========

GBS CC283 has been widespread in SE Asia for over two decades in humans, and at least a decade in aquatic animals, but is rare outside SE Asia. As ST283 was found in the first year of every available collection of invasive GBS from SE Asia (Lao PDR, Thailand, Vietnam and Singapore), it may have been present even earlier. The estimated date of emergence, of 1985, predates the first known human cases and is roughly contemporaneous with the start of the expansion of aquaculture in SE Asia, based on reports from the Food and Agricultural Organisation of the United Nations \[[@pntd.0007421.ref056]\]. The 2015 human outbreak of GBS in Singapore was fish borne; although our current study did not address transmission or consumption patterns, human cases in other parts of Southeast Asia could potentially also be fish-borne, especially since consumption of undercooked aquaculture foods is common in Asia, as shown by the high rate of trematode infections \[[@pntd.0007421.ref057], [@pntd.0007421.ref058]\]. It is not known to what extent healthy fish may carry ST283, or what the infectious dose for humans may be, but ST283 was isolated from healthy looking fish in Malaysian farms and in Singapore ports and markets, so fish sold for human consumption could be the source of exposure even if visibly diseased fish were excluded from sale. The high proportions of ST283 in invasive human GBS collections suggest that if human ST283 is acquired from aquaculture, or another undetermined food source, then invasive GBS is primarily a foodborne infection in Thailand and Lao PDR, and largely foodborne in Vietnam. The invasive human data from these three countries are all from areas bordering the Mekong River, separated by up to 1,500 km. If the ST283 data is not representative of each whole country, it may represent an epidemiology peculiar to the ecosystem surrounding the river.

Aquaculture and one health {#sec030}
--------------------------

The phylogenetic tree shows that where human and fish GBS were collected concurrently, they are intermingled, whereas separate clustering of isolates of human and piscine origin tends to reflect collections being from different countries or time periods. We found numerous examples of pairs of human and piscine GBS distinguished by few to zero SNPs. Of note, while strict use of SNP distances is not fully reliable for determination of transmission, many of the examples have SNP distances smaller than 21 and are monophyletic in the ST283 tree. These criteria both support a possible "match", as recently defined \[[@pntd.0007421.ref059]\], although any SNP number cutoff is arbitrary, and the fewer the SNPs, the more likely is the linkage. GBS from the Singapore outbreak cluster with GBS from multiple countries, suggesting multiple sources contributed to the outbreak; although this may not seem surprising, given that Singapore imports fish from multiple countries, it leaves an unanswered question as to why imports from across the region suddenly, and simultaneously, were associated with increased GBS incidence in 2015. Of note, the 2015/2016 El Niño broke warming records in the central Pacific \[[@pntd.0007421.ref060]\], and higher temperatures are associated with increased GBS load \[[@pntd.0007421.ref061]\] and outbreaks of streptococcosis \[[@pntd.0007421.ref062]\] in fish, and increased human cases \[[@pntd.0007421.ref024]\].

GBS is known to affect fish species other than tilapia, including farmed and wild freshwater and marine fishes \[[@pntd.0007421.ref063], [@pntd.0007421.ref064]\]. Our detection of CC283 in fish farms was from tilapia, which were preferentially sampled because they are the most commercially important species affected by GBS in SE Asia, but the involvement of other species remains uncertain. Although the outbreak in Singapore was linked to several freshwater fish species, it is unclear whether the fish left the farms as carriers, or were contaminated during transport and handling. Cross contamination and post-contamination bacterial amplification at ambient temperature might explain the data, from Singapore, that showed that ST283 was recovered from 1% of fish taken from ports, but 28% of fish from markets \[[@pntd.0007421.ref010]\].

Fish account for up to 37% of protein consumed in SE Asia \[[@pntd.0007421.ref065]\], and Lao PDR and Thailand are among the top ten global tilapia producers \[[@pntd.0007421.ref056]\]. Tilapia were introduced to SE Asia in the 1940s with repeated importations of fish fry: if these fry were contaminated at source, this could explain the finding of ST283 across the region. Alternatively, ST283 may have evolved regionally, from other animals, and acquired the ability to infect fishes through lateral gene transfer of virulence elements from other piscine GBS, followed by onward transmission \[[@pntd.0007421.ref066]\]. Interestingly, GBS transmission from a tilapia hatchery into a new farm was reported in Malaysia where, in an effort to prevent economic losses due to streptococcosis, farms were advised to source fry from disease-free hatcheries \[[@pntd.0007421.ref067]\]; GBS from these outbreaks in 2007/2008 were included in our study and are ST283. However, a survey of invasive GBS isolated from 13 adults in Malaysia in 2010 did not find Serotype III GBS \[[@pntd.0007421.ref068]\]; this was a small sample and the predominant Malay population do not habitually eat raw fish, perhaps explaining this lack of human ST283.

There are a limited number of GBS clades that cause streptococcosis in fish, and ST283 and its variants are the only known serotype III GBS that naturally affect tilapia \[[@pntd.0007421.ref009]\]. This suggests that ST283 caused streptococcosis in multiple sites across Thailand between 2003/2011, when 12% to 56% of GBS isolated from diseased tilapia were serotype III \[[@pntd.0007421.ref069]--[@pntd.0007421.ref071]\]; furthermore, some were serotype III-4, and the only described example of serotype III-4 is ST283 \[[@pntd.0007421.ref007]\]. Interestingly, untyped GBS were reported as an emerging cause of septic arthritis, in humans, in Thailand between 1990 and 2010 \[[@pntd.0007421.ref072]\].

The three predominant clades of GBS in fish are associated with different serotypes: ST7 with serotype Ia, CC552 (including ST260) with serotype Ib, and ST283 with serotype III \[[@pntd.0007421.ref009]\]. Commercially available vaccines cover serotype Ib only, or serotypes Ia and III, but without cross-protection between serotypes. Strain confirmation is therefore recommended prior to use of vaccination, but diagnostic infrastructure is very limited in SE Asia, so fish farmers rely on antimicrobial treatment rather than on strain typing and vaccination. Antibiotics, including tetracyclines, are commonly used in controlling streptococcosis \[[@pntd.0007421.ref073]\]. Concentrations above maximum limits have been found in fish sold for consumption in Vietnam \[[@pntd.0007421.ref074]\], which may explain the high prevalence of tetracycline resistance genes found in ST283 from Vietnam. However, although similar antimicrobial use is reported in Thailand \[[@pntd.0007421.ref075]\], tetracycline resistance amongst human ST283 from Thailand disappeared after 2012: this loss was also seen in human ST283 from Lao PDR and Singapore, possibly through three separate resistance gene loss events ([Fig 2](#pntd.0007421.g002){ref-type="fig"}). In contrast, tetracycline resistance was 88% amongst 712 non-ST283 invasive human GBS, isolated in Singapore from 2001 to 2018 (T. Barkham, personal communication).

Potential global dissemination {#sec031}
------------------------------

Whereas fish-associated CC7 and CC552 are geographically widespread, our literature review, up to December 2017, showed that CC283 had only been confirmed from SE Asia, although serotype III GBS had been described as the cause of disease outbreaks that occurred in 2016 in tilapia farms in Brazil \[[@pntd.0007421.ref076]\]. Serotype III GBS isolates from Brazilian tilapia farms have now been confirmed as ST283, while this paper was in revision, and cluster with ST283 from Asia \[[@pntd.0007421.ref025], [@pntd.0007421.ref077]\]. This observation, in combination with import records of live Nile tilapia from Singapore to Brazil in 2014, suggests that ST283 may have been introduced into South America from Southeast Asia, underlining the potential threat of Southeast Asian ST283 to expanding aquaculture worldwide. The global trade in tilapia has previously been reported to account for the dissemination of GBS CC552 (ST260/261), from Israel to Australia, Africa (Ghana), Asia (China) and America (USA) \[[@pntd.0007421.ref078]\] but CC552 does not affect homeothermic species. If ST283 spreads in a similar manner, human disease may also occur beyond SE Asia.

Limitations and further studies {#sec032}
-------------------------------

Our data is limited by the use of existing GBS collections. More extensive and systematic geographical sampling and transmission studies are needed. Future studies might specifically address the prevalence of GBS in fish hatcheries, animal feed, humans (including healthy carriers), and both healthy and diseased aquatic and non-aquatic animals, as well as the evolutionary origin and virulence mechanisms of ST283. Pathogenicity factors that explain the virulence of ST283 in humans are yet to be described, although a recent report identified the *bceR* gene as important for antibiotic resistance, biofilm formation, and lethality in a mouse model of infection \[[@pntd.0007421.ref079]\]. Humans have previously been implicated as the source of GBS that caused a mass die-off of fish \[[@pntd.0007421.ref080]\], so the interplay between aquaculture and human waste, which is commonly discharged into fish ponds and rivers in Asia, in part to recycle nutrients in the food production system, might also be studied.

Summary {#sec033}
-------

GBS ST283 is widespread in Southeast Asia, where it has been causing disease for over 25 years. Human ST283 is almost absent outside Asia but accounts for over 70% of invasive human GBS in collections from Thailand and Lao PDR. As both are significant producers of tilapia, and consumption of raw fish is common in these countries, we hypothesise that their ST283 infections are acquired from fish, as in the Singapore outbreak. ST283 has been detected in healthy and diseased farmed fish in SE Asia and is estimated to have emerged in 1985, corresponding with the beginning of the rapid expansion of the aquaculture industry. ST283 causing tilapia deaths in Brazil since 2016 are thought to have been introduced with tilapia imported from Asia in 2014. Evolutionary origins and routes of transmission within and between host species need further study, as it is possible that ST283 is transmitted from humans to aquaculture, or that there is another common source. Creation of cross-border collaborations in human and animal health are needed to complete the epidemiological picture, may lead to improved human and fish welfare, and may contribute to safer economic development in affected countries.

Disclaimer {#sec034}
----------

The findings and conclusions expressed by authors contributing in this manuscript do not necessarily reflect the official position of the Centers for Disease Control and Prevention or the institutions with which the authors are affiliated.

Supporting information {#sec035}
======================

###### Human group B *Streptococcus* (GBS) sequence type (ST) ST283 in Asia: demographics and data sources.

This table shows details for each human GBS collection, showing proportion of ST283, sample origins, dates and types, and age breakdown, foci of infection and gender.

(DOCX)

###### 

Click here for additional data file.

###### Glossary of scientific names of fish.

(DOCX)

###### 

Click here for additional data file.

###### Details of known single locus variants, and one double locus variant, of group B *Streptococcus* sequence type (ST) 283.

(DOCX)

###### 

Click here for additional data file.

###### Global studies reporting group B streptococcus (GBS) Multi Locus Sequence Typing data, in which clonal complex (CC) 283 was not found.

The absence of CC283 in this data, of over 4,000 human and over 1,300 animal GBS, demonstrate how unusual CC283 is. The data also show that Asia is not well represented in this group, with only one study found, from China.

(DOCX)

###### 

Click here for additional data file.

###### Single Nucleotide Polymorphism (SNP) distances between pairs of group B Streptococcus clonal complex (CC) 283, from humans and animals.

Examples were selected to show minimum and maximum SNP differences found between pairs within and between countries and hosts.

(DOCX)

###### 

Click here for additional data file.

The authors are grateful for the support of clinical and laboratory staff involved in supplying bacteria and data. Potential conflicts of interest: TB and SLC are named applicants on a patent for the ST283-specific PCR test used in this study.

[^1]: I have read the journal\'s policy and the authors of this manuscript have the following competing interests: SLC and TB are named applicants on a patent for the ST83-specific PCR test used in this study.
