Yersinia is a Gram-negative bacteria that includes serious pathogens such as the Yersinia pestis, which causes plague, Yersinia pseudotuberculosis, Yersinia enterocolitica. The remaining species are generally considered non-pathogenic to humans, although there is evidence that at least some of these species can cause occasional infections using distinct mechanisms from the more pathogenic species. With the advances in sequencing technologies, many genomes of Yersinia have been sequenced. However, there is currently no specialized platform to hold the rapidly-growing Yersinia genomic data and to provide analysis tools particularly for comparative analyses, which are required to provide improved insights into their biology, evolution and pathogenicity.
Twenty percent of genes that encode for hypothetical proteins from Klebsiella pneumoniae MGH78578 were identified, leading to KPN00728 and KPN00729 after bioinformatics analysis. Both open reading frames showed high sequence homology to Succinate dehydrogenase Chain C (SdhC) and D (SdhD) from Escherichia coli. Recently, KPN00729 was assigned as SdhD. KPN00728 thus remains of particular interest as no annotated genes from the complete genome sequence encode for SdhC. We discovered KPN00728 has a missing region with conserved residues important for ubiquinone (UQ) and heme group binding. Structure and function prediction of KPN00728 coupled with secondary structure analysis and transmembrane topology showed KPN00728 adopts SDH-(subunit C)-like structure. To further probe its functionality, UQ was docked on the built model (consisting KPN00728 and KPN00729) and formation of hydrogen bonds between UQ and Ser27, Arg31 (KPN00728) and Tyr84 (KPN00729) further reinforces and supports that KPN00728 is indeed SDH. This is the first report on the structural and function prediction of KPN00728 of K. pneumoniae MGH78578 as SdhC.
Pseudomonas aeruginosa has been implicated in a wide range of post-operation wound and lung infections. A wide range of acquired resistance and virulence markers indicate surviving strategy of P. aeruginosa. Complete-genome analysis has been identified as efficient approach towards understanding the pathogenicity of this organism. This study was designed to sequence the entire genome of P. aeruginosa UY1PSABAL and UY1PSABAL2; determine drug-resistance profiles and virulence factors of the isolates; assess factors that contribute toward stability of the genomes; and thereafter determine evolutionary relationships between the strains and other isolates from similar sources. The genomes of the MDR P. aeruginosa UY1PSABAL and UY1PSABAL2 were sequenced on the Illumina Miseq platform. The raw sequenced reads were assessed for quality using FastQC v.0.11.5 and filtered for low quality reads and adapter regions using Trimmomatic v.0.36. The de novo genome assembly was made with SPAdes v.3.13 and annotated using Prokka v.2.1.1 annotation pipeline; Rapid Annotation using Subsytems Technology (RAST) server v.2.0; and PATRIC annotation tool v.3.6.2. Antimicrobial resistance genes and virulence determinants were searched through the functional annotation data generated from Prokka, RAST and PATRIC annotation pipelines; In addition to ResFinder and Comprehensive Antibiotic Resistance Database (CARD) which were employed to determine resistance genes. The PHAge Search Tool Enhanced Release (PHASTER) web server was used for the rapid identification and annotation of prophage sequences within bacterial genome. Predictive secondary metabolites were identified with AntiSMASH v.5.0. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and cas genes regions were also investigated with the CRISPRone and CRISPRFinder server. The genome sizes of 7.0 and 6.4 Mb were determined for UY1PSABAL and UY1PSABAL2 strains with G+C contents of 66.1% and 66.48% respectively. β-lactamines resistance genes blaPAO, aminoglycoside phosphorylating enzymes genes aph(3')-IIb, fosfomycine resistance gene fosA, vancomycin vanW and tetracycline tetA were among identified resistance genes harboured in both isolates. UY1PSABAL bore additional aph(6)-Id, aph(3'')-Ib, ciprofloxacin-modifying enzyme crpP and ribosomal methylation enzyme rmtB. Both isolates were found harbouring virulence markers such as flagella and type IV pili; and also present various type III secretion systems such as exoA, exoS, exoU, exoT. Secondary metabolites such as pyochelin and pyoverdine with iron uptake activity were found within the genomes as well as quorum-sensing systems, and various fragments for prophages and insertion sequences. Only the UY1PSABAL2 contains CRISPR-Cas system. The phylogeny revealed a very close evolutionary relationship between UY1PSABAL and the similar strain isolated from Malaysia; the same trend was observed between UY1PSABAL2 and the strain from Chinese origin. Complete analyses of the entire genomes provide a wide range of information towards understanding pathogenicity of the pathogens in question.
Acinetobacter sp. strain GG2 is a quorum-sensing and quorum-quenching bacterium isolated from the ginger rhizosphere. It degrades a broad range of N-acylhomoserine lactone molecules via lactonase. The genome sequence of strain GG2 may provide insights on the regulation of quorum-sensing and quorum-quenching mechanisms in this bacterium.
Enterobacter sp. strain SST3 is an endophytic bacterium isolated from Saccharum spp. Here we present its annotated draft genome that may shed light on its role as a bacterial endophyte of sugarcane. To our knowledge, this is the first genome announcement of a sugarcane-associated bacterium from the genus Enterobacter.
Burkholderia pseudomallei, a soil-dwelling Gram-negative bacterium, is the causative agent of the endemic tropical disease melioidosis. Clinical manifestations of B. pseudomallei infection range from acute or chronic localized infection in a single organ to fulminant septicaemia in multiple organs. The diverse clinical manifestations are attributed to various factors, including the genome plasticity across B. pseudomallei strains. We previously characterized B. pseudomallei strains isolated in Malaysia and noted different levels of virulence in model hosts. We hypothesized that the difference in virulence might be a result of variance at the genome level. In this study, we sequenced and assembled four Malaysian clinical B. pseudomallei isolates, UKMR15, UKMPMC2000, UKMD286 and UKMH10. Phylogenomic analysis showed that Malaysian subclades emerged from the Asian subclade, suggesting that the Malaysian strains originated from the Asian region. Interestingly, the low-virulence strain, UKMH10, was the most distantly related compared to the other Malaysian isolates. Genomic island (GI) prediction analysis identified a new island of 23 kb, GI9c, which is present in B. pseudomallei and Burkholderia mallei, but not Burkholderia thailandensis. Genes encoding known B. pseudomallei virulence factors were present across all four genomes, but comparative analysis of the total gene content across the Malaysian strains identified 104 genes that are absent in UKMH10. We propose that these genes may encode novel virulence factors, which may explain the reduced virulence of this strain. Further investigation on the identity and role of these 104 proteins may aid in understanding B. pseudomallei pathogenicity to guide the design of new therapeutics for treating melioidosis.
Acinetobacter baumannii has emerged as an important nosocomial pathogen owing to its increasing resistance to most, if not all, antibiotics in clinical use. We recently reported the occurrence of extensively drug-resistant (XDR) A. baumannii isolates in a Malaysian tertiary hospital. The genome of one of these XDR isolates (A. baumannii AC12) was completely sequenced and comparative genome analyses were performed to elucidate the genetic basis of its antimicrobial resistance. The A. baumannii AC12 genome consists of a 3.8 Mbp circular chromosome and an 8731 bp cryptic plasmid, pAC12. It belongs to the ST195 lineage and is most closely related to A. baumannii BJAB0715 as well as other strains of the international clone III (IC-III) group. Two antibiotic resistance islands (RIs), designated AC12-RI1 and AC12-RI2, were found in the AC12 chromosome along with a 7 kb Tn1548::armA island conferring resistance to aminoglycosides and macrolides. The 22.8 kb AC12-RI1 interrupts the comM gene and harbours the carbapenem resistance gene blaOXA-23 flanked by ISAba1 within a Tn2006-like structure. AC12-RI1 also harbours resistance determinants for aminoglycosides, tetracyclines and sulphonamides. The 10.3 kb IS26-flanked AC12-RI2 is a derivative of AbGRI2-1, containing aphA1b and blaTEM genes (conferring aminoglycoside and β-lactam resistance, respectively). The presence of numerous genes mediating resistance to various antibiotics in novel RI structures as well as other genes encoding drug transporters and efflux pumps in A. baumannii AC12 most likely contributed to its XDR characteristics.
Mycobacterium abscessus is an emerging clinical pathogen commonly associated with non-tuberculous mycobacterial infections. We report herein the draft genome of M. abscessus strain M156.
This is a report of an annotated genome sequence of Mycobacterium tuberculosis MTBR1/09. The organism was isolated from a sputum sample from a male patient in Malaysia.
Here, we report of the annotated genome sequence of Mycobacterium tuberculosis MTB221/11. The organism was isolated from the cerebrospinal fluid of a patient in Malaysia.
Unlike environmental P. koreensis isolated from soil, which has been studied extensively for its role in promoting plant growth, pathogenic P. koreensis isolated from fish has been rarely reported. Therefore, we investigated and isolated the possible pathogen that is responsible for the diseased state of Tor tambroides. Herein, we reported the morphological and biochemical characteristics, as well as whole-genome sequences of a newly identified P. koreensis strain. We assembled a high-quality draft genome of P. koreensis CM-01 with a contig N50 value of 233,601 bp and 99.5% BUSCO completeness. The genome assembly of P. koreensis CM-01 is consists of 6,171,880 bp with a G+C content of 60.5%. Annotation of the genome identified 5538 protein-coding genes, 3 rRNA genes, 54 tRNAs, and no plasmids were found. Besides these, 39 interspersed repeat and 141 tandem repeat sequences, 6 prophages, 51 genomic islands, 94 insertion sequences, 4 clustered regularly interspaced short palindromic repeats, 5 antibiotic-resistant genes, and 150 virulence genes were also predicted in the P. koreensis CM-01 genome. Culture-based approach showed that CM-01 strain exhibited resistance against ampicillin, aztreonam, clindamycin, and cefoxitin with a calculated multiple antibiotic resistance (MAR) index value of 0.4. In addition, the assembled CM-01 genome was successfully annotated against the Cluster of Orthologous Groups of proteins database, Gene Ontology database, and Kyoto Encyclopedia of Genes and Genome pathway database. A comparative analysis of CM-01 with three representative strains of P. koreensis revealed that 92% of orthologous clusters were conserved among these four genomes, and only the CM-01 strain possesses unique elements related to pathogenicity and virulence. This study provides fundamental phenotypic and genomic information for the newly identified P. koreensis strain.
Mycobacterium brisbanense is a member of Mycobacterium fortuitum third biovariant complex, which includes rapidly growing Mycobacterium spp. that normally inhabit soil, dust and water, and can sometimes cause respiratory tract infections in humans. We present the first whole-genome analysis of M. brisbanense UM_WWY which was isolated from a 70-year-old Malaysian patient. Molecular phylogenetic analyses confirmed the identification of this strain as M. brisbanense and showed that it has an unusually large genome compared with related mycobacteria. The large genome size of M. brisbanense UM_WWY (~7.7Mbp) is consistent with further findings that this strain has a highly variable genome structure that contains many putative horizontally transferred genomic islands and prophage. Comparative analysis showed that M. brisbanense UM_WWY is the only Mycobacterium species that possesses a complete set of genes encoding enzymes involved in the urea cycle, suggesting that this soil bacterium is able to synthesize urea for use as plant fertilizers. It is likely that M. brisbanense UM_WWY is adapted to live in soil as its primary habitat since the genome contains many genes associated with nitrogen metabolism. Nevertheless, a large number of predicted virulence genes were identified in M. brisbanense UM_WWY that are mostly shared with well-studied mycobacterial pathogens such as Mycobacterium tuberculosis and Mycobacterium abscessus. These findings are consistent with the role of M. brisbanense as an opportunistic pathogen of humans. The whole-genome study of UM_WWY has provided the basis for future work of M. brisbanense.
Bacillus sp. is a Gram-positive bacterium that is commonly found in seawater. In this study, the genome of marine Bacillus sp. strain G3(2015) was sequenced using MiSeq. The fosfomycin resistant gene fosB was identified upon bacterial genome annotation.
An integral approach to decoding both culturable and uncultured microorganisms' metabolic activity involves the whole genome sequencing (WGS) of individual/complex microbial communities. WGS of culturable microbes, amplicon sequencing, metagenomics, and single-cell genome analysis are selective techniques integrating genetic information and biochemical mechanisms. These approaches transform microbial biotechnology into a quick and high-throughput culture-independent evaluation and exploit pollutant-degrading microbes. They are windows into enzyme regulatory bioremediation pathways (i.e., dehalogenase) and the complete bioremediation process of organohalide pollutants. While the genome sequencing technique is gaining the scientific community's interest, it is still in its infancy in the field of pollutant bioremediation. The techniques are becoming increasingly helpful in unraveling and predicting the enzyme structure and explore metabolic and biodegradation capabilities.
Here we report the draft genomes and annotation of four N-acyl homoserine lactone (AHL)-producing members from the family Sphingomonadaceae. Comparative genomic analyses of 62 Sphingomonadaceae genomes were performed to gain insights into the distribution of the canonical luxI/R-type quorum sensing (QS) network within this family. Forty genomes contained at least one luxR homolog while the genome of Sphingobium yanoikuyae B1 contained seven Open Reading Frames (ORFs) that have significant homology to that of luxR. Thirty-three genomes contained at least one luxI homolog while the genomes of Sphingobium sp. SYK6, Sphingobium japonicum, and Sphingobium lactosutens contained four luxI. Using phylogenetic analysis, the sphingomonad LuxR homologs formed five distinct clades with two minor clades located near the plant associated bacteria (PAB) LuxR solo clade. This work for the first time shows that 13 Sphingobium and one Sphingomonas genome(s) contain three convergently oriented genes composed of two tandem luxR genes proximal to one luxI (luxR-luxR-luxI). Interestingly, luxI solos were identified in two Sphingobium species and may represent species that contribute to AHL-based QS system by contributing AHL molecules but are unable to perceive AHLs as signals. This work provides the most comprehensive description of the luxI/R circuitry and genome-based taxonomical description of the available sphingomonad genomes to date indicating that the presence of luxR solos and luxI solos are not an uncommon feature in members of the Sphingomonadaceae family.
To facilitate the ongoing research of Vibrio spp., a dedicated platform for the Vibrio research community is needed to host the fast-growing amount of genomic data and facilitate the analysis of these data. We present VibrioBase, a useful resource platform, providing all basic features of a sequence database with the addition of unique analysis tools which could be valuable for the Vibrio research community. VibrioBase currently houses a total of 252 Vibrio genomes developed in a user-friendly manner and useful to enable the analysis of these genomic data, particularly in the field of comparative genomics. Besides general data browsing features, VibrioBase offers analysis tools such as BLAST interfaces and JBrowse genome browser. Other important features of this platform include our newly developed in-house tools, the pairwise genome comparison (PGC) tool, and pathogenomics profiling tool (PathoProT). The PGC tool is useful in the identification and comparative analysis of two genomes, whereas PathoProT is designed for comparative pathogenomics analysis of Vibrio strains. Both of these tools will enable researchers with little experience in bioinformatics to get meaningful information from Vibrio genomes with ease. We have tested the validity and suitability of these tools and features for use in the next-generation database development.
Social bacteria use chemical communication to coordinate and synchronize gene expression via the quorum-sensing (QS) regulatory pathway. In Pectobacterium, a causative agent of the blackleg and soft-rot diseases on potato plants and tubers, expression of the virulence factors is collectively controlled by the QS-signals N-acylhomoserine lactones (NAHLs). Several soil bacteria, such as the actinobacterium Rhodococcus erythropolis, are able to degrade NAHLs, hence quench the chemical communication and virulence of Pectobacterium. Here, next-generation sequencing was used to investigate structural and functional genomics of the NAHL-degrading R. erythropolis strain R138. The R. erythropolis R138 genome (6.7 Mbp) contained a single circular chromosome, one linear (250 kbp) and one circular (84 kbp) plasmid. Growth of R. erythropolis and P. atrosepticum was not altered in mixed-cultures as compared with monocultures on potato tuber slices. HiSeq-transcriptomics revealed that no R. erythropolis genes were differentially expressed when R. erythropolis was cultivated in the presence vs absence of the avirulent P. atrosepticum mutant expI, which is defective for QS-signal synthesis. By contrast 50 genes (<1% of the R. erythropolis genome) were differentially expressed when R. erythropolis was cultivated in the presence vs absence of the NAHL-producing virulent P. atrosepticum. Among them, quantitative real-time reverse-transcriptase-PCR confirmed that the expression of some alkyl-sulfatase genes decreased in the presence of a virulent P. atrosepticum, as well as deprivation of organic sulfur such as methionine, which is a key precursor in the synthesis of NAHL by P. atrosepticum.
Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.
The application of modern research tools has broadened our understanding of the chlamydiae and their role in disease. Chlamydial genome analysis showed the presence of genes for ATP and peptidoglycan synthesis, contradicting the common belief that chlamydiae lack the ability to produce these compounds. Phylogenetic tree analysis suggests that chlamydiae could have evolved from an intracellular existence in amoebae. Newly discovered obligate intracellular organisms with chlamydia-like life-cycles have been classified as chlamydiae by rRNA homology with existing chlamydial species. A proposed new classification adds three new families to the order Chlamydiales as well as creates two genera and nine species within the family Chlamydiaceae. Chlamydiae are incriminated in an increasingly large spectrum of diseases both in humans and in animals. The emergence of multi-drug resistant C. trachomatis strains forewarns therapeutic problems with this organism. While C. pneumoniae remains a significant respiratory pathogen, the role it plays in the pathogenesis of atherosclerosis and ischaemic heart disease awaits definition.
The ST131 multilocus sequence type (MLST) of Escherichia coli is a globally successful pathogen whose dissemination is increasing rates of antibiotic resistance. Numerous global surveys have demonstrated the pervasiveness of this clone; in some regions ST131 accounts for up to 30% of all E. coli isolates. However, many regions are underrepresented in these published surveys, including Africa, South America, and Asia. We collected consecutive bloodstream E. coli isolates from three countries in Southeast Asia; ST131 was the most common MLST type. As in other studies, the C2/H30Rx clade accounted for the majority of ST131 strains. Clinical risk factors were similar to other reported studies. However, we found that nearly all of the C2 strains in this study were closely related, forming what we denote the SEA-C2 clone. The SEA-C2 clone is enriched for strains from Asia, particularly Southeast Asia and Singapore. The SEA-C2 clone accounts for all of the excess resistance and virulence of ST131 relative to non-ST131 E. coli. The SEA-C2 strains appear to be locally circulating and dominant in Southeast Asia, despite the intuition that high international connectivity and travel would enable frequent opportunities for other strains to establish themselves.