Inferring gene function from evolutionary change in signatures of translation efficiency

The F plasmid conjutome: the repertoire of E. coli proteins translocated through an F-encoded type IV secretion system

Bacterial conjugation systems pose a major threat to human health through their widespread dissemination of mobile genetic elements (MGEs) carrying cargoes of antibiotic resistance genes. Using the Cre Recombinase Assay for Translocation (CRAfT), we recently reported that the IncFV pED208 conjugation system also translocates at least 16 plasmid-encoded proteins to recipient bacteria. Here, we deployed a high-throughput CRAfT screen to identify the repertoire of chromosomally encoded protein substrates of the pED208 system. We identified 32 substrates encoded by the Escherichia coli W3110 genome with functions associated with (i) DNA/nucleotide metabolism, (ii) stress tolerance/physiology, (iii) transcriptional regulation, or (iv) toxin inhibition. The respective gene deletions did not impact pED208 transfer proficiencies, nor did Group 1 (DNA/nucleotide metabolism) mutations detectably alter the SOS response elicited in new transconjugants upon acquisition of pED208. However, MC4100(pED208) donor cells intrinsically exhibit significantly higher SOS activation than plasmid-free MC4100 cells, and this plasmid carriage-induced stress response is further elevated in donor cells deleted of several Group 1 genes. Among 10 characterized substrates, we gained evidence of C-terminal or internal translocation signals that could function independently or synergistically for optimal protein transfer. Remarkably, nearly all tested proteins were also translocated through the IncN pKM101 and IncP RP4 conjugation systems. This repertoire of E. coli protein substrates, here termed the F plasmid "conjutome," is thus characterized by functions of potential benefit to new transconjugants, diverse TSs, and the capacity for promiscuous transfer through heterologous conjugation systems.

IMPORTANCE

Conjugation systems comprise a major subfamily of the type IV secretion systems (T4SSs) and are the progenitors of a second large T4SS subfamily dedicated to translocation of protein effectors. This study examined the capacity of conjugation machines to function as protein translocators. Using a high-throughput reporter screen, we determined that 32 chromosomally encoded proteins are delivered through an F plasmid conjugation system. The translocated proteins potentially enhance the establishment of the co-transferred F plasmid or mitigate mating-induced stresses. Translocation signals located C-terminally or internally conferred substrate recognition by the F system and, remarkably, many substrates also were translocated through heterologous conjugation systems. Our findings highlight the plasticity of conjugation systems in their capacities to co-translocate DNA and many protein substrates.

From structural design to delivery: mRNA therapeutics for cancer immunotherapy

mRNA therapeutics have emerged as powerful tools for cancer immunotherapy in accordance with their superiority in expressing all sequence-known proteins in vivo. In particular, with a small dosage of delivered mRNA, antigen-presenting cells (APCs) can synthesize mutant neo-antigens and multi-antigens and present epitopes to T lymphocytes to elicit antitumor effects. In addition, expressing receptors like chimeric antigen receptor (CAR), T-cell receptor (TCR), CD134, and immune-modulating factors including cytokines, interferons, and antibodies in specific cells can enhance immunological response against tumors. With the maturation of in vitro transcription (IVT) technology, large-scale and pure mRNA encoding specific proteins can be synthesized quickly. However, the clinical translation of mRNA-based anticancer strategies is restricted by delivering mRNA into target organs or cells and the inadequate endosomal escape efficiency of mRNA. Recently, there have been some advances in mRNA-based cancer immunotherapy, which can be roughly classified as modifications of the mRNA structure and the development of delivery systems, especially the lipid nanoparticle platforms. In this review, the latest strategies for overcoming the limitations of mRNA-based cancer immunotherapies and the recent advances in delivering mRNA into specific organs and cells are summarized. Challenges and opportunities for clinical applications of mRNA-based cancer immunotherapy are also discussed.

Epistasis decreases with increasing antibiotic pressure but not temperature

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Hfq protein and GcvB small RNA tailoring of oppA target mRNA to levels allowing translation activation by MicF small RNA in Escherichia coli

Traffic of molecules across the bacterial membrane mainly relies on porins and transporters, whose expression must adapt to environmental conditions. To ensure bacterial fitness, synthesis and assembly of functional porins and transporters are regulated through a plethora of mechanisms. Among them, small regulatory RNAs (sRNAs) are known to be powerful post-transcriptional regulators. In Escherichia coli, the MicF sRNA is known to regulate only four targets, a very narrow targetome for a sRNA responding to various stresses, such as membrane stress, osmotic shock, or thermal shock. Using an in vivo pull-down assay combined with high-throughput RNA sequencing, we sought to identify new targets of MicF to better understand its role in the maintenance of cellular homoeostasis. Here, we report the first positively regulated target of MicF, the oppA mRNA. The OppA protein is the periplasmic component of the Opp ATP-binding cassette (ABC) oligopeptide transporter and regulates the import of short peptides, some of them bactericides. Mechanistic studies suggest that oppA translation is activated by MicF through a mechanism of action involving facilitated access to a translation-enhancing region in oppA 5'UTR. Intriguingly, MicF activation of oppA translation depends on cross-regulation by negative trans-acting effectors, the GcvB sRNA and the RNA chaperone protein Hfq.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Translational remodeling by RNA-binding proteins and noncoding RNAs

Responsible for generating the proteome that controls phenotype, translation is the ultimate convergence point for myriad upstream signals that influence gene expression. System-wide adaptive translational reprogramming has recently emerged as a pillar of cellular adaptation. As classic regulators of mRNA stability and translation efficiency, foundational studies established the concept of collaboration and competition between RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) on individual mRNAs. Fresh conceptual innovations now highlight stress-activated, evolutionarily conserved RBP networks and ncRNAs that increase the translation efficiency of populations of transcripts encoding proteins that participate in a common cellular process. The discovery of post-transcriptional functions for long noncoding RNAs (lncRNAs) was particularly intriguing given their cell-type-specificity and historical definition as nuclear-functioning epigenetic regulators. The convergence of RBPs, lncRNAs, and microRNAs on functionally related mRNAs to enable adaptive protein synthesis is a newer biological paradigm that highlights their role as "translatome (protein output) remodelers" and reinvigorates the paradigm of "RNA operons." Together, these concepts modernize our understanding of cellular stress adaptation and strategies for therapeutic development. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Translation Regulation Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

OTTER, a new method quantifying absolute amounts of tRNAs

To maintain optimal proteome, both codon choice of each mRNA and supply of aminoacyl-tRNAs are two principal factors in translation. Recent reports have revealed that the amounts of tRNAs in cells are more dynamic than we had expected. High-throughput methods such as RNA-Seq and microarrays are versatile for comprehensive detection of changes in individual tRNA amounts, but they suffer from inability to assess signal production efficiencies of individual tRNA species. Thus, they are not the perfect choice to measure absolute amounts of tRNAs. Here, we introduce a novel method for this purpose, termed Oligonucleotide-directed Three-prime Terminal Extension of RNA (OTTER), which employs fluorescence-labeling at the 3'-terminus of a tRNA by optimized reverse primer extension and an assessment step of each labeling efficiency by northern blotting. Using this method, we quantified the absolute amounts of the 34 individual and 4 pairs of isoacceptor tRNAs out of the total 42 nuclear-encoded isoacceptors in the yeast Saccharomyces cerevisiae. We found that the amounts of tRNAs in log phase yeast cells grown in a rich glucose medium range from 0.030 to 0.73 pmol/µg RNA. The tRNA amounts seem to be altered at the isoacceptor level by a few folds in response to physiological growing conditions. The data obtained by OTTER are poorly correlated with those by simple RNA-Seq, marginally with those by microarrays and by microscale thermophoresis. However, the OTTER data showed good agreement with the data obtained by 2D-gel analysis of in vivo radiolabeled RNAs. Thus, OTTER is a suitable method for quantifying absolute amounts of tRNAs at the level of isoacceptor resolution.

Protein structure, amino acid composition and sequence determine proteome vulnerability to oxidation-induced damage

Oxidative stress alters cell viability, from microorganism irradiation sensitivity to human aging and neurodegeneration. Deleterious effects of protein carbonylation by reactive oxygen species (ROS) make understanding molecular properties determining ROS susceptibility essential. The radiation‐resistant bacterium Deinococcus radiodurans accumulates less carbonylation than sensitive organisms, making it a key model for deciphering properties governing oxidative stress resistance. We integrated shotgun redox proteomics, structural systems biology, and machine learning to resolve properties determining protein damage by γ‐irradiation in Escherichia coli and D. radiodurans at multiple scales. Local accessibility, charge, and lysine enrichment accurately predict ROS susceptibility. Lysine, methionine, and cysteine usage also contribute to ROS resistance of the D. radiodurans proteome. Our model predicts proteome maintenance machinery, and proteins protecting against ROS are more resistant in D. radiodurans. Our findings substantiate that protein‐intrinsic protection impacts oxidative stress resistance, identifying causal molecular properties.

Short- and Long-Term Transcriptomic Responses of Escherichia coli to Biocides: a Systems Analysis

The mechanisms of the bacterial response to biocides are poorly understood, despite their broad application. To identify the genetic basis and pathways implicated in the biocide stress response, we exposed Escherichia coli populations to 10 ubiquitous biocides. By comparing the transcriptional responses between a short-term exposure (30 min) and a long-term exposure (8 to 12 h) to biocide stress, we established the common gene and pathway clusters that are implicated in general and biocide-specific stress responses. Our analysis revealed a temporal choreography, starting from the upregulation of chaperones to the subsequent repression of motility and chemotaxis pathways and the induction of an anaerobic pool of enzymes and biofilm regulators. A systematic analysis of the transcriptional data identified a zur-regulated gene cluster to be highly active in the stress response against sodium hypochlorite and peracetic acid, presenting a link between the biocide stress response and zinc homeostasis. Susceptibility assays with knockout mutants further validated our findings and provide clear targets for downstream investigation of the implicated mechanisms of action.IMPORTANCE Antiseptics and disinfectant products are of great importance to control and eliminate pathogens, especially in settings such as hospitals and the food industry. Such products are widely distributed and frequently poorly regulated. Occasional outbreaks have been associated with microbes resistant to such compounds, and researchers have indicated potential cross-resistance with antibiotics. Despite that, there are many gaps in knowledge about the bacterial stress response and the mechanisms of microbial resistance to antiseptics and disinfectants. We investigated the stress response of the bacterium Escherichia coli to 10 common disinfectant and antiseptic chemicals to shed light on the potential mechanisms of tolerance to such compounds.

Mechanistic insights into UV-A mediated bacterial disinfection via endogenous photosensitizers

UV-A and visible light are thought to excite endogenous photosensitizers in microbes, thereby initiating complex chemical interactions that ultimately kill cells. Natural solar-based disinfection methods have been adapted into commercial lighting technologies with varying degrees of reported efficacy and associated safety hazards for human exposure. Here we utilize a narrow-spectrum UV-A LED prototype (currently in development for health care applications) to investigate the mechanism of bacterial photoinactivation using 365 nm light. Using a combination of reverse genetics and biochemical investigation, we report mechanistic evidence that 365nm light initiates a chain-reaction of superoxide-mediated damage via auto-excitation of vitamin-based electron carriers, specifically vitamin K2 menaquinones and the FAD flavoprotein in Complex II in the electron transport chain. We observe that photoinactivation is modifiable through supplementation of the environment to bypass cell damage. Lastly, we observe that bacteria forced into metabolic dormancy by desiccation become hypersensitized to the effects of UV-A light, thereby permitting photoinactivation at fluences that are significantly lower than the industry threshold for safe human exposure. In total, these results substantiate the mechanism and potential application of narrow- spectrum UV-A light for bacterial disinfection purposes.

Genes and Proteomes Associated With Increased Mutation Frequency and Multidrug Resistance of Naturally Occurring Mismatch Repair-Deficient Salmonella Hypermutators

The emergence of antibiotic-resistant Salmonella through mutations led to mismatch repair (MMR) deficiency that represents a potential hazard to public health. Here, four representative MMR-deficient Salmonella hypermutator strains and Salmonella Typhimurium LT2 were used to comprehensively reveal the influence of MMR deficiency on antibiotic resistance among Salmonella. Our results indicated that the mutation frequency ranged from 3.39 × 10^–4 to 5.46 × 10^–2 in the hypermutator. Mutation sites in MutS, MutL, MutT, and UvrD of the four hypermutators were all located in the essential and core functional regions. Mutation frequency of the hypermutator was most highly correlated with the extent of mutation in MutS. Mutations in MMR genes (mutS, mutT, mutL, and uvrD) were correlated with increased mutation in antibiotic resistance genes, and the extent of antibiotic resistance was significantly correlated with the number of mutation sites in MutL and in ParC. The number of mutation sites in MMR genes and antibiotic resistance genes exhibited a significant positive correlation with the number of antibiotics resisted and with expression levels of mutS, mutT, and mutL. Compared to Salmonella Typhimurium LT2, a total of 137 differentially expressed and 110 specifically expressed proteins were identified in the four hypermutators. Functional enrichment analysis indicated that the proteins significantly overexpressed in the hypermutators primarily associated with translation and stress response. Interaction network analysis revealed that the ribosome pathway might be a critical factor for high mutation frequency and multidrug resistance in MMR-deficient Salmonella hypermutators. These results help elucidate the mutational dynamics that lead to hypermutation, antibiotic resistance, and activation of stress response pathways in Salmonella.

Variation and selection on codon usage bias across an entire subphylum

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

Overlap of Promoter Recognition Specificity of Stress Response Sigma Factors SigD and SigH in Corynebacterium glutamicum ATCC 13032

Corynebacterium glutamicum ATCC 13032 harbors five sigma subunits of RNA polymerase belonging to Group IV, also called extracytoplasmic function (ECF) σ factors. These factors σ^C, σ^D, σ^E, σ^H, and σ^M are mostly involved in stress responses. The role of σ^D consists in the control of cell wall integrity. The σ^D regulon is involved in the synthesis of components of the mycomembrane which is part of the cell wall in C. glutamicum. RNA sequencing of the transcriptome from a strain overexpressing the sigD gene provided 29 potential σ^D-controlled genes and enabled us to precisely localize their transcriptional start sites. Analysis of the respective promoters by both in vitro transcription and the in vivo two-plasmid assay confirmed that transcription of 11 of the tested genes is directly σ^D-dependent. The key sequence elements of all these promoters were found to be identical or closely similar to the motifs -35 GTAAC^A/_G and -10 GAT. Surprisingly, nearly all of these σ^D-dependent promoters were also active to a much lower extent with σ^Hin vivo and one (Pcg0607) also in vitro, although the known highly conserved consensus sequence of the σ^H-dependent promoters is different (-35 GGAA^T/_C and -10 GTT). In addition to the activity of σ^H at the σ^D-controlled promoters, we discovered separated or overlapping σ^A- or σ^B-regulated or σ^H-regulated promoters within the upstream region of 8 genes of the σ^D-regulon. We found that phenol in the cultivation medium acts as a stress factor inducing expression of some σ^D-dependent genes. Computer modeling revealed that σ^H binds to the promoter DNA in a similar manner as σ^D to the analogous promoter elements. The homology models together with mutational analysis showed that the key amino acids, Ala 60 in σ^D and Lys 53 in σ^H, bind to the second nucleotide within the respective -10 promoter elements (GAT and GTT, respectively). The presented data obtained by integrating in vivo, in vitro and in silico approaches demonstrate that most of the σ^D-controlled genes also belong to the σ^H-regulon and are also transcribed from the overlapping or closely located housekeeping (σ^A-regulated) and/or general stress (σ^B-regulated) promoters.

Can Protein Expression Be Regulated by Modulation of tRNA Modification Profiles?

tRNAs are the central adaptor molecules in translation. Their decoding properties are influenced by post-transcriptional modifications, particularly in the critical anticodon-stem-loop (ASL) region. Synonymous codon choice, also called codon usage bias, affects both translation efficiency and accuracy, and ASL modifications play key roles in both of these processes. In combination with a handful of historical examples, recent studies integrating ribosome profiling, proteomics, codon-usage analyses, and modification quantifications show that levels of tRNA modifications can change under stress, during development, or under specific metabolic conditions and can modulate the expression of specific genes. Deconvoluting the different responses (global or specific) to tRNA modification deficiencies can be difficult because of pleiotropic effects, but, as more cases emerge, it does seem that tRNA modification changes could add another layer of regulation in the transfer of information from DNA to protein.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

The evolutionary signal in metagenome phyletic profiles predicts many gene functions

BACKGROUND

The function of many genes is still not known even in model organisms. An increasing availability of microbiome DNA sequencing data provides an opportunity to infer gene function in a systematic manner.

RESULTS

We evaluated if the evolutionary signal contained in metagenome phyletic profiles (MPP) is predictive of a broad array of gene functions. The MPPs are an encoding of environmental DNA sequencing data that consists of relative abundances of gene families across metagenomes. We find that such MPPs can accurately predict 826 Gene Ontology functional categories, while drawing on human gut microbiomes, ocean metagenomes, and DNA sequences from various other engineered and natural environments. Overall, in this task, the MPPs are highly accurate, and moreover they provide coverage for a set of Gene Ontology terms largely complementary to standard phylogenetic profiles, derived from fully sequenced genomes. We also find that metagenomes approximated from taxon relative abundance obtained via 16S rRNA gene sequencing may provide surprisingly useful predictive models. Crucially, the MPPs derived from different types of environments can infer distinct, non-overlapping sets of gene functions and therefore complement each other. Consistently, simulations on > 5000 metagenomes indicate that the amount of data is not in itself critical for maximizing predictive accuracy, while the diversity of sampled environments appears to be the critical factor for obtaining robust models.

CONCLUSIONS

In past work, metagenomics has provided invaluable insight into ecology of various habitats, into diversity of microbial life and also into human health and disease mechanisms. We propose that environmental DNA sequencing additionally constitutes a useful tool to predict biological roles of genes, yielding inferences out of reach for existing comparative genomics approaches.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Lost in Translation: Defects in Transfer RNA Modifications and Neurological Disorders

Transfer RNAs (tRNAs) are key molecules participating in protein synthesis. To augment their functionality they undergo extensive post-transcriptional modifications and, as such, are subject to regulation at multiple levels including transcription, transcript processing, localization and ribonucleoside base modification. Post-transcriptional enzyme-catalyzed modification of tRNA occurs at a number of base and sugar positions and influences specific anticodon-codon interactions and regulates translation, its efficiency and fidelity. This phenomenon of nucleoside modification is most remarkable and results in a rich structural diversity of tRNA of which over 100 modified nucleosides have been characterized. Most often these hypermodified nucleosides are found in the wobble position of tRNAs, where they play a direct role in codon recognition as well as in maintaining translational efficiency and fidelity, etc. Several recent studies have pointed to a link between defects in tRNA modifications and human diseases including neurological disorders. Therefore, defects in tRNA modifications in humans need intensive characterization at the enzymatic and mechanistic level in order to pave the way to understand how lack of such modifications are associated with neurological disorders with the ultimate goal of gaining insights into therapeutic interventions.

Codon optimization underpins generalist parasitism in fungi

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

Multiple Osmotic Stress Responses in Acidihalobacter prosperus Result in Tolerance to Chloride Ions

Extremely acidophilic microorganisms (pH optima for growth of ≤3) are utilized for the extraction of metals from sulfide minerals in the industrial biotechnology of “biomining.” A long term goal for biomining has been development of microbial consortia able to withstand increased chloride concentrations for use in regions where freshwater is scarce. However, when challenged by elevated salt, acidophiles experience both osmotic stress and an acidification of the cytoplasm due to a collapse of the inside positive membrane potential, leading to an influx of protons. In this study, we tested the ability of the halotolerant acidophile Acidihalobacter prosperus to grow and catalyze sulfide mineral dissolution in elevated concentrations of salt and identified chloride tolerance mechanisms in Ac. prosperus as well as the chloride susceptible species, Acidithiobacillus ferrooxidans. Ac. prosperus had optimum iron oxidation at 20 g L⁻¹ NaCl while At. ferrooxidans iron oxidation was inhibited in the presence of 6 g L⁻¹ NaCl. The tolerance to chloride in Ac. prosperus was consistent with electron microscopy, determination of cell viability, and bioleaching capability. The Ac. prosperus proteomic response to elevated chloride concentrations included the production of osmotic stress regulators that potentially induced production of the compatible solute, ectoine uptake protein, and increased iron oxidation resulting in heightened electron flow to drive proton export by the F₀F₁ ATPase. In contrast, At. ferrooxidans responded to low levels of Cl⁻ with a generalized stress response, decreased iron oxidation, and an increase in central carbon metabolism. One potential adaptation to high chloride in the Ac. prosperus Rus protein involved in ferrous iron oxidation was an increase in the negativity of the surface potential of Rus Form I (and Form II) that could help explain how it can be active under elevated chloride concentrations. These data have been used to create a model of chloride tolerance in the salt tolerant and susceptible species Ac. prosperus and At. ferrooxidans, respectively.

Dietary nitrogen alters codon bias and genome composition in parasitic microorganisms

BACKGROUND

Genomes are composed of long strings of nucleotide monomers (A, C, G and T) that are either scavenged from the organism's environment or built from metabolic precursors. The biosynthesis of each nucleotide differs in atomic requirements with different nucleotides requiring different quantities of nitrogen atoms. However, the impact of the relative availability of dietary nitrogen on genome composition and codon bias is poorly understood.

RESULTS

Here we show that differential nitrogen availability, due to differences in environment and dietary inputs, is a major determinant of genome nucleotide composition and synonymous codon use in both bacterial and eukaryotic microorganisms. Specifically, low nitrogen availability species use nucleotides that require fewer nitrogen atoms to encode the same genes compared to high nitrogen availability species. Furthermore, we provide a novel selection-mutation framework for the evaluation of the impact of metabolism on gene sequence evolution and show that it is possible to predict the metabolic inputs of related organisms from an analysis of the raw nucleotide sequence of their genes.

CONCLUSIONS

Taken together, these results reveal a previously hidden relationship between cellular metabolism and genome evolution and provide new insight into how genome sequence evolution can be influenced by adaptation to different diets and environments.

The landscape of microbial phenotypic traits and associated genes

Bacteria and Archaea display a variety of phenotypic traits and can adapt to diverse ecological niches. However, systematic annotation of prokaryotic phenotypes is lacking. We have therefore developed ProTraits, a resource containing ∼545 000 novel phenotype inferences, spanning 424 traits assigned to 3046 bacterial and archaeal species. These annotations were assigned by a computational pipeline that associates microbes with phenotypes by text-mining the scientific literature and the broader World Wide Web, while also being able to define novel concepts from unstructured text. Moreover, the ProTraits pipeline assigns phenotypes by drawing extensively on comparative genomics, capturing patterns in gene repertoires, codon usage biases, proteome composition and co-occurrence in metagenomes. Notably, we find that gene synteny is highly predictive of many phenotypes, and highlight examples of gene neighborhoods associated with spore-forming ability. A global analysis of trait interrelatedness outlined clusters in the microbial phenotype network, suggesting common genetic underpinnings. Our extended set of phenotype annotations allows detection of 57 088 high confidence gene-trait links, which recover many known associations involving sporulation, flagella, catalase activity, aerobicity, photosynthesis and other traits. Over 99% of the commonly occurring gene families are involved in genetic interactions conditional on at least one phenotype, suggesting that epistasis has a major role in shaping microbial gene content.

Extensive complementarity between gene function prediction methods

MOTIVATION

The number of sequenced genomes rises steadily but we still lack the knowledge about the biological roles of many genes. Automated function prediction (AFP) is thus a necessity. We hypothesized that AFP approaches that draw on distinct genome features may be useful for predicting different types of gene functions, motivating a systematic analysis of the benefits gained by obtaining and integrating such predictions.

RESULTS

Our pipeline amalgamates 5 133 543 genes from 2071 genomes in a single massive analysis that evaluates five established genomic AFP methodologies. While 1227 Gene Ontology (GO) terms yielded reliable predictions, the majority of these functions were accessible to only one or two of the methods. Moreover, different methods tend to assign a GO term to non-overlapping sets of genes. Thus, inferences made by diverse genomic AFP methods display a striking complementary, both gene-wise and function-wise. Because of this, a viable integration strategy is to rely on a single most-confident prediction per gene/function, rather than enforcing agreement across multiple AFP methods. Using an information-theoretic approach, we estimate that current databases contain 29.2 bits/gene of known Escherichia coli gene functions. This can be increased by up to 5.5 bits/gene using individual AFP methods or by 11 additional bits/gene upon integration, thereby providing a highly-ranking predictor on the Critical Assessment of Function Annotation 2 community benchmark. Availability of more sequenced genomes boosts the predictive accuracy of AFP approaches and also the benefit from integrating them.

AVAILABILITY AND IMPLEMENTATION

The individual and integrated GO predictions for the complete set of genes are available from http://gorbi.irb.hr/ CONTACT: fran.supek@irb.hrSupplementary information: Supplementary materials are available at Bioinformatics online.

Nucleoside modifications in the regulation of gene expression: focus on tRNA

Both, DNA and RNA nucleoside modifications contribute to the complex multi-level regulation of gene expression. Modified bases in tRNAs modulate protein translation rates in a highly dynamic manner. Synonymous codons, which differ by the third nucleoside in the triplet but code for the same amino acid, may be utilized at different rates according to codon-anticodon affinity. Nucleoside modifications in the tRNA anticodon loop can favor the interaction with selected codons by stabilizing specific base pairs. Similarly, weakening of base pairing can discriminate against binding to near-cognate codons. mRNAs enriched in favored codons are translated in higher rates constituting a fine-tuning mechanism for protein synthesis. This so-called codon bias establishes a basic protein level, but sometimes it is necessary to further adjust the production rate of a particular protein to actual requirements, brought by, e.g., stages in circadian rhythms, cell cycle progression or exposure to stress. Such an adjustment is realized by the dynamic change of tRNA modifications resulting in the preferential translation of mRNAs coding for example for stress proteins to facilitate cell survival. Furthermore, tRNAs contribute in an entirely different way to another, less specific stress response consisting in modification-dependent tRNA cleavage that contributes to the general down-regulation of protein synthesis. In this review, we summarize control functions of nucleoside modifications in gene regulation with a focus on recent findings on protein synthesis control by tRNA base modifications.

Transcriptomic analysis reveals how a lack of potassium ions increases Sulfolobus acidocaldarius sensitivity to pH changes

Extremely acidophilic microorganisms (optimum growth pH of ≤3) maintain a near neutral cytoplasmic pH via several homeostatic mechanisms, including an inside positive membrane potential created by potassium ions. Transcriptomic responses to pH stress in the thermoacidophilic archaeon, Sulfolobus acidocaldarius were investigated by growing cells without added sodium and/or potassium ions at both optimal and sub-optimal pH. Culturing the cells in the absence of added sodium or potassium ions resulted in a reduced growth rate compared to full-salt conditions as well as 43 and 75 significantly different RNA transcript ratios, respectively. Differentially expressed RNA transcripts during growth in the absence of added sodium ions included genes coding for permeases, a sodium/proline transporter and electron transport proteins. In contrast, culturing without added potassium ions resulted in higher RNA transcripts for similar genes as a lack of sodium ions plus genes related to spermidine that has a general role in response to stress and a decarboxylase that potentially consumes protons. The greatest RNA transcript response occurred when S. acidocaldarius cells were grown in the absence of potassium and/or sodium at a sub-optimal pH. These adaptations included those listed above plus osmoregulated glucans and mechanosensitive channels that have previously been shown to respond to osmotic stress. In addition, data analyses revealed two co-expressed IclR family transcriptional regulator genes with a previously unknown role in the S. acidocaldarius pH stress response. Our study provides additional evidence towards the importance of potassium in acidophile growth at acidic pH.

Pausing on Polyribosomes: Make Way for Elongation in Translational Control

Among the three phases of mRNA translation-initiation, elongation, and termination-initiation has traditionally been considered to be rate limiting and thus the focus of regulation. Emerging evidence, however, demonstrates that control of ribosome translocation (polypeptide elongation) can also be regulatory and indeed exerts a profound influence on development, neurologic disease, and cell stress. The correspondence of mRNA codon usage and the relative abundance of their cognate tRNAs is equally important for mediating the rate of polypeptide elongation. Here, we discuss recent results showing that ribosome pausing is a widely used mechanism for controlling translation and, as a result, biological transitions in health and disease.

Global translational impacts of the loss of the tRNA modification t6A in yeast

The universal tRNA modification t⁶A is found at position 37 of nearly all tRNAs decoding ANN codons. The absence of t⁶A₃₇ leads to severe growth defects in baker’s yeast, phenotypes similar to those caused by defects in mcm⁵s²U₃₄ synthesis. Mutants in mcm⁵s²U₃₄ can be suppressed by overexpression of tRNA^Lys_UUU, but we show t⁶A phenotypes could not be suppressed by expressing any individual ANN decoding tRNA, and t⁶A and mcm⁵s²U are not determinants for each other’s formation. Our results suggest that t⁶A deficiency, like mcm⁵s²U deficiency, leads to protein folding defects, and show that the absence of t⁶A led to stress sensitivities (heat, ethanol, salt) and sensitivity to TOR pathway inhibitors. Additionally, L-homoserine suppressed the slow growth phenotype seen in t⁶A-deficient strains, and proteins aggregates and Advanced Glycation End-products (AGEs) were increased in the mutants. The global consequences on translation caused by t⁶A absence were examined by ribosome profiling. Interestingly, the absence of t⁶A did not lead to global translation defects, but did increase translation initiation at upstream non-AUG codons and increased frame-shifting in specific genes. Analysis of codon occupancy rates suggests that one of the major roles of t⁶A is to homogenize the process of elongation by slowing the elongation rate at codons decoded by high abundance tRNAs and I₃₄:C₃ pairs while increasing the elongation rate of rare tRNAs and G₃₄:U₃ pairs. This work reveals that the consequences of t⁶A absence are complex and multilayered and has set the stage to elucidate the molecular basis of the observed phenotypes.

The Code of Silence: Widespread Associations Between Synonymous Codon Biases and Gene Function

Some mutations in gene coding regions exchange one synonymous codon for another, and thus do not alter the amino acid sequence of the encoded protein. Even though they are often called 'silent,' these mutations may exhibit a plethora of effects on the living cell. Therefore, they are often selected during evolution, causing synonymous codon usage biases in genomes. Comparative analyses of bacterial, archaeal, fungal, and human cancer genomes have found many links between a gene's biological role and the accrual of synonymous mutations during evolution. In particular, highly expressed genes in certain functional categories are enriched with optimal codons, which are decoded by the abundant tRNAs, thus enhancing the speed and accuracy of the translating ribosome. The set of genes exhibiting codon adaptation differs between genomes, and these differences show robust associations to organismal phenotypes. In addition to selection for translation efficiency, other distinct codon bias patterns have been found in: amino acid starvation genes, cyclically expressed genes, tissue-specific genes in animals and plants, oxidative stress response genes, cellular differentiation genes, and oncogenes. In addition, genomes of organisms harboring tRNA modifications exhibit particular codon preferences. The evolutionary trace of codon bias patterns across orthologous genes may be examined to learn about a gene's relevance to various phenotypes, or, more generally, its function in the cell.

Predictive analytics of environmental adaptability in multi-omic network models

Bacterial phenotypic traits and lifestyles in response to diverse environmental conditions depend on changes in the internal molecular environment. However, predicting bacterial adaptability is still difficult outside of laboratory controlled conditions. Many molecular levels can contribute to the adaptation to a changing environment: pathway structure, codon usage, metabolism. To measure adaptability to changing environmental conditions and over time, we develop a multi-omic model of Escherichia coli that accounts for metabolism, gene expression and codon usage at both transcription and translation levels. After the integration of multiple omics into the model, we propose a multiobjective optimization algorithm to find the allowable and optimal metabolic phenotypes through concurrent maximization or minimization of multiple metabolic markers. In the condition space, we propose Pareto hypervolume and spectral analysis as estimators of short term multi-omic (transcriptomic and metabolic) evolution, thus enabling comparative analysis of metabolic conditions. We therefore compare, evaluate and cluster different experimental conditions, models and bacterial strains according to their metabolic response in a multidimensional objective space, rather than in the original space of microarray data. We finally validate our methods on a phenomics dataset of growth conditions. Our framework, named METRADE, is freely available as a MATLAB toolbox.

Global Shifts in Genome and Proteome Composition Are Very Tightly Coupled

The amino acid composition (AAC) of proteomes differs greatly between microorganisms and is associated with the environmental niche they inhabit, suggesting that these changes may be adaptive. Similarly, the oligonucleotide composition of genomes varies and may confer advantages at the DNA/RNA level. These influences overlap in protein-coding sequences, making it difficult to gauge their relative contributions. We disentangle these effects by systematically evaluating the correspondence between intergenic nucleotide composition, where protein-level selection is absent, the AAC, and ecological parameters of 909 prokaryotes. We find that G + C content, the most frequently used measure of genomic composition, cannot capture diversity in AAC and across ecological contexts. However, di-/trinucleotide composition in intergenic DNA predicts amino acid frequencies of proteomes to the point where very little cross-species variability remains unexplained (91% of variance accounted for). Qualitatively similar results were obtained for 49 fungal genomes, where 80% of the variability in AAC could be explained by the composition of introns and intergenic regions. Upon factoring out oligonucleotide composition and phylogenetic inertia, the residual AAC is poorly predictive of the microbes' ecological preferences, in stark contrast with the original AAC. Moreover, highly expressed genes do not exhibit more prominent environment-related AAC signatures than lowly expressed genes, despite contributing more to the effective proteome. Thus, evolutionary shifts in overall AAC appear to occur almost exclusively through factors shaping the global oligonucleotide content of the genome. We discuss these results in light of contravening evidence from biophysical data and further reading frame-specific analyses that suggest that adaptation takes place at the protein level.

Codon optimality is a major determinant of mRNA stability

mRNA degradation represents a critical regulated step in gene expression. Although the major pathways in turnover have been identified, accounting for disparate half-lives has been elusive. We show that codon optimality is one feature that contributes greatly to mRNA stability. Genome-wide RNA decay analysis revealed that stable mRNAs are enriched in codons designated optimal, whereas unstable mRNAs contain predominately non-optimal codons. Substitution of optimal codons with synonymous, non-optimal codons results in dramatic mRNA destabilization, whereas the converse substitution significantly increases stability. Further, we demonstrate that codon optimality impacts ribosome translocation, connecting the processes of translation elongation and decay through codon optimality. Finally, we show that optimal codon content accounts for the similar stabilities observed in mRNAs encoding proteins with coordinated physiological function. This work demonstrates that codon optimization exists as a mechanism to finely tune levels of mRNAs and, ultimately, proteins.

Homology-independent metrics for comparative genomics

A mainstream procedure to analyze the wealth of genomic data available nowadays is the detection of homologous regions shared across genomes, followed by the extraction of biological information from the patterns of conservation and variation observed in such regions. Although of pivotal importance, comparative genomic procedures that rely on homology inference are obviously not applicable if no homologous regions are detectable. This fact excludes a considerable portion of “genomic dark matter” with no significant similarity — and, consequently, no inferred homology to any other known sequence — from several downstream comparative genomic methods. In this review we compile several sequence metrics that do not rely on homology inference and can be used to compare nucleotide sequences and extract biologically meaningful information from them. These metrics comprise several compositional parameters calculated from sequence data alone, such as GC content, dinucleotide odds ratio, and several codon bias metrics. They also share other interesting properties, such as pervasiveness (patterns persist on smaller scales) and phylogenetic signal. We also cite examples where these homology-independent metrics have been successfully applied to support several bioinformatics challenges, such as taxonomic classification of biological sequences without homology inference. They where also used to detect higher-order patterns of interactions in biological systems, ranging from detecting coevolutionary trends between the genomes of viruses and their hosts to characterization of gene pools of entire microbial communities. We argue that, if correctly understood and applied, homology-independent metrics can add important layers of biological information in comparative genomic studies without prior homology inference.