Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Three challenging cases of infections by multidrug-resistant Serratia marcescens in patients admitted to intensive care units

The occurrence of multidrug-resistant Serratia marcescens strains represents a serious public health threat. The purpose here is to report three cases of carbapenem-resistant S. marcescens infections with unfavorable clinical outcomes and provide a molecular description of the antibiotic resistance determinants at a genomic level. We performed bacterial identification by VITEK 2 and MALDI-TOF. The minimal inhibitory concentrations of antimicrobials were determined according to the Clinical and Laboratory Standards Institute guidelines, except for tigecycline, for which they were determined using Etest strips. Preliminary screening for the presence of carbapenemases was performed by ertapenem hydrolysis using MALDI-TOF MS. Whole-genome sequencing was provided to identify genes responsible for virulence and antimicrobial resistance. Here we report three challenging cases of S. marcescens that were resistant to the most commonly used antibiotics. Otherwise, we performed a genome description, which includes several genes involved in the resistance and virulence. These cases illustrate serious infection due to multidrug-resistant organisms and the complexity of treatment. Our results highlight the need to evaluate isolates regularly during long-term hospital stay to achieve optimal quality of clinical care and thus improve patient outcomes.

Campylobacter jejuni isolated from poultry meat in Brazil: in silico analysis and genomic features of two strains with different phenotypes of antimicrobial susceptibility

Campylobacter jejuni is the most common bacterial cause of foodborne diarrheal disease worldwide and is among the antimicrobial resistant "priority pathogens" that pose greatest threat to public health. The genomes of two C. jejuni isolated from poultry meat sold on the retail market in Southern Brazil phenotypically characterized as multidrug-resistant (CJ100) and susceptible (CJ104) were sequenced and analyzed by bioinformatic tools. The isolates CJ100 and CJ104 showed distinct multilocus sequence types (MLST). Comparative genomic analysis revealed a large number of single nucleotide polymorphisms, rearrangements, and inversions in both genomes, in addition to virulence factors, genomic islands, prophage sequences, and insertion sequences. A circular 103-kilobase megaplasmid carrying virulence factors was identified in the genome of CJ100, in addition to resistance mechanisms to aminoglycosides, beta-lactams, macrolides, quinolones, and tetracyclines. The molecular characterization of distinct phenotypes of foodborne C. jejuni and the discovery of a novel virulence megaplasmid provide useful data for pan-genome and large-scale studies to monitor the virulent C. jejuni in poultry meat is warranted.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain Xff49: a new isolate obtained from common beans in Southern Brazil

The genus Xanthomonas comprises Gram-negative bacteria, many of which are phytopathogens. Xanthomonas fuscans subsp. fuscans is one of the most devastating pathogens affecting the bean plant, resulting in the common bacterial blight of bean (CBB). The disease is mainly foliar and affects a wide variety of bean species, thus acting as the yield-limiting factor for the bean crop. Here, we report the whole-genome sequencing of a new strain of X. fuscans subsp. fuscans, named Xff49, isolated from the infected and symptomatic beans from Capão do Leão, Southern Brazil. The genetic analysis demonstrated the presence of single-nucleotide variants (SNVs) in this strain, potentially affecting the mobilome, cell mobility, and inorganic ion metabolism. In addition, the analysis resulted in the identification of a new plasmid similar to the pAX22 derived from Achromobacter denitrificans, which was named plX, along with plA and plC, previously reported in other strains of X. fuscans subsp. fuscans. Xff49 represents the first Brazilian genome of X. fuscans subsp. fuscans and might provide useful information applicable to the studies of phylogenetics, evolution, and pathogenomics, thereby allowing a better understanding of the genomic features present in the Brazilian strains.

Genome of Leptospira borgpetersenii strain 4E, a highly virulent isolate obtained from Mus musculus in southern Brazil

A previous study by our group reported the isolation and characterisation of Leptospira borgpetersenii serogroup Ballum strain 4E. This strain is of particular interest because it is highly virulent in the hamster model. In this study, we performed whole-genome shotgun genome sequencing of the strain using the SOLiD sequencing platform. By assembling and analysing the new genome, we were able to identify novel features that have been previously overlooked in genome annotations of other strains belonging to the same species.

Whole-genome sequencing of Leptospira interrogans from southern Brazil: genetic features of a highly virulent strain

BACKGROUND

Leptospirosis is the most widespread zoonotic disease. It is caused by infection with pathogenic Leptospira species, of which over 300 serovars have been described. The accurate identification of the causative Leptospira spp. is required to ascertain the pathogenic status of the local isolates.

OBJECTIVES

This study aimed to obtain the complete genome sequence of a virulent Leptospira interrogans strain isolated from southern Brazil and to describe its genetic features.

METHODS

The whole genome was sequenced by next-generation sequencing (Ion Torrent). The genome was assembled, scaffolded, annotated, and manually reviewed. Mutations were identified based on a variant calling analysis using the genome of L. interrogans strain Fiocruz L1-130 as a reference.

FINDINGS

The entire genome had an average GC content of 35%. The variant calling analysis identified 119 single nucleotide polymorphisms (SNPs), from which 30 led to a missense mutation. The structural analyses identified potential evidence of genomic inversions, translocations, and deletions in both the chromosomes.

MAIN CONCLUSIONS

The genome properties provide comprehensive information about the local isolates of Leptospira spp., and thereby, could facilitate the identification of new targets for the development of diagnostic kits and vaccines.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

Draft genome sequence of Bacillus velezensis 2A-2B strain: a rhizospheric inhabitant of Sporobolus airoides (Torr.) Torr., with antifungal activity against root rot causing phytopathogens

A Bacillus velezensis strain from the rhizosphere of Sporobolus airoides (Torr.) Torr., a grass in central-north México, was isolated during a biocontrol of phytopathogens scrutiny study. The 2A-2B strain exhibited at least 60% of growth inhibition of virulent isolates of phytopathogens causing root rot. These phytopathogens include Phytophthora capsici, Fusarium solani, Fusarium oxysporum and Rhizoctonia solani. Furthermore, the 2A-2B strain is an indolacetic acid producer, and a plant inducer of PR1, which is an induced systemic resistance related gene in chili pepper plantlets. Whole genome sequencing was performed to generate a draft genome assembly of 3.953 MB with 46.36% of GC content, and a N50 of 294,737. The genome contains 3713 protein coding genes and 89 RNA genes. Moreover, comparative genome analysis revealed that the 2A-2B strain had the greatest identity (98.4%) with Bacillus velezensis.

Complete genome sequence of native Bacillus cereus strains isolated from intestinal tract of the crab Ucides sp

Bacillus cereus is a gram positive bacterium with sporulation capacity. Here, we report the complete genome sequence of two native B. cereus strains (#25 and #29) isolated from intestinal tract of the crab Ucides sp. from Pacoti River in the State of Ceará, Brazil. The findings of this study might increase the molecular information for Bacillus strains. The data can be used in comparative analyses, origin and distribution, as well support for genetic engineering.