Eukaryotes are typically depicted as descendants of archaea, but their genomes are evolutionary chimeras with genes stemming from archaea and bacteria. Which prokaryotic heritage predominates? Here, we have clustered 19,050,992 protein sequences from 5,443 bacteria and 212 archaea with 3,420,731 protein sequences from 150 eukaryotes spanning six eukaryotic supergroups. By downsampling, we obtain estimates for the bacterial and archaeal proportions. Eukaryotic genomes possess a bacterial majority of genes. On average, the majority of bacterial genes is 56% overall, 53% in eukaryotes that never possessed plastids, and 61% in photosynthetic eukaryotic lineages, where the cyanobacterial ancestor of plastids contributed additional genes to the eukaryotic lineage. Intracellular parasites, which undergo reductive evolution in adaptation to the nutrient rich environment of the cells that they infect, relinquish bacterial genes for metabolic processes. Such adaptive gene loss is most pronounced in the human parasite Encephalitozoon intestinalis with 86% archaeal and 14% bacterial derived genes. The most bacterial eukaryote genome sampled is rice, with 67% bacterial and 33% archaeal genes. The functional dichotomy, initially described for yeast, of archaeal genes being involved in genetic information processing and bacterial genes being involved in metabolic processes is conserved across all eukaryotic supergroups.

Essential genes, those indispensable for the survival of an organism, play a key role in the emerging field, synthetic biology. Characterization of functions encoded by essential genes not only has important practical implications, such as in identifying antibiotic drug targets, but can also enhance our understanding of basic biology, such as functions needed to support cellular life. Enzymes are critical for almost all cellular activities. However, essential genes have not been systematically examined from the aspect of enzymes and the chemical reactions that they catalyze. Here, by comprehensively analyzing essential genes in 14 bacterial genomes in which large-scale gene essentiality screens have been performed, we found that enzymes are enriched in essential genes. Essential enzymes have overrepresented ligases (especially those forming carbon-oxygen bonds and carbon-nitrogen bonds), nucleotidyltransferases and phosphotransferases, while have underrepresented oxidoreductases. Furthermore, essential enzymes tend to associate with more gene ontology domains. These results, from the aspect of chemical reactions, provide further insights into the understanding of functions needed to support natural cellular life, as well as synthetic cells, and provide additional parameters that can be integrated into gene essentiality prediction algorithms.

Diverse soil-resident bacteria can contribute to plant growth and health, but the molecular mechanisms enabling them to effectively colonize their plant hosts remain poorly understood. We used randomly barcoded transposon mutagenesis sequencing (RB-TnSeq) in Pseudomonas simiae, a model root-colonizing bacterium, to establish a genome-wide map of bacterial genes required for colonization of the Arabidopsis thaliana root system. We identified 115 genes (2% of all P. simiae genes) with functions that are required for maximal competitive colonization of the root system. Among the genes we identified were some with obvious colonization-related roles in motility and carbon metabolism, as well as 44 other genes that had no or vague functional predictions. Independent validation assays of individual genes confirmed colonization functions for 20 of 22 (91%) cases tested. To further characterize genes identified by our screen, we compared the functional contributions of P. simiae genes to growth in 90 distinct in vitro conditions by RB-TnSeq, highlighting specific metabolic functions associated with root colonization genes. Our analysis of bacterial genes by sequence-driven saturation mutagenesis revealed a genome-wide map of the genetic determinants of plant root colonization and offers a starting point for targeted improvement of the colonization capabilities of plant-beneficial microbes.

Bacterial natural products (NPs) are still a major source of new drug leads. Polyketides (PKs) and non-ribosomal peptides (NRP) are two pharmaceutically important families of NPs and recent studies have revealed Antarctica to harbor endemic polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) genes, likely to be involved in the production of novel metabolites. Despite this, the diversity of secondary metabolites genes in Antarctica is still poorly explored. In this study, a computational bioprospection approach was employed to study the diversity and identity of PKS and NRPS genes to one of the most biodiverse areas in maritime Antarctica-Maxwell Bay. Amplicon sequencing of soil samples targeting ketosynthase (KS) and adenylation (AD) domains of PKS and NRPS genes, respectively, revealed abundant and unexplored chemical diversity in this peninsula. About 20% of AD domain sequences were only distantly related to characterized biosynthetic genes. Several PKS and NRPS genes were found to be closely associated to recently described metabolites including those from uncultured and candidate phyla. The combination of new approaches in computational biology and new culture-dependent and -independent strategies is thus critical for the recovery of the potential novel chemistry encoded in Antarctica microorganisms.

Bacterial-derived polyketide and non-ribosomal peptide natural products are crucial sources of therapeutics and yet little is known about the conditions that favor activation of natural product genes or the regulatory machinery controlling their transcription. Recent findings suggest that the σ54 system, which includes σ54-loaded RNA polymerase and transcriptional activators called enhancer binding proteins (EBPs), might be a common regulator of natural product genes. Here, we explored this idea by analyzing a selected group of putative σ54 promoters identified in Myxococcus xanthus natural product gene clusters. We show that mutations in putative σ54-RNA polymerase binding regions and in putative Nla28 EBP binding sites dramatically reduce in vivo promoter activities in growing and developing cells. We also show in vivo promoter activities are reduced in a nla28 mutant, that Nla28 binds to wild-type fragments of these promoters in vitro, and that in vitro binding is lost when the Nla28 binding sites are mutated. Together, our results indicate that M. xanthus uses σ54 promoters for transcription of at least some of its natural product genes. Interestingly, the vast majority of experimentally confirmed and putative σ54 promoters in M. xanthus natural product loci are located within genes and not in intergenic sequences.

In this work, we report likely recurrent horizontal (lateral) gene transfer events of genes encoding pore-forming toxins of the aerolysin family between species belonging to different kingdoms of life. Clustering based on pairwise similarity and phylogenetic analysis revealed several distinct aerolysin sequence groups, each containing proteins from multiple kingdoms of life. These results strongly support at least six independent transfer events between distantly related phyla in the evolutionary history of one protein family and discount selective retention of ancestral genes as a plausible explanation for this patchy phylogenetic distribution. We discuss the possible roles of these proteins and show evidence for a convergent new function in two extant species. We hypothesize that certain gene families are more likely to be maintained following horizontal gene transfer from commensal or pathogenic organism to its host if they 1) can function alone; and 2) are immediately beneficial for the ecology of the organism, as in the case of pore-forming toxins which can be utilized in multicellular organisms for defense and predation.

Haloarchaea, a major group of archaea, are able to metabolize sugars and to live in oxygenated salty environments. Their physiology and lifestyle strongly contrast with that of their archaeal ancestors. Amino acid optimizations, which lowered the isoelectric point of haloarchaeal proteins, and abundant lateral gene transfers from bacteria have been invoked to explain this deep evolutionary transition. We use network analyses to show that the evolution of novel genes exclusive to Haloarchaea also contributed to the evolution of this group.

Otitis media (OM) is a common polymicrobial infection of the middle ear in children under the age of 15 years. A widely used experimental strategy to analyse roles of specific phenotypes of bacterial pathogens of OM is to study changes in co-infection kinetics of bacterial populations in animal models when a wild-type bacterial strain is replaced by a specific isogenic mutant strain in the co-inoculating mixtures. As relationships between the OM bacterial pathogens within the host are regulated by many interlinked processes, connecting the changes in the co-infection kinetics to a bacterial phenotype can be challenging. We investigated middle ear co-infections in adult chinchillas (Chinchilla lanigera) by two major OM pathogens: non-typeable Haemophilus influenzae (NTHi) and Moraxella catarrhalis (Mcat), as well as isogenic mutant strains in each bacterial species. We analysed the infection kinetic data using Lotka-Volterra population dynamics, maximum entropy inference and Akaike information criteria-(AIC)-based model selection. We found that changes in relationships between the bacterial pathogens that were not anticipated in the design of the co-infection experiments involving mutant strains are common and were strong regulators of the co-infecting bacterial populations. The framework developed here allows for a systematic analysis of host-host variations of bacterial populations and small sizes of animal cohorts in co-infection experiments to quantify the role of specific mutant strains in changing the infection kinetics. Our combined approach can be used to analyse the functional footprint of mutant strains in regulating co-infection kinetics in models of experimental OM and other polymicrobial diseases.

Acinetobacter baumannii is emerging as a leading global multiple-antibiotic-resistant nosocomial pathogen. The identity of genes essential for pathogenesis in a mammalian host remains largely unknown. Using transposon-directed insertion-site sequencing (TraDIS), we identified A. baumannii genes involved in bacterial survival in a leukopenic mouse model of bloodstream infection. Mice were inoculated with a pooled transposon mutant library derived from 109,000 mutants, and TraDIS was used to map transposon insertion sites in the genomes of bacteria in the inoculum and of bacteria recovered from mouse spleens. Unique transposon insertion sites were mapped and used to calculate a fitness factor for every insertion site based on its relative abundance in the inoculum and postinfection libraries. Eighty-nine transposon insertion mutants that were underrepresented after experimental infection in mice compared to their presence in the inocula were delineated as candidates for further evaluation. Genetically defined mutants lacking feoB (ferrous iron import), ddc (d-ala-d-ala-carboxypeptidase), and pntB (pyridine nucleotide transhydrogenase subunit) exhibited a fitness defect during systemic infection resulting from bacteremia. In vitro, these mutants, as well as a fepA (ferric enterobactin receptor) mutant, are defective in survival in human serum and within macrophages and are hypersensitive to killing by antimicrobial peptides compared to the survival of the parental strain under these conditions. Our data demonstrate that FepA is involved in the uptake of exogenous enterobactin in A. baumannii. Genetic complementation rescues the phenotypes of mutants in assays that emulate conditions encountered during infection. In summary, we have determined novel A. baumannii fitness genes involved in the pathogenesis of mammalian infection. IMPORTANCE A. baumannii is a significant cause of bacterial bloodstream infection in humans. Since multiple antibiotic resistance is becoming more common among strains of A. baumannii, there is an urgent need to develop novel tools to treat infections caused by this dangerous pathogen. To develop knowledge-guided treatment approaches for A. baumannii, a thorough understanding of the mechanism by which this pathogen causes bloodstream infection is required. Here, using a mouse model of infection, we report the identification of A. baumannii genes that are critical for the ability of this pathogen to cause bloodstream infections. This study lays the foundation for future research on A. baumannii genes that can be targeted to develop novel therapeutics against this emerging human pathogen.

The lack of effective and well-tolerated therapies against antibiotic-resistant bacteria is a global public health problem leading to prolonged treatment and increased mortality. To improve the efficacy of existing antibiotic compounds, we introduce a new method for strategically inducing antibiotic hypersensitivity in pathogenic bacteria. Following the systematic verification that the AcrAB-TolC efflux system is one of the major determinants of the intrinsic antibiotic resistance levels in Escherichia coli, we have developed a short antisense oligomer designed to inhibit the expression of acrA and increase antibiotic susceptibility in E. coli. By employing this strategy, we can inhibit E. coli growth using 2- to 40-fold lower antibiotic doses, depending on the antibiotic compound utilized. The sensitizing effect of the antisense oligomer is highly specific to the targeted gene's sequence, which is conserved in several bacterial genera, and the oligomer does not have any detectable toxicity against human cells. Finally, we demonstrate that antisense oligomers improve the efficacy of antibiotic combinations, allowing the combined use of even antagonistic antibiotic pairs that are typically not favored due to their reduced activities.

The plant rhizosphere harbors a diverse population of microorganisms, including beneficial plant growth-promoting bacteria (PGPB), that colonize plant roots and enhance growth and productivity. In order to specifically define bacterial traits that contribute to this beneficial interaction, we used high-throughput transposon mutagenesis sequencing (TnSeq) in two model root-bacterium systems associated with Setaria viridis: Azoarcus olearius DQS4T and Herbaspirillum seropedicae SmR1. This approach identified ∼100 significant genes for each bacterium that appeared to confer a competitive advantage for root colonization. Most of the genes identified specifically in A. olearius encoded metabolism functions, whereas genes identified in H. seropedicae were motility related, suggesting that each strain requires unique functions for competitive root colonization. Genes were experimentally validated by site-directed mutagenesis, followed by inoculation of the mutated bacteria onto S. viridis roots individually, as well as in competition with the wild-type strain. The results identify key bacterial functions involved in iron uptake, polyhydroxybutyrate metabolism, and regulation of aromatic metabolism as important for root colonization. The hope is that by improving our understanding of the molecular mechanisms used by PGPB to colonize plants, we can increase the adoption of these bacteria in agriculture to improve the sustainability of modern cropping systems.IMPORTANCE There is growing interest in the use of associative, plant growth-promoting bacteria (PGPB) as biofertilizers to serve as a sustainable alternative for agriculture application. While a variety of mechanisms have been proposed to explain bacterial plant growth promotion, the molecular details of this process remain unclear. The current research supports the idea that PGPB use in agriculture will be promoted by gaining more knowledge as to how these bacteria colonize plants, promote growth, and do so consistently. Specifically, the research seeks to identify those bacterial genes involved in the ability of two, PGPB strains, Azoarcus olearius and Herbaspirillum seropedicae, to colonize the roots of the C4 model grass Setaria viridis. Applying a transposon mutagenesis (TnSeq) approach, we assigned phenotypes and function to genes that affect bacterial competitiveness during root colonization. The results suggest that each bacterial strain requires unique functions for root colonization but also suggests that a few, critical functions are needed by both bacteria, pointing to some common mechanisms. The hope is that such information can be exploited to improve the use and performance of PGPB in agriculture.

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Gene arrays and operons that encode functionally linked proteins form the most basic unit of transcriptional regulation in bacteria. Rules that govern the order and orientation of genes in these systems have been defined; however, these were based on a small set of genomes that may not be representative. The growing availability of large genomic datasets presents an opportunity to test these rules, to define the full range and diversity of these systems, and to understand their evolution. Here we present SLING, a tool to Search for LINked Genes by searching for a single functionally essential gene, along with its neighbours in a rule-defined proximity (https://github.com/ghoresh11/sling/wiki). Examining this subset of genes enables us to understand the basic diversity of these genetic systems in large datasets. We demonstrate the utility of SLING on a clinical collection of enteropathogenic Escherichia coli for two relevant operons: toxin antitoxin (TA) systems and RND efflux pumps. By examining the diversity of these systems, we gain insight on distinct classes of operons which present variable levels of prevalence and ability to be lost or gained. The importance of this analysis is not limited to TA systems and RND pumps, and can be expanded to understand the diversity of many other relevant gene arrays.

The population structure and diversity of Lactococcus lactis subsp. lactis, a major industrial bacterium involved in milk fermentation, was determined at both gene and genome level. Seventy-six lactococcal isolates of various origins were studied by different genotyping methods and thirty-six strains displaying unique macrorestriction fingerprints were analyzed by a new multilocus sequence typing (MLST) scheme. This gene-based analysis was compared to genomic characteristics determined by pulsed-field gel electrophoresis (PFGE).

This study was aimed at identifying differentially expressed genes (DEGs) in bacterial and fungal keratitis. The candidate genes can be selected and quantified to distinguish between causative agents of infectious keratitis to improve therapeutic outcomes.

Over the last years a great number of bacterial genomes were sequenced. Now one of the most important challenges of computational genomics is the functional annotation of nucleic acid sequences. In this study we presented the computational method and the annotation system for predicting biological functions using phylogenetic profiles. The phylogenetic profile of a gene was created by way of searching for similarities between the nucleotide sequence of the gene and 1204 reference genomes, with further estimation of the statistical significance of found similarities. The profiles of the genes with known functions were used for prediction of possible functions and functional groups for the new genes. We conducted the functional annotation for genes from 104 bacterial genomes and compared the functions predicted by our system with the already known functions. For the genes that have already been annotated, the known function matched the function we predicted in 63% of the time, and in 86% of the time the known function was found within the top five predicted functions. Besides, our system increased the share of annotated genes by 19%. The developed system may be used as an alternative or complementary system to the current annotation systems.

Antibiotic resistance spreads among bacteria through horizontal transfer of antibiotic resistance genes (ARGs). Here, we set out to determine predictive features of ARG transfer among bacterial clades. We use a statistical framework to identify putative horizontally transferred ARGs and the groups of bacteria that disseminate them. We identify 152 gene exchange networks containing 22,963 bacterial genomes. Analysis of ARG-surrounding sequences identify genes encoding putative mobilisation elements such as transposases and integrases that may be involved in gene transfer between genomes. Certain ARGs appear to be frequently mobilised by different mobile genetic elements. We characterise the phylogenetic reach of these mobilisation elements to predict the potential future dissemination of known ARGs. Using a separate database with 472,798 genomes from Streptococcaceae, Staphylococcaceae and Enterobacteriaceae, we confirm 34 of 94 predicted mobilisations. We explore transfer barriers beyond mobilisation and show experimentally that physiological constraints of the host can explain why specific genes are largely confined to Gram-negative bacteria although their mobile elements support dissemination to Gram-positive bacteria. Our approach may potentially enable better risk assessment of future resistance gene dissemination.

Because of a loss-of-function mutation in the GGTA1 gene, humans are unable to synthetize α1,3-Galactose (Gal) decorated glycans and develop high levels of circulating anti-α1,3-Galactose antibodies (anti-Gal Abs). Anti-Gal Abs have been identified as a major obstacle of organ xenotransplantation and play a role in several host-pathogen relationships including potential susceptibility to infection. Anti-Gal Abs are supposed to stem from immunization against the gut microbiota, an assumption derived from the observation that some pathogens display α1,3-Gal and that antibiotic treatment decreases the level of anti-Gal. However, there is little information to date concerning the microorganisms producing α1,3-Gal in the human gut microbiome. Here, available α1,3-Galactosyltransferase (GT) gene sequences from gut bacteria were selectively quantified for the first time in the gut microbiome shotgun sequences of 163 adult individuals from three published population-based metagenomics analyses. We showed that most of the gut microbiome of adult individuals contained a small set of bacteria bearing α1,3-GT genes. These bacteria belong mainly to the Enterobacteriaceae family, including Escherichia coli, but also to Pasteurellaceae genera, Haemophilus influenza and Lactobacillus species. α1,3-Gal antigens and α1,3-GT activity were detected in healthy stools of individuals exhibiting α1,3-GT bacterial gene sequences in their shotgun data.

Transcription start sites (TSSs) lying inside annotated genes, on the same or opposite strand, have been observed in diverse bacteria, but the function of these unexpected transcripts is unclear. Here, we use the metal-reducing bacterium Shewanella oneidensis MR-1 and its relatives to study the evolutionary conservation of unexpected TSSs. Using high-resolution tiling microarrays and 5'-end RNA sequencing, we identified 2,531 TSSs in S. oneidensis MR-1, of which 18% were located inside coding sequences (CDSs). Comparative transcriptome analysis with seven additional Shewanella species revealed that the majority (76%) of the TSSs within the upstream regions of annotated genes (gTSSs) were conserved. Thirty percent of the TSSs that were inside genes and on the sense strand (iTSSs) were also conserved. Sequence analysis around these iTSSs showed conserved promoter motifs, suggesting that many iTSS are under purifying selection. Furthermore, conserved iTSSs are enriched for regulatory motifs, suggesting that they are regulated, and they tend to eliminate polar effects, which confirms that they are functional. In contrast, the transcription of antisense TSSs located inside CDSs (aTSSs) was significantly less likely to be conserved (22%). However, aTSSs whose transcription was conserved often have conserved promoter motifs and drive the expression of nearby genes. Overall, our findings demonstrate that some internal TSSs are conserved and drive protein expression despite their unusual locations, but the majority are not conserved and may reflect noisy initiation of transcription rather than a biological function. Importance: The first step of gene expression is the initiation of transcription from promoters, which have been traditionally thought to be located upstream of genes. Recently, studies showed that in diverse bacteria, promoters are often located inside genes. It has not been clear if these unexpected promoters are important to the organism or if they result from transcriptional noise. Here, we identify and examine promoters in eight related bacterial species. Promoters that lie within genes on the sense strand are often conserved as locations and in their sequences. Furthermore, these promoters often affect the bacterium's growth. Thus, many of these unexpected promoters are likely functional. Fewer promoters that lie within genes on the antisense strand are conserved, but the conserved ones seem to drive the expression of nearby genes.

Although bioluminescent bacteria are the most abundant and widely distributed of all light-emitting organisms, the biological role and evolutionary history of bacterial luminescence are still shrouded in mystery. Bioluminescence has so far been observed in the genomes of three families of Gammaproteobacteria in the form of canonical lux operons that adopt the CDAB(F)E(G) gene order. LuxA and luxB encode the two subunits of bacterial luciferase responsible for light-emission. Our deep exploration of public marine environmental databases considerably expands this view by providing a catalog of new lux homolog sequences, including 401 previously unknown luciferase-related genes. It also reveals a broader diversity of the lux operon organization, which we observed in previously undescribed configurations such as CEDA, CAED and AxxCE. This expanded operon diversity provides clues for deciphering lux operon evolution and propagation within the bacterial domain. Leveraging quantitative tracking of marine bacterial genes afforded by planetary scale metagenomic sampling, our study also reveals that the novel lux genes and operons described herein are more abundant in the global ocean than the canonical CDAB(F)E(G) operon.