The degree of redundancy in metabolic genes is linked to mode of metabolism

An understanding of the factors favoring the maintenance of duplicate genes in microbial genomes is essential for developing models of microbial evolution. A genome-scale flux-balance analysis of the metabolic network of Saccharomyces cerevisiae has suggested that gene duplications primarily provide increased enzyme dosage to enhance metabolic flux because the incidence of gene duplications in essential genes is no higher than that in nonessential genes. Here, we used genome-scale metabolic models to analyze the extent of genetic and biochemical redundancy in prokaryotes that are either specialists, with one major mode of energy generation, or generalists, which have multiple metabolic strategies for conservation of energy. Surprisingly, the results suggest that generalists, such as Escherichia coli and Bacillus subtilis, are similar to the eukaryotic generalist, S. cerevisiae, in having a low percentage (<10%) of essential genes and few duplications of these essential genes, whereas metabolic specialists, such as Geobacter sulfurreducens and Methanosarcina barkeri, have a high percentage (>30%) of essential genes and a high degree of genetic redundancy in these genes compared to nonessential genes. Furthermore, the specialist organisms appear to rely more on gene duplications rather than alternative-but-equivalent metabolic pathways to provide resilience to gene loss. Generalists rely more on alternative pathways. Thus, the concept that the role of gene duplications is to boost enzymatic flux rather than provide metabolic resilience may not be universal. Rather, the degree of gene duplication in microorganisms may be linked to mode of metabolism and environmental niche.     

Bacterial Protein Interaction Networks: Connectivity is Ruled by Gene Conservation,Essentiality and Function

BackgroundProtein-protein interaction (PPI) networks are the backbone of all processes in living cells. In this work, we relate conservation, essentiality and functional repertoire of a gene to the connectivity k (i.e. the number of interactions, links) of the corresponding protein in the PPI network.MethodsOn a set of 42 bacterial genomes of different sizes, and with reasonably separated evolutionary trajectories, we investigate three issues: i) whether the distribution of connectivities changes between PPI subnetworks of essential and nonessential genes; ii) how gene conservation, measured both by the evolutionary retention index (ERI) and by evolutionary pressures, is related to the connectivity of the corresponding protein; iii) how PPI connectivities are modulated by evolutionary and functional relationships, as represented by the Clusters of Orthologous Genes (COGs).ResultsWe show that conservation, essentiality and functional specialisation of genes constrain the connectivity of the corresponding proteins in bacterial PPI networks. In particular, we isolated a core of highly connected proteins (connectivities k≥40), which is ubiquitous among the species considered here, though mostly visible in the degree distributions of bacteria with small genomes (less than 1000 genes).ConclusionThe genes that support this highly connected core are conserved, essential and, in most cases, belong to the COG cluster J, related to ribosomal functions and the processing of genetic information.     

Comparative analysis of essential genes and nonessential genes in Escherichia coli K12

Genes can be classified as essential or nonessential based on their indispensability for a living organism. Previous researches have suggested that essential genes evolve more slowly than nonessential genes and the impact of gene dispensability on a gene’s evolutionary rate is not as strong as expected. However, findings have not been consistent and evidence is controversial regarding the relationship between the gene indispensability and the rate of gene evolution. Understanding how different classes of genes evolve is essential for a full understanding of evolutionary biology, and may have medical relevance in the design of new antibacterial agents. We therefore performed an investigation into the properties of essential and nonessential genes. Analysis of evolutionary conservation, protein length distribution and amino acid usage between essential and nonessential genes in Escherichia coli K12 demonstrated that essential genes are relatively preserved throughout the bacterial kingdom when compared to nonessential genes. Furthermore, results show that essential genes, compared to nonessential genes, have a significantly higher proportion of large (>534 amino acids) and small proteins (<139 amino acids) relative to medium-sized proteins. The pattern of amino acids usage shows a similar trend for essential and nonessential genes, although some notable exceptions are observed. These findings help to clarify our understanding of the evolutionary mechanisms of essential and nonessential genes, relevant to the study of mutagenesis and possibly allowing prediction of gene properties in other poorly understood organisms.     

Patterns of evolutionary conservation of essential genes correlate with their compensability

Essential genes code for fundamental cellular functions required for the viability of an organism. For this reason, essential genes are often highly conserved across organisms. However, this is not always the case: orthologues of genes that are essential in one organism are sometimes not essential in other organisms or are absent from their genomes. This suggests that, in the course of evolution, essential genes can be rendered nonessential. How can a gene become non-essential? Here we used genetic manipulation to deplete the products of 26 different essential genes in Escherichia coli. This depletion results in a lethal phenotype, which could often be rescued by the overexpression of a non-homologous, non-essential gene, most likely through replacement of the essential function. We also show that, in a smaller number of cases, the essential genes can be fully deleted from the genome, suggesting that complete functional replacement is possible. Finally, we show that essential genes whose function can be replaced in the laboratory are more likely to be non-essential or not present in other taxa. These results are consistent with the notion that patterns of evolutionary conservation of essential genes are influenced by their compensability—that is, by how easily they can be functionally replaced, for example through increased expression of other genes.     

Two Rapidly Evolving Genes Contribute to Male Fitness in Drosophila

Purifying selection often results in conservation of gene sequence and function. The most functionally conserved genes are also thought to be among the most biologically essential. These observations have led to the use of sequence conservation as a proxy for functional conservation. Here we describe two genes that are exceptions to this pattern. We show that lack of sequence conservation among orthologs of CG15460 and CG15323—herein named jean-baptiste (jb) and karr, respectively—does not necessarily predict lack of functional conservation. These two Drosophila melanogaster genes are among the most rapidly evolving protein-coding genes in this species, being nearly as diverged from their D. yakuba orthologs as random sequences are. jb and karr are both expressed at an elevated level in larval males and adult testes, but they are not accessory gland proteins and their loss does not affect male fertility. Instead, knockdown of these genes in D. melanogaster via RNA interference caused male-biased viability defects. These viability effects occur prior to the third instar for jb and during late pupation for karr. We show that putative orthologs to jb and karr are also expressed strongly in the testes of other Drosophila species and have similar gene structure across species despite low levels of sequence conservation. While standard molecular evolution tests could not reject neutrality, other data hint at a role for natural selection. Together these data provide a clear case where a lack of sequence conservation does not imply a lack of conservation of expression or function.     

Extensive Identification of Bacterial Riboflavin Transporters and Their Distribution across Bacterial Species

Riboflavin, the precursor for the cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide, is an essential metabolite in all organisms. While the functions for de novo riboflavin biosynthesis and riboflavin import may coexist in bacteria, the extent of this co-occurrence is undetermined. The RibM, RibN, RfuABCD and the energy-coupling factor-RibU bacterial riboflavin transporters have been experimentally characterized. In addition, ImpX, RfnT and RibXY are proposed as riboflavin transporters based on positional clustering with riboflavin biosynthetic pathway (RBP) genes or conservation of the FMN riboswitch regulatory element. Here, we searched for the FMN riboswitch in bacterial genomes to identify genes encoding riboflavin transporters and assessed their distribution among bacteria. Two new putative riboflavin transporters were identified: RibZ in Clostridium and RibV in Mesoplasma florum. Trans-complementation of an Escherichia coli riboflavin auxotroph strain confirmed the riboflavin transport activity of RibZ from Clostridium difficile, RibXY from Chloroflexus aurantiacus, ImpX from Fusobacterium nucleatum and RfnT from Ochrobactrum anthropi. The analysis of the genomic distribution of all known bacterial riboflavin transporters revealed that most occur in species possessing the RBP and that some bacteria may even encode functional riboflavin transporters from two different families. Our results indicate that some species possess ancestral riboflavin transporters, while others possess transporters that appear to have evolved recently. Moreover, our data suggest that unidentified riboflavin transporters also exist. The present study doubles the number of experimentally characterized riboflavin transporters and suggests a specific, non-accessory role for these proteins in riboflavin-prototrophic bacteria.     

A genome-wide screen for essential yeast genes that affect telomere length maintenance

Telomeres are structures composed of repetitive DNA and proteins that protect the chromosomal ends in eukaryotic cells from fusion or degradation, thus contributing to genomic stability. Although telomere length varies between species, in all organisms studied telomere length appears to be controlled by a dynamic equilibrium between elongating mechanisms (mainly addition of repeats by the enzyme telomerase) and nucleases that shorten the telomeric sequences. Two previous studies have analyzed a collection of yeast deletion strains (deleted for nonessential genes) and found over 270 genes that affect telomere length (Telomere Length Maintenance or TLM genes). Here we complete the list of TLM by analyzing a collection of strains carrying hypomorphic alleles of most essential genes (DAmP collection). We identify 87 essential genes that affect telomere length in yeast. These genes interact with the nonessential TLM genes in a significant manner, and provide new insights on the mechanisms involved in telomere length maintenance. The newly identified genes span a variety of cellular processes, including protein degradation, pre-mRNA splicing and DNA replication.     

Essential Plasticity and Redundancy of Metabolism Unveiled by Synthetic Lethality Analysis

We unravel how functional plasticity and redundancy are essential mechanisms underlying the ability to survive of metabolic networks. We perform an exhaustive computational screening of synthetic lethal reaction pairs in Escherichia coli in a minimal medium and we find that synthetic lethal pairs divide in two different groups depending on whether the synthetic lethal interaction works as a backup or as a parallel use mechanism, the first corresponding to essential plasticity and the second to essential redundancy. In E. coli, the analysis of pathways entanglement through essential redundancy supports the view that synthetic lethality affects preferentially a single function or pathway. In contrast, essential plasticity, the dominant class, tends to be inter-pathway but strongly localized and unveils Cell Envelope Biosynthesis as an essential backup for Membrane Lipid Metabolism. When comparing E. coli and Mycoplasma pneumoniae, we find that the metabolic networks of the two organisms exhibit a large difference in the relative importance of plasticity and redundancy which is consistent with the conjecture that plasticity is a sophisticated mechanism that requires a complex organization. Finally, coessential reaction pairs are explored in different environmental conditions to uncover the interplay between the two mechanisms. We find that synthetic lethal interactions and their classification in plasticity and redundancy are basically insensitive to medium composition, and are highly conserved even when the environment is enriched with nonessential compounds or overconstrained to decrease maximum biomass formation.     

Genome-Scale Co-Expression Network Comparison across Escherichia coli and Salmonella enterica Serovar Typhimurium Reveals Significant Conservation at the Regulon Level of Local Regulators Despite Their Dissimilar Lifestyles

Availability of genome-wide gene expression datasets provides the opportunity to study gene expression across different organisms under a plethora of experimental conditions. In our previous work, we developed an algorithm called COMODO (COnserved MODules across Organisms) that identifies conserved expression modules between two species. In the present study, we expanded COMODO to detect the co-expression conservation across three organisms by adapting the statistics behind it. We applied COMODO to study expression conservation/divergence between Escherichia coli, Salmonella enterica, and Bacillus subtilis. We observed that some parts of the regulatory interaction networks were conserved between E. coli and S. enterica especially in the regulon of local regulators. However, such conservation was not observed between the regulatory interaction networks of B. subtilis and the two other species. We found co-expression conservation on a number of genes involved in quorum sensing, but almost no conservation for genes involved in pathogenicity across E. coli and S. enterica which could partially explain their different lifestyles. We concluded that despite their different lifestyles, no significant rewiring have occurred at the level of local regulons involved for instance, and notable conservation can be detected in signaling pathways and stress sensing in the phylogenetically close species S. enterica and E. coli. Moreover, conservation of local regulons seems to depend on the evolutionary time of divergence across species disappearing at larger distances as shown by the comparison with B. subtilis. Global regulons follow a different trend and show major rewiring even at the limited evolutionary distance that separates E. coli and S. enterica.     

Conserved synteny at the protein family level reveals genes underlying Shewanella species�� cold tolerance and predicts their novel phenotypes

Bacteria of the genus Shewanella can thrive in different environments and demonstrate significant variability in their metabolic and ecophysiological capabilities including cold and salt tolerance. Genomic characteristics underlying this variability across species are largely unknown. In this study, we address the problem by a comparison of the physiological, metabolic, and genomic characteristics of 19 sequenced Shewanella species. We have employed two novel approaches based on association of a phenotypic trait with the number of the trait-specific protein families (Pfam domains) and on the conservation of synteny (order in the genome) of the trait-related genes. Our first approach is top-down and involves experimental evaluation and quantification of the species?? cold tolerance followed by identification of the correlated Pfam domains and genes with a conserved synteny. The second, a bottom-up approach, predicts novel phenotypes of the species by calculating profiles of each Pfam domain among their genomes and following pair-wise correlation of the profiles and their network clustering. Using the first approach, we find a link between cold and salt tolerance of the species and the presence in the genome of a Na⁺/H⁺ antiporter gene cluster. Other cold-tolerance-related genes include peptidases, chemotaxis sensory transducer proteins, a cysteine exporter, and helicases. Using the bottom-up approach, we found several novel phenotypes in the newly sequenced Shewanella species, including degradation of aromatic compounds by an aerobic hybrid pathway in Shewanella woodyi, degradation of ethanolamine by Shewanella benthica, and propanediol degradation by Shewanella putrefaciens CN32 and Shewanella sp. W3-18-1.     

What are housekeeping genes?

The concept of “housekeeping gene” has been used for four decades but remains loosely defined. Housekeeping genes are commonly described as “essential for cellular existence regardless of their specific function in the tissue or organism”, and “stably expressed irrespective of tissue type, developmental stage, cell cycle state, or external signal”. However, experimental support for the tenet that gene essentiality is linked to stable expression across cell types, conditions, and organisms has been limited. Here we use genome-scale functional genomic screens together with bulk and single-cell sequencing technologies to test this link and optimize a quantitative and experimentally validated definition of housekeeping gene. Using the optimized definition, we identify, characterize, and provide as resources, housekeeping gene lists extracted from several human datasets, and 10 other animal species that include primates, chicken, and C. elegans. We find that stably expressed genes are not necessarily essential, and that the individual genes that are essential and stably expressed can considerably differ across organisms; yet the pathways enriched among these genes are conserved. Further, the level of conservation of housekeeping genes across the analyzed organisms captures their taxonomic groups, showing evolutionary relevance for our definition. Therefore, we present a quantitative and experimentally supported definition of housekeeping genes that can contribute to better understanding of their unique biological and evolutionary characteristics.     

The transferome of metabolic genes explored: analysis of the horizontal transfer of enzyme encoding genes in unicellular eukaryotes

Background  

Metabolic networks are responsible for many essential cellular processes, and exhibit a high level of evolutionary conservation from bacteria to eukaryotes. If genes encoding metabolic enzymes are horizontally transferred and are advantageous, they are likely to become fixed. Horizontal gene transfer (HGT) has played a key role in prokaryotic evolution and its importance in eukaryotes is increasingly evident. High levels of endosymbiotic gene transfer (EGT) accompanied the establishment of plastids and mitochondria, and more recent events have allowed further acquisition of bacterial genes. Here, we present the first comprehensive multi-species analysis of E/HGT of genes encoding metabolic enzymes from bacteria to unicellular eukaryotes.     

The Highly Conserved Bacterial RNase YbeY Is Essential in Vibrio cholerae,Playing a Critical Role in Virulence,Stress Regulation,and RNA Processing

YbeY, a highly conserved protein, is an RNase in E. coli and plays key roles in both processing of the critical 3′ end of 16 S rRNA and in 70 S ribosome quality control under stress. These central roles account for YbeY''s inclusion in the postulated minimal bacterial genome. However, YbeY is not essential in E. coli although loss of ybeY severely sensitizes it to multiple physiological stresses. Here, we show that YbeY is an essential endoribonuclease in Vibrio cholerae and is crucial for virulence, stress regulation, RNA processing and ribosome quality control, and is part of a core set of RNases essential in most representative pathogens. To understand its function, we analyzed the rRNA and ribosome profiles of a V. cholerae strain partially depleted for YbeY and other RNase mutants associated with 16 S rRNA processing; our results demonstrate that YbeY is also crucial for 16 S rRNA 3′ end maturation in V. cholerae and that its depletion impedes subunit assembly into 70 S ribosomes. YbeY''s importance to V. cholerae pathogenesis was demonstrated by the complete loss of mice colonization and biofilm formation, reduced cholera toxin production, and altered expression levels of virulence-associated small RNAs of a V. cholerae strain partially depleted for YbeY. Notably, the ybeY genes of several distantly related pathogens can fully complement an E. coli ΔybeY strain under various stress conditions, demonstrating the high conservation of YbeY''s activity in stress regulation. Taken together, this work provides the first comprehensive exploration of YbeY''s physiological role in a human pathogen, showing its conserved function across species in essential cellular processes.     

Natural Selection Shapes Maintenance of Orthologous sRNAs in Divergent Host-Restricted Bacterial Genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.     

Characteristics of Plant Essential Genes Allow for within- and between-Species Prediction of Lethal Mutant Phenotypes

Essential genes represent critical cellular components whose disruption results in lethality. Characteristics shared among essential genes have been uncovered in fungal and metazoan model systems. However, features associated with plant essential genes are largely unknown and the full set of essential genes remains to be discovered in any plant species. Here, we show that essential genes in Arabidopsis thaliana have distinct features useful for constructing within- and cross-species prediction models. Essential genes in A. thaliana are often single copy or derived from older duplications, highly and broadly expressed, slow evolving, and highly connected within molecular networks compared with genes with nonlethal mutant phenotypes. These gene features allowed the application of machine learning methods that predicted known lethal genes as well as an additional 1970 likely essential genes without documented phenotypes. Prediction models from A. thaliana could also be applied to predict Oryza sativa and Saccharomyces cerevisiae essential genes. Importantly, successful predictions drew upon many features, while any single feature was not sufficient. Our findings show that essential genes can be distinguished from genes with nonlethal phenotypes using features that are similar across kingdoms and indicate the possibility for translational application of our approach to species without extensive functional genomic and phenomic resources.     

Comparative Metabolic Systems Analysis of Pathogenic Burkholderia

Burkholderia cenocepacia and Burkholderia multivorans are opportunistic drug-resistant pathogens that account for the majority of Burkholderia cepacia complex infections in cystic fibrosis patients and also infect other immunocompromised individuals. While they share similar genetic compositions, B. cenocepacia and B. multivorans exhibit important differences in pathogenesis. We have developed reconciled genome-scale metabolic network reconstructions of B. cenocepacia J2315 and B. multivorans ATCC 17616 in parallel (designated iPY1537 and iJB1411, respectively) to compare metabolic abilities and contextualize genetic differences between species. The reconstructions capture the metabolic functions of the two species and give insight into similarities and differences in their virulence and growth capabilities. The two reconstructions have 1,437 reactions in common, and iPY1537 and iJB1411 have 67 and 36 metabolic reactions unique to each, respectively. After curating the extensive reservoir of metabolic genes in Burkholderia, we identified 6 genes essential to growth that are unique to iPY1513 and 13 genes uniquely essential to iJB1411. The reconstructions were refined and validated by comparing in silico growth predictions to in vitro growth capabilities of B. cenocepacia J2315, B. cenocepacia K56-2, and B. multivorans ATCC 17616 on 104 carbon sources. Overall, we identified functional pathways that indicate B. cenocepacia can produce a wider array of virulence factors compared to B. multivorans, which supports the clinical observation that B. cenocepacia is more virulent than B. multivorans. The reconciled reconstructions provide a framework for generating and testing hypotheses on the metabolic and virulence capabilities of these two related emerging pathogens.