Tool for rapid annotation of microbial SNPs (TRAMS): a simple program for rapid annotation of genomic variation in prokaryotes

Genome Analysis of Erwinia amylovora Strains Responsible for a Fire Blight Outbreak in Korea

Erwinia amylovora is a plant-pathogenic bacterium that causes fire blight disease in Rosaceae plants. Since fire blight is highly contagious and results in serious losses once introduced, it is regulated as a quarantine disease. Recently, for the first time in East Asia, fire blight emerged in Korea with strains of E. amylovora being isolated from lesions of infected trees. Five of those strains were selected and subjected to whole-genome shotgun sequencing. Each strain had two circular replicons, a 3.8-Mb chromosome and a 28-kb plasmid. The genome sequences were compared with those of other E. amylovora strains isolated from different hosts or geographical regions. Genome synteny was analyzed and sequence variations including nucleotide substitutions, inversions, insertions, and deletions were detected. Analysis of the population genomic structure revealed that the five strains form a distinct structural group. Phylogenomic analysis was performed to infer the evolutionary relationships among E. amylovora strains, which indicated that the Korean isolates, all descended from a common ancestor, are closely related to a lineage of North American strains. These results provide useful information for understanding the genomic dynamics of E. amylovora strains including those in Korea, developing genetic markers for surveillance of the pathogen or diagnosis of the disease, and eventually developing measures to eradicate it.

Genomic analyses of Staphylococcus aureus clonal complex 45 isolates does not distinguish nasal carriage from bacteraemia

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus. To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

Elucidation of drug resistance mutations in Mycobacterium tuberculosis isolates from North India by whole-genome sequencing

OBJECTIVES

Rapid diagnosis of drug-resistant tuberculosis (TB) is required for better patient management and treatment outcomes. Whole-genome sequencing (WGS) can be used to detect single nucleotide polymorphisms (SNPs) and deletions/insertions that are responsible for mostMycobacterium tuberculosis drug resistance. WGS is being performed at scale in high-income countries, but there are limited reports of its use in India.

METHODS

In this study, 33 clinicalM. tuberculosis isolates from the Mycobacterial Repository in Chandigarh underwent WGS. Phenotypic drug susceptibility testing was performed according to World Health Organization (WHO) recommendations. Four isolates were excluded from the analysis due to culture contamination or mislabelling during the study.

RESULTS

Among the remaining 29 isolates, 21 (72.4%) were multidrug-resistant TB (MDR-TB) and 1 (3.4%) was extensively-drug resistant TB (XDR-TB). The most common mutations observed for isoniazid, rifampicin, ofloxacin and kanamycin resistance werekatG(S315T), rpoB(S450L), gyrA(A90V) and rrs(A1401G), respectively. The isolates mainly belonged to lineages 2 and 3, with most MDR-TB among lineage 2 isolates.

CONCLUSION

WGS ofM. tuberculosis isolates allows the detection of drug resistance to all drugs in a single test and also provides insight into the evolution and drug-resistant TB.

Genome sequence analysis of 91 Salmonella Enteritidis isolates from mice caught on poultry farms in the mid 1990s

A total of 91 draft genome sequences were used to analyze isolates of Salmonella enterica serovar Enteritidis obtained from feral mice caught on poultry farms in Pennsylvania. One objective was to find mutations disrupting open reading frames (ORFs) and another was to determine if ORF-disruptive mutations were present in isolates obtained from other sources. A total of 83 mice were obtained between 1995-1998. Isolates separated into two genomic clades and 12 subgroups due to 742 mutations. Nineteen ORF-disruptive mutations were found, and in addition, bigA had exceptional heterogeneity requiring additional evaluation. The TRAMS algorithm detected only 6 ORF disruptions. The sefD mutation was the most frequently encountered mutation and it was prevalent in human, poultry, environmental and mouse isolates. These results confirm previous assessments of the mouse as a rich source of Salmonella enterica serovar Enteritidis that varies in genotype and phenotype.

Whole-genome sequencing of bloodstream Staphylococcus aureus isolates does not distinguish bacteraemia from endocarditis

Most Staphylococcus aureus isolates can cause invasive disease given the right circumstances, but it is unknown if some isolates are more likely to cause severe infections than others. S. aureus bloodstream isolates from 120 patients with definite infective endocarditis and 121 with S. aureus bacteraemia without infective endocarditis underwent whole-genome sequencing. Genome-wide association analysis was performed using a variety of bioinformatics approaches including SNP analysis, accessory genome analysis and k-mer based analysis. Core and accessory genome analyses found no association with either of the two clinical groups. In this study, the genome sequences of S. aureus bloodstream isolates did not discriminate between bacteraemia and infective endocarditis. Based on our study and the current literature, it is not convincing that a specific S. aureus genotype is clearly associated to infective endocarditis in patients with S. aureus bacteraemia.

Comparison of advanced whole genome sequence-based methods to distinguish strains of Salmonella enterica serovar Heidelberg involved in foodborne outbreaks in Québec

Salmonella enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. This serovar ranks second and third among serovars that cause human infections in Québec and Canada, respectively, and has been associated with severe infections. Traditional typing methods such as PFGE do not display adequate discrimination required to resolve outbreak investigations due to the low level of genetic diversity of isolates belonging to this serovar. This study evaluates the ability of four whole genome sequence (WGS)-based typing methods to differentiate among 145 S. Heidelberg strains involved in four distinct outbreak events and sporadic cases of salmonellosis that occurred in Québec between 2007 and 2016. Isolates from all outbreaks were indistinguishable by PFGE. The core genome single nucleotide variant (SNV), core genome multilocus sequence typing (MLST) and whole genome MLST approaches were highly discriminatory and separated outbreak strains into four distinct phylogenetic clusters that were concordant with the epidemiological data. The clustered regularly interspaced short palindromic repeats (CRISPR) typing method was less discriminatory. However, CRISPR typing may be used as a secondary method to differentiate isolates of S. Heidelberg that are genetically similar but epidemiologically unrelated to outbreak events. WGS-based typing methods provide a highly discriminatory alternative to PFGE for the laboratory investigation of foodborne outbreaks.

Effector Diversification Contributes to Xanthomonas oryzae pv. oryzae Phenotypic Adaptation in a Semi-Isolated Environment

Understanding the processes that shaped contemporary pathogen populations in agricultural landscapes is quite important to define appropriate management strategies and to support crop improvement efforts. Here, we took advantage of an historical record to examine the adaptation pathway of the rice pathogen Xanthomonas oryzae pv. oryzae (Xoo) in a semi-isolated environment represented in the Philippine archipelago. By comparing genomes of key Xoo groups we showed that modern populations derived from three Asian lineages. We also showed that diversification of virulence factors occurred within each lineage, most likely driven by host adaptation, and it was essential to shape contemporary pathogen races. This finding is particularly important because it expands our understanding of pathogen adaptation to modern agriculture.

A Primer on Infectious Disease Bacterial Genomics

The number of large-scale genomics projects is increasing due to the availability of affordable high-throughput sequencing (HTS) technologies. The use of HTS for bacterial infectious disease research is attractive because one whole-genome sequencing (WGS) run can replace multiple assays for bacterial typing, molecular epidemiology investigations, and more in-depth pathogenomic studies. The computational resources and bioinformatics expertise required to accommodate and analyze the large amounts of data pose new challenges for researchers embarking on genomics projects for the first time. Here, we present a comprehensive overview of a bacterial genomics projects from beginning to end, with a particular focus on the planning and computational requirements for HTS data, and provide a general understanding of the analytical concepts to develop a workflow that will meet the objectives and goals of HTS projects.

IonGAP: integrative bacterial genome analysis for Ion Torrent sequence data

We introduce IonGAP, a publicly available Web platform designed for the analysis of whole bacterial genomes using Ion Torrent sequence data. Besides assembly, it integrates a variety of comparative genomics, annotation and bacterial classification routines, based on the widely used FASTQ, BAM and SRA file formats. Benchmarking with different datasets evidenced that IonGAP is a fast, powerful and simple-to-use bioinformatics tool. By releasing this platform, we aim to translate low-cost bacterial genome analysis for microbiological prevention and control in healthcare, agroalimentary and pharmaceutical industry applications.

AVAILABILITY AND IMPLEMENTATION

IonGAP is hosted by the ITER's Teide-HPC supercomputer and is freely available on the Web for non-commercial use at http://iongap.hpc.iter.es.

CONTACT

mcolesan@ull.edu.es or cflores@ull.edu.es

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

FEATnotator: A tool for integrated annotation of sequence features and variation, facilitating interpretation in genomics experiments

As approaches are sought for more efficient and democratized uses of non-model and expanded model genomics references, ease of integration of genomic feature datasets is especially desirable in multidisciplinary research communities. Valuable conclusions are often missed or slowed when researchers refer experimental results to a single reference sequence that lacks integrated pan-genomic and multi-experiment data in accessible formats. Association of genomic positional information, such as results from an expansive variety of next-generation sequencing experiments, with annotated reference features such as genes or predicted protein binding sites, provides the context essential for conclusions and ongoing research. When the experimental system includes polymorphic genomic inputs, rapid calculation of gene structural and protein translational effects of sequence variation from the reference can be invaluable. Here we present FEATnotator, a lightweight, fast and easy to use open source software program that integrates and reports overlap and proximity in genomic information from any user-defined datasets including those from next generation sequencing applications. We illustrate use of the tool by summarizing whole genome sequence variation of a widely used natural isolate of Arabidopsis thaliana in the context of gene models of the reference accession. Previous discovery of a protein coding deletion influencing root development is replicated rapidly. Appropriate even in investigations of a single gene or genic regions such as QTL, comprehensive reports provided by FEATnotator better prepare researchers for interpretation of their experimental results. The tool is available for download at http://featnotator.sourceforge.net.

Comparative genomic analyses reveal a lack of a substantial signature of host adaptation in Rhodococcus equi ('Prescottella equi')

Rhodococcus equi ('Prescottella equi') is a pathogenic actinomycete primarily infecting horses but has emerged as an opportunistic human pathogen. We have sequenced the genome of the type strain of this species, R. equi strain C7(T) , and compared the genome with that of another foal isolate 103S and of a human isolate ATCC 33707. The R. equi strains are closely related to each other and yet distantly related to other rhodococci and Nocardia brasiliensis. The comparison of gene contents among R. equi strains revealed minor differences that could be associated with host adaptation from foals to humans, including the presence of a paa operon in the human isolate, which is potentially involved in pathogenesis.

A lack of genetic basis for biovar differentiation in clinically important Corynebacterium diphtheriae from whole genome sequencing

The differentiation of clinically important Corynebacterium diphtheriae into specific biovars is complex and phylogenetically unclear. Comparative genomic analyses of 17 strains indicate that the division of C. diphtheriae into different biovars does not correlate with the variation in the gene content in the relevant metabolic categories that are potentially involved in the biovar discrimination. The biochemical separation is also not supported by phylogenetic analyses, suggesting molecular methods of typing C. diphtheriae strains should be adopted much more widely.