Streptococcus suis has been well-recognized as a zoonotic pathogen worldwide, and the diversity and unpredictable adaptive potential of sporadic human strains represent a great risk to the public health. In this study, S. suis LSM178, isolated from a patient in contact with pigs and raw pork, was assessed as a hyper-virulent strain and interpreted for the virulence based on its genetic information. The strain was more invasive for Caco-2 cells than two other S. suis strains, SC19 and P1/7. Sequence analysis designated LSM178 with serotype 2 and a novel sequence type 1005. Phylogenetic analysis showed that LSM178 clustered with highly virulent strains including all human strains and epidemic strains. Compared with other strains, these S. suis have the most and the same virulent factors and a type I-89 K pathogenicity island. Further, groups of genes were identified to distinguish these highly virulent strains from other generally virulent strains, emphasizing the key roles of genes modeling transcription, cell barrier, replication, recombination and repair on virulence regulation. Additionally, LSM178 contains a novel prophage conducive potentially to pathogenicity.
Streptococcus suis (S. suis) is one of the most important swine pathogens leading to severe economic losses to the porcine industry worldwide. However, S. suis has emerged as a zoonotic agent, causing fever, septicemia, meningitis, arthritis and a variety of other symptoms in humans1. Since the first case of human S. suis infection reported in 19682, it has spread in more than 30 countries and regions, particularly the southeast Asian countries where the pathogen represents a significant public health concern3,4. Two serious outbreaks have occurred in China (Sichuan in 2005 and Jiangsu in 1998)5 and four in Thailand (Phayao in 2007, Chiang Mai and Lamphoon in 2008, Phetchabul in 2010 and Uttaradit in 2019)6. A recent study reported that, upon testing of raw pork and edible pig organs collected from 88 sales locations in central Thailand, the positive rate of S. suis was as high as 85.23% and the positive rate of serotype 2 was 17.05%7. In addition, 70.4% of S. suis isolates of serotypes 2 and 14 from slaughtered pigs revealed sequence types and pulsotypes identical to human isolates in Thailand, suggesting transmission from pigs to humans8.
Typing of S. suis strains is epidemiologically important to control the infection. The most commonly used method, serotyping, is not only used for identification and diagnosis of clinical of S. suis isolates, but has also been suggested to be indicative of pathogenicity6. Of the 35 serotypes (types 1–34 and 1/2) originally identified according to the antigenicity of capsular polysaccharide (CPS), six S. suis-like strains (serotypes 20, 22, 26, 32, 33, and 34) have been taxonomically removed from the S. suis species based on phylogenetic and sequence analyses9. Additionally, new variants with serotype Chz and novel cps loci were recently investigated, although their relationship to virulence potential remains unclear10,11,12,13. The prevalence of S. suis serotypes in countries and regions is different. For instance, isolates associated with pig disease were predominantly identified as serotypes 2 and 9 in Europe14, and serotypes 2 and 3 in North America15. However, serotype 2 is considered to be the most toxic and prevalent serotype causing both pig and human infection worldwide16,17, although other serotypes such as serotype 9 and 14 are of increasing concern18.
Besides serotyping, genetic classification by multiple sequence locus typing (MLST)19 has become increasingly important because of its higher resolution for investigating strain evolution and delineating the relationship between subtype and pathogenicity. For instance, while serotype 2 ST1 strains present high zoonotic potential worldwide, ST7 from serotypes 2 and 14 is common in China6,20. And, in Thailand for human infections, ST104 are almost exclusively predominant in serotype 2 and the majority of serotype 14 isolates are ST1056. So far as to April 8, 2021, 2,808 STs have been recorded in the S. suis MLST database, showing that S. suis is constantly evolving. Though only a small number of STs have been found to be mainly responsible for human infections, the increasing diversity brings new risks and challenges, such as ST658 isolated in China. The strain originates from the ST1 complex and showed more mortality rates (90%) than the international reference virulent strain P1/7 (ST1) (70%) in a murine model21,22.
There is no doubt that the repertoire of S. suis virulence determinants play a role in human infection, since S. suis was suggested to be a cause of community-acquired disease23. It is difficult to make any clear distinction about the virulence factors belonging exclusively to pigs versus humans. In a recent study, no defined genomic differences between human strains and pig strains were suggested, although human disease isolates are limited to a single virulent population whose origin nevertheless coincided with the first intensification of pig production24. Another study found that virulent strains could not be identified only by the presence of proposed virulence factors, making the definition of virulence factors ambiguous25. Virulence factors play key roles in many aspects. One reason why serotype 2 strains show high zoonotic potential and virulence in humans is their better adherence to human intestinal epithelial cells6,26,27. Deletion of factors involved in adhesion to host cells greatly attenuated virulence28,29,30,31,32. The ability to escape immune clearance is also necessary for strain survival, dissemination and pathogenesis. For example, CPS has been shown to enhance bacterial resistance against killing by host phagocytes33.
In this study, a human S. suis LSM178 with serotype 2 and a novel ST1005, which causes fever, nausea and general malaise, was isolated and comprehensively assessed as a hyper-virulent strain based on toxicity tests, a zebrafish model of infection and Caco-2 cell toxicity and adhesion/invasion assays and genomic analysis. With genome sequencing, genetic features including virulence factors, pathogenicity islands (PAIs), prophages and core virulence genes were characterized to understand the pathogenic potential of the strain. The analysis suggested that genome plasticity contributes to virulence evolution of S. suis, and specific elements predict the virulence change and adaptation of S. suis to humans.
The newly isolated strain LSM178 was compared for virulence in a zebrafish intraperitoneal infection model with two other well characterised S. suis isolates, SC19 and P1/7. Before 10 h, all three infection groups (LSM178, SC19 and P1/7) of infected zebrafish showed no signs of disease except for a slight decrease in swimming ability. After that, pathological changes appeared in challenged zebrafish such as systemic hemorrhage or abdominal hemorrhage, abdominal swelling and dyspnea. The survival rate did not show significant difference (P = 0.5287) and stabilized at 10% for three groups (Fig. 1A). However, the 50% and 90% lethal time by LSM178 (24–28 h and 34–38 h) appeared always shorter than that by SC19 (26–32 h and 38–48 h) and P1/7 (30–36 h and 36–40 h) in each of 3 independent repeated tests (Fig. 1A and Supplementary Fig. 1). The main symptoms of the death caused by LSM178 infection are serious abdominal swelling, blood spots and ecchymosis (Fig. 1B), which were similar to those caused by SC19 and P1/7.
Human whole blood resistance assay, cytotoxic assay and cell adhesion and invasion test
In a human blood resistance assay, LSM178 survived similarly to P1/7 (P = 0.7948), but seemed grow better than SC19 (P = 0.2905) (Fig. 2A). Although there was no significance in cytotoxity to Caco-2 among the three strains, LSM178 appeared more cytotoxic than SC19 (P = 0.3244) and P1/7 (P = 0.1492) and was almost twice as toxic as P1/7 strain (Fig. 2B). Further, LSM178 showed significantly more invasion of Caco-2 cells than P1/7 (P = 0.0090), but similar to SC19 (P = 0.8454) (Fig. 2C). Both LSM178 (P = 0.0099) and P1/7 (P = 0.0079) had significantly lower adhesion to Caco-2 cells than SC19 (Fig. 2D).
The LSM178 genome (accession number, CP047248) consists of a single circular chromosome of 2,115,437 base pairs (bp) with average GC content of 41.19% and 2,065 ORFs (Fig. 3). Sequence analysis designated LSM178 with serotype 2 and a novel ST (1005). With goeBURST analysis, ST1005 was shown to be an individual ST (Supplementary Fig. 2).
Analysis of antimicrobial susceptibility profiles
The scanning results showed that LSM178 possesses 10 potential specific antibiotic-resistant genes. Only two of these corresponded to a resistance phenotype in antimicrobial susceptibility tests, a tetracycline resistant gene (tetM, 04,995) conferring the tetracycline resistance and a MSL methylase (04,940) encoding gene resulting in both erythromycin and clindamycin resistance (Table 1). Other potential antibiotic-resistance genes include those encoding five penicillin-binding proteins (PBPs, PBP1b (00,760), PBP1a (02,025), PBP1a (09,665), PBP2b (03,000) and PBP2x (08,260)), an Aminoglycoside 6-adenylyltansferase (ANT(6)-Ia, 04,805) and two primary targets of quinolone (DNA gyrase (GyrA, 04,275) and topoisomerase IV (ParC, 05,695)). These are ineffective since LSM178 was sensitive to penicillin, streptomycin and levofloxacin (Table 1).
Phylogenetic tree analysis
Using the genome sequence of LSM178 and 52 published S. suis complete genomes, a phylogenetic tree was generated with the 51,520 core-genome SNPs from the nonrecombinant regions (Fig. 4). Within the 19 strains clustered with LSM178 (LSM178 clade), 7 out of 12 pig isolates and 6 out of 7 human strains were from China. Among the 53 strains, most ST1 (6/9) and ST7 (9/12) were clustered in LSM178 clade. It was interesting that in LSM178 clade, the STs of human isolates (GZ1, LSM178, LSM102, 05ZYH33, SC84, BM407 and 98HAH33 with ST1, ST1005, ST658, ST945, ST7, ST1 and ST890 respectively) are more diverse than those of pig strains (just including ST1 and ST7).
Analysis of prophages and toxin-antitoxin (TA)
Two prophages, pha17801 (positions 1,417,367–1,458,272) and pha17802 (positions 1,871,779–1,926,900), were found in LSM178. The pha17801, with 50 ORFs and GC content of 39.36%, contains two integrases (positions 1,437,316 − 1,437,858 / 1,446,215–1,447,357) and is highly similar to the partial genome of Streptococcus phage 20,617 (RefSeq accession NC_023503). The prophage included two overlapping GIs (positions 1,421,901–1,447,702 and positions 1,437,316–1,447,702). These two GIs were shared only by 18 S. suis (LSM178, 05ZYH33, 98HAH33, A7, BM407, CS100322, GZ1, JS14, LSM102, P1/7, S10, SC070731, SC19, SC84, SS12, SS2-1, T15 and ZY05719) (Fig. 5A). Excluding avirulent T1534, these strains almost constituted the LSM178 clade (Fig. 4).
In pha17802, the overlapped GI (positions 1,866,902–1,904,229) was characteristic to LSM178 in comparison with other 52 S. suis (Fig. 5A). The left part (positions 1,871,862–1,892,563) of the GI is only highly similar to (identity of 90%) a region of S. suis 6407, whereas the right part (positions 1,892,563–1,903,867) together with sequence of positions 903,868–1,906,268 is not homologous with any S. suis but with Streptococcus agalactiae Sag37 (identity of 91%) (Fig. 5B). The prophage comprises 55,121 bp with an average GC content of 41.87%, containing a total of 68 ORFs including prophage core component genes encoding lysin, tail, head, recombinase, capsid, portal, integrase, portal, and Cro/Cl-type repressor. The absence of excisionase confirmed the non-plaques on several strains containing LSM178 (data not shown). The short directly repeated sequences (5′GGTTTCAATTTTA3′) located the prophage between 09,360 (adenylate kinase) and 09,655 (preprotein translocase subunit SecE). The sequence of positions 1,866,902–1,906,268 contained 23 function-known genes along with 29 encoding hypothetical proteins and the unique sequence (positions 1,892,563–1,906,268) involved 29 genes where 21 encode unknown proteins. The TA of HicA (09,375)-HicB (09,370) is the only candidate virulence elements in this phage. Other TAs were listed in Supplementary Table 1.
Comparative analysis of virulence factors
Examining the 96 virulence factors among 53 S. suis genomes, LSM178 is one of the 8 strains (LSM178, ZY05719, 98HAH33, SC84, 05ZYH33, SS2-1, SC19 and LSM102) containing the most (95) virulence factors (Supplementary Table 2). The 95 factors contain 15 linked to immune evasion or systemic infection, including SalKR, NisKR, Epf, Fhb, IgA1, IdeS. suis, MRP, Sly, Nudp, SsnA, EndA, ScpA and SsadS. Excluding GZ1 and BM407 with 87 and 93 virulence factors respectively, these 8 strains are all human and/or epidemic strains. The Rgg is the only virulence factor missing in the 8 S. suis and it is found only in 5 other S. suis (D12, 0061, 1081, CZ130302 and HN105). The common characteristics of virulence factors in the 8 highly virulent strains could explain the highly pathogenic phenotype of LSM178.
An 84 kb PAI (84 K-PAI) (positions 926,718–1,010,299) was found to be highly similar to type I-89 K-PAI of SC19 and both of them are also similar to ICESsu05SC260 belonging to ICESa2603 family (Fig. 5C). The direct repeat sequences indicated the location of the 84 K-PAI just downstream of rplL (04,665). The PAI encodes characterized virulence factors such as SalKR, NisKR and several type IV secretion systems (Vir D4, Vir B6, Vir B1 and Vir B4). The presence of integrase, transposase, excisionase and helicase support activity of transposition and propagation of 84 K-PAI. The major difference of LSM178 84 K-PAI from SC19 89 K-PAI is the absence genes encoding PqqD family protein, asparagine synthetase and 4 ABC transporter units, and addition of genes encoding transposase and MSL methylase as described. Apart from the methylase which may cause epigenetic changes associated with virulence, the encoded proteins are not defined as virulence factors and should not contribute to virulence.
The 89 K-PAI has evolved to be diverse22. The type I-89 K-PAI was first identified in SS2-1 (diseased pig, 1998) and 98HAH33 (human, 1998), close relatives of SC84 (human, 2005), ZY05719 (pig, 2005) and SC19 (pig, 2005) (Fig. 5A and Fig. 4). From 2005, the PAI was reported in almost all human strains (SC84, 2005; 05ZYH33, 2005; LSM102, 2014 and LSM178, 2016) apart from GZ1 (2005), but only in two (ZY0571, 2005 and SC19, 2005) rather than the other 10 pig strains. The human strains not from China (861,180/Netherlands/1986 and BM407/Vietnam/2004) do not contain type I-89 K-PAI (Figs. 4 and 5A). The distribution of the type I-89 K-PAI among the strains coincides interestingly with the combination of the most (95) virulent factors (Fig. 5A and Supplementary Table 2). Since LSM178 shares these important virulence markers with epidemic strains and almost all human strains, these strains were referred here as the highly virulent strains (HVS) at least for serotype 2 strains from China, those are LSM178, ZY05719, 98HAH33, SC84, 05ZYH33, SS2-1, SC19 and LSM102.
Genomic comparative analysis
Comparison of virulence factors between strains in LSM178 clade and closely related avirulent T15 identified 87 virulence factors which were shared by T15 and the LSM178 clade strains, except A7, P1/7, S10 and GZ1 (Supplementary Table 2). These four strains had 87 virulence factors, but they share Epf, NadR, RevS and SBP2 instead of Trag, VirB1, VirB4 and VirD4 as in T15. However, there were no virulence factors belonging exclusively to avirulent T15 or all LSM178 clade strains. To globally identify the characteristic virulence factors, the difference was checked between strains in LSM178 clade and closely related avirulent T15. Referring to T15, the numbers of unique genes of each virulent strain were between 318 and 582, and 224 shared genes were found (Fig. 6A). HVSs possess more core genes (332) than other virulent strains (249, called generally virulent strain here) (Fig. 6B and C). Addition of any one generally virulent strain except BM407 (Vietnam) resulted sharp decreased number of core high virulence genes (Fig. 6D). However, little change was observed when any one highly virulent strain was removed. This supports HVSs as a seperate group in serotype 2 strains. The characteristic genes from various categories encode a wide range of functions (Fig. 6E–G). Compared with generally virulent strains, HVSs show predominant increase on genes with functions of transcription, cell barrier, replication, recombination, repair and mobile elements (Fig. 6F and G). GIs of LSM178 occupy 10.4% of the genome (Fig. 6H). While there are 3.7% of core virulence genes in GIs of LSM178, 23.5% of core high virulence genes locate in the GIs (Fig. 6H).
Generally, S. suis infections in humans are restricted to workers in close contact with pigs or swine byproducts. However, in southeast Asia, the bacterium has been reported to affect the general population23. The sporadic human S. suis, particular those with serotype 2, should get attention in epidemiological monitoring because of their unpredictable adaptive potential, as demonstrated by LSM178, a serotype 2 strain with a novel ST1005. LSM178 was more efficiently invasive to Caco-2 cells compared with P1/7 and SC19. However, no significant virulence was detected in cytotoxic activity and zebrafish challenge assays which have been used as a model to evaluate the virulence of S. suis35. Indeed, virulence of S. suis could not be intensively evaluated and compared in other models36. Probably, virulent strains have their own outstanding virulence aspects, which can balance the weak ones and eventually lead to a similar pathogenicity at least under the inoculation. For instance, P1/7 and SC19 were prominent in biofilm formation (Supplementary Fig. 3) and adhesion respectively.
Recent studies suggested that quinolones, beta-lactams, florfenicol and trimethoprim/sulfamethoxazole could still effectively treat clinical S. suis infection37,38,39. However, strains resistant to widely used effective beta-lactams have been increasely reported40,41. LSM178 showed the sensitivity at least to beta-lactams and quinolones, consistent with the features of resistance genes. The PBPs of LSM178 are entirely the same as those of sensitive A7, but habor substitutions throughout the sequence compared with the resistant R61 (Supplementary Fig. 4A–E)40. Several mutations in the quinolone resistance-determining region (QRDR) of both GyrA and ParC could increase the resistance to fluoroquinolone42. However, no amino acid changes were revealed in the QRDRs between LSM178 and four quinolone-sensitive strains (A7, BM407, P1/7 and SC84)40 (Supplementary Fig. 5). Additionally, it is interesting that the P1/7 without ANT(6)-Ia confered resistance against streptomycin, whereas ANT(6)-Ia containing LSM178 could not38. However, ANT(6)-Ia of LSM178 is only 49.5% identical to the functional homologue (UniProtKB-P12055 (AADK_STAAU)) from Staphylococcus aureus43 (Supplementary Fig. 6).
Several factors are deemed to be important for the pathogenesis, such as CPS, Fbps, enolases, dipeptidylpeptidase DppIV and SrtA25, which were all found in LSM178. It has been demonstrated that S. suis serotype 2 virulent strains are able to exacerbate inflammatory activation scavenging bacteria. All the 15 known anti-immunity factors22 exist in LSM178 and would modulate the immune responses improving its survival at the inflammation area. Although it is difficult to define a virulent strain only by proposed virulence genes, strains isolated from healthy animal are found with less virulence factors (less than 87), such as WUSS351, 0061, 1081, HA1003, TL13, 05HAS68, YB51, LS9N, SRD478, 90–1330 and T15. These S. suis lack at least one of the marker virulence genes (sly, mrp, epf and cps2). Notably, four of those strains contain less than 70 virulent factors (0061, 1081, HA1003 and LS9N). On the contrary, the HVSs contain consistently the most (95/96) virulent factors. However, it should be realized that the combination of various virulence factors may cause pathogenicity despite the number of the virulence factors. For instance, the CZ130302 contains 69 virulence factors but shown to be a virulent strain11, and the strains other than T15 isolated from healthy animal are not necessarily the avirulent strains.
It is very clear that intermediately pathogenic strain could evolve to become a highly pathogenic and then epidemic strain4. It should be noted that the evolutionary relationships could be lack of reliably based on the raw sequence due to the extensive genetic recombination in S. suis24. With the nonrecombinant regions, phylogenetic analysis should show the reliable evolutionary relationships1. In LSM178 clade, all 19 strains possessed the serotype 2. However, while 7 human strains have 6 STs, 12 pig strains were just ST1 (4 strains) or ST7 (8 strains). It may imply that while S. suis has been purified to relatively stable genotype in pig, multiple evolutionary directions are in progress to be epidemic during their adaption to human, at least for serotype 2 strains in China.
Type I-89 K-PAI is specific to highly pathogenic S. suis linked to Chinese epidemics22 and could horizontally transfer among strains44. Humans can carry avirulent S. suis without clinical signs3,45. However, the 89 K-PAI of human virulent strain should not be obtained from pig strain owning to the transfer, since all strains with 89 K-PAI presented an extremely short evolutionary distance from each other. It suggested that human pathogenic infection was due to the interspecies transmission of a swine-origin strain. The gradual increase of 89 K-PAIs in human strains and decrease in pig strains may suggest that the 89 K-PAI is one of the markers adapting pig strain to human. On the contrary, the existence of PAI would not be conducive to the adaptation of strains in pigs, that might explain the reason why the pig strains lose it more and more rapidly. Under this consideration, the combination of 95 virulence factors may also be regarded as indicative one of the potential for the adaptation to human.
In published S. suis genomes, many prophages were held as remnants46,47. A few intact prophages was described and one of them was reported to be induced to lyse S. suis48. The two prophages in LSM178 could not be induced to form plaques on several S. suis strains. The reason might be that improper induction method was adopted or that the phages produced are defective on infection even for the almost intact pha1780248. Alternatively, the lysis spectrum of the phages is very narrow, or no phages were induced at all, which could be supported by the facts that pha17801 contains only integrase and pha17802 does not include excisionase. Genomic comparative analysis suggested that the pha17802 may integrate two elements horizontally transferred from S. suis 6407 and S. agalactiae Sag37 respectively. The unique gene fragment from S. agalactiae Sag37 has not been found in any other S. suis genome published so far. The empty target site might accommodate potentially unpredictable integration of other mobile genetic elements. Prophages make up a platform for the dissemination of virulence determinants between intra- and inter-species, contributing to the evolution of pathogenic bacteria30,49. For pha17802, except the hicAB located in the arm homologous to 6407, no known genes contributing to virulence were identified. There was no evidence demonstrating the contribution of the two transferring elements to virulence. However, the rarity of pha17802 suggested that it might provide an advantage under certain circumstances. At least, the lysogenic state would increase the survival in the environment by resisting to infection of similar viruses. Furthermore, it may increase the fitness of the bacteria by modulating host metabolism48.
While the core virulence genes were identified through the comparation between avirulent T15 and virulent strains in LSM178 clade, they should be responsible for the enhanced pathogenicity. Moreover, some genes could be used to make a distinction between HVSs and other general virulent strains. Particularly, it should emphasize the roles of genes with functions of transcription, cell barrier, replication, recombination and repair in virulence enhancement, since their number increased obviously in HVSs core genes. Thus, the importance of the mobilizable elements is beyond all doubt, because it is the carrier responsible for genetic differences50. In fact, GIs make up a considerable part of the genome for LSM178 (10.4%). And, higher proportion of virulence-enhancing genes are dispersed in GIs, for example that GIs of LSM178 are colonized with 3.7% of core general virulence genes but with 23.5% of core high virulence genes. However, there are less GIs in LSM178 (number of 23 and total length of 21, 5087 bp) than T15 (29 and 31, 4245 bp) (Supplementary Fig. 7B). These suggested that fusion of specific PAIs increased the virulence. The non-existence of CRISPRs defending against foreign invading elements should be one of the reasons interpreting the rich GIs51. In addition, assay showed that there are more core avirulence genes (535) (Supplementary Fig. 7A) than core virulence genes (224). Probably, both the gain of virulence genes and loss of avirulence genes contribute to the increased virulence. Other potential difference, such as SNPs and patches of insertion and deletion, deserve also to be explored and should not be ignored. These differences are small but numerous and they would definitely interpret the change of strain virulence through just affecting the genes expression or protein activity52.
Materials and methods
Strains and antimicrobial susceptibility testing
S. suis LSM178 was isolated at 2016 from a patient with clinical symptoms of fever, nausea, and general malaise. The patient had been in contact with pigs and handled raw pork before admission. The S. suis SC19 and S. suis P1/7 were stored in our laboratory. Antimicrobial susceptibility was tested by E-test (AB Biodisk, Sweden) with Streptococcus pneumoniae ATCC49619 as a control. All protocols was approved by committee of State Key Laboratory of Agriculture Microbiolgy and the ethics committee of Huazhong Agricultural University. Experiments were performed in ABSL 3 laboratory, Huazhong Agricultural University.
Human whole blood resistance assay
Blood assays were conducted according to an approval issued by the Medical Ethics Committee of the Huazhong Agriculture University (Wuhan, China). Strain suspension (100 μL, 5.0 × 107 CFU/ml) was transferred into 900 μL of fresh human whole blood followed by incubation at 37℃. Samples were withdrawn every 1 h and diluted to incubate on TSA solid medium (containing 10% fresh FBS) at 37 ℃ overnight to count colony-forming units (CFU). Growth index (%) = (CFUat a certain time point—CFUoriginal inoculum) / CFUoriginal inoculum × 100%.
Zebrafish were fed as previously described53. Inoculum was collected at the end of the logarithm period, cleaned twice with phosphate buffered saline (PBS), and adjusted to the appropriate dose (2.5 × 109 CFU/ml)35. Adult zebrafish were infected by intraperitoneal inoculation with 20 μL of bacterial solution per tail. Each group contained 10 zebrafish. The symptoms of zebrafish were recorded every 2 h. Euthanasia of zebrafish was conducted at 60 h after challenge using Tris-buffered tricaine at a concentration of 320 μg/ml.
Biofilm formation assay
Strains (20 μL, 5.0 × 107 CFU/ml) were inoculated into 2 mL TSB medium (containing 10% fresh FBS) and cultured in a 24-well cell plate at 37℃ for 3 days, and the un-inoculated medium was used as the control. After that, the strains were washed twice with sterilized PBS, and fixed with 500 μL methanol for 30 min to attach the strains to the wall. Then, methanol was removed and the plate was air-dried at room temperature. Next, 500 μL of 0.1% crystal violet dye solution was added in and removed out until 30 min later. After drying at 56℃, 500 μL of 33% acetic acid solution was added and placed on a shaker for 30 min to release the crystal violet bound to the biofilm. The released solution (200 μL) from each well was measured at a wavelength of 600 nm.
Strains at log phase was used in the expriments. For cytotoxic assay, Caco-2 cells (1 × 104)54 in 96 well plate was used to detect cytotoxicity of the strains (2 × 105 CFU) with lactate dehydrogenase kit (Beyotime, Beijing, China). The percentage of cytotoxicity was calculated referring to the protocol of the kit: cytotoxicity (%) = (LDH release from infected cells –spontaneous release of LDH from uninfected cells) / (maximum LDH release from cell lysate–spontaneous release of LDH from uninfected cells) × 100%.
For cell adhesion and invasion, single layer Caco-2 cells in the 24 well culture plate were inoculated with 500 μL bacterial suspension (1 × 106 CFU) for 3 h. After washing to remove unadhesive strain, cells were then treated by trypsin digestion for 2 h. In invasion, extracellular bacteria were treated with gentamicin (100 μg/ml) and penicillin G (5 μg/ml) before trypsin treatment. The digested cells were lysed using 1% saponin and the lysis was inoculated on THB plate. The rate of adhesion (Ra) and invasion (Ri) was expressed as (CFU determined from plate / CFU original inoculum) × 100%. The relative invasion rate was expressed as Ri strain / Ri P1/7 × 100%.
LSM178 culture in exponential growth phase was induced by mitomycin C (500 ng/ml) (Sigma, St. Louis, USA) for 5–15 min. The culture (100 μL) was mixed with 3 ml TSA (45 ℃) containing 10% fresh FBS to prepare sandwich plaque assay at 37 ℃. Plaque formation was observed after 12 h.
The genomic DNA was extracted using a DNA extraction kit (TaKaRa DNAiso, TaKaRa Biotechnology Co., Ltd., Dalian, China). The genome of LMS178 was sequenced using combined plaforms of Illumina Miniseq and PacBio sequel. Illumina Miniseq generated 4,234,826 reads producing a total of 926,438,020 bp with Q30 of 76.81%. PacBio sequel generated 234,808 reads and a total sequence length of 1,538,801,765 bp with N50 of 9,236 bp. The data from Illumina Miniseq and PacBio sequel were assembled by A5-Miseq v2015052255 and CANU56 respectively. After correction of the results using Pilon57, the complete circular genome was constructed.
Local virulence factor database of S. suis was established (Supplementary Table 2) and used to detect the potential virulence genes in genome. Of the 96 virulence factors, 84 are from the 2016 study by Willemse et al58. Another 12 putative virulence factors were identifed by systematic search using the term ‘Streptococcus suis virulence factor’ on July 27th 2017 in NCBI PubMed. TAs were predicted by Rasta59 and TAfinder60. Antibiotic resistance genes were predicted with the comprehensive antibiotic resistance database (CARD) with default settings61. The ST was determined using the MLST typing scheme (https://pubmlst.org/ssuis/). Using BLAST (evalue ≤ 1e-10) 62, serotyping strategy was executed based on the levels (≥ 80%) of both identity and coverage between WZY amino acid sequence of the LSM178 and 33 standard strains with known serotypes63. The sequence of CpsK was used to discriminate the serotype 2 (W161) from 1/2 (C161) for all strains with serotype 2 or 1/264,65. ST complexes were analysed by goeBURST61 program (http://goeburst.phyloviz.net). Prediction of gene islands (GI) and prophages were performed using IslandViewer 466 and PHAST67 respectively. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) were predicted by CRISPR recognition tool (CRT)68. Open reading frames (ORFs) were predicted with Glimmer3.069. tRNA and rRNA were predicted with the Aragorn70 and RNAmmer71 embedded in Prokka72. The proteins in prophage and PAI were integrally annotated by databases of NR73, eggNOG74, KEGG75, Swiss-Prot76 and GO77. Genome mapping with information was generated by CGView78. To identify the unique regions in the genome of LSM178, the whole-genome sequence was used as a reference to compare with that of the other 52 strains using BLAST Ring Image Generator (BRIG)79. The 84 K-PAI in LSM178, 89 K-PAI in SC 19 and ICESsu05SC260 were compared to determine their similarity using BLAST embedded in Easyfig80. Similarly, the characteristics of the prophages in LSM178 were evaluated. The differential genes between two strains were extracted with Roary81 to create Venn diagrams and their COG functions were classified. Circos82 was conducted to map the distribution of GIs and interesting genes in genomes of LSM178 and T15. The.
Complete genomic sequences of 52 S. suis sequences were downloaded from NCBI database (https://www.ncbi.nlm.nih.gov/genome/genomes/199) (Supplementary Table 3). The chromosomal sequences were aligned using Parsnp83 generating core-genome single-nucleotide polymorphisms (SNPs). Phylogenetic tree based on core SNP of 53 complete S. suis genome (including LSM178) was constructed using maximum-likelihood phylogenetic trees by FastTree embedded in the Parsnp. In sequence alignment, recombinant regions were filtered using Gubbins v2.4.084. The bootstrap value was set at 1000 times. The phylogenetic tree was displayed using the online website iTOL (http://itol.embl.de/).
Independent determinations were performed in dupicate or triplicate and experiment was repeated at least 3 times in each group. The significance was analyzed with unpaired student’s test using GraphPad Prism 5. The bar represent the mean ± standard error of mean (sem). The P < 0.05 and P < 0.01 were represented as * and ** respectively.
This study was carried out in compliance with the ARRIVE guidelines. The study was approved by ethics committee of Huazhong Agricultural University and all experiments were performed in accordance with guidelines of State Key Laboratory of Agriculture Microbiology. The informed consent was obtained from all participants and/or their legal guardians.
Statement on Guidelines for Human
All procedures performed in studies involving human participants were in accordance with the ethical standards of national research committee and with the 1975 Helsinki Declaration (or its later amendments).
Statement on Guidelines for Animal
All procedures performed in studies involving animal were in accordance with the ethical standards of national institutional guidelines on the care and use of animals.
The whole-genome sequence of the LSM178 was deposited in NCBI Genbank (Accession Number CP047248). Accession numbers for raw data from Illumina Miniseq and PacBio sequel are SRR15853891 and SRR15943338 respectively.
Dong, X. et al. The global emergence of a novel Streptococcus suis clade associated with human infections. EMBO Mol. Med. 13, e13810. https://doi.org/10.15252/emmm.202013810 (2021).
Perch, B., Kristjansen, P. & Skadhauge, K. Group R streptococci pathogenic for man. Two cases of meningitis and one fatal case of sepsis. Acta Pathol. Microbiol. Scand. 74, 69–76 (1968).
Gottschalk, M., Xu, J., Calzas, C. & Segura, M. Streptococcus suis: A new emerging or an old neglected zoonotic pathogen?. Future Microbiol. 5, 371–391. https://doi.org/10.2217/fmb.10.2 (2010).
Wertheim, H. F., Nghia, H. D., Taylor, W. & Schultsz, C. Streptococcus suis: An emerging human pathogen. Clin. Infect. Dis. 48, 617–625. https://doi.org/10.1086/596763 (2009).
Ye, C. et al. Spread of Streptococcus suis sequence type 7 China. Emerg. Infect. Dis. 14, 787–791. https://doi.org/10.3201/eid1405.070437 (2008).
Segura, M. et al. Update on Streptococcus suis Research and Prevention in the Era of Antimicrobial Restriction: 4th International Workshop on S. suis. Pathogens https://doi.org/10.3390/pathogens9050374 (2020).
Boonyong, N., Kaewmongkol, S., Khunbutsri, D., Satchasataporn, K. & Meekhanon, N. Contamination of Streptococcus suis in pork and edible pig organs in central Thailand. Vet. World 12, 165–169. https://doi.org/10.14202/vetworld.2019.165-169 (2019).
Kerdsin, A. et al. Genotypic comparison between streptococcus suis isolated from pigs and humans in Thailand. Pathogens 9, 50. https://doi.org/10.3390/pathogens9010050 (2020).
Okura, M. et al. Current taxonomical situation of streptococcus suis. Pathogens https://doi.org/10.3390/pathogens5030045 (2016).
Huang, J. et al. Identification of six novel capsular polysaccharide loci (NCL) from Streptococcus suis multidrug resistant non-typeable strains and the pathogenic characteristic of strains carrying new NCLs. Transbound. Emerg. Dis. 66, 995–1003. https://doi.org/10.1111/tbed.13123 (2019).
Pan, Z. et al. Novel variant serotype of streptococcus suis isolated from piglets with meningitis. Appl. Environ. Microbiol. 81, 976–985. https://doi.org/10.1128/AEM.02962-14 (2015).
Zheng, H. et al. Eight novel capsular polysaccharide synthesis gene loci identified in nontypeable streptococcus suis isolates. Appl. Environ. Microbiol. 81, 4111–4119. https://doi.org/10.1128/AEM.00315-15 (2015).
Qiu, X., Bai, X., Lan, R., Zheng, H. & Xu, J. Novel capsular polysaccharide loci and new diagnostic tools for high-throughput capsular gene typing in streptococcus suis. Appl. Environ. Microbiol. 82, 7102–7112. https://doi.org/10.1128/AEM.02102-16 (2016).
Wisselink, H. J., Smith, H. E., Stockhofe-Zurwieden, N., Peperkamp, K. & Vecht, U. Distribution of capsular types and production of muramidase-released protein (MRP) and extracellular factor (EF) of Streptococcus suis strains isolated from diseased pigs in seven European countries. Vet. Microbiol. 74, 237–248. https://doi.org/10.1016/s0378-1135(00)00188-7 (2000).
Goyette-Desjardins, G., Auger, J. P., Xu, J., Segura, M. & Gottschalk, M. Streptococcus suis, an important pig pathogen and emerging zoonotic agent-an update on the worldwide distribution based on serotyping and sequence typing. Emerg. Microbes. Infect. 3, e45. https://doi.org/10.1038/emi.2014.45 (2014).
Ye, C. et al. Streptococcus suis sequence type 7 outbreak, Sichuan China. Emerg. Infect. Dis. 12, 1203–1208. https://doi.org/10.3201/eid1708.060232 (2006).
Higgins, R., Gottschalk, M., Boudreau, M., Lebrun, A. & Henrichsen, J. Description of six new capsular types (29–34) of Streptococcus suis. J. Vet. Diagn. Invest. 7, 405–406. https://doi.org/10.1177/104063879500700322 (1995).
Kerdsin, A. et al. Emergence of Streptococcus suis serotype 9 infection in humans. J. Microbiol. Immunol. Infect. 50, 545–546. https://doi.org/10.1016/j.jmii.2015.06.011 (2017).
Auger, J. P., Fittipaldi, N., Benoit-Biancamano, M. O., Segura, M. & Gottschalk, M. Virulence studies of different sequence types and geographical origins of streptococcus suis serotype 2 in a mouse model of infection. Pathogens https://doi.org/10.3390/pathogens5030048 (2016).
Yu, H. et al. Human Streptococcus suis outbreak, Sichuan China. Emerg. Infect. Dis. 12, 914–920. https://doi.org/10.3201/eid1206.051194 (2006).
Chen, C. et al. A glimpse of streptococcal toxic shock syndrome from comparative genomics of S. suis 2 Chinese isolates. PLoS ONE 2, e315. https://doi.org/10.1371/journal.pone.0000315 (2007).
Zhou, Y. et al. Predominance of Streptococcus suis ST1 and ST7 in human cases in China, and detection of a novel sequence type, ST658. Virulence 8, 1031–1035. https://doi.org/10.1080/21505594.2016.1243193 (2017).
Hui, A. C. et al. Bacterial meningitis in Hong Kong: 10-years’ experience. Clin. Neurol. Neurosurg. 107, 366–370. https://doi.org/10.1016/j.clineuro.2004.10.006 (2005).
Weinert, L. A. et al. Genomic signatures of human and animal disease in the zoonotic pathogen Streptococcus suis. Nat. Commun. 6, 6740. https://doi.org/10.1038/ncomms7740 (2015).
Berthelot-Herault, F., Gottschalk, M., Morvan, H. & Kobisch, M. Dilemma of virulence of Streptococcus suis: Canadian isolate 89–1591 characterized as a virulent strain using a standardized experimental model in pigs. Can. J. Vet. Res. 69, 236–240 (2005).
del Campo Sepulveda, E. M., Altman, E., Kobisch, M., D’Allaire, S. & Gottschalk, M. Detection of antibodies against Streptococcus suis capsular type 2 using a purified capsular polysaccharide antigen-based indirect ELISA. Vet. Microbiol. 52, 113–125. https://doi.org/10.1016/0378-1135(96)00056-9 (1996).
Benga, L., Goethe, R., Rohde, M. & Valentin-Weigand, P. Non-encapsulated strains reveal novel insights in invasion and survival of Streptococcus suis in epithelial cells. Cell Microbiol. 6, 867–881. https://doi.org/10.1111/j.1462-5822.2004.00409.x (2004).
Zhang, Y. et al. SssP1, a streptococcus suis fimbria-like protein transported by the SecY2/A2 system, contributes to bacterial virulence. Appl. Environ. Microbiol. https://doi.org/10.1128/AEM.01385-18 (2018).
Li, M. et al. The type II histidine triad protein HtpsC is a novel adhesion with the involvement of Streptococcus suis virulence. Virulence 6, 631–641. https://doi.org/10.1080/21505594.2015.1056971 (2015).
Ni, H. et al. Inactivation of the htpsA gene affects capsule development and pathogenicity of Streptococcus suis. Virulence 11, 927–940. https://doi.org/10.1080/21505594.2020.1792080 (2020).
Yao, X. et al. The chromosomal SezAT toxin-antitoxin system promotes the maintenance of the SsPI-1 pathogenicity island in epidemic Streptococcus suis. Mol. Microbiol. 98, 243–257. https://doi.org/10.1111/mmi.13116 (2015).
Liu, F. et al. Binding of fibronectin to SsPepO facilitates the development of streptococcus suis meningitis. J. Infect. Dis. 217, 973–982. https://doi.org/10.1093/infdis/jix523 (2018).
Chabot-Roy, G., Willson, P., Segura, M., Lacouture, S. & Gottschalk, M. Phagocytosis and killing of Streptococcus suis by porcine neutrophils. Microb. Pathog. 41, 21–32. https://doi.org/10.1016/j.micpath.2006.04.001 (2006).
Touil, F., Higgins, R. & Nadeau, M. Isolation of Streptococcus suis from diseased pigs in Canada. Vet. Microbiol. 17, 171–177. https://doi.org/10.1016/0378-1135(88)90008-9 (1988).
Wu, Z., Zhang, W., Lu, Y. & Lu, C. Transcriptome profiling of zebrafish infected with Streptococcus suis. Microb. Pathog. 48, 178–187. https://doi.org/10.1016/j.micpath.2010.02.007 (2010).
Zhao, Y. et al. Role of a type IV-like secretion system of Streptococcus suis 2 in the development of streptococcal toxic shock syndrome. J. Infect. Dis. 204, 274–281. https://doi.org/10.1093/infdis/jir261 (2011).
Huang, J., Shang, K., Kashif, J. & Wang, L. Genetic diversity of Streptococcus suis isolated from three pig farms of China obtained by acquiring antibiotic resistance genes. J. Sci. Food Agric. 95, 1454–1460. https://doi.org/10.1002/jsfa.6841 (2015).
Jones, M. E. et al. Prevalence of gyrA, gyrB, parC, and parE mutations in clinical isolates of Streptococcus pneumoniae with decreased susceptibilities to different fluoroquinolones and originating from Worldwide Surveillance Studies during the 1997–1998 respiratory season. Antimicrob. Agents. Chemother. 44, 462–466. https://doi.org/10.1128/AAC.44.2.462-466.2000 (2000).
Fittipaldi, N., Segura, M., Grenier, D. & Gottschalk, M. Virulence factors involved in the pathogenesis of the infection caused by the swine pathogen and zoonotic agent Streptococcus suis. Future Microbiol. 7, 259–279. https://doi.org/10.2217/fmb.11.149 (2012).
Hu, P. et al. Comparative genomics study of multi-drug-resistance mechanisms in the antibiotic-resistant Streptococcus suis R61 strain. PLoS ONE 6, e24988. https://doi.org/10.1371/journal.pone.0024988 (2011).
Moneret-Vautrin, A. Letter to the editor regarding the article entitled:"idiopathic capillary leak syndrome". Rev. Med. Int. 31, 180. https://doi.org/10.1016/j.revmed.2009.10.431 (2010).
Wu, Z. et al. Probing genomic diversity and evolution of Streptococcus suis serotype 2 by NimbleGen tiling arrays. BMC Genom. 12, 219. https://doi.org/10.1186/1471-2164-12-219 (2011).
Tang, F., Bossers, A., Harders, F., Lu, C. & Smith, H. Comparative genomic analysis of twelve Streptococcus suis (pro)phages. Genomics 101, 336–344. https://doi.org/10.1016/j.ygeno.2013.04.005 (2013).
Ma, Y. L. & Lu, C. P. Isolation and identification of a bacteriophage capable of infecting Streptococcus suis type 2 strains. Vet. Microbiol. 132, 340–347. https://doi.org/10.1016/j.vetmic.2008.05.013 (2008).
Nghia, H. D. et al. Risk factors of Streptococcus suis infection in Vietnam A case-control study. PLoS ONE 6, e17604. https://doi.org/10.1371/journal.pone.0017604 (2011).
Segura, M., Fittipaldi, N., Calzas, C. & Gottschalk, M. Critical streptococcus suis virulence factors: Are they all really critical?. Trends Microbiol. 25, 585–599. https://doi.org/10.1016/j.tim.2017.02.005 (2017).
Waldor, M. K. & Friedman, D. I. Phage regulatory circuits and virulence gene expression. Curr. Opin. Microbiol. 8, 459–465. https://doi.org/10.1016/j.mib.2005.06.001 (2005).
Fujisawa, H. & Minagawa, T. DNA packaging by double stranded DNA phages. Uirusu 36, 185–193. https://doi.org/10.2222/jsv.36.185 (1986).
Wu, Z. et al. Comparative genomic analysis shows that Streptococcus suis meningitis isolate SC070731 contains a unique 105K genomic island. Gene 535, 156–164. https://doi.org/10.1016/j.gene.2013.11.044 (2014).
Frost, L. S., Leplae, R., Summers, A. O. & Toussaint, A. Mobile genetic elements: The agents of open source evolution. Nat. Rev. Microbiol. 3, 722–732. https://doi.org/10.1038/nrmicro1235 (2005).
Brouns, S. J. et al. Small CRISPR RNAs guide antiviral defense in prokaryotes. Science 321, 960–964. https://doi.org/10.1126/science.1159689 (2008).
Pisarenko, S. V. et al. Global evolution and phylogeography of Brucella melitensis strains. BMC Genom. 19, 353. https://doi.org/10.1186/s12864-018-4762-2 (2018).
Neely, M. N., Pfeifer, J. D. & Caparon, M. Streptococcus-zebrafish model of bacterial pathogenesis. Infect. Immun. 70, 3904–3914. https://doi.org/10.1128/IAI.70.7.3904-3914.2002 (2002).
Ferrando, M. L. et al. Host-pathogen interaction at the intestinal mucosa correlates with zoonotic potential of streptococcus suis. J. Infect. Dis. 212, 95–105. https://doi.org/10.1093/infdis/jiu813 (2015).
Coil, D., Jospin, G. & Darling, A. E. A5-miseq: An updated pipeline to assemble microbial genomes from Illumina MiSeq data. Bioinformatics 31, 587–589. https://doi.org/10.1093/bioinformatics/btu661 (2015).
Koren, S. et al. Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. https://doi.org/10.1101/gr.215087.116 (2017).
Walker, B. J. et al. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963. https://doi.org/10.1371/journal.pone.0112963 (2014).
Willemse, N. et al. An emerging zoonotic clone in the Netherlands provides clues to virulence and zoonotic potential of Streptococcus suis. Sci. Rep. 6, 28984. https://doi.org/10.1038/srep28984 (2016).
Sevin, E. W. & Barloy-Hubler, F. RASTA-Bacteria: A web-based tool for identifying toxin-antitoxin loci in prokaryotes. Genome Biol. 8, R155. https://doi.org/10.1186/gb-2007-8-8-r155 (2007).
Shao, Y. et al. TADB: A web-based resource for Type 2 toxin-antitoxin loci in bacteria and archaea. Nucleic. Acids Res. 39, D606-611. https://doi.org/10.1093/nar/gkq908 (2011).
Alcock, B. P. et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 48, D517–D525. https://doi.org/10.1093/nar/gkz935 (2020).
McGinnis, S. & Madden, T. L. BLAST: At the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 32, W20-25. https://doi.org/10.1093/nar/gkh435 (2004).
Wang, J. P., Du, P. H., Bai, X. & Zheng, H. Relationship between wzy genes and serotype of Streptococcus suis. Disease Surveillance https://doi.org/10.3784/j.issn.1003-9961.2018.12.006 (2018).
Roy, D. et al. A single amino acid polymorphism in the glycosyltransferase CpsK defines four Streptococcus suis serotypes. Sci. Rep. 7, 4066. https://doi.org/10.1038/s41598-017-04403-3 (2017).
Athey, T. B. et al. Determining Streptococcus suis serotype from short-read whole-genome sequencing data. BMC Microbiol. 16, 162. https://doi.org/10.1186/s12866-016-0782-8 (2016).
Bertelli, C. et al. IslandViewer 4: Expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 45, W30–W35. https://doi.org/10.1093/nar/gkx343 (2017).
Zhou, Y., Liang, Y., Lynch, K. H., Dennis, J. J. & Wishart, D. S. PHAST: A fast phage search tool. Nucleic Acids Res 39, W347-352. https://doi.org/10.1093/nar/gkr485 (2011).
Bland, C. et al. CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinform. 8, 209. https://doi.org/10.1186/1471-2105-8-209 (2007).
Delcher, A. L., Bratke, K. A., Powers, E. C. & Salzberg, S. L. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 23, 673–679. https://doi.org/10.1093/bioinformatics/btm009 (2007).
Laslett, D. & Canback, B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 32, 11–16. https://doi.org/10.1093/nar/gkh152 (2004).
Lagesen, K. et al. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. https://doi.org/10.1093/nar/gkm160 (2007).
Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. https://doi.org/10.1093/bioinformatics/btu153 (2014).
O’Leary, N. A. et al. Reference sequence (RefSeq) database at NCBI: Current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 44, D733-745. https://doi.org/10.1093/nar/gkv1189 (2016).
Huerta-Cepas, J. et al. eggNOG 50: A hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. https://doi.org/10.1093/nar/gky1085 (2019).
Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y. & Hattori, M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 32, D277-280. https://doi.org/10.1093/nar/gkh063 (2004).
Bairoch, A. & Apweiler, R. The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999. Nucleic Acids Res. 27, 49–54. https://doi.org/10.1093/nar/27.1.49 (1999).
Ashburner, M. et al. Gene ontology: Tool for the unification of biology The Gene Ontology Consortium. Nat. Genet. 25, 25–29. https://doi.org/10.1038/75556 (2000).
Stothard, P. & Wishart, D. S. Circular genome visualization and exploration using CGView. Bioinformatics 21, 537–539. https://doi.org/10.1093/bioinformatics/bti054 (2005).
Alikhan, N. F., Petty, N. K., Ben Zakour, N. L. & Beatson, S. A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 12, 402. https://doi.org/10.1186/1471-2164-12-402 (2011).
Sullivan, M. J., Petty, N. K. & Beatson, S. A. Easyfig: A genome comparison visualizer. Bioinformatics 27, 1009–1010. https://doi.org/10.1093/bioinformatics/btr039 (2011).
Page, A. J. et al. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. https://doi.org/10.1093/bioinformatics/btv421 (2015).
Krzywinski, M. et al. Circos: An information aesthetic for comparative genomics. Genome Res. 19, 1639–1645. https://doi.org/10.1101/gr.092759.109 (2009).
Treangen, T. J., Ondov, B. D., Koren, S. & Phillippy, A. M. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 15, 524. https://doi.org/10.1186/s13059-014-0524-x (2014).
Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43, e15. https://doi.org/10.1093/nar/gku1196 (2015).
This work was supported by the National Key R&D Program of China (2017YFD0500204), National Natural Science Foundation of China (32072323, 31772083, 32073022), Special fund for Technology Innovation of Hubei Province (2019AHB074), Teaching Research Programs of Hubei Province (2018323), and the National Innovation and Entrepreneurship Training Program for Undergraduates (202110504076).
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hu, Y., Fu, S., Zou, G. et al. Genome analysis provides insight into hyper-virulence of Streptococcus suis LSM178, a human strain with a novel sequence type 1005. Sci Rep 11, 23919 (2021). https://doi.org/10.1038/s41598-021-03370-0