Comparative evaluation of Nanopore polishing tools for microbial genome assembly and polishing strategies for downstream analysis

Assembling high-quality microbial genomes using only cost-effective Nanopore long-read systems such as Flongle is important to accelerate research on the microbial genome and the most critical point for this is the polishing process. In this study, we performed an evaluation based on BUSCO and Prokka gene prediction in terms of microbial genome assembly for eight state-of-the-art Nanopore polishing tools and combinations available. In the evaluation of individual tools, Homopolish, PEPPER, and Medaka demonstrated better results than others. In combination polishing, the second round Homopolish, and the PEPPER × medaka combination also showed better results than others. However, individual tools and combinations have specific limitations on usage and results. Depending on the target organism and the purpose of the downstream research, it is confirmed that there remain some difficulties in perfectly replacing the hybrid polishing carried out by the addition of a short-read. Nevertheless, through continuous improvement of the protein pores, related base-calling algorithms, and polishing tools based on improved error models, a high-quality microbial genome can be achieved using only Nanopore reads without the production of additional short-read data. The polishing strategy proposed in this study is expected to provide useful information for assembling the microbial genome using only Nanopore reads depending on the target microorganism and the purpose of the research.

Advances in Accurate Microbial Genome-Editing CRISPR Technologies

Previous studies have modified microbial genomes by introducing gene cassettes containing selectable markers and homologous DNA fragments. However, this requires several steps including homologous recombination and excision of unnecessary DNA regions, such as selectable markers from the modified genome. Further, genomic manipulation often leaves scars and traces that interfere with downstream iterative genome engineering. A decade ago, the CRISPR/Cas system (also known as the bacterial adaptive immune system) revolutionized genome editing technology. Among the various CRISPR nucleases of numerous bacteria and archaea, the Cas9 and Cas12a (Cpf1) systems have been largely adopted for genome editing in all living organisms due to their simplicity, as they consist of a single polypeptide nuclease with a target-recognizing RNA. However, accurate and fine-tuned genome editing remains challenging due to mismatch tolerance and protospacer adjacent motif (PAM)-dependent target recognition. Therefore, this review describes how to overcome the aforementioned hurdles, which especially affect genome editing in higher organisms. Additionally, the biological significance of CRISPR-mediated microbial genome editing is discussed, and future research and development directions are also proposed.

MicrobeAnnotator: a user-friendly, comprehensive functional annotation pipeline for microbial genomes

BACKGROUND

High-throughput sequencing has increased the number of available microbial genomes recovered from isolates, single cells, and metagenomes. Accordingly, fast and comprehensive functional gene annotation pipelines are needed to analyze and compare these genomes. Although several approaches exist for genome annotation, these are typically not designed for easy incorporation into analysis pipelines, do not combine results from different annotation databases or offer easy-to-use summaries of metabolic reconstructions, and typically require large amounts of computing power for high-throughput analysis not available to the average user.

RESULTS

Here, we introduce MicrobeAnnotator, a fully automated, easy-to-use pipeline for the comprehensive functional annotation of microbial genomes that combines results from several reference protein databases and returns the matching annotations together with key metadata such as the interlinked identifiers of matching reference proteins from multiple databases [KEGG Orthology (KO), Enzyme Commission (E.C.), Gene Ontology (GO), Pfam, and InterPro]. Further, the functional annotations are summarized into Kyoto Encyclopedia of Genes and Genomes (KEGG) modules as part of a graphical output (heatmap) that allows the user to quickly detect differences among (multiple) query genomes and cluster the genomes based on their metabolic similarity. MicrobeAnnotator is implemented in Python 3 and is freely available under an open-source Artistic License 2.0 from https://github.com/cruizperez/MicrobeAnnotator .

CONCLUSIONS

We demonstrated the capabilities of MicrobeAnnotator by annotating 100 Escherichia coli and 78 environmental Candidate Phyla Radiation (CPR) bacterial genomes and comparing the results to those of other popular tools. We showed that the use of multiple annotation databases allows MicrobeAnnotator to recover more annotations per genome compared to faster tools that use reduced databases and is computationally efficient for use in personal computers. The output of MicrobeAnnotator can be easily incorporated into other analysis pipelines while the results of other annotation tools can be seemingly incorporated into MicrobeAnnotator to generate summary plots.

MetaGeneHunt for protein domain annotation in short-read metagenomes

The annotation of short-reads metagenomes is an essential process to understand the functional potential of sequenced microbial communities. Annotation techniques based solely on the identification of local matches tend to confound local sequence similarity and overall protein homology and thus don't mirror the complex multidomain architecture and the shuffling of functional domains in many protein families. Here, we present MetaGeneHunt to identify specific protein domains and to normalize the hit-counts based on the domain length. We used MetaGeneHunt to investigate the potential for carbohydrate processing in the mouse gastrointestinal tract. We sampled, sequenced, and analyzed the microbial communities associated with the bolus in the stomach, intestine, cecum, and colon of five captive mice. Focusing on Glycoside Hydrolases (GHs) we found that, across samples, 58.3% of the 4,726,023 short-read sequences matching with a GH domain-containing protein were located outside the domain of interest. Next, before comparing the samples, the counts of localized hits matching the domains of interest were normalized to account for the corresponding domain length. Microbial communities in the intestine and cecum displayed characteristic GH profiles matching distinct microbial assemblages. Conversely, the stomach and colon were associated with structurally and functionally more diverse and variable microbial communities. Across samples, despite fluctuations, changes in the functional potential for carbohydrate processing correlated with changes in community composition. Overall MetaGeneHunt is a new way to quickly and precisely identify discrete protein domains in sequenced metagenomes processed with MG-RAST. In addition, using the sister program "GeneHunt" to create custom Reference Annotation Table, MetaGeneHunt provides an unprecedented way to (re)investigate the precise distribution of any protein domain in short-reads metagenomes.

Autometa: automated extraction of microbial genomes from individual shotgun metagenomes

Shotgun metagenomics is a powerful, high-resolution technique enabling the study of microbial communities in situ. However, species-level resolution is only achieved after a process of 'binning' where contigs predicted to originate from the same genome are clustered. Such culture-independent sequencing frequently unearths novel microbes, and so various methods have been devised for reference-free binning. As novel microbiomes of increasing complexity are explored, sometimes associated with non-model hosts, robust automated binning methods are required. Existing methods struggle with eukaryotic contamination and cannot handle highly complex single metagenomes. We therefore developed an automated binning pipeline, termed 'Autometa', to address these issues. This command-line application integrates sequence homology, nucleotide composition, coverage and the presence of single-copy marker genes to separate microbial genomes from non-model host genomes and other eukaryotic contaminants, before deconvoluting individual genomes from single metagenomes. The method is able to effectively separate over 1000 genomes from a metagenome, allowing the study of previously intractably complex environments at the level of single species. Autometa is freely available at https://bitbucket.org/jason_c_kwan/autometa and as a docker image at https://hub.docker.com/r/jasonkwan/autometa under the GNU Affero General Public License 3 (AGPL 3).

Real-Time Analysis and Visualization of Pathogen Sequence Data

The rapid development of sequencing technologies has to led to an explosion of pathogen sequence data, which are increasingly collected as part of routine surveillance or clinical diagnostics. In public health, sequence data are used to reconstruct the evolution of pathogens, to anticipate future spread, and to target interventions. In clinical settings, whole-genome sequencing can identify pathogens at the strain level, can be used to predict phenotypes such as drug resistance and virulence, and can inform treatment by linking closely related cases. While sequencing has become cheaper, the analysis of sequence data has become an important bottleneck. Deriving interpretable and actionable results for a large variety of pathogens, each with its own complexity, from continuously updated data is a daunting task that requires flexible bioinformatic workflows and dissemination platforms. Here, we review recent developments in real-time analyses of pathogen sequence data, with a particular focus on the visualization and integration of sequence and phenotype data.

Ortholog Identification and Comparative Analysis of Microbial Genomes Using MBGD and RECOG

Comparative genomics is becoming an essential approach for identification of genes associated with a specific function or phenotype. Here, we introduce the microbial genome database for comparative analysis (MBGD), which is a comprehensive ortholog database among the microbial genomes available so far. MBGD contains several precomputed ortholog tables including the standard ortholog table covering the entire taxonomic range and taxon-specific ortholog tables for various major taxa. In addition, MBGD allows the users to create an ortholog table within any specified set of genomes through dynamic calculations. In particular, MBGD has a "My MBGD" mode where users can upload their original genome sequences and incorporate them into orthology analysis. The created ortholog table can serve as the basis for various comparative analyses. Here, we describe the use of MBGD and briefly explain how to utilize the orthology information during comparative genome analysis in combination with the stand-alone comparative genomics software RECOG, focusing on the application to comparison of closely related microbial genomes.

Intestinal microbiome landscaping: insight in community assemblage and implications for microbial modulation strategies

High individuality, large complexity and limited understanding of the mechanisms underlying human intestinal microbiome function remain the major challenges for designing beneficial modulation strategies. Exemplified by the analysis of intestinal bacteria in a thousand Western adults, we discuss key concepts of the human intestinal microbiome landscape, i.e. the compositional and functional 'core', the presence of community types and the existence of alternative stable states. Genomic investigation of core taxa revealed functional redundancy, which is expected to stabilize the ecosystem, as well as taxa with specialized functions that have the potential to shape the microbiome landscape. The contrast between Prevotella- and Bacteroides-dominated systems has been well described. However, less known is the effect of not so abundant bacteria, for example, Dialister spp. that have been proposed to exhibit distinct bistable dynamics. Studies employing time-series analysis have highlighted the dynamical variation in the microbiome landscape with and without the effect of defined perturbations, such as the use of antibiotics or dietary changes. We incorporate ecosystem-level observations of the human intestinal microbiota and its keystone species to suggest avenues for designing microbiome modulation strategies to improve host health.

Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system

The subterranean world hosts up to one-fifth of all biomass, including microbial communities that drive transformations central to Earth's biogeochemical cycles. However, little is known about how complex microbial communities in such environments are structured, and how inter-organism interactions shape ecosystem function. Here we apply terabase-scale cultivation-independent metagenomics to aquifer sediments and groundwater, and reconstruct 2,540 draft-quality, near-complete and complete strain-resolved genomes that represent the majority of known bacterial phyla as well as 47 newly discovered phylum-level lineages. Metabolic analyses spanning this vast phylogenetic diversity and representing up to 36% of organisms detected in the system are used to document the distribution of pathways in coexisting organisms. Consistent with prior findings indicating metabolic handoffs in simple consortia, we find that few organisms within the community can conduct multiple sequential redox transformations. As environmental conditions change, different assemblages of organisms are selected for, altering linkages among the major biogeochemical cycles.

The Age of Effectors: Genome-Based Discovery and Applications

Microbial pathogens cause devastating diseases on economically and ecologically important plant species, threatening global food security, and causing billions of dollars of losses annually. During the infection process, pathogens secrete so-called effectors that support host colonization, often by deregulating host immune responses. Over the last decades, much of the research on molecular plant-microbe interactions has focused on the identification and functional characterization of such effectors. The increasing availability of sequenced plant pathogen genomes has enabled genomics-based discovery of effector candidates. Nevertheless, identification of full plant pathogen effector repertoires is often hampered by erroneous gene annotation and the localization effector genes in genomic regions that are notoriously difficult to assemble. Here, we argue that recent advances in genome sequencing technologies, genome assembly, gene annotation, as well as effector identification methods hold promise to disclose complete and correct effector repertoires. This allows to exploit complete effector repertoires, and knowledge of their diversity within pathogen populations, to develop durable and sustainable resistance breeding strategies, disease control, and management of plant pathogens.

Disease Management in the Genomics Era-Summaries of Focus Issue Papers

The genomics revolution has contributed enormously to research and disease management applications in plant pathology. This development has rapidly increased our understanding of the molecular mechanisms underpinning pathogenesis and resistance, contributed novel markers for rapid pathogen detection and diagnosis, and offered further insights into the genetics of pathogen populations on a larger scale. The availability of whole genome resources coupled with next-generation sequencing (NGS) technologies has helped fuel genomics-based approaches to improve disease resistance in crops. NGS technologies have accelerated the pace at which whole plant and pathogen genomes have become available, and made possible the metagenomic analysis of plant-associated microbial communities. Furthermore, NGS technologies can now be applied routinely and cost effectively to rapidly generate plant and/or pathogen genome or transcriptome marker sequences associated with virulence phenotypes in the pathogen or resistance phenotypes in the plant, potentially leading to improvements in plant disease management. In some systems, investments in plant and pathogen genomics have led to immediate, tangible benefits. This focus issue covers some of the systems. The articles in this focus issue range from overall perspective articles to research articles describing specific genomics applications for detection and control of diseases caused by nematode, viral, bacterial, fungal, and oomycete pathogens. The following are representative short summaries of the articles that appear in this Focus Issue .

Practical Approaches for Detecting Selection in Microbial Genomes

Microbial genome evolution is shaped by a variety of selective pressures. Understanding how these processes occur can help to address important problems in microbiology by explaining observed differences in phenotypes, including virulence and resistance to antibiotics. Greater access to whole-genome sequencing provides microbiologists with the opportunity to perform large-scale analyses of selection in novel settings, such as within individual hosts. This tutorial aims to guide researchers through the fundamentals underpinning popular methods for measuring selection in pathogens. These methods are transferable to a wide variety of organisms, and the exercises provided are designed for researchers with any level of programming experience.

Whole-genome sequencing in outbreak analysis

In addition to the ever-present concern of medical professionals about epidemics of infectious diseases, the relative ease of access and low cost of obtaining, producing, and disseminating pathogenic organisms or biological toxins mean that bioterrorism activity should also be considered when facing a disease outbreak. Utilization of whole-genome sequencing (WGS) in outbreak analysis facilitates the rapid and accurate identification of virulence factors of the pathogen and can be used to identify the path of disease transmission within a population and provide information on the probable source. Molecular tools such as WGS are being refined and advanced at a rapid pace to provide robust and higher-resolution methods for identifying, comparing, and classifying pathogenic organisms. If these methods of pathogen characterization are properly applied, they will enable an improved public health response whether a disease outbreak was initiated by natural events or by accidental or deliberate human activity. The current application of next-generation sequencing (NGS) technology to microbial WGS and microbial forensics is reviewed.

Enediynes: Exploration of microbial genomics to discover new anticancer drug leads

The enediyne natural products have been explored for their phenomenal cytotoxicity. The development of enediynes into anticancer drugs has been successfully achieved through the utilization of polymer- and antibody-drug conjugates (ADCs) as drug delivery systems. An increasing inventory of enediynes would benefit current application of ADCs in many oncology programs. Innovations in expanding the enediyne inventory should take advantage of the current knowledge of enediyne biosynthesis and post-genomics technologies. Bioinformatics analysis of microbial genomes reveals that enediynes are underexplored, in particular from Actinomycetales. This digest highlights the emerging opportunities to explore microbial genomics for the discovery of novel enediyne natural products.

SNP-associations and phenotype predictions from hundreds of microbial genomes without genome alignments

SNP-association studies are a starting point for identifying genes that may be responsible for specific phenotypes, such as disease traits. The vast bulk of tools for SNP-association studies are directed toward SNPs in the human genome, and I am unaware of any tools designed specifically for such studies in bacterial or viral genomes. The PPFS (Predict Phenotypes From SNPs) package described here is an add-on to kSNP, a program that can identify SNPs in a data set of hundreds of microbial genomes. PPFS identifies those SNPs that are non-randomly associated with a phenotype based on the χ² probability, then uses those diagnostic SNPs for two distinct, but related, purposes: (1) to predict the phenotypes of strains whose phenotypes are unknown, and (2) to identify those diagnostic SNPs that are most likely to be causally related to the phenotype. In the example illustrated here, from a set of 68 E. coli genomes, for 67 of which the pathogenicity phenotype was known, there were 418,500 SNPs. Using the phenotypes of 36 of those strains, PPFS identified 207 diagnostic SNPs. The diagnostic SNPs predicted the phenotypes of all of the genomes with 97% accuracy. It then identified 97 SNPs whose probability of being causally related to the pathogenic phenotype was >0.999. In a second example, from a set of 116 E. coli genome sequences, using the phenotypes of 65 strains PPFS identified 101 SNPs that predicted the source host (human or non-human) with 90% accuracy.

[Advances in microbial genome reduction and modification]

Microbial genome reduction and modification are important strategies for constructing cellular chassis used for synthetic biology. This article summarized the essential genes and the methods to identify them in microorganisms, compared various strategies for microbial genome reduction, and analyzed the characteristics of some microorganisms with the minimized genome. This review shows the important role of genome reduction in constructing cellular chassis.

Amino acid usage is asymmetrically biased in AT- and GC-rich microbial genomes

INTRODUCTION

Genomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates.

RESULTS

We found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB.

CONCLUSION

Genomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study.