AssessORF: combining evolutionary conservation and proteomics to assess prokaryotic gene predictions

High-throughput transposon mutagenesis in the family Enterobacteriaceae reveals core essential genes and rapid turnover of essentiality

The Enterobacteriaceae are a scientifically and medically important clade of bacteria, containing the model organism Escherichia coli, as well as major human pathogens including Salmonella enterica and Klebsiella pneumoniae. Essential gene sets have been determined for several members of the Enterobacteriaceae, with the Keio E. coli single-gene deletion library often regarded as a gold standard. However, it remains unclear how gene essentiality varies between related strains and species. To investigate this, we have assembled a collection of 13 sequenced high-density transposon mutant libraries from five genera within the Enterobacteriaceae. We first assess several gene essentiality prediction approaches, investigate the effects of transposon density on essentiality prediction, and identify biases in transposon insertion sequencing data. Based on these investigations, we develop a new classifier for gene essentiality. Using this new classifier, we define a core essential genome in the Enterobacteriaceae of 201 universally essential genes. Despite the presence of a large cohort of variably essential genes, we find an absence of evidence for genus-specific essential genes. A clear example of this sporadic essentiality is given by the set of genes regulating the σE extracytoplasmic stress response, which appears to have independently acquired essentiality multiple times in the Enterobacteriaceae. Finally, we compare our essential gene sets to the natural experiment of gene loss in obligate insect endosymbionts that have emerged from within the Enterobacteriaceae. This isolates a remarkably small set of genes absolutely required for survival and identifies several instances of essential stress responses masked by redundancy in free-living bacteria.IMPORTANCEThe essential genome, that is the set of genes absolutely required to sustain life, is a core concept in genetics. Essential genes in bacteria serve as drug targets, put constraints on the engineering of biological chassis for technological or industrial purposes, and are key to constructing synthetic life. Despite decades of study, relatively little is known about how gene essentiality varies across related bacteria. In this study, we have collected gene essentiality data for 13 bacteria related to the model organism Escherichia coli, including several human pathogens, and investigated the conservation of essentiality. We find that approximately a third of the genes essential in any particular strain are non-essential in another related strain. Surprisingly, we do not find evidence for essential genes unique to specific genera; rather it appears a substantial fraction of the essential genome rapidly gains or loses essentiality during evolution. This suggests that essentiality is not an immutable characteristic but depends crucially on the genomic context. We illustrate this through a comparison of our essential genes in free-living bacteria to genes conserved in 34 insect endosymbionts with naturally reduced genomes, finding several cases where genes generally regarded as being important for specific stress responses appear to have become essential in endosymbionts due to a loss of functional redundancy in the genome.

A Practical Guide for the Quality Evaluation of Fluobodies/Chromobodies

BACKGROUND

Fluorescent proteins (FPs) are pivotal reagents for flow cytometry analysis or fluorescent microscopy. A new generation of immunoreagents (fluobodies/chromobodies) has been developed by fusing recombinant nanobodies to FPs.

METHODS

We analyzed the quality of such biomolecules by a combination of gel filtration and SDS-PAGE to identify artefacts due to aggregation or material degradation.

RESULTS

In the SDS-PAGE run, unexpected bands corresponding to separate fluobodies were evidenced and characterized as either degradation products or artefacts that systematically resulted in the presence of specific FPs and some experimental conditions. The elimination of N-terminal methionine from FPs did not impair the appearance of FP fragments, whereas the stability and migration characteristics of some FP constructs were strongly affected by heating in loading buffer, which is a step samples undergo before electrophoretic separation.

CONCLUSIONS

In this work, we provide explanations for some odd results observed during the quality control of fluobodies and summarize practical suggestions for the choice of the most convenient FPs to fuse to antibody fragments.

Long-term application of organic fertilizer prompting the dispersal of antibiotic resistance genes and their health risks in the soil plastisphere

Microplastic (MP) pollution is a rapidly growing global environmental concern that has led to the emergence of a new environmental compartment, the plastisphere, which is a hotspot for the accumulation of antibiotic resistance genes (ARGs) and human bacterial pathogens (HBPs). However, studies on the effects of long-term organic fertilizer application on the dispersal of ARGs and virulence factor genes (VFGs) in the plastisphere of farmland soil have been limited. Here, we performed a field culture experiment by burying nylon bags filled with MPs in paddy soil that had been treated with different fertilizers for over 30 years to explore the changes of ARGs and VFGs in soil plastisphere. Our results show that the soil plastisphere amplified the ARG and VFG pollution caused by organic fertilization by 1.5 and 1.4 times, respectively. And it also led to a 2.7-fold increase in the risk of horizontal gene transfer. Meanwhile, the plastisphere tended to promote deterministic process in the community assembly of HBPs, with an increase of 1.4 times. Network analysis found a significant correlation between ARGs, VFGs, and bacteria in plastisphere. Correlation analysis highlight the important role of mobile genetic elements (MGEs) and bacterial communities in shaping the abundance of ARGs and VFGs, respectively. Our findings provide new insights into the health risk associated with the soil plastisphere due ARGs and VFGs derived from organic fertilizers.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

INRI-seq enables global cell-free analysis of translation initiation and off-target effects of antisense inhibitors

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

KEMET - a python tool for KEGG Module evaluation and microbial genome annotation expansion

•

KEMET allows evaluating and expanding microbial genome annotations.

•

KEMET can identify unannotated genes having partial sequence truncations.

•

Annotation expansion can improve contextual draft genome-scale metabolic model reconstruction.

Background

The rapid accumulation of sequencing data from metagenomic studies is enabling the generation of huge collections of microbial genomes, with new challenges for mapping their functional potential. In particular, metagenome-assembled genomes are typically incomplete and harbor partial gene sequences that can limit their annotation from traditional tools. New scalable solutions are thus needed to facilitate the evaluation of functional potential in microbial genomes.

Methods

To resolve annotation gaps in microbial genomes, we developed KEMET, an open-source Python library devised for the analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) functional units. KEMET focuses on the in-depth analysis of metabolic reaction networks to identify missing orthologs through hidden Markov model profiles.

Results

We evaluate the potential of KEMET for expanding functional annotations by simulating the effect of assembly issues on real gene sequences and showing that our approach can identify missing KEGG orthologs. Additionally, we show that recovered gene annotations can sensibly increase the quality of draft genome-scale metabolic models obtained from metagenome-assembled genomes, in some cases reaching the accuracy of models generated from complete genomes.

Conclusions

KEMET therefore allows expanding genome annotations by targeted searches for orthologous sequences, enabling a better qualitative and quantitative assessment of metabolic capabilities in novel microbial organisms.

FindNonCoding: rapid and simple detection of non-coding RNAs in genomes

SUMMARY

Non-coding RNAs are often neglected during genome annotation due to their difficulty of detection relative to protein coding genes. FindNonCoding takes a pattern mining approach to capture the essential sequence motifs and hairpin loops representing a non-coding RNA family and quickly identify matches in genomes. FindNonCoding was designed for ease of use and accurately finds non-coding RNAs with a low false discovery rate.

AVAILABILITY AND IMPLEMENTATION

FindNonCoding is implemented within the DECIPHER package (v2.19.3) for R (v4.1) available from Bioconductor. Pre-trained models of common non-coding RNA families are included for bacteria, archaea and eukarya.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.