Abstract:Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%-90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943-0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including "hypothetical proteins" was accurately predicted (F1 = 0.902 [0.898-0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876-0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data. IMPORTANCE Having the ability to predict the protein-encoding gene content of a genome is important for assessing genome quality, binning genomes from shotgun metagenomic assemblies, and assessing risk due to the presence of antimicrobial resistance and other virulence genes. In this study, we built a set of binary classifiers for predicting the presence or absence of variable genes occurring in 10%-90% of all publicly available E. coli genomes. Overall, the results show that a large portion of the E. coli variable gene content can be predicted with high accuracy, including genes with functions relating to horizontal gene transfer. This study offers a strategy for predicting gene content using limited input sequence data.

Genomic law guided gene prediction in fungi and metazoans.

How to overcome constraints imposed by microsporidian genome features to ensure gene prediction?

Re-Annotation of Protein-Coding Genes in the Genome of Saccharomyces Cerevisiae Based on Support Vector Machines

A New Computational Strategy for Predicting Essential Genes

A Brief Review of Computational Gene Prediction Methods

Function Prediction For Hypothetical Proteins In Yeast Saccharomyces Cerevisiae Using Multiple Sources Of High-Throughput Data

A Kind of Method on Predicting the Number of Genes in Saccharomyces Cerevisiae Genome

Using a Euclid distance discriminant method to find protein coding genes in the yeast genome.

A Critical Assessment of Mus Musculus Gene Function Prediction Using Integrated Genomic Evidence

MED: a New Non-Supervised Gene Prediction Algorithm for Bacterial and Archaeal Genomes

Predicting variable gene content in Escherichia coli using conserved genes

Accurate Prediction of Human Essential Genes Using Only Nucleotide Composition and Association Information

Prediction of Protein-Coding Small ORFs in Multi-Species Using Integrated Sequence-Derived Features and the Random Forest Model.

Conserved Codon Composition of Ribosomal Protein Coding Genes in Escherichia Coli, Mycobacterium Tuberculosis and Saccharomyces Cerevisiae: Lessons from Supervised Machine Learning in Functional Genomics

Intra-genomic genes-to-genes correlation enables genome representation

SCGPred： A Score-based Method for Gene Structure Prediction by Combining Multiple Sources of Evidence

GeneAlign: a Coding Exon Prediction Tool Based on Phylogenetical Comparisons.

Computational Methods for the Prediction of Microbial Essential Genes

Well-Characterized Sequence Features of Eukaryote Genomes and Implications for Ab Initio Gene Prediction

ProtiGeno: a prokaryotic short gene finder using protein language models

Exploring The Optimal Strategy To Predict Essential Genes In Microbes