Human oropharynx as natural reservoir of Streptobacillus hongkongensis

Recently, we reported the isolation of Streptobacillus hongkongensis sp. nov. from patients with quinsy or septic arthritis. In this study, we developed a PCR sequencing test after sulfamethoxazole/trimethoprim and nalidixic acid enrichment for detection of S. hongkongensis. During a three-month study period, among the throat swabs from 132 patients with acute pharyngitis and 264 controls, PCR and DNA sequencing confirmed that S. hongkongensis and S. hongkongensis-like bacteria were detected in 16 patients and 29 control samples, respectively. Among these 45 positive samples, five different sequence variants were detected. Phylogenetic analysis based on the 16S rRNA gene showed that sequence variant 1 was clustered with S. hongkongensis HKU33T/HKU34 with high bootstrap support; while the other four sequence variants formed another distinct cluster. When compared with the 16S rRNA gene of S. hongkongensis HKU33T, the five sequence variants possessed 97.5–100% sequence identities. Among sequence variants 2–5, their sequences showed ≥99.5% nucleotide identities to each other. Forty-two individuals (93.3%) only harbored one sequence variant. We showed that the human oropharynx is a reservoir of S. hongkongensis, but the bacterium is not associated with acute pharyngitis. Another undescribed novel Streptobacillus species is probably also residing in the human oropharynx.

Streptobacillus is one of the four genera within the family Leptotrichiaceae. It was first isolated from the blood of a patient who suffered from a rat bite in 1914 1 . Since the establishment of the genus Streptobacillus in 1925, this genus had contained only one single species, Streptobacillus moniliformis, for almost 90 years 2 . S. moniliformis is the causative agent of streptobacillary rat-bite fever 3 . It is also associated with amnionitis 4 , bacteremia 5-7 , brain abscess 8 , cutaneous and subcutaneous abscesses 9,10 , endocarditis 11,12 , female genital tract abscess 13 , palmoplantar pustulosis 14 , septic arthritis [15][16][17] , spinal epidural abscess 18 , splenic abscess 19 , spondylodiscitis with psoas abscess 20 , and synovitis 21 . Microbiologically, S. moniliformis is a Gram-negative facultative anaerobic bacillus that grows in chains. It is naturally harbored in the oral cavity 22 and upper respiratory tract 23,24 of rats.
In 2014, we reported the isolation of a novel Streptobacillus species, named S. hongkongensis, from the pus of a patient with quinsy and the elbow joint fluid of another patient with tophaceous gout and left elbow septic arthritis 25 . Both strains of S. hongkongensis were resistant to cotrimoxazole and nalidixic acid as determined by disk diffusion test (unpublished data). Since S. moniliformis colonizes the oral cavity of rats while S. hongkongensis could be isolated from peritonsillar abscess pus of human, we hypothesize that the oropharynx of human may be the natural reservoir of this bacterium or S. hongkongensis may be associated with acute pharyngitis. To test these hypotheses, we first developed an in-house molecular test for the detection of S. hongkongensis on throat swabs of patients with acute pharyngitis. Preliminary study showed that S. hongkongensis or S. hongkongensis-like sequences were detected in around 10% of 100 patient samples collected (unpublished data), indicating the presence of this bacterium in human oropharynx. In this molecular epidemiological study, in order to further test for Primer specificity. PCR of the partial 16S rRNA gene of S. hongkongensis strains HKU33 T and HKU34 using the primer pair LPW21593/LPW21594 yielded DNA fragments of about 700 bp with strong signal. For S. moniliformis CCUG 13453 T , Sneathia sanguinegens CCUG 41628 T , and "Sneathia amnii" CCUG 52976, no PCR product of expected size was detected (Fig. 1).

Molecular detection of S. hongkongensis and S. hongkongensis-like bacteria.
During the three-month period, 132 throat swab samples from 132 patients [male:female = 7:4, age (median, range) = 18.5 years, 0-92 years] with acute pharyngitis were sent to our clinical microbiology laboratory and 264 throat swab samples from 264 controls [male:female = 85:47, age (median, range) = 23 years, 0-95 years], including 132 healthy individuals and 132 out-patients not on antibiotics, without acute pharyngitis were included. Among the 132 patient and 264 control samples collected, PCR and DNA sequencing confirmed that S. hongkongensis and S. hongkongensis-like bacteria were detected in 16 patient and 29 control samples (Table 1).

Sequence analysis and phylogenetic characterization.
Among all of the 45 positive samples, five different sequence variants were detected. Phylogenetic analysis based on the 16S rRNA gene showed that sequence variant 1 was clustered with S. hongkongensis HKU33 T and HKU34 with high bootstrap support; while the other four sequence variants formed another distinct cluster (Fig. 2). When compared with the 16S rRNA gene of S. hongkongensis HKU33 T , sequence variants 1, 2, 3, 4, and 5 possessed 100%, 97.5%, 97.7%, 97.5%, and 97.7% sequence identities, respectively. Among sequence variants 2, 3, 4, and 5, their sequences showed ≥ 99.5% nucleotide identities. Forty-two individuals (93.3%) only harbored one sequence variant (variants 1, 2 and 4). However, in one patient and one control sample, two sequence variants (variants 2 and 4) were detected while in another control sample all the five sequence variants were detected (variants 1-5).

Discussion
In this study, we showed that human oropharynx is a reservoir of S. hongkongensis. Since numerous bacterial species reside in the human oropharynx, we developed a PCR test after sulfamethoxazole/trimethoprim and nalidixic acid enrichment, as the two S. hongkongensis strains we isolated previously, HKU33 T and HKU34, were resistant to these antibiotics. These antibiotics were also widely used as enrichment supplements for the detection of Streptobacillus-like bacteria previously 22,26 . Using this molecular method for detection, overall, S. hongkongensis (sequence variant 1), which possessed a 100% 16S rRNA gene sequence identity with those of the two previously isolated S. hongkongensis strains HKU33 T and HKU34 25 , was found in 2.0% of the throat swabs from all subjects tested. There was no significant difference between the detection rates for patients and controls, indicating that the bacterium is not associated with acute pharyngitis. All positive samples (sequence variant 1) were collected from children under the age of 14 years (p < 0.002 by Fisher's exact test) ( Table 1) without any sex predilection, suggesting children are the major population carrying the bacterium. On the other hand, it is notable that the two strains of S. hongkongensis which we previously isolated were recovered from the peritonsillar abscess pus of a 38-year-old man and the joint fluid of a 64-year-old man, respectively 25 , indicating that diseases caused by S. hongkongensis can occur in adults. Isolation and identification of more strains of S. hongkongensis from infective sites would reveal a more detail epidemiology and disease spectrum of this bacterium.
Another undescribed novel Streptobacillus species is likely also residing in the human oropharynx. In the present study, in addition to S. hongkongensis (sequence variant 1), four other 16S rRNA gene sequence variants (sequence variants 2-5), which possessed ≥ 99.5% nucleotide identities among themselves, 97.5-97.7% nucleotide identities to that of S. hongkongensis, and ≤ 95.9% nucleotide identities to that of the type strains of other Streptobacillus species, were detected. Since these four sequence variants form a unique cluster distinct from S. hongkongensis (sequence variant 1) as well as other known Streptobacillus species, it is likely that this represent another novel Streptobacillus species (Fig. 2). It is notable that the 16S rRNA gene sequences of S. felis and S. More novel Streptobacillus species could be discovered from other animals. The genus Streptobacillus had remained monotypic for almost a century 2 . In the last two decades, the use of 16S rRNA gene sequencing has led to the discovery of numerous novel bacteria 29 . In the past two years, four novel Streptobacillus species were described, including S. hongkongensis, S. felis, S. notomylis, and S. ratti. S. felis was isolated from the lung of a cat with acute suppurative to fibrinous, focally necrotizing bronchopneumonia with multifocal desquamation of type II pneumocytes and alveolar macrophages 28 ; whereas S. notomylis from the heart of an Australian spinifex hopping mouse (Notomys alexis) with septicemia as well as the oral swabs of rats (Rattus rattus) 22,27,30,31 ; and S. ratti from the oral swab of a rat 32 . In two recent animal microbiome projects, other potentially novel Streptobacillus species with ≥ 91.5% partial 16S rRNA gene sequence identities to known Streptobacillus species were also found in the oral cavities of a Cape ground squirrel (Xerus inauris) (accession number: HM590422, unpublished data) and dogs (Canis familiaris) (accession numbers: EU082091 and JN713542) 33,34 . In human, the two anatomical sites with the largest number of novel bacteria discovered in the recent two decades are the gastrointestinal tract and oral cavity 29 . For example, Streptococcus sinensis, a bacterium discovered in 2002 35 and associated with infective endocarditis in patients from different parts of the world [35][36][37][38] , was finally also found to be part of human oral flora 39 . This indicates that there should be yet numerous novel bacterial species remained undescribed from the oral cavity. We speculate that the oral cavities and oropharynx of other animals might be the reservoirs of their own Streptobacillus species. We anticipate that more novel Streptobacillus species will be described and more potentially novel Streptobacillus species will be revealed by molecular epidemiology and microbiome studies in the next decade.

Materials and Methods
Antimicrobial susceptibility. The in vitro antimicrobial susceptibilities of S. hongkongensis against sulfamethoxazole/trimethoprim (range: 1.5625/0.3125 to 400/80 μg/ml) and nalidixic acid (range: 0.15625 to 40 μg/ml) were determined to test for the enrichment condition. Briefly, drug powders of sulfamethoxazole (Fluka, Switzerland), trimethoprim (Fluka), and nalidixic acid (Sigma-Aldrich, St. Louis, MO, USA) were dissolved in sterile dimethyl sulfoxide (Sigma-Aldrich) for the preparation of stock solutions. The stock solutions were then added into 5 ml of BHI broth (Oxoid, UK) to give a series of culture media supplemented with antimicrobial drugs at concentrations within the test range. Approximately 2.5 × 10 6 cfu of S. hongkongensis HKU33 T was then added into each of the test medium (final cell density = 5 × 10 5 cfu/ml) and the inoculated broths were incubated under an aerobic environment supplemented with 5% CO 2 at 37 °C for 3 days. The presence of growth at each drug concentration was then recorded. The test was also performed using another S. hongkongensis reference strain, HKU34.

Primer design for PCR detection. A pair of primers [LPW21593 (5′-TAGGCGGTTAAA
CAAGTCAGG-3′ ) and LPW21594 (5′-TGAGATTCGCTCCATATCGCTATTTCG-3′ )] (Invitrogen, Carlsbad, CA, USA) were designed by selecting the conserved regions of the 16S rRNA genes of S. hongkongensis HKU33 T and HKU34 (Fig. 3). The specificity of this primer pair was tested using the genomic DNA extracted from strains HKU33 T and HKU34, as well as the closely related reference strains S. moniliformis CCUG 13453 T , Sneathia sanguinegens CCUG 41628 T , and "Sneathia amnii" CCUG 52976, which were obtained from Culture Collection, University of Göteborg (CCUG), Sweden. Bacterial DNA was extracted using the alkaline lysis method. Briefly, cells from a single colony were resuspended in 200 μl of phosphate-buffered saline (PBS) (Oxoid), and 800 μl of NaOH (0.05 mol/L) (Sigma-Aldrich) was added into the bacterial suspensions. The cell suspension mixtures were then heated at 60 °C for 45 min. The lysed cell suspensions were then each neutralized by an addition of 240 μl of Tris-hydrochloric acid (pH 7.0) (Sigma-Aldrich). The neutralized suspensions were then centrifuged at 16,100 rcf for 5 min to remove the cell debris and the total bacterial DNA contained in the supernatant was collected and stored at − 20 °C until use.
PCR amplification of the partial 16S rRNA gene was performed using the primer pair LPW21593/LPW21594. To determine if S. hongkongensis is associated with acute pharyngitis, controls without acute pharyngitis were recruited from healthy volunteers and out-patients. Subjects were excluded if they had immunocompromising conditions including primary immunodeficiency, systemic or inhalational steroid, other immunosuppressive drugs, chemotherapy, diabetes mellitus, chronic organ failure, thalassemia major with regular transfusions and/or iron overload, and hematopoietic stem cell or solid organ transplantation. Subjects who received antibiotics in the preceding 2 weeks prior to throat swabs collection were also excluded. All throat swabs were processed immediately after their arrival at the laboratory.
Bacterial culture. Throat swabs were suspended in BHI broth supplemented with 400/80 μg/ml of sulfamethoxazole/trimethoprim and 40 μg/ml of nalidixic acid. The inoculated broths were then incubated under aerobic environment supplemented with 5% CO 2 at 37 °C for 3-5 days.
DNA extraction, PCR, cloning, and DNA sequencing. In order to extract total bacterial DNA from the broth cultures, cells in the broth suspensions were collected by centrifugation at 16,100 rcf for 10 min. Then, the cells were resuspended in 200 μl of phosphate-buffered saline (PBS), and DNA was extracted using the alkaline lysis method as described above. PCR of the partial 16S rRNA gene was performed using the primer pair LPW21593/LPW21594 as described above. DNA from S. hongkonensis HKU33 T and autoclaved distilled water were used as the positive and negative controls in each run of PCR, respectively. Preparation of PCR master mix, addition of DNA samples into the reaction tubes and post-PCR steps, including gel electrophoresis, purification of PCR products and DNA sequencing, were performed in separated rooms to avoid possible contamination. After agarose gel electrophoresis, the PCR products were purified using the QIAquick Gel Extraction Kit (QIAGEN) according to the manufacturer's protocol. Both strands of the PCR products were sequenced twice with an ABI Prism 3700 DNA Analyzer (Applied Biosystems), using the respective PCR primers. The sequencing electropherograms obtained were viewed using Chromas Lite 2.1.1 (Technelysium, Australia).
When ambiguous peaks were observed in the sequencing electropherograms, the PCR products obtained were cloned into plasmids for another trial of sequencing, which was performed according to our previous publication 40 . Briefly, freshly-prepared gel-purified PCR products were cloned into pCRII-TOPO vector using the TOPO TA Cloning Kit (Invitrogen) according to the manufacturer's protocol. The TA-ligated plasmids were then transformed into Escherichia coli DH5α (TaKaRa Bio, Japan) by electroporation. The electroporated cells were then grown on LB agar, Lennox (Difco, BD Diagnostic Systems, Sparks, MD, USA) with kanamycin (50 μg/ml) (Sigma-Aldrich) for the selection of positive transformants, as well as isopropyl β -D-1-thiogalactopyranoside (IPTG) (40 μg/ ml) (Sigma-Aldrich) and 5-bromo-4-chloro-3-indolyl-beta-D-galactopyranoside (X-gal) (100 μg/ml) (Promega, Madison, WI, USA) for blue-white screening. White colonies were then selected and grown in LB broth, Lennox (Difco) with kanamycin (50 μg/ml) overnight and the plasmids in the bacterial cells were extracted using the QIAprep Spin Miniprep Kit (QIAGEN) according to the manufacturer's protocol. The purified plasmids were then sequenced using the PCR primers directly. The sequencing electropherograms obtained were viewed using Chromas Lite 2.1.1.
Comparative sequence identity analyses and phylogenetic analyses. The sequences obtained from the cloned PCR products were analyzed by pairwise alignment, with the optimal GLOBAL alignment parameters, using BioEdit 7.2.0 41 . The sequences of the PCR products were also compared with sequences of closely related species from the DDBJ/ENA/GenBank databases by multiple sequence alignment using MUSCLE 3.8 42 and were then end-trimmed. Poorly aligned or divergent regions of the aligned, end-trimmed DNA sequences were removed using Gblocks 0.91b 43 with relaxed parameters, and 653 nucleotide positions were included in the final sequence alignment dataset. Tests for substitution models and phylogenetic tree construction using the maximum likelihood method were performed using MEGA 6.06 44 . Nucleotide sequence accession numbers. The 16S rRNA gene sequences of the five S. hongkongensis sequence variants detected have been deposited in the DDBJ/ENA/GenBank databases with the accession numbers LC097058-LC097062.