Search

The genotoxin colibactin is a secondary metabolite produced by the polyketide synthase (pks) island harbored by extraintestinal pathogenic E. coli (ExPEC) and other members of the Enterobacteriaceae that has been increasingly reported to have critical implications in human health. The present study entails a high-throughput whole-genome comparison and phylogenetic analysis of such pathogenic E. coli isolates to gain insights into the patterns of distribution, horizontal transmission, and evolution of the island. For the current study, 23 pks-positive ExPEC genomes were newly sequenced, and their virulome and resistome profiles indicated a preponderance of virulence encoding genes and a reduced number of genes for antimicrobial resistance. In addition, 4,090 E. coli genomes from the public domain were also analyzed for large-scale screening for pks-positive genomes, out of which a total of 530 pks-positive genomes were studied to understand the subtype-based distribution pattern(s). The pks island showed a significant association with the B2 phylogroup (82.2%) and a high prevalence in sequence type 73 (ST73; n = 179) and ST95 (n = 110) and the O6:H1 (n = 110) serotype. Maximum-likelihood (ML) phylogeny of the core genome and intergenic regions (IGRs) of the ST95 model data set, which was selected because it had both pks-positive and pks-negative genomes, displayed clustering in relation to their carriage of the pks island. Prevalence patterns of genes encoding RM systems in the pks-positive and pks-negative genomes were also analyzed to determine their potential role in pks island acquisition and the maintenance capability of the genomes. Further, the maximum-likelihood phylogeny based on the core genome and pks island sequences from 247 genomes with an intact pks island demonstrated horizontal gene transfer of the island across sequence types and serotypes, with few exceptions. This study vitally contributes to understanding of the lineages and subtypes that have a higher propensity to harbor the pks island-encoded genotoxin with possible clinical implications. IMPORTANCE Extraintestinal pathologies caused by highly virulent strains of E. coli amount to clinical implications with high morbidity and mortality rates. Pathogenic E. coli strains are evolving with the horizontal acquisition of mobile genetic elements, including pathogenicity islands such as the pks island, which produces the genotoxin colibactin, resulting in severe clinical outcomes, including colorectal cancer progression. The current study encompasses high-throughput comparative genomics and phylogenetic analyses to address the questions pertaining to the acquisition and evolution pattern of the genomic island in different E. coli subtypes. It is crucial to gain insights into the distribution, transfer, and maintenance of pathogenic islands, as they harbor multiple virulence genes involved in pathogenesis and clinical implications of the infection.

Shigellosis, caused by Shigella boydii type 1, is understudied and underreported. For 3 years, GEMS study identified 5.4% of all Shigella as S. boydii. We showed the prevalent serotypes of S. boydii in Bangladesh and phage-based diagnosis of S. boydii type 1, a rapid and low-cost approach. Previously typed 793 clinical S. boydii strains were used for serotype distribution. Twenty-eight environmental water samples were collected for isolation of Shigella phages. Forty-eight serotypes of Shigella and other enteric bacteria were used for testing the susceptibility to phage MK-13. Electron microscopy, restriction enzyme analysis, whole genome sequencing (WGS), and annotation were performed for extensive characterization. S. boydii type 1 is the second most prevalent serotype among 20 serotypes of S. boydii in Bangladesh. We isolated a novel phage, MK-13, which specifically lyses S. boydii type 1, but doesn’t lyse other 47 serotypes of Shigella or other enteric bacteria tested. The phage belongs to the Myoviridae family and distinct from other phages indicated by electron microscopy and restriction enzyme analysis, respectively. MK-13 genome consists of 158 kbp of circularly permuted double-stranded DNA with G + C content of 49.45%, and encodes 211 open reading frames including four tRNA-coding regions. The genome has 98% identity with previously reported phage, ΦSboM-AG3, reported to have a broader host range infecting most of the S. boydii and other species of Shigella tested. To our knowledge, MK-13 is the first phage reported to be used as a diagnostic marker to detect S. boydii type 1, especially in remote settings with limited laboratory infrastructure.

Colibactin, a genotoxin, encoded by the pks pathogenicity island of Escherichia coli belonging to the B2 phylogroup has been increasingly reported as a critical determinant of bacterial pathogenicity. The present study was designed to detect the pks pathogenicity island in extraintestinal pathogenic E. coli (ExPEC) cultured from clinical sources in India. Of 462 pathogenic E. coli analyzed, the pks genomic island was detected in 35 (7.6 %) isolates, which predominantly belonged to pathogenic phylogroup B2 (97%), and harbored virulence genes such as fimH, sfaD/E, and usp. Biofilm formation assay revealed 21 of 35 pks-carrying isolates to be strong (SBF > 1.0), 10 strains to be moderate (SBF = 0.5-1.0), and 4 strains to be weak (SBF < 0.5) biofilm formers. All of the pks-carrying isolates proved resistant against bactericidal activity of human serum and were clonal, as confirmed by PCR based on the Enterobacterial repetitive intergenic consensus (ERIC-PCR). Assays carried out to detect antimicrobial susceptibility revealed 31% of the isolates to be multidrug resistant with 37% carrying ESBL and 25% blaCTX-M-15. The observed prevalence of multidrug resistant and colibactin producing characteristics among pathogenic E. coli belonging to phylogenetic group B2 underscores urgent need for broader surveillance in order to understand and prevent transmission of this extra-intestinal pathogen in community and hospital settings.

Pyrazinamide (PZA) is an important first-line anti-tuberculosis drug, resistance to which occurs primarily due to mutations in pncA (Rv2043c) that encodes the pyrazinamidase enzyme responsible for conversion of pro-drug PZA into its active form. Previous studies have reported numerous resistance-conferring mutations distributed across the entire length of pncA without any hotspot regions. As different lineages of Mycobacterium tuberculosis display a strong geographic association, we sought to understand whether the genetic background influenced the distribution of mutations in pncA. We analyzed the whole genome sequence data of 1,480 clinical isolates representing four major M. tuberculosis lineages to identify the distribution of mutations in the complete operon (Rv2044c-pncA-Rv2042c) and its upstream promoter region. We observed a non-overlapping pattern of mutations among various lineages and identified a lineage 3-specific frame-shift deletion in gene Rv2044c upstream of pncA that disrupted the stop codon and led to its fusion with pncA. This resulted in the addition of a novel domain of unknown function (DUF2784) to the pyrazinamidase enzyme. The variant molecule was computationally modelled and physico-chemical parameters determined to ascertain stability. Although the functional impact of this mutation remains unknown, its lineage specific nature highlights the importance of genetic background and warrants further study.

The discordant prevalence of Helicobacter pylori and its related diseases, for a long time, fostered certain enigmatic situations observed in the countries of the southern world. Variation in H. pylori infection rates and disease outcomes among different populations in multi-ethnic Malaysia provides a unique opportunity to understand dynamics of host-pathogen interaction and genome evolution. In this study, we extensively analyzed and compared genomes of 27 Malaysian H. pylori isolates and identified three major phylogeographic lineages: hspEastAsia, hpEurope and hpSouthIndia. The analysis of the virulence genes within the core genome, however, revealed a comparable pathogenic potential of the strains. In addition, we identified four genes limited to strains of East-Asian lineage. Our analyses identified a few strain-specific genes encoding restriction modification systems and outlined 311 core genes possibly under differential evolutionary constraints, among the strains representing different ethnic groups. The cagA and vacA genes also showed variations in accordance with the host genetic background of the strains. Moreover, restriction modification genes were found to be significantly enriched in East-Asian strains. An understanding of these variations in the genome content would provide significant insights into various adaptive and host modulation strategies harnessed by H. pylori to effectively persist in a host-specific manner.

Tuberculosis (TB) is the deadly infectious disease challenging the public health globally and its impact is further aggravated by co-infection with HIV and the emergence of drug resistant strains of Mycobacterium tuberculosis. In this study, we attempted to characterise the Rv2004c encoded protein, a member of DosR regulon, for its role in drug resistance. In silico docking analysis revealed that Rv2004c binds with streptomycin (SM). Phosphotransferase assay demonstrated that Rv2004c possibly mediates SM resistance through the aminoglycoside phosphotransferase activity. Further, E. coli expressing Rv2004c conferred resistance to 100μM of SM in liquid broth cultures indicating a mild aminoglycoside phosphotransferase activity of Rv2004c. Moreover, we investigated the role of MSMEG_3942 (an orthologous gene of Rv2004c) encoded protein in intracellular survival, its effect on in-vitro growth and its expression in different stress conditions by over expressing it in Mycobacterium smegmatis (M. smegmatis). MSMEG_3942 overexpressing recombinant M. smegmatis strains grew faster in acidic medium and also showed higher bacillary counts in infected macrophages when compared to M. smegmatis transformed with vector alone. Our results are likely to contribute to the better understanding of the involvement of Rv2004c in partial drug resistance, intracellular survival and adaptation of bacilli to stress conditions.

A wide variety of genome sequencing platforms have emerged in the recent past. High-throughput platforms like Illumina and 454 are essentially adaptations of the shotgun approach generating millions of fragmented single or paired sequencing reads. To reconstruct whole genomes, the reads have to be assembled into contigs, which often require further downstream processing. The contigs can be directly ordered according to a reference, scaffolded based on paired read information, or assembled using a combination of the two approaches. While the reference-based approach appears to mask strain-specific information, scaffolding based on paired-end information suffers when repetitive elements longer than the size of the sequencing reads are present in the genome. Sequencing technologies that produce long reads can solve the problems associated with repetitive elements but are not necessarily easily available to researchers. The most common high-throughput technology currently used is the Illumina short read platform. To improve upon the shortcomings associated with the construction of draft genomes with Illumina paired-end sequencing, we developed Contig-Layout-Authenticator (CLA). The CLA pipeline can scaffold reference-sorted contigs based on paired reads, resulting in better assembled genomes. Moreover, CLA also hints at probable misassemblies and contaminations, for the users to cross-check before constructing the consensus draft. The CLA pipeline was designed and trained extensively on various bacterial genome datasets for the ordering and scaffolding of large repetitive contigs. The tool has been validated and compared favorably with other widely-used scaffolding and ordering tools using both simulated and real sequence datasets. CLA is a user friendly tool that requires a single command line input to generate ordered scaffolds.

Typhoid fever poses significant burden on healthcare systems in Southeast Asia and other endemic countries. Several epidemiological and genomic studies have attributed pseudogenisation to be the major driving force for the evolution of Salmonella Typhi although its real potential remains elusive. In the present study, we analyzed genomes of S. Typhi from different parts of Southeast Asia and Oceania, comprising of isolates from outbreak, sporadic and carrier cases. The genomes showed high genetic relatedness with limited opportunity for gene acquisition as evident from pan-genome structure. Given that pseudogenisation is an active process in S. Typhi, we further investigated core and pan-genome profiles of functional and pseudogenes separately. We observed a decline in core functional gene content and a significant increase in accessory pseudogene content. Upon functional classification, genes encoding metabolic functions formed a major constituent of pseudogenes as well as core functional gene clusters with SNPs. Further, an in-depth analysis of accessory pseudogene content revealed the existence of heterogeneous complements of functional and pseudogenes among the strains. In addition, these polymorphic genes were also enriched in metabolism related functions. Thus, the study highlights the existence of heterogeneous strains in a population with varying metabolic potential and that S. Typhi possibly resorts to metabolic fine tuning for its adaptation.

We report whole-genome sequences of two clinical isolates of Mycobacterium tuberculosis isolated from patients in Odisha, India. The sequence analysis revealed that these isolates are of an ancestral type and might represent some of the “pristine” isolates in India that have not admixed with other lineages.

Salmonella enterica serovar Typhi is a human pathogen that causes typhoid fever predominantly in developing countries. In this article, we describe the whole genome sequence of the S. Typhi strain CR0044 isolated from a typhoid fever carrier in Kelantan, Malaysia. These data will further enhance the understanding of its host persistence and adaptive mechanism.

Many of the developing countries of the Southeast Asian region are significantly affected by endemic typhoid fever, possibly as a result of marginal living standards. It is an important public health problem in countries such as Papua New Guinea, which is geographically close to some of the foci of endemicity in Asia. The severity of the disease varies in different regions, and this may be attributable to genetic diversity among the native strains. Genome sequence data on strains from different countries are needed to clearly understand their genetic makeup and virulence potential. We describe the genomes of twoSalmonellaTyphi isolates from patients with fatal and nonfatal cases of typhoid fever in Papua New Guinea. We discuss in brief the underlying sequencing methodology, assembly, genome statistics, and important features of the two draft genomes, which form an essential step in our functional molecular infection epidemiology program centering on typhoid fever. The comparative genomics of these and other isolates would enable us to identify genetic rearrangements and mechanisms responsible for endemicity and the differential severity of pathogenic salmonellae in Papua New Guinea and elsewhere.

Salmonella entericaserovar Typhi is the causative agent of typhoid fever, which causes nearly 21.7 million illnesses and 217,000 deaths globally. Herein, we describe the whole-genome sequence of theSalmonellaTyphi strain ST0208, isolated from a sporadic case of typhoid fever in Kuala Lumpur, Malaysia. The whole-genome sequence and comparative genomics allow an in-depth understanding of the genetic diversity, and its link to pathogenicity and evolutionary dynamics, of this highly clonal pathogen that is endemic to Malaysia.

Among enteric pathogens, Salmonella enterica serovar Typhi is responsible for the largest number of food-borne outbreaks and fatalities. The ability of the pathogen to cause systemic infection for extended durations leads to a high cost of disease control. Chronic carriers play important roles in the evolution of Salmonella Typhi; therefore, identification and in-depth characterization of isolates from clinical cases and carriers, especially those from zones of endemicity where the pathogen has not been extensively studied, are necessary. Here, we describe the genome sequence of the highly virulent Salmonella Typhi strain BL196/05 isolated during the outbreak of typhoid in Kelantan, Malaysia, in 2005. The whole-genome sequence and comparative genomics of this strain should enable us to understand the virulence mechanisms and evolutionary dynamics of this pathogen in Malaysia and elsewhere.

The bacterial genus Paracoccus is comprised of metabolically versatile organisms having diverse degradative capabilities and potential industrial and environmental applications for bioremediation in particular. We report a de novo -assembled sequence and annotation of the genome of a novel isolate of Paracoccus denitrificans originally sourced from coal mine tailings in India. The isolate was capable of utilizing N , N -dimethylformamide (DMF) as a source of carbon and nitrogen and therefore holds potential for bioremediation and mineralization of industrial pollutants. The genome sequence and biological circuitry revealed thereupon will be invaluable in understanding the metabolic capabilities, functioning, and evolution of this important bacterial organism.

Uropathogenic Escherichia coli (UPEC) causes serious infections in people at risk and has a significant environmental prevalence due to contamination by human and animal excreta. In developing countries, UPEC assumes importance in certain dwellings because of poor community/personal hygiene and exposure to contaminated water or soil. We report the complete genome sequence of E. coli strain NA114 from India, a UPEC strain with a multidrug resistance phenotype and the capacity to produce extended-spectrum beta-lactamase. The genome sequence and comparative genomics emanating from it will be significant in under-standing the genetic makeup of diverse UPEC strains and in boosting the development of new diagnostics/vaccines.

The recent outbreak of the Escherichia coli infection in north Germany reportedly has killed 35 people and sickened about 4000. Infections were reported in 2001 and 2002 in Germany and Central Africa respectively. In order to dissect all possible connections among the three infection/outbreak strains and to extract the most virulent genes responsible for the augmented pathogenicity of the present German strain (2011), whole genome sequence data of the two isolates (E. coli TY_2482 & LB226692) were analyzed. We report that the German outbreak strain isolated in 2011 is a seasoned clone of the 2002 isolate sequenced from Central Africa (55989) and perhaps a direct descendent of the German EAEC strain reported in 2001 infection cases.

The prevalence of different H. pylori genotypes in various geographical regions indicates region-specific adaptations during the course of evolution. Complete genomes of H. pylori from countries with high infection burdens, such as India, have not yet been described. Herein we present genome sequences of two H. pylori strains, NAB47 and NAD1, from India. In this report, we briefly mention the sequencing and finishing approaches, genome assembly with downstream statistics, and important features of the two draft genomes, including their phylogenetic status. We believe that these genome sequences and the comparative genomics emanating thereupon will help us to clearly understand the ancestry and biology of the Indian H. pylori genotypes, and this will be helpful in solving the so-called Indian enigma, by which high infection rates do not corroborate the minuscule number of serious outcomes observed, including gastric cancer.

The diverse clinical outcomes of colonization by Helicobacter pylori reflect the need to understand the genomic rearrangements enabling the bacterium to adapt to host niches and exhibit varied colonization/virulence potential. We describe the genome sequences of the two serial isolates, H. pylori 2017 and 2018 (the chronological subclones of H. pylori 908), cultured in 2003 from the antrum and corpus, respectively, of an African patient who suffered from recrudescent duodenal ulcer disease. When compared with the genome of the parent strain, 908 (isolated from the antrum of the same patient in 1994), the genome sequences revealed genomic alterations relevant to virulence optimization or host-specific adaptation.

Vibrio parahaemolyticus, an important human pathogen, is associated with gastroenteritis and transmitted through partially cooked seafood. It has become a major concern in the production and trade of marine food products. The prevalence of potentially virulent and pathogenic V. parahaemolyticus in raw seafood is of public health significance. Here we describe the genome sequence of a V. parahaemolyticus isolate of crustacean origin which was cultured from prawns in 2008 in Selangor, Malaysia (isolate PCV08-7). The next generation sequencing and analysis revealed that the genome of isolate PCV08-7 has closest similarity to that of V. parahaemolyticus RIMD2210633. However, there are certain unique features of the PCV08-7 genome such as the absence of TDH-related hemolysin (TRH), and the presence of HU-alpha insertion. The genome of isolate PCV08-7 encodes a thermostable direct hemolysin (TDH), an important virulence factor that classifies PCV08-7 isolate to be a serovariant of O3:K6 strain. Apart from these, we observed that there is certain pattern of genetic rearrangements that makes V. parahaemolyticus PCV08-7 a non-pandemic clone. We present detailed genome statistics and important genetic features of this bacterium and discuss how its survival, adaptation and virulence in marine and terrestrial hosts can be understood through the genomic blueprint and that the availability of genome sequence entailing this important Malaysian isolate would likely enhance our understanding of the epidemiology, evolution and transmission of foodborne Vibrios in Malaysia and elsewhere.