Search

Abstract Background Diabetic foot infections (DFIs) are the most common complications of diabetic foot ulcers (DFUs), and a significant cause of lower extremity amputation. In this study we used whole genome sequencing to characterize the clonal composition, virulence and resistance genetic determinants of 58 Staphylococcus/Mammaliicoccus spp. isolates from contralateral healthy skin and DFU from 44 hospitalized patients. Results S. aureus (n = 32) and S. epidermidis (n = 10) isolates were recovered from both DFUs and healthy skin, whereas, S. haemolyticus (n = 8), M. sciuri (n = 1), S. hominis (n = 1) and S. simulans (n = 3) were recovered exclusively from healthy skin. In contrast, S. caprae (n = 2) and S. saprophyticus (n = 1) were recovered only from DFUs. Among S. aureus isolates, MRSA were present with high prevalence (27/32, 84.4%), 18 of which (66.7%) were from DFUs and 9 (33.3%) from healthy skin. In contrast, the coagulase-negative Staphylococcus (CoNS)/Mammaliicoccus isolates (n = 26), in particular S. epidermidis and S. haemolyticus were more prevalent in healthy skin, (10/26, 38.5%) and (8/26, 30.8%), respectively. MLST, spa and SCCmec typing classified the 32 S. aureus isolates into 6 STs, ST672, ST80, ST241, ST1, ST97, ST291 and 4 unknown STs (STNF); 8 spa types, t044, t037, t3841, t1247, t127, t639, t937 and t9432 and 2 SCCmec types, type IV and type III(A). Among CoNS, the S. epidermidis isolates belonged to ST54, ST35 and ST640. S. haemolyticus belonged to ST3, ST25, ST29, ST1 and ST56. The sole M. sciuri isolate was found to carry an SCCmec type III(A). A wide range of virulence genes and antimicrobial resistance genes were found among our isolates, with varying distribution between species or STs. The pan-genome analysis revealed a highly clonal population of Staphylococcus isolates, particularly among S. aureus isolates. Interestingly, the majority of S. aureus isolates including MRSA, recovered from the healthy skin and DFUs of the same patient belonged to the same clone and exhibited similar virulence/resistance genotype. Conclusions Our study provides clinically relevant information on the population profile, virulence and antibiotic resistance of Staphylococcus/Mammaliicoccus spp. in DFIs, which could serve as a basis for further studies on these as well as other groups of pathogens associated with DFIs.

Epidemiologic and genomic investigation of SARS-CoV-2 infections in members of Japan’s national wrestling team after participation in international tournaments in 2021 revealed multiple lineages of SARS-CoV-2 not reported in Japan. The attack rate among wrestlers was high. Results suggest possible transmission during matches. We recommend early case detection and response practices.

Abstract The rapid spread of the SARS-CoV-2 Variant of Concern (VOC) Gamma in Amazonas during early 2021 fueled a second large COVID-19 epidemic wave and raised concern about the potential role of reinfections. Very few cases of reinfection associated with the VOC Gamma have been reported to date, and their potential impact on clinical, immunological, and virological parameters remains largely unexplored. Here we describe 25 cases of SARS-CoV-2 reinfection in Brazil. SARS-CoV-2 genomic analysis confirmed that individuals were primo-infected with distinct viral lineages between March and December 2020 (B.1.1, B.1.1.28, B.1.1.33, B.1.195, and P.2) and reinfected with the VOC Gamma between 3 to 12 months after primo-infection. We found a similar mean cycle threshold (Ct) value and limited intra-host viral diversity in both primo-infection and reinfection samples. Sera of 14 patients tested 10–75 days after reinfection displayed detectable neutralizing antibodies (NAb) titers against SARS-CoV-2 variants that circulated before (B.1.*), during (Gamma), and after (Delta and Omicron) the second epidemic wave in Brazil. All individuals had milder or no symptoms after reinfection, and none required hospitalization. These findings demonstrate that individuals reinfected with the VOC Gamma may display relatively high RNA viral loads at the upper respiratory tract after reinfection, thus contributing to onward viral transmissions. Despite this, our study points to a low overall risk of severe Gamma reinfections, supporting that the abrupt increase in hospital admissions and deaths observed in Amazonas and other Brazilian states during the Gamma wave was mostly driven by primary infections. Our findings also indicate that most individuals analyzed developed a high anti-SARS-CoV-2 NAb response after reinfection that may provide some protection against reinfection or disease by different SARS-CoV-2 variants.

ABSTRACTClostridium tetani is the etiological agent of the disease tetanus, resulting in various symptoms, including neck tension, tongue-twisting, a decline in walking ability, and death in severe cases. To date, publicly available genome sequences of C. tetani are limited. C. tetani was isolated from soil in Kumamoto Prefecture, Japan, between 2020 and 2022. Subsequently, bacteriological, biochemical, genetic, and immunological tests were conducted. Among the 46 soil sampling sites, 151 C. tetani strains were isolated from 33 sites, with a 71.7% positivity rate. Comparative genomic analysis subdivided clade 1 into clades 1-1, 1-2, and 1-3. The enzyme-linked immunosorbent assay revealed that partial clade 1-3 strains produced over sevenfold more tetanus neurotoxin (TeNT) than the reference strain (first isolate of C. tetani KHSU-154301-001). RNA-Seq analysis supported this elevated tetanus toxin tetX gene expression. Additionally, tetX-negative strains were isolated from the clade, indicating that TeNT may not be essential for ensuring the survival of their offspring without its representative toxin-causing tetanus. Five strains acquired the tetracycline resistance gene tet(M) in the tetX plasmid, resulting in significantly reduced susceptibility to tetracycline. In conclusion, although we investigated C. tetani isolates from a specific region of the Kumamoto Prefecture, we demonstrated that these isolates exhibit variable genomic features, toxin productivity, and substantial acquired antimicrobial resistance. These findings indicate that multiple C. tetani lineages have successfully adapted to the soil in Japan, as evidenced by the additional genome information and toxin productivity assays in this study.IMPORTANCEC. tetani is a spore-forming, anaerobic bacterium that produces a toxin causing muscle stiffness and paralysis. Tetanus is preventable with the toxoid vaccine, but it remains a significant public health threat in regions with low vaccine coverage. However, there are relatively few isolates and limited genomic information available worldwide. In Japan, about 100 cases are reported each year, but there have been no nationwide surveys of isolates, and no genomic information from Japanese isolates has been published. In our study, we analyzed the genomes of 151 strains from a limited survey of soil in Kumamoto, Japan. Our findings revealed a high degree of genetic diversity, and we also identified a subset of strains that produced significantly more toxin, which provides new insights into the pathogenesis of tetanus. Our findings lay the foundation for future studies to investigate the distribution and evolution of C. tetani in Japan and neighboring countries.

ABSTRACT Wastewater can be useful in monitoring the spread of antimicrobial resistance (AMR) within a hospital. The abundance of antibiotic resistance genes (ARGs) in hospital effluent was assessed using metagenomic sequencing (mDNA-seq) and hybrid capture (xHYB). mDNA-seq analysis and subsequent xHYB targeted enrichment were conducted on two effluent samples per month from November 2018 to May 2021. Reads per kilobase per million (RPKM) values were calculated for all 1,272 ARGs in the constructed database. The monthly numbers of patients with presumed extended-spectrum β-lactamase (ESBL)-producing and metallo-β-lactamase (MBL)-producing bacteria, methicillin-resistant Staphylococcus aureus (MRSA), and vancomycin-resistant enterococci (VRE) were compared with the monthly RPKM values of blaCTX-M, blaIMP, mecA, vanA, and vanB by xHYB. The average RPKM value for all ARGs detected by xHYB was significantly higher than that of mDNA-seq (665, 225, and 328, respectively, and P < 0.05). The average number of patients with ESBL producers and RPKM values of blaCTX-M-1 genes in 2020 were significantly higher than that in 2019 (17 and 13 patients per month and 921 vs 232 per month, respectively, both P < 0.05). The average numbers of patients with MBL-producers, MRSA, and VRE were 1, 28, and 0 per month, respectively, while the average RPKM values of blaIMP, mecA, vanA, and vanB were 6,163, 6, 0, and 126 per month, respectively. Monitoring ARGs in hospital effluent using xHYB was found to be more useful than conventional mDNA-seq in detecting ARGs including blaCTX-M, blaIMP, and vanB, which are important for infection control. IMPORTANCE Environmental ARGs play a crucial role in the emergence and spread of AMR that constitutes a significant global health threat. One major source of ARGs is effluent from healthcare facilities, where patients are frequently administered antimicrobials. Culture-independent methods, including metagenomics, can detect environmental ARGs carried by non-culturable bacteria and extracellular ARGs. mDNA-seq is one of the most comprehensive methods for environmental ARG surveillance; however, its sensitivity is insufficient for wastewater surveillance. This study demonstrates that xHYB appropriately monitors ARGs in hospital effluent for sensitive identification of nosocomial AMR dissemination. Correlations were observed between the numbers of inpatients with antibiotic-resistant bacteria and the ARG RPKM values in hospital effluent over time. ARG surveillance in hospital effluent using the highly sensitive and specific xHYB method could improve our understanding of the emergence and spread of AMR within a hospital.

PurposeAcute abdominal infections can be fatal if the causative organism (s) are misidentified. The spread of antimicrobial-resistant bacteria has become a serious problem worldwide, making antibiotic selection extremely difficult. Using quantitative metagenomic analysis, we evaluated a commercial multiplex polymerase chain reaction (PCR) system (FilmArray™, bioMérieux, Marcy-l’Étoile, France) for the rapid identification of causative bacteria.MethodsThe cases of 10 patients with acute abdominal infections were enrolled in this retrospective study. There were six cases of perforated peritonitis and four cases of intraabdominal abscess. Fluid collected from the acute surgical abdominal infections were examined.ResultsAll specimens tested positive for microorganisms in culture, and six involved two or more microorganisms. Using the multiplex PCR system, nine of ten specimens were found to involve at least one microorganism. One specimen was not included in the multiplex PCR system panel. Nineteen of 21 microorganisms (90.5%) isolated by culture were detected by the multiplex PCR system. Microorganisms with very small numbers of reads (19 reads) were detectable.ConclusionThis multiplex PCR system showed a high detection rate for causative microorganisms in ascites and intraabdominal abscesses. This system may be suitable as an affordable rapid identification system for causative bacteria in these cases.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection among pets owned by coronavirus disease 2019 (COVID-19) patients has been reported around the world. However, how often the animals are exposed to SARS-CoV-2 by their owners is still unclear. We have collected swab samples from COVID-19 patients’ pets and performed real-time RT-PCR to detect the viral genome. In total, 8 of 53 dogs (15.1%) and 5 of 34 cats (14.7%) tested positive for the SARS-CoV-2 N gene. The result of a virus neutralization (VN) test also showed VN antibodies in four cats and six dogs. Our results indicate that the virus often passed from infected owners to their pets, which then excreted the virus despite having no or mild clinical signs.

Mumps virus (MuV) is the etiological agent of mumps, a disease characterized by painful swelling of the parotid glands and often accompanied by severe complications. To understand the molecular mechanism of MuV infection, a functional analysis of the involved host factors is required. However, little is known about the host factors involved in MuV infection, especially those involved in the late stage of infection. Here, we identified 638 host proteins that have close proximity to MuV glycoproteins, which are a major component of the viral particles, by proximity labeling and examined comprehensive protein-protein interaction networks of the host proteins. From siRNA screening and immunoprecipitation results, we found that a SNARE subfamily protein, USE1, bound specifically to the MuV fusion (F) protein and was important for MuV propagation. In addition, USE1 plays a role in complete N-linked glycosylation and expression of the MuV F protein.

A 72-year-old patient was admitted to the intensive care unit due to acute respiratory distress syndrome caused by COVID-19. On day 20, the patient experienced shock. The electrocardiogram showed ST segment elevation in leads V3–V6 and severe left ventricular dysfunction with an ejection fraction of 35%–40%. The left ventricle showed basal hypokinesis and apical akinesis, while the creatine kinase level was normal, indicating Takotsubo cardiomyopathy. On day 24, the patient died of multiple organ failure. In post-mortem biopsy, SARS-CoV-2 antigen was detected in cardiomyocytes by immunostaining. Moreover, SARS-CoV-2 RNA was detected in heart tissue. We need to further analyse the direct link between SARS-CoV-2 and cardiomyocytes.

Abstract Borrelia miyamotoi, a member of the tick-borne relapsing fever spirochetes, shows a serum-resistant phenotype in vitro. This ability of B. miyamotoi may contribute to bacterial evasion of the host innate immune system. To investigate the molecular mechanism of serum-resistance, we constructed a membrane protein-encoding gene library of B. miyamotoi using Borrelia garinii strain HT59G, which shows a transformable and serum-susceptible phenotype. By screening the library, we found that bom1093 and bom1515 of B. miyamotoi provided a serum-resistant phenotype to the recipient B. garinii. These B. miyamotoi genes are predicted to encode P35-like antigen genes and are conserved among relapsing fever borreliae. Functional analysis revealed that BOM1093 bound to serum vitronectin and that the C-terminal region of BOM1093 was involved in the vitronectin-binding property. Importantly, the B. garinii transformant was not serum-resistant when the C terminus-truncated BOM1093 was expressed. We also observed that the depletion of vitronectin from human serum enhances the bactericidal activity of BOM1093 expressing B. garinii, and the survival rate of BOM1093 expressing B. garinii in vitronectin-depleted serum is enhanced by the addition of purified vitronectin. Our data suggests that B. miyamotoi utilize BOM1093-mediated binding to vitronectin as a mechanism of serum resistance.

Objectives: The Escherichia coli O25-ST131 clone is responsible for global dissemination of the blaCTX-M gene. However, the prevalence of this clone in the digestive tract, devoid of antimicrobial selection, and its molecular epidemiology remain unclear. In this study, we examined the origin of blaCTX-M-positive E. coli O25-ST131 and its distribution. Methods: We separately sequenced the chromosomal and plasmid genomes of 50 E. coli O25 isolates obtained from faecal samples of patients with diarrhoea in Japan. Results: Although 36 (72%) of 50 E. coli O25 isolates were ST131, only 6 harboured blaCTX-M. According to the fimH and ybbW sequences and fluoroquinolone susceptibility, H30R1 isolates were dominant (27/36; 75%) and possessed IncFII-FIA-FIB with FAB formula subtype F1:A2:B20 plasmids at a high frequency (24/27; 89%). The F1:A2:B20 plasmids possessed more resistance genes such as blaTEM-1, aminoglycoside resistance genes and trimethoprim/sulfamethoxazole resistance genes compared with non-F1:A2:B20 plasmids. In contrast, only one blaCTX-M-14 gene was located on the F1:A2:B20 plasmids, whereas the other three were located on IncFII (F4:A-:B-) (n = 1) and IncZ (n = 2) plasmids. Two H30Rx-ST131 isolates harboured blaCTX-M-15: one was on the chromosome and the other on the IncFIA-R plasmid. The stability and conjugation ability of the F1:A2:B20 plasmids were compared with those of non-F1:A2:B20 plasmids, which revealed higher stability but lower conjugative ability. Conclusions: These results suggest that E. coli H30R1-ST131 is a multidrug-resistant clone containing several resistance genes in the F1:A2:B20 plasmid, which were widely distributed before the acquisition of blaCTX-M.

We report a severe case of Chromobacterium haemolyticum pneumonia associated with near-drowning and detail the investigation of the pathogen and river water. Our genomic and environmental investigation demonstrated that river water in a temperate region can be a source of C. haemolyticum causing human infections.

An autopsy of a patient in Japan with coronavirus disease indicated pneumonia lung pathology, manifested as diffuse alveolar damage. We detected severe acute respiratory syndrome coronavirus 2 antigen in alveolar epithelial cells and macrophages. Coronavirus disease is essentially a lower respiratory tract disease characterized by direct viral injury of alveolar epithelial cells.

Hemotropic mycoplasmas are common pathogens in animals, but it remains unclear what role these pathogens play in human infections. We report clinical and biologic characterization of Candidatus Mycoplasma haemohominis infection in a 42-year-old man in Japan. The patient had severe hemophagocytic syndrome 1 month after an accidental needlestick injury. Metagenomic deep sequencing identified Candidatus M. haemohominis and determined its draft genome for an isolate from serum of the patient. A high copy number of the Candidatus M. haemohominis genome was detected in serum and bone marrow samples. Electron microscopy examination showed morphologic characteristics of Candidatus M. haemohominis. Levofloxacin monotherapy induced resistance caused by a gyrase A gene mutation in the quinolone resistance–determining region, but a combination treatment with moxifloxacin and minocycline was effective. We identified Candidatus M. haemohominis in a patient who had life-threatening symptoms related to multiple organ infection. Human infection with this mycoplasma might occur more frequently than has been generally recognized.

The coronavirus disease 2019 (COVID-19) has caused a serious disease burden and poses a tremendous public health challenge worldwide. Here, we report a comprehensive epidemiological and genomic analysis of SARS-CoV-2 from 63 patients in Niigata City, a medium-sized Japanese city, during the early phase of the pandemic, between February and May 2020. Among the 63 patients, 32 (51%) were female, with a mean (±standard deviation) age of 47.9 ± 22.3 years. Fever (65%, 41/63), malaise (51%, 32/63), and cough (35%, 22/63) were the most common clinical symptoms. The median Ct value after the onset of symptoms lowered within 9 days at 20.9 cycles (interquartile range, 17–26 cycles), but after 10 days, the median Ct value exceeded 30 cycles (p < 0.001). Of the 63 cases, 27 were distributed in the first epidemic wave and 33 in the second, and between the two waves, three cases from abroad were identified. The first wave was epidemiologically characterized by a single cluster related to indoor sports activity spread in closed settings, which included mixing indoors with families, relatives, and colleagues. The second wave showed more epidemiologically diversified events, with most index cases not related to each other. Almost all secondary cases were infected by droplets or aerosols from closed indoor settings, but at least two cases in the first wave were suspected to be contact infections. Results of the genomic analysis identified two possible clusters in Niigata City, the first of which was attributed to clade S (19B by Nexstrain clade) with a monophyletic group derived from the Wuhan prototype strain but that of the second wave was polyphyletic suggesting multiple introductions, and the clade was changed to GR (20B), which mainly spread in Europe in early 2020. These findings depict characteristics of SARS-CoV-2 transmission in the early stages in local community settings during February to May 2020 in Japan, and this integrated approach of epidemiological and genomic analysis may provide valuable information for public health policy decision-making for successful containment of chains of infection.

To characterize environmental antimicrobial resistance (AMR) in urban areas, extended-spectrum β-lactamase- (ESBL)/carbapenemase-producing bacteria (EPB/CPB, respectively) from urban wastewater treatment plant effluents in Tokyo were isolated on CHROMagar ESBL plate. Complete genome sequence analysis, including plasmids, indicated that 126 CTX-M-positive isolates (31%) were identified among the 404 obtained isolates. The CTX-M-9 group was predominant (n = 65, 52%), followed by the CTX-M-1 group (n = 44, 35%). Comparative genome analysis revealed that CTX-M-27-positive E. coli O16:H5-ST131-fimH41 exhibited a stable genome structure and clonal-global dissemination. Plasmidome network analysis revealed that 304 complete plasmid sequences among 85 isolates were grouped into 14 incompatibility (Inc) network communities (Co1 to Co14). Co10 consisted of primarily IncFIA/IncFIB plasmids harboring blaCTX-M in E. coli, whereas Co12 consisted primarily of IncFIA(HI1)/Inc FIB(K) plasmids harboring blaCTX-M, blaKPC, and blaGES in Klebsiella spp. Co11 was markedly located around Co10 and Co12. Co11 exhibited blaCTX-M, blaKPC, and blaNDM, and was mainly detected in E. coli and Klebsiella spp. from human and animal sources, suggesting a mutual role of Co11 in horizontal gene transfer between E. coli and Klebsiella spp. This comprehensive resistome analysis uncovers the mode of relational transfer among bacterial species, highlighting the potential source of AMR burden on public health in urban communities.

Salmonella enterica subsp. enterica serovar Typhimurium sequence type 34 (ST34) and its monophasic variant (Salmonella 4,[5],12:i:-) are among the most frequently isolated clones from both humans and animals worldwide. Our previous study demonstrated that Salmonella Typhimurium/4,[5],12:i:- strains isolated in Japan could be classified into nine clades and that clade 9 consisted of ST34 strains. In Japan, ST34/clade 9 was first found in the 1990s and has become predominant among food animals in recent years. In the present study, we analyzed the whole genome-based phylogenetic relationships and temporal information of 214 Salmonella Typhimurium/4,[5],12:i:- ST34/clade 9 strains isolated from 1998 to 2017 in Japan. The 214 strains were classified into two sublineages: the newly identified clade 9–2 diverged from clade 9 in the early 2000s and has predominated in recent years. Clonally expanding subclades in clades 9–1 or 9–2 lacked Gifsy-1 or HP1 prophages, respectively, and some strains in these subclades acquired plasmids encoding antimicrobial resistance genes. Additional genome reduction around the fljB gene encoding the phase 2-H antigen was generated by an IS26-mediated deletion adjacent to the transposon in clade 9–2. Although most of the clade 9 strains were isolated from cattle in Japan, the clonally expanding subclades in clade 9–2 (i.e., all and 24% strains of subclades 9–2a and 9–2b, respectively) were isolated from swine. The spread of clade 9 in recent years among food animals in Japan was responsible for the emergence of multiple host-adapted sublineages involving the clonally expanding subclades generated by mobile genetic element-mediated microevolution.

In Japan’s Kanto region, the number of Salmonella enterica serovar Chester infections increased temporarily between 2014 and 2016. Concurrently with this temporal increase in the Kanto region, S. Chester isolates belonging to one clonal group were causing repetitive outbreaks in Europe. A recent study reported that the European outbreaks were associated with travelers who had been exposed to contaminated food in Morocco, possibly seafood. Because Japan imports a large amount of seafood from Morocco, we aimed to establish whether the temporal increase in S. Chester infections in the Kanto region was associated with imported Moroccan seafood. Short sequence reads from the whole-genome sequencing of 47 S. Chester isolates from people in the Kanto region (2014–2016), and the additional genome sequences from 58 isolates from the European outbreaks, were analyzed. The reads were compared with the complete genome sequence from a S. Chester reference strain, and 347 single nucleotide polymorphisms (SNPs) were identified. These SNPs were used in this study. Cluster and Bayesian cluster analyses showed that the Japanese and European isolates fell into two different clusters. Therefore, ΦPT and IAS values were calculated to evaluate genetic differences between these clusters. The results revealed that the Japanese and European isolates were genetically distinct populations. Our root-to-tip analysis showed that the Japanese isolates originating from one clone had accumulated mutations, suggesting that an emergence of this organism occurred. A minimum spanning tree analysis demonstrated no correlation between genetic and geographical distances in the Japanese isolates, suggesting that the emergence of the serovar in the Kanto region did not involve person-to-person contact; rather, it occurred through food consumption. The dN/dS ratio indicated that the Japanese strain has evolved under positive selection pressure. Generally, a population of bacterial clones in a reservoir faces negative selection pressure. Therefore, the Japanese strain must have existed outside of any reservoir during its emergence. In conclusion, S. Chester isolates originating from one clone probably emerged in the Kanto region via the consumption of contaminated foods other than imported Moroccan seafood. The emerging strain may have not established a reservoir for survival in the food supply chain resulting in its disappearance after 2017.

We conducted molecular typing of a Corynebacterium ulcerans isolate from a woman who died in Japan in 2016. Genomic DNA modification might have affected the isolate’s ribotyping profile. Multilocus sequence typing results (sequence type 337) were more accurate. Whole-genome sequencing had greater ability to discriminate lineages at high resolution.

Abstract Background Torque teno virus-induced aseptic meningitis has not been documented, although torque teno virus infections still remain under consideration for etiological agents. This study identified a torque teno virus sequence using next generation sequencing and immunoglobulin M response to the torque teno virus antigen, therefore, that would be a comprehensive diagnosis for torque teno virus infection. Case presentation A 2-month-old Japanese boy was brought to our hospital because he was irritable, drowsy, and lethargic. He was admitted based on his test results which indicated the possibility of septic meningitis. He was started on treatment with high-dose antibiotics and steroids. On the third day of hospitalization, he became afebrile with improvement in his general status and was discharged on the sixth day. He had no developmental problems for up to 1 year after discharge. Metagenomic ribonucleic acid-Seq pathogen detection using next generation sequencing of a sample of his cerebrospinal fluid, which was collected at admission, revealed three short reads homologous to those in torque teno virus out of a total of 1,708,516 reads. This finding indicated that our patient was positive compared to the torque teno virus-negative cerebrospinal fluid samples (controls) from 13 other patients. The torque teno virus has been shown to have a whole genome sequence of 2810 nt by polymerase chain reaction. We prepared a recombinant GP2 antigen from torque teno virus and used it to study our patient’s anti-torque teno virus immune response. An anti-GP2 serum immunoglobulin M response was detected, providing further supportive evidence of torque teno virus infection. Conclusions This case speculates that torque teno virus-induced aseptic meningitis has a good course. New technologies like next generation sequencing can help in the identification of such cases, and an accumulation of future cases is expected.

Abstract Streptomyces rochei 7434AN4 produces two structurally unrelated polyketide antibiotics, lankacidin and lankamycin, and carries three linear plasmids, pSLA2-L (211 kb), -M (113 kb), and -S (18 kb), whose nucleotide sequences were previously reported. The complete nucleotide sequence of the S. rochei chromosome has now been determined using the long-read PacBio RS-II sequencing together with short-read Illumina Genome Analyzer IIx sequencing and Roche 454 pyrosequencing techniques. The assembled sequence revealed an 8,364,802-bp linear chromosome with a high G + C content of 71.7% and 7,568 protein-coding ORFs. Thus, the gross genome size of S. rochei 7434AN4 was confirmed to be 8,706,406 bp including the three linear plasmids. Consistent with our previous study, a tap-tpg gene pair, which is essential for the maintenance of a linear topology of Streptomyces genomes, was not found on the chromosome. Remarkably, the S. rochei chromosome contains seven ribosomal RNA (rrn) operons (16S-23S-5S), although Streptomyces species generally contain six rrn operons. Based on 2ndFind and antiSMASH platforms, the S. rochei chromosome harbors at least 35 secondary metabolite biosynthetic gene clusters, including those for the 28-membered polyene macrolide pentamycin and the azoxyalkene compound KA57-A.

Summary: Subtilase cytotoxin (SubAB) is a virulence factor produced by locus of enterocyte effacement-negative Shiga-toxigenic Escherichia coli strains. The toxin recognizes sialoglycans for entry and cleaves an endoplasmic reticulum chaperon, binding immunoglobulin protein, to cause cell death. However, no systematic screening has yet been performed to identify critical host factors. Here, we performed a genome-wide CRISPR/Cas9 knockout screen for SubAB-induced cell death and identified various sialoglycan-related and membrane-trafficking genes. Analysis of glycan-deficient cells demonstrated that not only N-glycans but also O-glycans serve as SubAB receptors. In addition, SLC39A9, which is a predicted zinc transporter, as well as KDELRs and JTB, were required for SubAB to induce maximal cell death. Disruption of the SLC39A9 gene markedly reduced both complex-type N-glycans and core 1 O-glycans, and the O-glycan reduction was attributed to the reduction of core 1 synthase (C1GalT1). These results provide insights into the post-transcriptional regulation of glycosyltransferases by SLC39A9, as well as sialoglycan species as SubAB receptors. : Biological Sciences; Biochemistry; Molecular Biology; Microbiology; Cell Biology Subject Areas: Biological Sciences, Biochemistry, Molecular Biology, Microbiology, Cell Biology

Summary: Glycosphingolipids (GSLs) are produced by various GSL-synthesizing enzymes, but post-translational regulation of these enzymes is incompletely understood. To address this knowledge disparity, we focused on biosynthesis of globotriaosylceramide (Gb3), the Shiga toxin (STx) receptor, and performed a genome-wide CRISPR/CAS9 knockout screen in HeLa cells using STx1-mediated cytotoxicity. We identified various genes including sphingolipid-related genes and membrane-trafficking genes. In addition, we found two proteins, LAPTM4A and TM9SF2, for which physiological roles remain elusive. Disruption of either LAPTM4A or TM9SF2 genes reduced Gb3 biosynthesis, resulting in accumulation of its precursor, lactosylceramide. Loss of LAPTM4A decreased endogenous Gb3 synthase activity in a post-transcriptional mechanism, whereas loss of TM9SF2 did not affect Gb3 synthase activity but instead disrupted localization of Gb3 synthase. Furthermore, the Gb3-regulating activity of TM9SF2 was conserved in the TM9SF family. These results provide mechanistic insight into the post-translational regulation of the activity and localization of Gb3 synthase. : Molecular Biology; Molecular Mechanism of Behavior; Cell Biology; Functional Aspects of Cell Biology Subject Areas: Molecular Biology, Molecular Mechanism of Behavior, Cell Biology, Functional Aspects of Cell Biology

The emergence and spread of antimicrobial resistance (AMR) has become a persistent problem globally. In this study, an ozone treatment facility was established for an advanced hospital wastewater treatment in a core hospital facility in an urban area in Japan to evaluate the inactivation of antimicrobial-resistant bacteria and antimicrobials. Metagenomic DNA-seq analysis and the isolation of potential extended-spectrum β-lactamase (ESBL)-producing bacteria suggested that ozone exposure for at least 20 min is required for the adequate inactivation of DNA and ESBL-producing bacteria. Escherichia coli and Klebsiella species were markedly susceptible to 20-min ozone exposure, whereas Raoultella ornithinolytica and Pseudomonas putida were isolated even after an 80-min exposure. These ozone-resistant bacteria might play a pivotal role as AMR reservoirs in the environment. Nine antimicrobials (ampicillin, cefdinir, cefpodoxime, ciprofloxacin, levofloxacin, clarithromycin, chlortetracycline, minocycline, and vancomycin) were detected at 373 ng/L to 27 μg/L in the hospital wastewater, and these were removed (96–100% removal) after a 40-min treatment. These results facilitate a comprehensive understanding of the AMR risk posed by hospital wastewater and provides insights for devising strategies to eliminate or mitigate the burden of antimicrobial-resistant bacteria and the flow of antimicrobials into the environment. To the best of our knowledge, this is the first report on the implementation of a batch-type, plant-scale ozone treatment system in a hospital facility to execute and evaluate the inactivation of drug-resistant bacteria and antimicrobials.

AbstractThe worldwide eruption of coronavirus disease 2019 (COVID-19) that began in Wuhan, China in late 2019 reached 10 million cases by late June 2020. In order to understand the epidemiological landscape of the COVID-19 pandemic, many studies have attempted to elucidate phylogenetic relationships between collected viral genome sequences using haplotype networks. However, currently available applications for network visualization are not suited to understand the COVID-19 epidemic spatiotemporally due to functional limitations that motivated us to develop Haplotype Explorer, an intuitive tool for visualizing and exploring haplotype networks. Haplotype Explorer enables to dissect epidemiological consequences via interactive node filters and provides the perspective on infectious disease dynamics depend on regions and time, such as introduction, outbreak, expansion, and containment. Here, we demonstrate the effectiveness of Haplotype Explorer by showing features and an example of visualization. The demo using severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes are available at https://github.com/TKSjp/HaplotypeExplorer/blob/master/Example/

ABSTRACT After the first case of coronavirus disease 2019 (COVID-19) in Japan on 15 January 2020, multiple nationwide COVID-19 clusters were identified by the end of February. The Japanese government focused on mitigating the emerging COVID-19 clusters by conducting active nationwide epidemiological surveillance. However, an increasing number of cases continued to appear until early April 2020, many with unclear infection routes and no recent history of travel outside Japan. We aimed to evaluate the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome sequences from the COVID-19 cases that appeared until early April 2020 and to characterize their genealogical networks in order to demonstrate possible routes of spread in Japan. Nasopharyngeal specimens were collected from patients, and reverse transcription-quantitative PCR tests for SARS-CoV-2 were performed. Positive RNA samples were subjected to whole-genome sequencing, and a haplotype network analysis was performed. Some of the primary clusters identified during January and February 2020 in Japan descended directly from the Wuhan-Hu-1-related isolates from China and other distinct clusters. Clusters were almost contained until mid-March; the haplotype network analysis demonstrated that the COVID-19 cases from late March through early April may have created an additional large cluster related to the outbreak in Europe, leading to additional spread within Japan. In conclusion, genome surveillance has suggested that there were at least two distinct SARS-CoV-2 introductions into Japan from China and other countries. IMPORTANCE This study aimed to evaluate the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome sequences from COVID-19 cases and to characterize their genealogical networks to demonstrate possible routes of spread in Japan. We found that there were at least two distinct SARS-CoV-2 introductions into Japan, initially from China and subsequently from other countries, including Europe. Our findings can help understand how SARS-CoV-2 entered Japan and contribute to increased knowledge of SARS-CoV-2 in Asia and its association with implemented stay-at-home/shelter-in-place/self-restraint/lockdown measures. This study suggested that it is necessary to formulate a more efficient containment strategy using real-time genome surveillance to support epidemiological field investigations in order to highlight potential infection linkages and mitigate the next wave of COVID-19 in Japan.

Nationwide increases in Mycoplasma pneumoniae pneumonia cases in Japan were reported in 2011, 2012, 2015, and 2016. In this study, we isolated 554 M. pneumoniae strains in 4 areas in Japan (Kanagawa, Okayama, Osaka, and Saitama) between 2006 and 2019, and performed genotyping analysis. More than 80% of the strains isolated in 2011 and 2012 harbored type 1 p1 adhesin gene; however, strains harboring type 2 or its variant p1 gene increased in 2015 and 2016 and dominated after 2017. These findings suggested that a shift in the prevalent genotype of M. pneumoniae clinical strains occurred recently in Japan. More than 90% of the type 1 strains isolated after 2010 harbored macrolide-resistance mutations in their 23S rRNA gene, whereas most type 2 lineage strains had no such mutations. Consequently, the increase in type 2 lineage strains in Japan has reduced the macrolide resistance rate of clinical M. pneumoniae strains. During this analysis, we also identified M. pneumoniae strains carrying a novel variant type 1 p1 gene, and we classified it as type 1b. We then sequenced the genomes of 81 selected M. pneumoniae strains that we collected between 1976 and 2017 in Japan, and compared them with 156 M. pneumoniae genomes deposited in public databases to provide insights into the interpretation of M. pneumoniae genotyping methods, including p1 typing, multiple-locus variable-number tandem repeat analysis (MLVA), multi-locus sequence typing (MLST), and typing by 8 single-nucleotide polymorphism markers (SNP-8). As expected, p1 typing, MLST, and SNP-8 results exhibited good correlation with whole-genome SNP analysis results in terms of phylogenetic relationships; however, MLVA typing results were less comparable to those of the other methods. MLVA may be useful for the discrimination of strains derived from a single outbreak within a limited area; however, is not reliable for classification of strains collected from distantly separated areas at different time points. This study showed the usefulness of genome-based comparison of M. pneumoniae for molecular epidemiology. Genome sequencing of more strains will improve our understanding of global propagation routes of this pathogen and evolutionary aspects of M. pneumoniae strains.

Japan has reported 26 cases of coronavirus disease 2019 (COVID-19) linked to cruise tours on the River Nile in Egypt between March 5 and 15, 2020. Here, we characterized the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome of isolates from 10 travelers who returned from Egypt and from patients possibly associated with these travelers. We performed haplotype network analysis of SARS-CoV-2 isolates using genome-wide single-nucleotide variations. Our analysis identified two potential Egypt-related clusters from these imported cases, and these clusters were related to globally detected viruses in different countries.

JC polyomavirus (JCPyV, JCV) causes progressive multifocal leukoencephalopathy (PML) in immunocompromised hosts. JCPyV replicates in oligodendrocytes within the brain tissue of patients with PML. The JCPyV genome encodes a microRNA (miRNA) in the region encoding the large T antigen. JCPyV-encoded miRNA (miR-J1) has been detected in the tissue and cerebrospinal fluid samples of patients with PML; however, there are no reports describing the localization of polyomavirus-encoded miRNA in histological samples of patients with virus-associated diseases. In the present study, we detected high miR-J1 expression in the nuclei of JCPyV-infected cells in PML tissue samples via in situ hybridization. Additionally, in situ hybridization also revealed the expression of BK polyomavirus (BKPyV, BKV)-encoded miRNA in lesions of BKPyV-associated nephropathy. In situ hybridization for miR-J1-5p and -3p showed positive signals in 24/25 (96%) of PML tissues that were positive for JCPyV by immunohistochemistry. Higher copy numbers of miR-J1 were detected in PML tissues than in non-PML tissues by real-time reverse transcription PCR. Next generation sequencing showed that miR-J1-5p, a mature miRNA of primary miRNA, was predominant in the lesions compared with miR-J1-3p, another mature miRNA. Deletion or mutation of miR-J1 in recombinant JCPyV promoted the production of JCPyV-encoded proteins in cells transfected with JCPyV DNA, suggesting that polyomavirus-encoded miRNA may have a repressive role in viral replication in PML tissues. In situ hybridization for viral miRNA may be a useful diagnostic tool for PML.

Since December 2019, the coronavirus disease 2019 (COVID-19) caused by a novel coronavirus SARS-CoV-2 has rapidly spread to almost every nation in the world. Soon after the pandemic was recognized by epidemiologists, a group of biologists comprising the ARTIC Network, has devised a multiplexed polymerase chain reaction (PCR) protocol and primer set for targeted whole-genome amplification of SARS-CoV-2. The ARTIC primer set amplifies 98 amplicons, which are separated only in two PCRs, across a nearly entire viral genome. The original primer set and protocol showed a fairly small amplification bias when clinical samples with relatively high viral loads were used. However, as sample's viral load become low, rapid decrease in abundances of several amplicons were seen. In this report, we will show that dimer formations between some primers are the major cause of coverage bias in the multiplex PCR. Based on this, we propose 12 alternative primers in total in the ARTIC primer set that were predicted to be involved in 14 primer interactions. The resulting primer set, version N1 (NIID-1), exhibits improved overall coverage compared to the ARTIC Network's original (V1) and modified (V3) primer set.

In insects, the voltage-gated sodium channel (VGSC) is the primary target site of pyrethroid insecticides. Various amino acid substitutions in the VGSC protein, which are selected under insecticide pressure, are known to confer insecticide resistance. In the genome, the VGSC gene consists of more than 30 exons sparsely distributed across a large genomic region, which often exceeds 100 kbp. Due to this complex genomic structure, it is often challenging to genotype full coding nucleotide sequences (CDSs) of VGSC from individual genomic DNA (gDNA). In this study, we designed biotinylated oligonucleotide probes from CDSs of VGSC of Asian tiger mosquito, Aedes albopictus. The probe set effectively concentrated (>80,000-fold) all targeted regions of gene VGSC from pooled barcoded Illumina libraries each constructed from individual A. albopictus gDNAs. The probe set also captured all orthologous VGSC CDSs, except some tiny exons, from the gDNA of other Culicinae mosquitos, A. aegypti and Culex pipiens complex, with comparable efficiency as a result of the high nucleotide-level conservation of VGSC. To improve efficiency of the downstream bioinformatic process, we developed an automated pipeline-MoNaS (Mosquito Na+ channel mutation Search)-which calls amino acid substitutions in the VGSC from NGS reads and compares those to known resistance mutations. The proposed method and our bioinformatic tool should facilitate the discovery of novel amino acid variants conferring insecticide resistance on VGSC and population genetic studies on resistance alleles (with respect to the origin, selection, and migration etc.) in both clinically and agriculturally important insect pests.

During the 2016–17 winter season in Japan, human norovirus GII.P16-GII.2 strains (2016 strains) caused large outbreaks of acute gastroenteritis. Phylogenetic analyses suggested that the 2016 strains derived from the GII.2 strains detected during 2010–12. Immunochromatography between 2016 strains and the pre-2016 GII.2 strains showed similar reactivity.

Abstract Background Streptococcus pyogenes (group A Streptococcus [GAS]) is a major human pathogen that causes a wide spectrum of clinical manifestations. Although invasive GAS (iGAS) infections are relatively uncommon, emm3/ST15 GAS is a highly virulent, invasive, and pathogenic strain. Global molecular epidemiology analysis has suggested that the frequency of emm3 GAS has been recently increasing. Case presentation A 14-year-old patient was diagnosed with streptococcal toxic shock syndrome and severe pneumonia, impaired renal function, and rhabdomyolysis. GAS was isolated from a culture of endotracheal aspirates and designated as KS030. Comparative genome analysis suggested that KS030 is classified as emm3 (emm-type) and ST15 (multilocus sequencing typing [MLST]), which is similar to iGAS isolates identified in the UK (2013) and Switzerland (2015). Conclusions We conclude that the global dissemination of emm3/ST15 GAS strain has the potential to cause invasive disease.

Shiga toxin (STx) or Vero toxin is a virulence factor produced by enterohemorrhagic Escherichia coli. The toxin binds to the glycosphingolipid globotriaosylceramide (Gb3) for its entry, and causes cell death by inhibiting ribosome function. Previously, we performed a loss-of-function screen in HeLa cells using a human CRISPR knockout (KO) library and identified various host genes required for STx-induced cell death. To determine whether this library targeted to the human genome is applicable to non-human primate cells and to identify previously unrecognized factors crucial for STx-induced cell death, we herein performed a similar screen in the African green monkey kidney-derived Vero C1008 subline. Many genes relevant to metabolic enzymes and membrane trafficking were enriched, although the number of enriched genes was less than that obtained in the screening for HeLa cells. Of note, several genes that had not been enriched in the previous screening were enriched: one of these genes was SYS1, which encodes a multi-spanning membrane protein in the Golgi apparatus. In SYS1 KO Vero cells, expression of Gb3 and sphingomyelin was decreased, while that of glucosylceramide and lactosylceramide was increased. In addition, loss of SYS1 inhibited the biosynthesis of protein glycans, deformed the Golgi apparatus, and perturbed the localization of trans-Golgi network protein (TGN) 46. These results indicate that the human CRISPR KO library is applicable to Vero cell lines, and SYS1 has a widespread effect on glycan biosynthesis via regulation of intra-Golgi and endosome–TGN retrograde transports.

Abstract A multiple DNA inversion system, the shufflon, exists in incompatibility (Inc) I1 and I2 plasmids. The shufflon generates variants of the PilV protein, a minor component of the thin pilus. The shufflon is one of the most difficult regions for de novo genome assembly because of its structural diversity even in an isolated bacterial clone. We determined complete genome sequences, including those of IncI2 plasmids carrying mcr-1, of three Escherichia coli strains using single-molecule, real-time (SMRT) sequencing and Illumina sequencing. The sequences assembled using only SMRT sequencing contained misassembled regions in the shufflon. A hybrid analysis using SMRT and Illumina sequencing resolved the misassembled region and revealed that the three IncI2 plasmids, excluding the shufflon region, were highly conserved. Moreover, the abundance ratio of whole-shufflon structures could be determined by quantitative structural variation analysis of the SMRT data, suggesting that a remarkable heterogeneity of whole-shufflon structural variations exists in IncI2 plasmids. These findings indicate that remarkable rearrangement regions should be validated using both long-read and short-read sequencing data and that the structural variation of PilV in the shufflon might be closely related to phenotypic heterogeneity of plasmid-mediated transconjugation involved in horizontal gene transfer even in bacterial clonal populations.

ABSTRACT We recently detected a novel variant of an IMP-type metallo-β-lactamase gene (blaIMP-68) from meropenem-resistant but imipenem-susceptible Klebsiella pneumoniae TA6363 isolated in Tokyo, Japan. blaIMP-68 encodes a Ser262Gly point mutant of IMP-11, and transformation experiments showed that blaIMP-68 increased the MIC of carbapenems in recipient strains, whereas the MIC of imipenem was not greatly increased relative to that of other carbapenems, including meropenem. Kinetics experiments showed that IMP-68 imipenem-hydrolyzing activity was lower than that for other carbapenems, suggesting that the antimicrobial susceptibility profile of TA6363 originated from IMP-68 substrate specificity. Whole-genome sequencing showed that blaIMP-68 is harbored by the class 1 integron located on the IncL/M plasmid pTMTA63632 (88,953 bp), which was transferable via conjugation. The presence of plasmid-borne blaIMP-68 is notable, because it conferred antimicrobial resistance to carbapenems, except for imipenem, on Enterobacteriaceae and will likely affect treatment plans using antibacterial agents in clinical settings. IMPORTANCE IMP-type metallo-β-lactamases comprise one group of the “Big 5” carbapenemases. Here, a novel blaIMP-68 gene encoding IMP-68 (harboring a Ser262Gly point mutant of IMP-11) was discovered from meropenem-resistant but imipenem-susceptible Klebsiella pneumoniae TA6363. The Ser262Gly substitution was previously identified as important for substrate specificity according to a study of other IMP variants, including IMP-6. We confirmed that IMP-68 exhibited weaker imipenem-hydrolyzing activity than that for other carbapenems, demonstrating that the antimicrobial susceptibility profile of TA6363 originated from IMP-68 substrate specificity, with this likely to affect treatment strategies using antibacterial agents in clinical settings. Notably, the carbapenem resistance conferred by IMP-68 was undetectable based on the MIC of imipenem as a carbapenem representative, which demonstrates a comparable antimicrobial susceptibility profile to IMP-6-producing Enterobacteriaceae that previously spread in Japan due to lack of awareness of its existence.

Vancomycin (VAN)-intermediate-resistant Staphylococcus aureus (VISA) is continually isolated globally, with a systematic review suggesting a prevalence of 2% in all blood culture samples. Most VISA strains exhibit common characteristics, such as a thickened cell wall, reduced autolysis, and attenuated virulence. Here, based on multi-omics approaches, we have characterized clinical VISA isolates obtained through prolonged antimicrobial treatment in a single patient. All VISA isolates were isogenic, based on multi-locus sequence typing (MLST) ST5, SCCmec type II (2A), and spa type t17639. Core-genome single nucleotide variations (SNVs) found among thirteen isolates during the patient’s hospitalization, indicated clonality, but not notable genetic features of the VISA phenotype. We determined the complete genome sequence of VAN-susceptible strain KG-03 (minimum inhibitory concentration [MIC] 0.5 μg/mL) and two VISA strains, KG-18 and KG-22 (MIC 8.0 and 4.0 μg/mL, respectively). Comparative genome analysis showed remarkable strain-specific IS256 insertions. RNA-Seq transcriptome analysis revealed IS256-mediated overexpression of the walKR two-component system in VISA KG-18, possibly leading to modulation of cell wall integrity (lytM and sceD) and surface charge (mprF and dltABCD). In addition, secretome analysis indicated that cell wall-anchored proteins (Protein A, SasG, and SdrD) were significantly decreased. KG-18 and KG-22 exhibit thickened cell wall, and are relatively resistant to lysostaphin, which cleaves a staphylococcus-unique pentaglycine chain in the peptidoglycan. We conclude that KG-18 achieved reduced susceptibility to VAN by IS256-mediated WalKR overexpression, leading to a markedly thickened cell wall for trapping free VAN molecules with redundant D-Ala-D-Ala targets. In addition, a positively charged surface with lysyl-phosphatidylglycerol and depolarization of wall teichoic acid could contribute to inhibiting cationic daptomycin and VAN antimicrobial activity. Comparative omics approaches in this study strongly suggest that fully complete and annotated genome sequences will be indispensable for characterizing overall VISA phenotype.

The regulation of paramyxovirus RNA synthesis by host proteins is poorly understood. Here, we identified a novel regulation mechanism of paramyxovirus RNA synthesis by the Hsp90 co-chaperone R2TP complex. We showed that the R2TP complex interacted with the paramyxovirus polymerase L protein and that silencing of the R2TP complex led to uncontrolled upregulation of mumps virus (MuV) gene transcription but not genome replication. Regulation by the R2TP complex was critical for MuV replication and evasion of host innate immune responses. The R2TP complex also regulated measles virus (MeV) RNA synthesis, but its function was inhibitory and not beneficial to MeV, as MeV evaded host innate immune responses in the absence of the R2TP complex. The identification of the R2TP complex as a critical host factor sheds new light on the regulation of paramyxovirus RNA synthesis.

Candida auris is an invasive and multidrug-resistant ascomycetous yeast that is under global surveillance. All clinical cases of C. auris infection diagnosed from 1997 to 2019 in Japan were non-invasive and sporadic otitis media cases. In the present study, we performed whole-genome sequencing of seven C. auris strains isolated from patients with otitis media in Japan, all of which belonged to clade II. Comparative genome analysis using the high-quality draft genome sequences JCM 15448T revealed that single nucleotide variations (SNVs), clade-specific accessory genes, and copy number variations (CNVs) were identified in each C. auris clade. A total of 61 genes involved in cell wall and stress response-related functions was absent in clade II, and the pattern of conserved CNVs in each clade was more stable in clade II than in other clades. Our data suggest that the genomic structural diversity is stable in C. auris isolated from each biogeographic location, and Japanese strains isolated from patients with otitis media might belong to an ancestral type of C. auris. One Japanese strain, TWCC 58362, with reduced susceptibility to fluconazole, exhibited no mutation in ergosterol biosynthesis-related genes (ERG). However, TWCC 58362-specific variations, including SNVs, indels, and CNVs were detected, suggesting that gene duplication events in C. auris might contribute to antifungal drug resistance. Taken together, we demonstrated that genomic structural variations in C. auris could correlate to geographical dissemination, epidemiology, lesions in the host, and antifungal resistance.

ObjectivesAutomated online software tools that analyse whole genome sequencing (WGS) data without the need for bioinformatics expertise can motivate the implementation of WGS-based molecular drug susceptibility testing (DST) in routine diagnostic settings for tuberculosis (TB). Pyrazinamide (PZA) is a key drug for current and future TB treatment regimens; however, it was reported that predictive power for PZA resistance by the available tools is low. Therefore, this low predictive power may make users hesitant to use the tools. This study aimed to elucidate why and to uncover the real performance of the tools when taking into account their variation calling lists (manual inspection), not just their automated reporting system (default setting) that was evaluated by previous studies.MethodsWGS data from 191 datasets comprising 108 PZA-resistant and 83 susceptible strains were used to evaluate the potential performance of the available online tools (TB Profiler, TGS-TB, PhyResSE, and CASTB) for predicting phenotypic PZA resistance.ResultsWhen taking into consideration the variation calling lists, 73 variants in total (47 non-synonymous mutations and 26 indels) in pncA were detected by TGS-TB and PhyResSE, covering all mutations for the 108 PZA-resistant strains. The 73 variants were confirmed by Sanger sequencing. TB Profiler also detected all but three complete loss, two large deletion at the 3'-end, and one relatively large insertion of pncA. On the other hand, many of the 73 variants were lacking in the automated reporting systems except by TGS-TB; of these variants, CASTB detected only 20. By applying the 'non-wild type sequence' approach for predicting PZA resistance, accuracy of the results significantly improved compared with that of the automated results obtained by each tool.ConclusionUsers can obtain more accurate predictions for PZA resistance than previously reported by manually checking the results and applying the 'non-wild type sequence' approach.

A novel species of carbapenemase-producing Enterobacteriaceae (CPE) was isolated from a patient diagnosed with sigmoid colon diverticulitis. At first, laboratory testing suggested it was Klebsiella oxytoca or Pantoea sp.; however, a complete genome sequence of the isolate, MRY16-398, revealed that it could be novel species, most similar to [Kluyvera] intestini, of which taxonomic nomenclature is still under discussion. Orthologous conserved gene analysis among 42 related bacterial strains indicated that MRY16-398 was classified as the newly proposed genus Metakosakonia. Further, MRY16-398 was found to harbor the blaIMP-6 gene-positive class 1 integron (In722) in plasmid pMRY16-398_2 (IncN replicon, 47.4 kb in size). This finding implies that rare and opportunistic bacteria could be potential infectious agents. In conclusion, our results highlight the need for continuous monitoring for CPE even in nonpathogenic bacteria in the nosocomial environment.

A 60-year-old woman experienced fever, headache, rash, and altered vision after returning to Japan from India. Testing detected elevated antibody titers to spotted fever group rickettsia; PCR on blood yielded positive results for the rickettsial outer membrane protein A gene. We isolated a unique rickettsial agent and performed a full-genome analysis.

Human metapneumovirus (HMPV) is a major etiological agent of acute respiratory infections in humans. HMPV has been circulating worldwide for more than six decades and is currently divided into five agreed-upon subtypes: A1, A2a, A2b, B1, and B2. Recently, the novel HMPV subtypes A2c, A2b1, and A2b2 have been proposed. However, the phylogenetic and evolutionary relationships between these recently proposed HMPV subtypes are unclear. Here, we report a genome-wide phylogenetic and evolutionary analysis of 161 HMPV strains, including unique HMPV subtype A2b strains with a 180- or 111-nucleotide duplication in the G gene (nt-dup). Our data demonstrate that the HMPV A2b subtype contains two distinct subtypes, A2b1 and A2b2, and that the HMPV subtypes A2c and A2b2 may be different names for the same subtype. HMPV A2b strains with a nt-dup also belong to subtype A2b2. Molecular evolutionary analyses indicate that subtypes A2b1 and A2b2 diverged from subtype A2b around a decade after the subtype A2 was divided into the subtypes A2a and A2b. These data support the A2b1 and A2b2 subtypes proposed in 2012 and are essential for the unified classification of HMPV subtype A2 strains, which is important for future HMPV surveillance and epidemiological studies.

ABSTRACT Antimicrobial resistance genes (ARGs) and the bacteria that harbor them are widely distributed in the environment, especially in surface water, sewage treatment plant effluent, soil, and animal waste. In this study, we isolated a KPC-2-producing Klebsiella pneumoniae strain (GSU10-3) from a sampling site in Tokyo Bay, Japan, near a wastewater treatment plant (WWTP) and determined its complete genome sequence. Strain GSU10-3 is resistant to most β-lactam antibiotics and other antimicrobial agents (quinolones and aminoglycosides). This strain is classified as sequence type 11 (ST11), and a core genome phylogenetic analysis indicated that strain GSU10-3 is closely related to KPC-2-positive Chinese clinical isolates from 2011 to 2017 and is clearly distinct from strains isolated from the European Union (EU), United States, and other Asian countries. Strain GSU10-3 harbors four plasmids, including a blaKPC-2-positive plasmid, pGSU10-3-3 (66.2 kb), which is smaller than other blaKPC-2-positive plasmids and notably carries dual replicons (IncFII [pHN7A8] and IncN). Such downsizing and the presence of dual replicons may promote its maintenance and stable replication, contributing to its broad host range with low fitness costs. A second plasmid, pGSU10-3-1 (159.0 kb), an IncA/C2 replicon, carries a class 1 integron (containing intI1, dfrA12, aadA2, qacEΔ1, and sul1) with a high degree of similarity to a broad-host-range plasmid present in the family Enterobacteriaceae. The plasmid pGSU10-3-2 (134.8 kb), an IncFII(K) replicon, carries the IS26-mediated ARGs [aac(6')Ib-cr , blaOXA-1, catB4 (truncated), and aac(3)-IId], tet(A), and a copper/arsenate resistance locus. GSU10-3 is the first nonclinical KPC-2-producing environmental Enterobacteriaceae isolate from Japan for which the whole genome has been sequenced. IMPORTANCE We isolated and determined the complete genome sequence of a KPC-2-producing K. pneumoniae strain from a sampling site in Tokyo Bay, Japan, near a wastewater treatment plant (WWTP). In Japan, the KPC type has been very rarely detected, while IMP is the most predominant type of carbapenemase in clinical carbapenemase-producing Enterobacteriaceae (CPE) isolates. Although laboratory testing thus far suggested that Japan may be virtually free of KPC-producing Enterobacteriaceae, we have detected it from effluent from a WWTP. Antimicrobial resistance (AMR) monitoring of WWTP effluent may contribute to the early detection of future AMR bacterial dissemination in clinical settings and communities; indeed, it will help illuminate the whole picture in which environmental contamination through WWTP effluent plays a part.

Linezolid (LZD) has become one of the most important antimicrobial agents for infections caused by gram-positive bacteria, including those caused by Enterococcus species. LZD-resistant (LR) genetic features include mutations in 23S rRNA/ribosomal proteins, a plasmid-borne 23S rRNA methyltransferase gene cfr, and ribosomal protection genes (optrA and poxtA). Recently, a cfr gene variant, cfr(B), was identified in a Tn6218-like transposon (Tn) in a Clostridioides difficile isolate. Here, we isolated an LR Enterococcus faecalis clinical isolate, KUB3006, from a urine specimen of a patient with urinary tract infection during hospitalization in 2017. Comparative and whole-genome analyses were performed to characterize the genetic features and overall antimicrobial resistance genes in E. faecalis isolate KUB3006. Complete genome sequencing of KUB3006 revealed that it carried cfr(B) on a chromosomal Tn6218-like element. Surprisingly, this Tn6218-like element was almost (99%) identical to that of C. difficile Ox3196, which was isolated from a human in the UK in 2012, and to that of Enterococcus faecium 5_Efcm_HA-NL, which was isolated from a human in the Netherlands in 2012. An additional oxazolidinone and phenicol resistance gene, optrA, was also identified on a plasmid. KUB3006 is sequence type (ST) 729, suggesting that it is a minor ST that has not been reported previously and is unlikely to be a high-risk E. faecalis lineage. In summary, LR E. faecalis KUB3006 possesses a notable Tn6218-like-borne cfr(B) and a plasmid-borne optrA. This finding raises further concerns regarding the potential declining effectiveness of LZD treatment in the future.

[This corrects the article DOI: 10.1371/journal.ppat.1006882.].

The disease state of amebiasis, caused by Entamoeba histolytica, varies from asymptomatic to severe manifestations that include dysentery and extraintestinal abscesses. The virulence factors of the pathogen, and host defense mechanisms, contribute to the outcomes of infection; however, the underlying genetic factors, which affect clinical outcomes, remain to be fully elucidated. To identify these genetic factors in E. histolytica, we used Illumina next-generation sequencing to conduct a comparative genomic analysis of two clinical isolates obtained from diarrheal and asymptomatic patients (strains KU50 and KU27, respectively). By mapping KU50 and KU27 reads to the genome of a reference HM-1:IMSS strain, we identified two genes (EHI_089440 and EHI_176590) that were absent in strain KU27. In KU27, a single AIG1 (avrRpt2-induced gene 1) family gene (EHI_176590) was found to be deleted, from a tandem array of three AIG1 genes, by homologous recombination between the two flanking genes. Overexpression of the EHI_176590 gene, in strain HM-1:IMSS cl6, resulted in increased formation of cell-surface protrusions and enhanced adhesion to human erythrocytes. The EHI_176590 gene was detected by PCR in 56% of stool samples from symptomatic patients infected with E. histolytica, but only in 15% of stool samples from asymptomatic individuals. This suggests that the presence of the EHI_176590 gene is correlated with the outcomes of infection. Taken together, these data strongly indicate that the AIG1 family protein plays a pivotal role in E. histolytica virulence via regulation of host cell adhesion. Our in-vivo experiments, using a hamster liver abscess model, showed that overexpression or gene silencing of EHI_176590 reduced and increased liver abscess formation, respectively. This suggests that the AIG1 genes may have contrasting roles in virulence depending on the genetic background of the parasite and host environment.

Carbapenemase-producing Enterobacteriaceae (CPE) are a global concern because these bacteria are resistant to almost all β-lactams. Horizontal interspecies gene transfer via plasmid conjugation has increased the global dissemination of CPE. Recently, an Enterobacteriaceae strain positive for carbapenemase gene but showing a carbapenem-susceptible phenotype was identified, suggesting that these susceptible strains may be challenging to detect solely via antimicrobial susceptibility tests without molecular analysis. Here, we isolated a blaIMP-6 carbapenemase-gene positive but imipenem- and meropenem-susceptible Escherichia coli (ISMS-E) strain A56-1S (imipenem and meropenem minimum inhibitory concentration, ≤ 0.125 mg/L), from a human urine specimen in Japan. A56-1S was carbapenemase negative by the Carba NP test, suggesting that IMP-6 production was low or undetectable. Thus, to characterize the mechanism of this phenotype, a meropenem-resistant E. coli A56-1R strain was obtained using meropenem-selection. A56-1R was positive for carbapenemase production by the Carba NP test, and blaIMP-6 transcription in A56-1R was 53-fold higher than in A56-1S, indicating that blaIMP-6 in A56-1S is negatively regulated at the transcriptional level. Comparative genomic analysis between the two strains revealed that the alleviation of restriction of DNA (ardK) gene encoding a putative transcription factor is disrupted by the IS26 insertion in A56-1R. A cotransformation assay of ardK and the regulatory element upstream of blaIMP-6 showed repression of blaIMP-6 transcription, indicating that ArdK negatively modulates blaIMP-6 transcription. ArdK binding and affinity assays demonstrated that ArdK directly binds to the regulatory element upstream of blaIMP-6 with dissociation constant values comparable to those of general transcription factors. The IMP-6 carbapenemase showed low hydrolytic activity against imipenem, resulting in an imipenem-susceptible and meropenem-resistant (ISMR) phenotype (previously reported as a stealth phenotype). However, the low expression of IMP-6 in the A56-1S strain could be a typical characteristic of ISMS-E due to gene repression, indicating that conventional antimicrobial susceptibility tests might be unable to detect such strains even when using both imipenem and meropenem. Bacteria that exhibit the ISMS phenotype could play a potential role as undetectable reservoirs and might facilitate gene transfer to other organisms while avoiding detection.

In the 2016/2017 winter season in Japan, HuNoV GII.P16-GII.2 strains (2016 strains) emerged and caused large outbreaks of acute gastroenteritis. To better understand the outbreaks, we examined the molecular evolution of the VP1 gene and RdRp region in 2016 strains from patients by studying their time-scale evolutionary phylogeny, positive/negative selection, conformational epitopes, and phylodynamics. The time-scale phylogeny suggested that the common ancestors of the 2016 strains VP1 gene and RdRp region diverged in 2006 and 1999, respectively, and that the 2016 strain was the progeny of a pre-2016 GII.2. The evolutionary rates of the VP1 gene and RdRp region were around 10-3 substitutions/site/year. Amino acid substitutions (position 341) in an epitope in the P2 domain of 2016 strains were not found in pre-2016 GII.2 strains. Bayesian skyline plot analyses showed that the effective population size of the VP1 gene in GII.2 strains was almost constant for those 50 years, although the number of patients with NoV GII.2 increased in 2016. The 2016 strain may be involved in future outbreaks in Japan and elsewhere.

Human norovirus (HuNoV) is a leading cause of viral gastroenteritis worldwide, of which GII.4 is the most predominant genotype. Unlike other genotypes, GII.4 has created various variants that escaped from previously acquired immunity of the host and caused repeated epidemics. However, the molecular evolutionary differences among all GII.4 variants, including recently discovered strains, have not been elucidated. Thus, we conducted a series of bioinformatic analyses using numerous, globally collected, full-length GII.4 major capsid (VP1) gene sequences (466 strains) to compare the evolutionary patterns among GII.4 variants. The time-scaled phylogenetic tree constructed using the Bayesian Markov chain Monte Carlo (MCMC) method showed that the common ancestor of the GII.4 VP1 gene diverged from GII.20 in 1840. The GII.4 genotype emerged in 1932, and then formed seven clusters including 14 known variants after 1980. The evolutionary rate of GII.4 strains was estimated to be 7.68 × 10−3 substitutions/site/year. The evolutionary rates probably differed among variants as well as domains [protruding 1 (P1), shell, and P2 domains]. The Osaka 2007 variant strains probably contained more nucleotide substitutions than any other variant. Few conformational epitopes were located in the shell and P1 domains, although most were contained in the P2 domain, which, as previously established, is associated with attachment to host factors and antigenicity. We found that positive selection sites for the whole GII.4 genotype existed in the shell and P1 domains, while Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants were under positive selection in the P2 domain. Amino acid substitutions overlapped with putative epitopes or were located around the epitopes in the P2 domain. The effective population sizes of the present strains increased stepwise for Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants. These results suggest that HuNoV GII.4 rapidly evolved in a few decades, created various variants, and altered its evolutionary rate and antigenicity.