Differential codon adaptation between dsDNA and ssDNA phages in Escherichia coli

Differential Selection for Translation Efficiency Shapes Translation Machineries in Bacterial Species

Different bacterial species have dramatically different generation times, from 20-30 min in Escherichia coli to about two weeks in Mycobacterium leprae. The translation machinery in a cell needs to synthesize all proteins for a new cell in each generation. The three subprocesses of translation, i.e., initiation, elongation, and termination, are expected to be under stronger selection pressure to optimize in short-generation bacteria (SGB) such as Vibrio natriegens than in the long-generation Mycobacterium leprae. The initiation efficiency depends on the start codon decoded by the initiation tRNA, the optimal Shine-Dalgarno (SD) decoded by the anti-SD (aSD) sequence on small subunit rRNA, and the secondary structure that may embed the initiation signals and prevent them from being decoded. The elongation efficiency depends on the tRNA pool and codon usage. The termination efficiency in bacteria depends mainly on the nature of the stop codon and the nucleotide immediately downstream of the stop codon. By contrasting SGB with long-generation bacteria (LGB), we predict (1) SGB to have more ribosome RNA operons to produce ribosomes, and more tRNA genes for carrying amino acids to ribosomes, (2) SGB to have a higher percentage of genes using AUG as the start codon and UAA as the stop codon than LGB, (3) SGB to exhibit better codon and anticodon adaptation than LGB, and (4) SGB to have a weaker secondary structure near the translation initiation signals than LGB. These differences between SGB and LGB should be more pronounced in highly expressed genes than the rest of the genes. We present empirical evidence in support of these predictions.

Systematic analysis of prophages carried by Porphyromonas gingivalis

To systematically investigate the prophages carrying in Porphyromonas gingivalis (P. gingivalis) strains, analyze potential antibiotic resistance genes (ARGs) and virulence genes in these prophages. We collected 90 whole genome sequences of P. gingivalis from NCBI and utilized the Prophage Hunter online software to predict prophages; Comprehensive antibiotic research database (CARD) and virulence factors database (VFDB) were adopted to analyze the ARGs and virulence factors (VFs) carried by the prophages. Sixty-nine prophages were identified among 24/90 P. gingivalis strains, including 17 active prophages (18.9%) and 52 ambiguous prophages (57.8%). The proportion of prophages carried by each P. gingivalis genome ranged from 0.5% to 6.7%. A total of 188 antibiotic resistance genes belonging to 25 phenotypes and 46 different families with six mechanisms of antibiotic resistance were identified in the 17 active prophages. Three active prophages encoded 4 virulence genes belonging to type III and type VI secretion systems. The potential hosts of these virulence genes included Escherichia coli, Shigella sonnei, Salmonella typhi, and Klebsiella pneumoniae. In conclusion, 26.7% P. gingivalis strains carry prophages, while the proportion of prophage genes in the P. gingivalis genome is relatively low. In addition, approximately 39.7% of the P. gingivalis prophage genes have ARGs identified, mainly against streptogramin, peptides, and aminoglycosides. Only a few prophages carry virulence genes. Prophages may play an important role in the acquisition, dissemination of antibiotic resistance genes, and pathogenicity evolution in P. gingivalis.

A Novel Alteromonas Phage Lineage with a Broad Host Range and Small Burst Size

Alteromonas is an opportunistic marine bacterium that persists in the global ocean and has important ecological significance. However, current knowledge about the diversity and ecology of alterophages (phages that infect Alteromonas) is lacking. Here, three similar phages infecting Alteromonas macleodii ATCC 27126^T were isolated and physiologically characterized. Transmission electron microscopy revealed Siphoviridae morphology, with an oblate icosahedral head and a long noncontractile tail. Notably, these members displayed a small burst size (15–19 plaque-forming units/cell) yet an extensively broad host spectrum when tested on 175 Alteromonas strains. Such unique infection kinetics are potentially associated with discrepancies in codon usage bias from the host tRNA inventory. Phylogenetic analysis indicated that the three phages are closely evolutionarily related; they clustered at the species level and represent a novel genus. Three auxiliary metabolic genes with roles in nucleotide metabolism and putative biofilm dispersal were found in these phage genomes, which revealed important biogeochemical significance of these alterophages in marine ecosystems. Our isolation and characterization of these novel phages expand the current understanding of alterophage diversity, evolution, and phage–host interactions.

IMPORTANCE The marine bacterium Alteromonas is prevalent in the global ocean with crucial ecological significance; however, little is known about the diversity and evolution of its bacteriophages that profoundly affect the bacterial communities. Our study characterized a novel genus of three newly isolated Alteromonas phages that exhibited a distinct infection strategy of broad host spectrum and small burst size. This strategy is likely a consequence of the viral trade-off between virulence and lysis profiles during phage–host coevolution, and our work provides new insight into viral evolution and infection strategies.

GPR120 modulates epileptic seizure and neuroinflammation mediated by NLRP3 inflammasome

Background

The complex pathophysiology of epilepsy hampers the development of effective treatments. Although more than ten kinds of anti-seizures drugs (ASDs) have good effects on seizure control worldwide, about 30% of patients still display pharmacoresistance against ASDs. Neuroinflammation seems to play a crucial role in disease progression. G protein-coupled receptor 120 (GPR120) has been shown to negatively regulate inflammation and apoptosis. However, the role of GPR120 in epilepsy remains unclear. In this study, we aimed to explore the mechanism of GPR120 in epilepsy.

Methods

Male adult C57BL/6 mice were intracranially injected with kainic acid (KA) to establish epilepsy model, and the adeno associated virus (AAV) was administered intracranially at 3 weeks before KA injection. VX765 was administered by intragastric administration at 30 min before KA induced and an equal dose administrated twice a day (10 a.m. and 4 p.m.) lasting 7 days until the mice were killed. Western blot analysis, immunofluorescence staining, video monitoring of seizure, LFP recording, Nissl staining were performed.

Results

GPR120 was increased in both the hippocampus and cortex in the KA-induced model with temporal lobe epilepsy (TLE), and both were most highly expressed at 7 days after KA injection. Overexpression of GPR120 significantly alleviated epileptic activity, reduced neuronal death after status epilepticus (SE), downregulated the expression of IL-1β, IL-6, IL-18, and pyrin domain-containing protein 3 (NLRP3) inflammasome, whereas knockdown GPR120 showed the opposite effect. The effects of GPR120 knockdown were reversed by VX765 inhibition cysteinyl aspartate specific proteinase-1 (Caspase-1).

Conclusion

GPR120 modulates epileptic seizure activity and affects neuronal survival in KA-induced mouse model of temporal lobe epilepsy. Furthermore, GPR120 regulated neuroinflammation in epileptic animals through NLRP3/Caspase-1/IL-1β signaling pathway.

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3'untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.

Nucleotide Composition and Codon Usage Across Viruses and Their Respective Hosts

The genetic material of the three domains of life (Bacteria, Archaea, and Eukaryota) is always double-stranded DNA, and their GC content (molar content of guanine plus cytosine) varies between ≈ 13% and ≈ 75%. Nucleotide composition is the simplest way of characterizing genomes. Despite this simplicity, it has several implications. Indeed, it is the main factor that determines, among other features, dinucleotide frequencies, repeated short DNA sequences, and codon and amino acid usage. Which forces drive this strong variation is still a matter of controversy. For rather obvious reasons, most of the studies concerning this huge variation and its consequences, have been done in free-living organisms. However, no recent comprehensive study of all known viruses has been done (that is, concerning all available sequences). Viruses, by far the most abundant biological entities on Earth, are the causative agents of many diseases. An overview of these entities is important also because their genetic material is not always double-stranded DNA: indeed, certain viruses have as genetic material single-stranded DNA, double-stranded RNA, single-stranded RNA, and/or retro-transcribing. Therefore, one may wonder if what we have learned about the evolution of GC content and its implications in prokaryotes and eukaryotes also applies to viruses. In this contribution, we attempt to describe compositional properties of ∼ 10,000 viral species: base composition (globally and according to Baltimore classification), correlations among non-coding regions and the three codon positions, and the relationship of the nucleotide frequencies and codon usage of viruses with the same feature of their hosts. This allowed us to determine how the base composition of phages strongly correlate with the value of their respective hosts, while eukaryotic viruses do not (with fungi and protists as exceptions). Finally, we discuss some of these results concerning codon usage: reinforcing previous results, we found that phages and hosts exhibit moderate to high correlations, while for eukaryotes and their viruses the correlations are weak or do not exist.

Rapid discovery of novel prophages using biological feature engineering and machine learning

Prophages are phages that are integrated into bacterial genomes and which are key to understanding many aspects of bacterial biology. Their extreme diversity means they are challenging to detect using sequence similarity, yet this remains the paradigm and thus many phages remain unidentified. We present a novel, fast and generalizing machine learning method based on feature space to facilitate novel prophage discovery. To validate the approach, we reanalyzed publicly available marine viromes and single-cell genomes using our feature-based approaches and found consistently more phages than were detected using current state-of-the-art tools while being notably faster. This demonstrates that our approach significantly enhances bacteriophage discovery and thus provides a new starting point for exploring new biologies.

Controlled expression of lysis gene E by a mutant of the promoter pL of the thermo-inducible λcI857-pL system

AIMS

To identify a lambda promoter pL mutant that could extend the thermal stability of the thermo-inducible λcI857-pR/pL system and to evaluate the effects of the modified system for the controlled expression of lysis gene E during the production of bacterial ghosts (BGs).

METHODS AND RESULTS

The promoter pL mutant was identified by random mutagenesis and site-directed mutagenesis. The results showed that a T → 35C mutation in the pL promoter was responsible for the phenotype alteration. Under the same induction conditions, the lysis rates of the modified lytic system on Escherichia coli and Salmonella enteritidis were significantly lower than that of the control, while the lysis rates of Escherichia coli with the thermo-inducible lytic system were significantly higher than that of S. enteritidis with the corresponding plasmid (P < 0·05).

CONCLUSIONS

Increasing the heat stability of the thermo-inducible lytic systems decreased lysis efficiency during the production of BGs. There exist differences in the lysis efficiency of thermo-inducible lytic systems between different bacterial strains.

SIGNIFICANCE AND IMPACT OF THE STUDY

These findings enrich current knowledge about modifications to thermo-inducible systems and provide a reference for the application of these modified systems for the production of BGs and controlled gene expression in bacteria.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Improving Phylogenetic Signals of Mitochondrial Genes Using a New Method of Codon Degeneration

Recovering deep phylogeny is challenging with animal mitochondrial genes because of their rapid evolution. Codon degeneration decreases the phylogenetic noise and bias by aiming to achieve two objectives: (1) alleviate the bias associated with nucleotide composition, which may lead to homoplasy and long-branch attraction, and (2) reduce differences in the phylogenetic results between nucleotide-based and amino acid (AA)-based analyses. The discrepancy between nucleotide-based analysis and AA-based analysis is partially caused by some synonymous codons that differ more from each other at the nucleotide level than from some nonsynonymous codons, e.g., Leu codon TTR in the standard genetic code is more similar to Phe codon TTY than to synonymous CTN codons. Thus, nucleotide similarity conflicts with AA similarity. There are many such examples involving other codon families in various mitochondrial genetic codes. Proper codon degeneration will make synonymous codons more similar to each other at the nucleotide level than they are to nonsynonymous codons. Here, I illustrate a "principled" codon degeneration method that achieves these objectives. The method was applied to resolving the mammalian basal lineage and phylogenetic position of rheas among ratites. The codon degeneration method was implemented in the user-friendly and freely available DAMBE software for all known genetic codes (genetic codes 1 to 33).

Viral adaption of staphylococcal phage: A genome-based analysis of the selective preference based on codon usage Bias

Given the high therapeutic value of the staphylococcal phage, the genome co-evolution of the phage and the host has gained great attention. Though the genome-wide AT richness in staphylococcal phages has been well-studied with nucleotide usage bias, here we proved that host factor, lifestyle and taxonomy are also important factors in understanding the phage nucleotide usages bias using information entropy formula. Such correlation is especially prominent when it comes to the synonymous codon usages of staphylococcal phages, despite the overall scattered codon usage pattern represented by principal component analysis. This strong relationship is explained by nucleotide skew which testified that the usage biases of nucleotide at different codon positions are acting on synonymous codons. Therefore, our study reveals a hidden relationship of genome evolution with host limitation and phagic phenotype, providing new insight into phage genome evolution at genetic level.

Dissimilation of synonymous codon usage bias in virus-host coevolution due to translational selection

Eighteen of the 20 amino acids are each encoded by more than one synonymous codon. Due to differential transfer RNA supply within the cell, synonymous codons are not used with equal frequency, a phenomenon termed codon usage bias (CUB). Previous studies have demonstrated that CUB of endogenous genes trans-regulates the translational efficiency of other genes. We hypothesized similar effects for CUB of exogenous genes on host translation, and tested it in the case of viral infection, a common form of naturally occurring exogenous gene translation. We analysed public Ribo-Seq datasets from virus-infected yeast and human cells and showed that virus CUB trans-regulated tRNA availability, and therefore the relative decoding time of codons. Manipulative experiments in yeast using 37 synonymous fluorescent proteins confirmed that an exogenous gene with CUB more similar to that of the host would apply decreased translational load on the host per unit of expression, whereas expression of the exogenous gene was elevated. The combination of these two effects was that exogenous genes with CUB overly similar to that of the host severely impeded host translation. Finally, using a manually curated list of viruses and natural and symptomatic hosts, we found that virus CUB tended to be more similar to that of symptomatic hosts than that of natural hosts, supporting a general deleterious effect of excessive CUB similarity between virus and host. Our work revealed repulsion between virus and host CUBs when they are overly similar, a previously unrecognized complexity in the coevolution of virus and host.

An improved estimation of tRNA expression to better elucidate the coevolution between tRNA abundance and codon usage in bacteria

The degree to which codon usage can be explained by tRNA abundance in bacterial species is often inadequate, partly because differential tRNA abundance is often approximated by tRNA copy numbers. To better understand the coevolution between tRNA abundance and codon usage, we provide a better estimate of tRNA abundance by profiling tRNA mapped reads (tRNA tpm) using publicly available RNA Sequencing data. To emphasize the feasibility of our approach, we demonstrate that tRNA tpm is consistent with tRNA abundances derived from RNA fingerprinting experiments in Escherichia coli, Bacillus subtilis, and Salmonella enterica. Furthermore, we do not observe an appreciable reduction in tRNA sequencing efficiency due to post-transcriptional methylations in the seven bacteria studied. To determine optimal codons, we calculate codon usage in highly and lowly expressed genes determined by protein per transcript. We found that tRNA tpm is sensitive to identify more translationally optimal codons than gene copy number and early tRNA fingerprinting abundances. Additionally, tRNA tpm improves the predictive power of tRNA adaptation index over codon preference. Our results suggest that dependence of codon usage on tRNA availability is not always associated with species growth-rate. Conversely, tRNA availability is better optimized to codon usage in fast-growing than slow-growing species.

Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Viruses are known to be highly dependent on the host translation machinery for their protein synthesis. However, tRNA genes are occasionally identified in such organisms, and in addition, few of them harbor tRNA gene clusters comprising dozens of genes. Recently, tRNA gene clusters have been shown to occur among the three domains of life. In such a scenario, the viruses could play a role in the dispersion of such structures among these organisms. Thus, in order to reveal the prevalence of tRNA genes as well as tRNA gene clusters in viruses, we performed an unbiased large-scale genome survey. Interestingly, tRNA genes were predicted in ssDNA (single-stranded DNA) and ssRNA (single-stranded RNA) viruses as well in many other dsDNA viruses of families from Caudovirales order. In the latter group, tRNA gene clusters composed of 15 to 37 tRNA genes were characterized, mainly in bacteriophages, enlarging the occurrence of such structures within viruses. These bacteriophages were from hosts that encompass five phyla and 34 genera. This in-silico study presents the current global scenario of tRNA genes and their organization in virus genomes, contributing and opening questions to be explored in further studies concerning the role of the translation apparatus in these organisms.

Genomic characterization of three novel Basilisk-like phages infecting Bacillus anthracis

Background

In the present study, we sequenced the complete genomes of three novel bacteriophages v_B-Bak1, v_B-Bak6, v_B-Bak10 previously isolated from historical anthrax burial sites in the South Caucasus country of Georgia. We report here major trends in the molecular evolution of these phages, which we designate as “Basilisk-Like-Phages” (BLPs), and illustrate patterns in their evolution, genomic plasticity and core genome architecture.

Results

Comparative whole genome sequence analysis revealed a close evolutionary relationship between our phages and two unclassified Bacillus cereus group phages, phage Basilisk, a broad host range phage (Grose JH et al., J Vir. 2014;88(20):11846-11860) and phage PBC4, a highly host-restricted phage and close relative of Basilisk (Na H. et al. FEMS Microbiol. letters. 2016;363(12)). Genome comparisons of phages v_B-Bak1, v_B-Bak6, and v_B-Bak10 revealed significant similarity in sequence, gene content, and synteny with both Basilisk and PBC4. Transmission electron microscopy (TEM) confirmed the three phages belong to the Siphoviridae family. In contrast to the broad host range of phage Basilisk and the single-strain specificity of PBC4, our three phages displayed host specificity for Bacillus anthracis. Bacillus species including Bacillus cereus, Bacillus subtilis, Bacillus anthracoides, and Bacillus megaterium were refractory to infection.

Conclusions

Data reported here provide further insight into the shared genomic architecture, host range specificity, and molecular evolution of these rare B. cereus group phages. To date, the three phages represent the only known close relatives of the Basilisk and PBC4 phages and their shared genetic attributes and unique host specificity for B. anthracis provides additional insight into candidate host range determinants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5056-4) contains supplementary material, which is available to authorized users.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

The Evolution of Molecular Compatibility between Bacteriophage ΦX174 and its Host

Viruses rely upon their hosts for biosynthesis of viral RNA, DNA and protein. This dependency frequently engenders strong selection for virus genome compatibility with potential hosts, appropriate gene regulation and expression necessary for a successful infection. While bioinformatic studies have shown strong correlations between codon usage in viral and host genomes, the selective factors by which this compatibility evolves remain a matter of conjecture. Engineered to include codons with a lesser usage and/or tRNA abundance within the host, three different attenuated strains of the bacterial virus ФX174 were created and propagated via serial transfers. Molecular sequence data indicate that biosynthetic compatibility was recovered rapidly. Extensive computational simulations were performed to assess the role of mutational biases as well as selection for translational efficiency in the engineered phage. Using bacteriophage as a model system, we can begin to unravel the evolutionary processes shaping codon compatibility between viruses and their host.

Bacteriophage evolution differs by host, lifestyle and genome

Bacteriophages play key roles in microbial evolution^1,2, marine nutrient cycling³ and human disease⁴. Phages are genetically diverse, and their genome architectures are characteristically mosaic, driven by horizontal gene transfer with other phages and host genomes⁵. As a consequence, phage evolution is complex and their genomes are composed of genes with distinct and varied evolutionary histories^6,7. However, there are conflicting perspectives on the roles of mosaicism and the extent to which it generates a spectrum of genome diversity⁸ or genetically discrete populations^9,10. Here, we show that bacteriophages evolve within two general evolutionary modes that differ in the extent of horizontal gene transfer by an order of magnitude. Temperate phages distribute into high and low gene flux modes, whereas lytic phages share only the lower gene flux mode. The evolutionary modes are also a function of the bacterial host and different proportions of temperate and lytic phages are distributed in either mode depending on the host phylum. Groups of genetically related phages fall into either the high or low gene flux modes, suggesting there are genetic as well as ecological drivers of horizontal gene transfer rates. Consequently, genome mosaicism varies depending on the host, lifestyle and genetic constitution of phages.

Isolation and Characterization of a Shewanella Phage-Host System from the Gut of the Tunicate, Ciona intestinalis

Outnumbering all other biological entities on earth, bacteriophages (phages) play critical roles in structuring microbial communities through bacterial infection and subsequent lysis, as well as through horizontal gene transfer. While numerous studies have examined the effects of phages on free-living bacterial cells, much less is known regarding the role of phage infection in host-associated biofilms, which help to stabilize adherent microbial communities. Here we report the cultivation and characterization of a novel strain of Shewanella fidelis from the gut of the marine tunicate Ciona intestinalis, inducible prophages from the S. fidelis genome, and a strain-specific lytic phage recovered from surrounding seawater. In vitro biofilm assays demonstrated that lytic phage infection affects biofilm formation in a process likely influenced by the accumulation and integration of the extracellular DNA released during cell lysis, similar to the mechanism that has been previously shown for prophage induction.

The Role of +4U as an Extended Translation Termination Signal in Bacteria

Termination efficiency of stop codons depends on the first 3' flanking (+4) base in bacteria and eukaryotes. In both Escherichia coli and Saccharomyces cerevisiae, termination read-through is reduced in the presence of +4U; however, the molecular mechanism underlying +4U function is poorly understood. Here, we perform comparative genomics analysis on 25 bacterial species (covering Actinobacteria, Bacteriodetes, Cyanobacteria, Deinococcus-Thermus, Firmicutes, Proteobacteria, and Spirochaetae) with bioinformatics approaches to examine the influence of +4U in bacterial translation termination by contrasting highly- and lowly-expressed genes (HEGs and LEGs, respectively). We estimated gene expression using the recently formulated Index of Translation Elongation, ITE, and identified stop codon near-cognate transfer RNAs (tRNAs) from well-annotated genomes. We show that +4U was consistently overrepresented in UAA-ending HEGs relative to LEGs. The result is consistent with the interpretation that +4U enhances termination mainly for UAA. Usage of +4U decreases in GC-rich species where most stop codons are UGA and UAG, with few UAA-ending genes, which is expected if UAA usage in HEGs drives up +4U usage. In HEGs, +4U usage increases significantly with abundance of UAA nc_tRNAs (near-cognate tRNAs that decode codons differing from UAA by a single nucleotide), particularly those with a mismatch at the first stop codon site. UAA is always the preferred stop codon in HEGs, and our results suggest that UAAU is the most efficient translation termination signal in bacteria.

Bioinformatics and Drug Discovery

Bioinformatic analysis can not only accelerate drug target identification and drug candidate screening and refinement, but also facilitate characterization of side effects and predict drug resistance. High-throughput data such as genomic, epigenetic, genome architecture, cistromic, transcriptomic, proteomic, and ribosome profiling data have all made significant contribution to mechanismbased drug discovery and drug repurposing. Accumulation of protein and RNA structures, as well as development of homology modeling and protein structure simulation, coupled with large structure databases of small molecules and metabolites, paved the way for more realistic protein-ligand docking experiments and more informative virtual screening. I present the conceptual framework that drives the collection of these high-throughput data, summarize the utility and potential of mining these data in drug discovery, outline a few inherent limitations in data and software mining these data, point out news ways to refine analysis of these diverse types of data, and highlight commonly used software and databases relevant to drug discovery.

Evolutionary interpretations of mycobacteriophage biodiversity and host-range through the analysis of codon usage bias

In an genomics course sponsored by the Howard Hughes Medical Institute (HHMI), undergraduate students have isolated and sequenced the genomes of more than 1,150 mycobacteriophages, creating the largest database of sequenced bacteriophages able to infect a single host, Mycobacterium smegmatis, a soil bacterium. Genomic analysis indicates that these mycobacteriophages can be grouped into 26 clusters based on genetic similarity. These clusters span a continuum of genetic diversity, with extensive genomic mosaicism among phages in different clusters. However, little is known regarding the primary hosts of these mycobacteriophages in their natural habitats, nor of their broader host ranges. As such, it is possible that the primary host of many newly isolated mycobacteriophages is not M. smegmatis, but instead a range of closely related bacterial species. However, determining mycobacteriophage host range presents difficulties associated with mycobacterial cultivability, pathogenicity and growth. Another way to gain insight into mycobacteriophage host range and ecology is through bioinformatic analysis of their genomic sequences. To this end, we examined the correlations between the codon usage biases of 199 different mycobacteriophages and those of several fully sequenced mycobacterial species in order to gain insight into the natural host range of these mycobacteriophages. We find that UPGMA clustering tends to match, but not consistently, clustering by shared nucleotide sequence identify. In addition, analysis of GC content, tRNA usage and correlations between mycobacteriophage and mycobacterial codon usage bias suggests that the preferred host of many clustered mycobacteriophages is not M. smegmatis but other, as yet unknown, members of the mycobacteria complex or closely allied bacterial species.

PhyPA: Phylogenetic method with pairwise sequence alignment outperforms likelihood methods in phylogenetics involving highly diverged sequences

While pairwise sequence alignment (PSA) by dynamic programming is guaranteed to generate one of the optimal alignments, multiple sequence alignment (MSA) of highly divergent sequences often results in poorly aligned sequences, plaguing all subsequent phylogenetic analysis. One way to avoid this problem is to use only PSA to reconstruct phylogenetic trees, which can only be done with distance-based methods. I compared the accuracy of this new computational approach (named PhyPA for phylogenetics by pairwise alignment) against the maximum likelihood method using MSA (the ML+MSA approach), based on nucleotide, amino acid and codon sequences simulated with different topologies and tree lengths. I present a surprising discovery that the fast PhyPA method consistently outperforms the slow ML+MSA approach for highly diverged sequences even when all optimization options were turned on for the ML+MSA approach. Only when sequences are not highly diverged (i.e., when a reliable MSA can be obtained) does the ML+MSA approach outperforms PhyPA. The true topologies are always recovered by ML with the true alignment from the simulation. However, with MSA derived from alignment programs such as MAFFT or MUSCLE, the recovered topology consistently has higher likelihood than that for the true topology. Thus, the failure to recover the true topology by the ML+MSA is not because of insufficient search of tree space, but by the distortion of phylogenetic signal by MSA methods. I have implemented in DAMBE PhyPA and two approaches making use of multi-gene data sets to derive phylogenetic support for subtrees equivalent to resampling techniques such as bootstrapping and jackknifing.

Coevolution between Stop Codon Usage and Release Factors in Bacterial Species

Three stop codons in bacteria represent different translation termination signals, and their usage is expected to depend on their differences in translation termination efficiency, mutation bias, and relative abundance of release factors (RF1 decoding UAA and UAG, and RF2 decoding UAA and UGA). In 14 bacterial species (covering Proteobacteria, Firmicutes, Cyanobacteria, Actinobacteria and Spirochetes) with cellular RF1 and RF2 quantified, UAA is consistently over-represented in highly expressed genes (HEGs) relative to lowly expressed genes (LEGs), whereas UGA usage is the opposite even in species where RF2 is far more abundant than RF1. UGA usage relative to UAG increases significantly with P_RF2 [=RF2/(RF1 + RF2)] as expected from adaptation between stop codons and their decoders. P_RF2 is > 0.5 over a wide range of AT content (measured by P_AT3 as the proportion of AT at third codon sites), but decreases rapidly toward zero at the high range of P_AT3. This explains why bacterial lineages with high P_AT3 often have UGA reassigned because of low RF2. There is no indication that UAG is a minor stop codon in bacteria as claimed in a recent publication. The claim is invalid because of the failure to apply the two key criteria in identifying a minor codon: (1) it is least preferred by HEGs (or most preferred by LEGs) and (2) it corresponds to the least abundant decoder. Our results suggest a more plausible explanation for why UAA usage increases, and UGA usage decreases, with P_AT3, but UAG usage remains low over the entire P_AT3 range.

Escherichia coli and Staphylococcus phages: effect of translation initiation efficiency on differential codon adaptation mediated by virulent and temperate lifestyles

Rapid biosynthesis is key to the success of bacteria and viruses. Highly expressed genes in bacteria exhibit a strong codon bias corresponding to the differential availability of tRNAs. However, a large clade of lambdoid coliphages exhibits relatively poor codon adaptation to the host translation machinery, in contrast to other coliphages that exhibit strong codon adaptation to the host. Three possible explanations were previously proposed but dismissed: (1) the phage-borne tRNA genes that reduce the dependence of phage translation on host tRNAs, (2) lack of time needed for evolving codon adaptation due to recent host switching, and (3) strong strand asymmetry with biased mutation disrupting codon adaptation. Here, we examined the possibility that phages with relatively poor codon adaptation have poor translation initiation which would weaken the selection on codon adaptation. We measured translation initiation by: (1) the strength and position of the Shine-Dalgarno (SD) sequence, and (2) the stability of the secondary structure of sequences flanking the SD and start codon known to affect accessibility of the SD sequence and start codon. Phage genes with strong codon adaptation had significantly stronger SD sequences than those with poor codon adaptation. The former also had significantly weaker secondary structure in sequences flanking the SD sequence and start codon than the latter. Thus, lambdoid phages do not exhibit strong codon adaptation because they have relatively inefficient translation initiation and would benefit little from increased elongation efficiency. We also provided evidence suggesting that phage lifestyle (virulent versus temperate) affected selection intensity on the efficiency of translation initiation and elongation.

A major controversy in codon-anticodon adaptation resolved by a new codon usage index

Two alternative hypotheses attribute different benefits to codon-anticodon adaptation. The first assumes that protein production is rate limited by both initiation and elongation and that codon-anticodon adaptation would result in higher elongation efficiency and more efficient and accurate protein production, especially for highly expressed genes. The second claims that protein production is rate limited only by initiation efficiency but that improved codon adaptation and, consequently, increased elongation efficiency have the benefit of increasing ribosomal availability for global translation. To test these hypotheses, a recent study engineered a synthetic library of 154 genes, all encoding the same protein but differing in degrees of codon adaptation, to quantify the effect of differential codon adaptation on protein production in Escherichia coli. The surprising conclusion that “codon bias did not correlate with gene expression” and that “translation initiation, not elongation, is rate-limiting for gene expression” contradicts the conclusion reached by many other empirical studies. In this paper, I resolve the contradiction by reanalyzing the data from the 154 sequences. I demonstrate that translation elongation accounts for about 17% of total variation in protein production and that the previous conclusion is due to the use of a codon adaptation index (CAI) that does not account for the mutation bias in characterizing codon adaptation. The effect of translation elongation becomes undetectable only when translation initiation is unrealistically slow. A new index of translation elongation I_TE is formulated to facilitate studies on the efficiency and evolution of the translation machinery.