Characterization of the gut microbiota of Kawasaki disease patients by metagenomic analysis

Kawasaki disease (KD) is an acute febrile illness of early childhood. Previous reports have suggested that genetic disease susceptibility factors, together with a triggering infectious agent, could be involved in KD pathogenesis; however, the precise etiology of this disease remains unknown. Additionally, previous culture-based studies have suggested a possible role of intestinal microbiota in KD pathogenesis. In this study, we performed metagenomic analysis to comprehensively assess the longitudinal variation in the intestinal microbiota of twenty-eight KD patients. Several notable bacterial genera were commonly extracted during the acute phase, whereas a relative increase in the number of Ruminococcus bacteria was observed during the non-acute phase of KD. The metagenomic analysis results based on bacterial species classification suggested that the number of sequencing reads with similarity to five Streptococcus spp. (S. pneumonia, pseudopneumoniae, oralis, gordonii, and sanguinis), in addition to patient-derived Streptococcus isolates, markedly increased during the acute phase in most patients. Streptococci include a variety of pathogenic bacteria and probiotic bacteria that promote human health; therefore, this further species discrimination could comprehensively illuminate the KD-associated microbiota. The findings of this study suggest that KD-related Streptococci might be involved in the pathogenesis of this disease

Corynebacterium ulcerans 0102 carries the gene encoding diphtheria toxin on a prophage different from the C. diphtheriae NCTC 13129 prophage

BACKGROUND: Corynebacterium ulcerans can cause a diphtheria-like illness, especially when the bacterium is lysogenized with a tox gene-carrying bacteriophage that produces diphtheria toxin. Acquisition of toxigenicity upon phage lysogenization is a common feature of C. ulcerans and C. diphtheriae. However, because of a lack of C. ulcerans genome information, a detailed comparison of prophages has not been possible between these two clinically important and closely related bacterial species. RESULTS: We determined the whole genome sequence of the toxigenic C. ulcerans 0102 isolated in Japan. The genomic sequence showed a striking similarity with that of Corynebacterium pseudotuberculosis and, to a lesser extent, with that of C. diphtheriae. The 0102 genome contained three distinct prophages. One of these, ΦCULC0102-I, was a tox-positive prophage containing genes in the same structural order as for tox-positive C. diphtheriae prophages. However, the primary structures of the individual genes involved in the phage machinery showed little homology between the two counterparts. CONCLUSION: Taken together, these results suggest that the tox-positive prophage in this strain of C. ulcerans has a distinct origin from that of C. diphtheriae NCTC 13129

Molecular characterization of the sequences of the 16S-23S rDNA internal spacer region (ISR) from isolates of Taylorella asinigenitalis

<p>Abstract</p> <p>Background</p> <p>Sequence information on the 16S-23S rDNA internal spacer region (ISR) exhibits a large degree of sequence and length variation at both the genus and species levels. A primer pair for the amplification of 16S-23S rDNA ISR generated three amplicons for each of isolates of <it>Taylorella asinigenitalis </it>(UCD-1^T, UK-1 and UK-2).</p> <p>Findings</p> <p>Following TA cloning and sequencing, the three isolates of <it>T. asinigenitalis </it>were demonstrated to possess three ISR units of different lengths. Although the three corresponding ISRs (A, B and C) were identified to be identical to each other (UK-1 and UK-2 isolates), the ISRs shared approximately 95.3–98.9% nucleotide sequence similarities between the UCD-1^Tand UK-1/-2 isolates. A typical order of two intercistronic tRNA genes (5'-tRNA^Ile-tRNA^Ala-3') with the different nucleotide spacers [44 through 51 base pairs (bp)] in length was identified among the isolates. The consensus sequences of the antiterminators of <b>boxB </b>and <b>boxA </b>were also identified in all ISRs. Thus, three ISRs were identified for each isolate, and therefore, at least three distinctly different ribosomal RNA operons were suggested to occur in the genome of <it>T. asinigenitalis</it>. This was also confirmed by Southern hybridization procedure.</p> <p>Conclusion</p> <p>The present study represents a dendrogram constructed based on the nucleotide sequence data of 16S-23S rDNA ISR for <it>T. asinigenitalis</it>, which may aid in the phylogenetic positioning of <it>T. asinigenitalis </it>within the genus <it>Taylorella</it>, and in the molecular discrimination of <it>T. asinigenitalis</it>.</p

Complete Genome Sequence of blaIMP–6-Positive Metakosakonia sp. MRY16-398 Isolate From the Ascites of a Diverticulitis Patient

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

Complete Genome Sequence and Characterization of Linezolid-Resistant Enterococcus faecalis Clinical Isolate KUB3006 Carrying a cfr(B)-Transposon on Its Chromosome and optrA-Plasmid

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

Characterization of Quasispecies of Pandemic 2009 Influenza A Virus (A/H1N1/2009) by De Novo Sequencing Using a Next-Generation DNA Sequencer

Pandemic 2009 influenza A virus (A/H1N1/2009) has emerged globally. In this study, we performed a comprehensive detection of potential pathogens by de novo sequencing using a next-generation DNA sequencer on total RNAs extracted from an autopsy lung of a patient who died of viral pneumonia with A/H1N1/2009. Among a total of 9.4×106 40-mer short reads, more than 98% appeared to be human, while 0.85% were identified as A/H1N1/2009 (A/Nagano/RC1-L/2009(H1N1)). Suspected bacterial reads such as Streptococcus pneumoniae and other oral bacteria flora were very low at 0.005%, and a significant bacterial infection was not histologically observed. De novo assembly and read mapping analysis of A/Nagano/RC1-L/2009(H1N1) showed more than ×200 coverage on average, and revealed nucleotide heterogeneity on hemagglutinin as quasispecies, specifically at two amino acids (Gly172Glu and Gly239Asn of HA) located on the Sa and Ca2 antigenic sites, respectively. Gly239 and Asn239 on antigenic site Ca2 appeared to be minor amino acids compared with the highly distributed Asp239 in H1N1 HAs. This study demonstrated that de novo sequencing can comprehensively detect pathogens, and such in-depth investigation facilitates the identification of influenza A viral heterogeneity. To better characterize the A/H1N1/2009 virus, unbiased comprehensive techniques will be indispensable for the primary investigations of emerging infectious diseases

Comparative genome analysis of VSP-II and SNPs reveals heterogenic variation in contemporary strains of Vibrio cholerae O1 isolated from cholera patients in Kolkata, India

Cholera is an acute diarrheal disease and a major public health problem in many developing countries in Asia, Africa, and Latin America. Since the Bay of Bengal is considered the epicenter for the seventh cholera pandemic, it is important to understand the genetic dynamism of Vibrio cholerae from Kolkata, as a representative of the Bengal region. We analyzed whole genome sequence data of V. cholerae O1 isolated from cholera patients in Kolkata, India, from 2007 to 2014 and identified the heterogeneous genomic region in these strains. In addition, we carried out a phylogenetic analysis based on the whole genome single nucleotide polymorphisms to determine the genetic lineage of strains in Kolkata. This analysis revealed the heterogeneity of the Vibrio seventh pandemic island (VSP)-II in Kolkata strains. The ctxB genotype was also heterogeneous and was highly related to VSP-II types. In addition, phylogenetic analysis revealed the shifts in predominant strains in Kolkata. Two distinct lineages, 1 and 2, were found between 2007 and 2010. However, the proportion changed markedly in 2010 and lineage 2 strains were predominant thereafter. Lineage 2 can be divided into four sublineages, I, II, III and IV. The results of this study indicate that lineages 1 and 2-I were concurrently prevalent between 2007 and 2009, and lineage 2-III observed in 2010, followed by the predominance of lineage 2-IV in 2011 and continued until 2014. Our findings demonstrate that the epidemic of cholera in Kolkata was caused by several distinct strains that have been constantly changing within the genetic lineages of V. cholerae O1 in recent years

The genome sequence of Streptomyces rochei 7434AN4, which carries a linear chromosome and three characteristic linear plasmids

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.This work was supported by Grants-in-Aid for Scientific Research on Innovative Areas (23108515, 25108718 and 17H05446 to K.A.) from Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT), Grants-in-Aid for Scientific Research (B) (16H04917 to K.A.) from the Japan Society for the Promotion of Science (JSPS), and the Sasakawa Scientific Research Grant from the Japan Science Society to Y.N. This work was partly supported by a JSPS A3 Foresight Program. A.A.F. and R.M. were supported by the Indonesia Endowment Fund for Education (LPDP). Sequencing analysis using an Illumina GAIIx sequencer was supported by the Grant in Aid for Scientific Research on Innovative Areas (22108010 to J.I.) from MEXT

Complete Sequencing of the blaNDM-1-Positive IncA/C Plasmid from Escherichia coli ST38 Isolate Suggests a Possible Origin from Plant Pathogens

The complete sequence of the plasmid pNDM-1_Dok01 carrying New Delhi metallo-β-lactamase (NDM-1) was determined by whole genome shotgun sequencing using Escherichia coli strain NDM-1_Dok01 (multilocus sequence typing type: ST38) and the transconjugant E. coli DH10B. The plasmid is an IncA/C incompatibility type composed of 225 predicted coding sequences in 195.5 kb and partially shares a sequence with blaCMY-2-positive IncA/C plasmids such as E. coli AR060302 pAR060302 (166.5 kb) and Salmonella enterica serovar Newport pSN254 (176.4 kb). The blaNDM-1 gene in pNDM-1_Dok01 is terminally flanked by two IS903 elements that are distinct from those of the other characterized NDM-1 plasmids, suggesting that the blaNDM-1 gene has been broadly transposed, together with various mobile elements, as a cassette gene. The chaperonin groES and groEL genes were identified in the blaNDM-1-related composite transposon, and phylogenetic analysis and guanine-cytosine content (GC) percentage showed similarities to the homologs of plant pathogens such as Pseudoxanthomonas and Xanthomonas spp., implying that plant pathogens are the potential source of the blaNDM-1 gene. The complete sequence of pNDM-1_Dok01 suggests that the blaNDM-1 gene was acquired by a novel composite transposon on an extensively disseminated IncA/C plasmid and transferred to the E. coli ST38 isolate

Characterization of Antimicrobial Resistance Dissemination across Plasmid Communities Classified by Network Analysis

The global clustering of gene families through network analysis has been demonstrated in whole genome, plasmid, and microbiome analyses. In this study, we carried out a plasmidome network analysis of all available complete bacterial plasmids to determine plasmid associations. A blastp clustering search at 100% aa identity cut-off and sharing at least one gene between plasmids, followed by a multilevel community network analysis revealed that a surprisingly large number of the plasmids were connected by one largest connected component (LCC), with dozens of community sub-groupings. The LCC consisted mainly of Bacilli and Gammaproteobacteria plasmids. Intriguingly, horizontal gene transfer (HGT) was noted between different phyla (i.e., Staphylococcus and Pasteurellaceae), suggesting that Pasteurellaceae can acquire antimicrobial resistance (AMR) genes from closely contacting Staphylococcus spp., which produce the external supplement of V-factor (NAD). Such community network analysis facilitate displaying possible recent HGTs like a class 1 integron, str and tet resistance markers between communities. Furthermore, the distribution of the Inc replicon type and AMR genes, such as the extended-spectrum ß-lactamase (ESBL) CTX-M or the carbapenemases KPC NDM-1, implies that such genes generally circulate within limited communities belonging to typical bacterial genera. Thus, plasmidome network analysis provides a remarkable discriminatory power for plasmid-related HGT and evolution