Selection on start codons in prokaryotes and potential compensatory nucleotide substitutions

A Practical Guide for the Quality Evaluation of Fluobodies/Chromobodies

BACKGROUND

Fluorescent proteins (FPs) are pivotal reagents for flow cytometry analysis or fluorescent microscopy. A new generation of immunoreagents (fluobodies/chromobodies) has been developed by fusing recombinant nanobodies to FPs.

METHODS

We analyzed the quality of such biomolecules by a combination of gel filtration and SDS-PAGE to identify artefacts due to aggregation or material degradation.

RESULTS

In the SDS-PAGE run, unexpected bands corresponding to separate fluobodies were evidenced and characterized as either degradation products or artefacts that systematically resulted in the presence of specific FPs and some experimental conditions. The elimination of N-terminal methionine from FPs did not impair the appearance of FP fragments, whereas the stability and migration characteristics of some FP constructs were strongly affected by heating in loading buffer, which is a step samples undergo before electrophoretic separation.

CONCLUSIONS

In this work, we provide explanations for some odd results observed during the quality control of fluobodies and summarize practical suggestions for the choice of the most convenient FPs to fuse to antibody fragments.

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.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

The pAblo·pCasso self-curing vector toolset for unconstrained cytidine and adenine base-editing in Gram-negative bacteria

A synthetic biology toolkit, exploiting clustered regularly interspaced short palindromic repeats (CRISPR) and modified CRISPR-associated protein (Cas) base-editors, was developed for genome engineering in Gram-negative bacteria. Both a cytidine base-editor (CBE) and an adenine base-editor (ABE) have been optimized for precise single-nucleotide modification of plasmid and genome targets. CBE comprises a cytidine deaminase conjugated to a Cas9 nickase from Streptococcus pyogenes (SpnCas9), resulting in C→T (or G→A) substitutions. Conversely, ABE consists of an adenine deaminase fused to SpnCas9 for A→G (or T→C) editing. Several nucleotide substitutions were achieved using these plasmid-borne base-editing systems and a novel protospacer adjacent motif (PAM)-relaxed SpnCas9 (SpRY) variant. Base-editing was validated in Pseudomonas putida and other Gram-negative bacteria by inserting premature STOP codons into target genes, thereby inactivating both fluorescent proteins and metabolic (antibiotic-resistance) functions. The functional knockouts obtained by engineering STOP codons via CBE were reverted to the wild-type genotype using ABE. Additionally, a series of induction-responsive vectors was developed to facilitate the curing of the base-editing platform in a single cultivation step, simplifying complex strain engineering programs without relying on homologous recombination and yielding plasmid-free, modified bacterial cells.

StORF-Reporter: finding genes between genes

Large regions of prokaryotic genomes are currently without any annotation, in part due to well-established limitations of annotation tools. For example, it is routine for genes using alternative start codons to be misreported or completely omitted. Therefore, we present StORF-Reporter, a tool that takes an annotated genome and returns regions that may contain missing CDS genes from unannotated regions. StORF-Reporter consists of two parts. The first begins with the extraction of unannotated regions from an annotated genome. Next, Stop-ORFs (StORFs) are identified in these unannotated regions. StORFs are open reading frames that are delimited by stop codons and thus can capture those genes most often missing in genome annotations. We show this methodology recovers genes missing from canonical genome annotations. We inspect the results of the genomes of model organisms, the pangenome of Escherichia coli, and a set of 5109 prokaryotic genomes of 247 genera from the Ensembl Bacteria database. StORF-Reporter extended the core, soft-core and accessory gene collections, identified novel gene families and extended families into additional genera. The high levels of sequence conservation observed between genera suggest that many of these StORFs are likely to be functional genes that should now be considered for inclusion in canonical annotations.

Maximizing Heterologous Expression of Engineered Type I Polyketide Synthases: Investigating Codon Optimization Strategies

Type I polyketide synthases (T1PKSs) hold enormous potential as a rational production platform for the biosynthesis of specialty chemicals. However, despite great progress in this field, the heterologous expression of PKSs remains a major challenge. One of the first measures to improve heterologous gene expression can be codon optimization. Although controversial, choosing the wrong codon optimization strategy can have detrimental effects on the protein and product levels. In this study, we analyzed 11 different codon variants of an engineered T1PKS and investigated in a systematic approach their influence on heterologous expression in Corynebacterium glutamicum, Escherichia coli, and Pseudomonas putida. Our best performing codon variants exhibited a minimum 50-fold increase in PKS protein levels, which also enabled the production of an unnatural polyketide in each of these hosts. Furthermore, we developed a free online tool (https://basebuddy.lbl.gov) that offers transparent and highly customizable codon optimization with up-to-date codon usage tables. In this work, we not only highlight the significance of codon optimization but also establish the groundwork for the high-throughput assembly and characterization of PKS pathways in alternative hosts.

Chloroplasts evolved an additional layer of translational regulation based on non-AUG start codons for proteins with different turnover rates

Chloroplasts have evolved from photosynthetic cyanobacteria-like progenitors through endosymbiosis. The chloroplasts of present-day land plants have their own transcription and translation systems that show several similarities with prokaryotic organisms. A remarkable feature of the chloroplast translation system is the use of non-AUG start codons in the protein synthesis of certain genes that are evolutionarily conserved from Algae to angiosperms. However, the biological significance of such use of non-AUG codons is not fully understood. The present study was undertaken to unravel the significance of non-AUG start codons in vivo using the chloroplast genetic engineering approach. For this purpose, stable transplastomic tobacco plants expressing a reporter gene i.e. uidA (GUS) under four different start codons (AUG/UUG/GUG/CUG) were generated and β-glucuronidase (GUS) expression was compared. To investigate further the role of promoter sequences proximal to the start codon, uidA was expressed under two different chloroplast gene promoters psbA and psbC that use AUG and a non-AUG (GUG) start codons, respectively, and also showed significant differences in the DNA sequence surrounding the start codon. Further, to delineate the role of RNA editing that creates AUG start codon by editing non-AUG codons, if any, which is another important feature of the chloroplast transcription and translation system, transcripts were sequenced. In addition, a proteomic approach was used to identify the translation initiation site(s) of GUS and the N-terminal amino acid encoded when expressed under different non-AUG start codons. The results showed that chloroplasts use non-AUG start codons in combination with the translation initiation site as an additional layer of gene regulation to over-express proteins that are required at high levels due to their high rates of turnover.

Heterologous Gene Regulation in Clostridia: Rationally Designed Gene Regulation for Industrial and Medical Applications

Several species from the Clostridium genus show promise as industrial solvent producers and cancer therapeutic delivery vehicles. Previous development of shuttle plasmids and genome editing tools has aided the study of these species and enabled their exploitation in industrial and medical applications. Nevertheless, the precise control of gene expression is still hindered by the limited range of characterized promoters. To address this, libraries of promoters (native and synthetic), 5' UTRs, and alternative start codons were constructed. These constructs were tested in Escherichia coli K-12, Clostridium sporogenes NCIMB 10696, and Clostridium butyricum DSM 10702, using β-glucuronidase (gusA) as a gene reporter. Promoter activity was corroborated using a second gene reporter, nitroreductase (nmeNTR) from Neisseria meningitides. A strong correlation was observed between the two reporters. In C. sporogenes and C. butyricum, respectively, changes in GusA activity between the weakest and strongest expressing levels were 129-fold and 78-fold. Similar results were obtained with the nmeNTR. Using the GusA reporter, translation initiation from six alternative (non-AUG) start codons was measured in E. coli, C. sporogenes, and C. butyricum. Clearly, species-specific differences between clostridia and E. coli in translation initiation were observed, and the performance of the start codons was influenced by the upstream 5' UTR sequence. These results highlight a new opportunity for gene control in recombinant clostridia. To demonstrate the value of these results, expression of the sacB gene from Bacillus subtilis was optimized for use as a novel negative selection marker in C. butyricum. In summary, these results indicate improvements in the understanding of heterologous gene regulation in Clostridium species and E. coli cloning strains. This new knowledge can be utilized for rationally designed gene regulation in Clostridium-mediated industrial and medical applications, as well as fundamental research into the biology of Clostridium species.

Comparative Analysis of Diverse Acetyltransferase-Type Toxin-Antitoxin Loci in Klebsiella pneumoniae

Toxin-antitoxin (TA) modules containing a Gcn5-related N-acetyltransferase (GNAT) toxin domain regulate bacterial physiology under adverse environmental stresses. Multiple GNAT-ribbon-helix-helix domain (RHH) TA loci have been identified in single bacterial genomes. However, their diversity and interactions are still obscure. Our previous analysis showed that the GNAT toxin of Klebsiella pneumoniae, KacT, introduces antibiotic tolerance and the toxicity of GNAT is neutralized by KacA, an RHH antitoxin. We here present a phylogenetic analysis of GNAT toxins of more than 1,000 GNAT-RHH pairs detected in completely sequenced K. pneumoniae genomes, revealing that the GNAT toxins are diverse and grouped into four distinct clades. Overexpression of GNAT toxins representative of each of the four clades halts the cell growth of K. pneumoniae, while the coexpression of the cognate RHH antitoxin neutralizes GNAT toxicity. We also identify point mutations that inactivate the GNAT toxins. Moreover, we observe a cross-interaction between GNAT-RHH pairs encoded by different replicons, where a chromosomal toxin (KacT2) can be neutralized by its cognate RHH antitoxin (KacA2 on a chromosome) and another antitoxin (KacA3 on a plasmid). Finally, statistical analysis of the distribution of GNAT-RHH loci in K. pneumoniae strains shows pronounced deviation from random distribution within the same clades. Moreover, we also obtain statistically significant correlations between different clades, which we discuss in terms of the experimental results.

IMPORTANCE Elucidating the roles of multifaceted GNAT-RHH TA loci is essential for understanding how these TAs interact among themselves. Recently, the reaction mechanisms and structures of several GNAT-RHH pairs have been reported. While bacterial strains can carry multiple GNAT-RHH loci with diverse origins, studies on the possible cross-interactions of these TA pairs are still limited. Here, we find that 1,000 predicted GNAT toxins of K. pneumoniae can be grouped into four distinct clades. The distributions of TA loci within these clades in K. pneumoniae strains are highly nonrandom, with the presence of a single locus of each clade per strain being highly overrepresented. Moreover, the toxicity of a GNAT toxin encoded by a chromosome was alleviated by a noncognate RHH antitoxin on a plasmid. These results might yield a profound understanding of the widespread GNAT-RHH TA pairs and the cross-interactions between noncognate TA pairs located on different replicons.

Did Homocysteine Take Part in the Start of the Synthesis of Peptides on the Early Earth?

Unlike its shorter analog, cysteine, and its methylated derivative, methionine, homocysteine is not today a proteinogenic amino acid. However, this thiol containing amino acid is capable of forming an activated species intramolecularly. Its thiolactone could have made it an interesting molecular building block at the origin of life on Earth. Here we study the cyclization of homocysteine in water and show theoretically and experimentally that in an acidic medium the proportion of thiolactone is significant. This thiolactone easily reacts with amino acids to form dipeptides. We envision that these reactions may help interpret why a methionine residue is introduced at the start of all protein synthesis.

Collaborative Mining of Whole Genome Sequences for Intelligent HIV-1 Sub-Strain(s) Discovery

BACKGROUND

Effective global antiretroviral vaccines and therapeutic strategies depend on the diversity, evolution, and epidemiology of their various strains as well as their transmission and pathogenesis. Most viral disease-causing particles are clustered into a taxonomy of subtypes to suggest pointers toward nucleotide-specific vaccines or therapeutic applications of clinical significance sufficient for sequence-specific diagnosis and homologous viral studies. These are very useful to formulate predictors to induce cross-resistance to some retroviral control drugs being used across study areas.

OBJECTIVE

This research proposed a collaborative framework of hybridized (Machine Learning and Natural Language Processing) techniques to discover hidden genome patterns and feature predictors, for HIV-1 genome sequences mining.

METHOD

630 human HIV-1 genome sequences above 8500 bps were excavated from the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov) for 21 countries across different continents, Antarctica exempt. These sequences were transformed and learned using a self-organizing map (SOM). To discriminate emerging/new sub-strain(s), the HIV-1 reference genome was included as part of the input isolates/samples during the training. After training the SOM, component planes defining pattern clusters of the input datasets were generated, for cognitive knowledge mining and subsequent labelling of the datasets. Additional genome features including dinucleotide transmission recurrences, codon recurrences, and mutation recurrences, were finally extracted from the raw genomes to construct output classification targets for supervised learning.

RESULTS

SOM training explains the inherent pattern diversity of HIV-1 genomes as well as inter- and intra-country transmissions in which mobility might play an active role, as corroborated by literature. Nine sub-strains were discovered after disassembling the SOM correlation hunting matrix space attributed to disparate clusters. Cognitive knowledge mining separated similar pattern clusters bounded by a certain degree of correlation range, discovered by the SOM. A Kruskal-Wallis rank-sum test and Wilcoxon rank-sum test showed statistically significant variations in dinucleotide, codon, and mutation patterns.

CONCLUSION

Results of the discovered sub-strains and response clusters visualizations corroborate existing literature, with significant haplotype variations. The proposed framework would assist in the development of decision support systems for easy contact tracing, infectious disease surveillance, and studying the progressive evolution of the reference HIV-1 genome.

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.

A roadmap for metagenomic enzyme discovery

Covering: up to 2021Metagenomics has yielded massive amounts of sequencing data offering a glimpse into the biosynthetic potential of the uncultivated microbial majority. While genome-resolved information about microbial communities from nearly every environment on earth is now available, the ability to accurately predict biocatalytic functions directly from sequencing data remains challenging. Compared to primary metabolic pathways, enzymes involved in secondary metabolism often catalyze specialized reactions with diverse substrates, making these pathways rich resources for the discovery of new enzymology. To date, functional insights gained from studies on environmental DNA (eDNA) have largely relied on PCR- or activity-based screening of eDNA fragments cloned in fosmid or cosmid libraries. As an alternative, shotgun metagenomics holds underexplored potential for the discovery of new enzymes directly from eDNA by avoiding common biases introduced through PCR- or activity-guided functional metagenomics workflows. However, inferring new enzyme functions directly from eDNA is similar to searching for a 'needle in a haystack' without direct links between genotype and phenotype. The goal of this review is to provide a roadmap to navigate shotgun metagenomic sequencing data and identify new candidate biosynthetic enzymes. We cover both computational and experimental strategies to mine metagenomes and explore protein sequence space with a spotlight on natural product biosynthesis. Specifically, we compare in silico methods for enzyme discovery including phylogenetics, sequence similarity networks, genomic context, 3D structure-based approaches, and machine learning techniques. We also discuss various experimental strategies to test computational predictions including heterologous expression and screening. Finally, we provide an outlook for future directions in the field with an emphasis on meta-omics, single-cell genomics, cell-free expression systems, and sequence-independent methods.

A tripartite carbohydrate-binding module to functionalize cellulose nanocrystals

The development of protein and microorganism engineering have led to rising expectations of biotechnology in the design of emerging biomaterials, putatively of high interest to reduce our dependence on fossil carbon resources. In this way, cellulose, a renewable carbon based polysaccharide and derived products, displays unique properties used in many industrial applications. Although the functionalization of cellulose is common, it is however limited in terms of number and type of functions. In this work, a Carbohydrate-Binding Module (CBM) was used as a central core to provide a versatile strategy to bring a large diversity of functions to cellulose surfaces. CBM3a from Clostridium thermocellum, which has a high affinity for crystalline cellulose, was flanked through linkers with a streptavidin domain and an azide group introduced through a non-canonical amino acid. Each of these two extra domains was effectively produced and functionalized with a variety of biological and chemical molecules. Structural properties of the resulting tripartite chimeric protein were investigated using molecular modelling approaches, and its potential for the multi-functionalization of cellulose was confirmed experimentally. As a proof of concept, we show that cellulose can be labelled with a fluorescent version of the tripartite protein grafted to magnetic beads and captured using a magnet.

The Contribution of the Predicted Sorting Platform Component HrcQ to Type III Secretion in Xanthomonas campestris pv. vesicatoria Depends on an Internal Translation Start Site

Pathogenicity of the Gram-negative bacterium Xanthomonas campestris pv. vesicatoria depends on a type III secretion (T3S) system which translocates effector proteins into plant cells. T3S systems are conserved in plant- and animal-pathogenic bacteria and consist of at least nine structural core components, which are designated Sct (secretion and cellular translocation) in animal-pathogenic bacteria. Sct proteins are involved in the assembly of the membrane-spanning secretion apparatus which is associated with an extracellular needle structure and a cytoplasmic sorting platform. Components of the sorting platform include the ATPase SctN, its regulator SctL, and pod-like structures at the periphery of the sorting platform consisting of SctQ proteins. Members of the SctQ family form a complex with the C-terminal protein domain, SctQ_C, which is translated as separate protein and likely acts either as a structural component of the sorting platform or as a chaperone for SctQ. The sorting platform has been intensively studied in animal-pathogenic bacteria but has not yet been visualized in plant pathogens. We previously showed that the SctQ homolog HrcQ from X. campestris pv. vesicatoria assembles into complexes which associate with the T3S system and interact with components of the ATPase complex. Here, we report the presence of an internal alternative translation start site in hrcQ leading to the separate synthesis of the C-terminal protein region (HrcQ_C). The analysis of genomic hrcQ mutants showed that HrcQ_C is essential for pathogenicity and T3S. Increased expression levels of hrcQ or the T3S genes, however, compensated the lack of HrcQ_C. Interaction studies and protein analyses suggest that HrcQ_C forms a complex with HrcQ and promotes HrcQ stability. Furthermore, HrcQ_C colocalizes with HrcQ as was shown by fluorescence microscopy, suggesting that it is part of the predicted cytoplasmic sorting platform. In agreement with this finding, HrcQ_C interacts with the inner membrane ring protein HrcD and the SctK-like linker protein HrpB4 which contributes to the docking of the HrcQ complex to the membrane-spanning T3S apparatus. Taken together, our data suggest that HrcQ_C acts as a chaperone for HrcQ and as a structural component of the predicted sorting platform.

Sulfur Amino Acids: From Prebiotic Chemistry to Biology and Vice Versa

Two sulfur-containing amino acids are included in the list of the 20 classical protein amino acids. A methionine residue is introduced at the start of the synthesis of all current proteins. Cysteine, thanks to its thiol function, plays an essential role in a very large number of catalytic sites. Here we present what is known about the prebiotic synthesis of these two amino acids and homocysteine, and we discuss their introduction into primitive peptides and more elaborate proteins.

1 Introduction

2 Sulfur Sources

3 Prebiotic Synthesis of Cysteine

4 Prebiotic Synthesis of Methionine

5 Homocysteine and Its Thiolactone

6 Methionine and Cystine in Proteins

7 Prebiotic Scenarios Using Sulfur Amino Acids

8 Introduction of Cys and Met in the Genetic Code

9 Conclusion

What Do We See in Spectra?: Assignment of High-Intensity Peaks of Cutibacterium and Staphylococcus Spectra of MALDI-TOF Mass Spectrometry by Interspecies Comparative Proteogenomics

Matrix-assisted laser-desorption/ionization time-of-flight (MALDI–TOF) mass spectrometry is a widely used and reliable technology to identify microbial species and subspecies. The current methodology is based on spectral fingerprinting, analyzing protein peaks, most of which are yet to be characterized. In order to deepen the understanding of these peaks and to develop a more reasonable identification workflow, we applied proteogenomic approaches to assign the high-intensity peaks of MALDI–TOF spectra of two bacterial genera. First, the 3–22 kD proteomes of 5 Cutibacterium strains were profiled by UPLC–MS/MS, and the amino acid sequences were refined by referring to their genome in the public database. Then, the sequences were converted to m/z (x-axis) values based on their molecular masses. When the interspecies comparison of calculated m/z values was well-fitted to the observed peaks, the peak assignments for the five Cutibacterium species were confirmed. Second, the peak assignments for six Staphylococcus species were performed by using the above result for Cutibacterium and referring to ribosomal subunit proteins coded on the S10-spc-alpha operon (the S10-GERMS method), a previous proteomics report by Becher et al., and comprehensive genome analysis. We successfully assigned 13 out of 15 peaks for the Cutibacterium species and 11 out of 13 peaks for the Staphylococcus species. DNA-binding protein HU, the CsbD-like protein, and 50S ribosomal protein L7/L12 were observed in common. The commonality suggests they consist of high-intensity peaks in the MALDI spectra of other bacterial species. Our workflow may lead to the development of a more accurate species identification database of MALDI–TOF mass spectrometry based on genome data.

Analysis of Stop Codons within Prokaryotic Protein-Coding Genes Suggests Frequent Readthrough Events

Nonsense mutations turn a coding (sense) codon into an in-frame stop codon that is assumed to result in a truncated protein product. Thus, nonsense substitutions are the hallmark of pseudogenes and are used to identify them. Here we show that in-frame stop codons within bacterial protein-coding genes are widespread. Their evolutionary conservation suggests that many of them are not pseudogenes, since they maintain dN/dS values (ratios of substitution rates at non-synonymous and synonymous sites) significantly lower than 1 (this is a signature of purifying selection in protein-coding regions). We also found that double substitutions in codons—where an intermediate step is a nonsense substitution—show a higher rate of evolution compared to null models, indicating that a stop codon was introduced and then changed back to sense via positive selection. This further supports the notion that nonsense substitutions in bacteria are relatively common and do not necessarily cause pseudogenization. In-frame stop codons may be an important mechanism of regulation: Such codons are likely to cause a substantial decrease of protein expression levels.

Genetic Determinants Enabling Medium-Dependent Adaptation to Nafcillin in Methicillin-Resistant Staphylococcus aureus

Antimicrobial susceptibility testing standards driving clinical decision-making have centered around the use of cation-adjusted Mueller-Hinton broth (CA-MHB) as the medium with the notion of supporting bacterial growth, without consideration of recapitulating the in vivo environment. However, it is increasingly recognized that various medium conditions have tremendous influence on antimicrobial activity, which in turn may have major implications on the ability of in vitro susceptibility assays to predict antibiotic activity in vivo. To elucidate differential growth optimization and antibiotic resistance mechanisms, adaptive laboratory evolution was performed in the presence or absence of the antibiotic nafcillin with methicillin-resistant Staphylococcus aureus (MRSA) TCH1516 in either (i) CA-MHB, a traditional bacteriological nutritionally rich medium, or (ii) Roswell Park Memorial Institute (RPMI), a medium more reflective of the in vivo host environment. Medium adaptation analysis showed an increase in growth rate in RPMI, but not CA-MHB, with mutations in apt, adenine phosphoribosyltransferase, and the manganese transporter subunit, mntA, occurring reproducibly in parallel replicate evolutions. The medium-adapted strains showed no virulence attenuation. Continuous exposure of medium-adapted strains to increasing concentrations of nafcillin led to medium-specific evolutionary strategies. Key reproducibly occurring mutations were specific for nafcillin adaptation in each medium type and did not confer resistance in the other medium environment. Only the vraRST operon, a regulator of membrane- and cell wall-related genes, showed mutations in both CA-MHB- and RPMI-evolved strains. Collectively, these results demonstrate the medium-specific genetic adaptive responses of MRSA and establish adaptive laboratory evolution as a platform to study clinically relevant resistance mechanisms.IMPORTANCE The ability of pathogens such as Staphylococcus aureus to evolve resistance to antibiotics used in the treatment of infections has been an important concern in the last decades. Resistant acquisition usually translates into treatment failure and puts patients at risk of unfavorable outcomes. Furthermore, the laboratory testing of antibiotic resistance does not account for the different environment the bacteria experiences within the human body, leading to results that do not translate into the clinic. In this study, we forced methicillin-resistant S. aureus to develop nafcillin resistance in two different environments, a laboratory environment and a physiologically more relevant environment. This allowed us to identify genetic changes that led to nafcillin resistance under both conditions. We concluded that not only does the environment dictate the evolutionary strategy of S. aureus to nafcillin but also that the evolutionary strategy is specific to that given environment.

Non-AUG start codons: Expanding and regulating the small and alternative ORFeome

Recent ribosome profiling and proteomic studies have revealed the presence of thousands of novel coding sequences, referred to as small open reading frames (sORFs), in prokaryotic and eukaryotic genomes. These genes have defied discovery via traditional genomic tools not only because they tend to be shorter than standard gene annotation length cutoffs, but also because they are, as a class, enriched in sequence properties previously assumed to be unusual, including non-AUG start codons. In this review, we summarize what is currently known about the incidence, efficiency, and mechanism of non-AUG start codon usage in prokaryotes and eukaryotes, and provide examples of regulatory and functional sORFs that initiate at non-AUG codons. While only a handful of non-AUG-initiated novel genes have been characterized in detail to date, their participation in important biological processes suggests that an improved understanding of this class of genes is needed.

Crossing fitness valleys via double substitutions within codons

BACKGROUND

Single nucleotide substitutions in protein-coding genes can be divided into synonymous (S), with little fitness effect, and non-synonymous (N) ones that alter amino acids and thus generally have a greater effect. Most of the N substitutions are affected by purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases potentially could alleviate the deleterious effect of single substitutions, making them subject to positive selection. To elucidate the effects of selection on double substitutions in all codons, it is critical to differentiate selection from mutational biases.

RESULTS

We addressed the evolutionary regimes of within-codon double substitutions in 37 groups of closely related prokaryotic genomes from diverse phyla by comparing the fractions of double substitutions within codons to those of the equivalent double S substitutions in adjacent codons. Under the assumption that substitutions occur one at a time, all within-codon double substitutions can be represented as "ancestral-intermediate-final" sequences (where "intermediate" refers to the first single substitution and "final" refers to the second substitution) and can be partitioned into four classes: (1) SS, S intermediate-S final; (2) SN, S intermediate-N final; (3) NS, N intermediate-S final; and (4) NN, N intermediate-N final. We found that the selective pressure on the second substitution markedly differs among these classes of double substitutions. Analogous to single S (synonymous) substitutions, SS double substitutions evolve neutrally, whereas analogous to single N (non-synonymous) substitutions, SN double substitutions are subject to purifying selection. In contrast, NS show positive selection on the second step because the original amino acid is recovered. The NN double substitutions are heterogeneous and can be subject to either purifying or positive selection, or evolve neutrally, depending on the amino acid similarity between the final or intermediate and the ancestral states.

CONCLUSIONS

The results of the present, comprehensive analysis of the evolutionary landscape of within-codon double substitutions reaffirm the largely conservative regime of protein evolution. However, the second step of a double substitution can be subject to positive selection when the first step is deleterious. Such positive selection can result in frequent crossing of valleys on the fitness landscape.

Structure and functional implications of WYL domain-containing bacterial DNA damage response regulator PafBC

In mycobacteria, transcriptional activator PafBC is responsible for upregulating the majority of genes induced by DNA damage. Understanding the mechanism of PafBC activation is impeded by a lack of structural information on this transcription factor that contains a widespread, but poorly understood WYL domain frequently encountered in bacterial transcription factors. Here, we determine the crystal structure of Arthrobacter aurescens PafBC. The protein consists of two modules, each harboring an N-terminal helix-turn-helix DNA-binding domain followed by a central WYL and a C-terminal extension (WCX) domain. The WYL domains exhibit Sm-folds, while the WCX domains adopt ferredoxin-like folds, both characteristic for RNA-binding proteins. Our results suggest a mechanism of regulation in which WYL domain-containing transcription factors may be activated by binding RNA or other nucleic acid molecules. Using an in vivo mutational screen in Mycobacterium smegmatis, we identify potential co-activator binding sites on PafBC.

In eubacteria, unlike eukaryotes, there is no evidence for selection favouring fail-safe 3' additional stop codons

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense ‘codons’ downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3’ sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates.

In all organisms, gene expression is error-prone. One such error, translational read-through, occurs where the primary stop codon of an expressed gene is missed by the translational machinery. Failure to terminate is likely to be costly, hence genomes are under selection to prevent this from happening. One proposed error-proofing strategy involves in-frame proximal additional stop codons (ASCs) which may act as a ‘fail-safe’ mechanism by providing another opportunity for translation to terminate. There is evidence for ASC enrichment in several eukaryotes. We extend this evidence showing it to be common but not universal in single celled eukaryotes. However, the situation in bacteria is poorly understood, despite bacteria having high read-through rates. Here, we test the fail-safe hypothesis within a broad range of bacteria. To our surprise, we find that not only are ASCs not enriched, but they may even be selected against. This provides evidence for an unusual circumstance where eukaryotes and prokaryotes could solve the same problem the same way but don’t. What are we to make of this? We suggest that if read-through is the problem, ASCs are not necessarily the expected solution. Owing to the absence of nonsense-mediated decay, a process that makes gene truncation in eukaryotes less viable, we propose bacteria may rescue a leaky stop by mutation that creates a new stop upstream. Alternatively, raised read-through rates in some particular conditions in eukaryotes might explain the difference.

The plant pathogen Pectobacterium atrosepticum contains a functional formate hydrogenlyase-2 complex

Pectobacterium atrosepticum SCRI1043 is a phytopathogenic Gram-negative enterobacterium. Genomic analysis has identified that genes required for both respiration and fermentation are expressed under anaerobic conditions. One set of anaerobically expressed genes is predicted to encode an important but poorly understood membrane-bound enzyme termed formate hydrogenlyase-2 (FHL-2), which has fascinating evolutionary links to the mitochondrial NADH dehydrogenase (Complex I). In this work, molecular genetic and biochemical approaches were taken to establish that FHL-2 is fully functional in P. atrosepticum and is the major source of molecular hydrogen gas generated by this bacterium. The FHL-2 complex was shown to comprise a rare example of an active [NiFe]-hydrogenase-4 (Hyd-4) isoenzyme, itself linked to an unusual selenium-free formate dehydrogenase in the final complex. In addition, further genetic dissection of the genes encoding the predicted membrane arm of FHL-2 established surprisingly that the majority of genes encoding this domain are not required for physiological hydrogen production activity. Overall, this study presents P. atrosepticum as a new model bacterial system for understanding anaerobic formate and hydrogen metabolism in general, and FHL-2 function and structure in particular.

gFACs: Gene Filtering, Analysis, and Conversion to Unify Genome Annotations Across Alignment and Gene Prediction Frameworks

Published genomes frequently contain erroneous gene models that represent issues associated with identification of open reading frames, start sites, splice sites, and related structural features. The source of these inconsistencies is often traced back to integration across text file formats designed to describe long read alignments and predicted gene structures. In addition, the majority of gene prediction frameworks do not provide robust downstream filtering to remove problematic gene annotations, nor do they represent these annotations in a format consistent with current file standards. These frameworks also lack consideration for functional attributes, such as the presence or absence of protein domains that can be used for gene model validation. To provide oversight to the increasing number of published genome annotations, we present a software package, the Gene Filtering, Analysis, and Conversion (gFACs), to filter, analyze, and convert predicted gene models and alignments. The software operates across a wide range of alignment, analysis, and gene prediction files with a flexible framework for defining gene models with reliable structural and functional attributes. gFACs supports common downstream applications, including genome browsers, and generates extensive details on the filtering process, including distributions that can be visualized to further assess the proposed gene space. gFACs is freely available and implemented in Perl with support from BioPerl libraries at https://gitlab.com/PlantGenomicsLab/gFACs.

Persister Cells - a Plausible Outcome of Neutral Coevolutionary Drift

The phenomenon of bacterial persistence - a non-inherited antibiotic tolerance in a minute fraction of the bacterial population, was observed more than 70 years ago. Nowadays, it is suggested that "persister cells" undergo an alternative scenario of the cell cycle; however, pathways involved in its emergence are still not identified. We present a mathematically grounded scenario of such possibility. We have determined that population drift in the space of multiple neutrally coupled mutations, which we called "neutrally coupled co-evolution" (NCCE), leads to increased dynamic complexity of bacterial populations via appearance of cells capable of carrying out a single cell cycle in two or more alternative ways and that universal properties of the coupled transcription-translation system underlie this phenotypic multiplicity. According to our hypothesis, modern persister cells have derived from such cells and regulatory mechanisms that govern the consolidation of this phenomenon represented the trigger. We assume that the described type of neutrally coupled co-evolution could play an important role in the origin of extremophiles, both in bacteria and archaea.

Purifying and positive selection in the evolution of stop codons

Modes of evolution of stop codons in protein-coding genes, especially the conservation of UAA, have been debated for many years. We reconstructed the evolution of stop codons in 40 groups of closely related prokaryotic and eukaryotic genomes. The results indicate that the UAA codons are maintained by purifying selection in all domains of life. In contrast, positive selection appears to drive switches from UAG to other stop codons in prokaryotes but not in eukaryotes. Changes in stop codons are significantly associated with increased substitution frequency immediately downstream of the stop. These positions are otherwise more strongly conserved in evolution compared to sites farther downstream, suggesting that such substitutions are compensatory. Although GC content has a major impact on stop codon frequencies, its contribution to the decreased frequency of UAA differs between bacteria and archaea, presumably, due to differences in their translation termination mechanisms.

Genetic Analysis of the Major Capsid Protein of the Archaeal Fusellovirus SSV1: Mutational Flexibility and Conformational Change

Viruses with spindle or lemon-shaped virions are rare in the world of viruses, but are common in viruses of archaeal extremophiles, possibly due to the extreme conditions in which they thrive. However, the structural and genetic basis for the unique spindle shape is unknown. The best-studied spindle-shaped virus, Sulfolobus Spindle-shaped Virus 1 (SSV1), is composed mostly of the major capsid protein VP1. Similar to many other viruses, proteolytic cleavage of VP1 is thought to be critical for virion formation. Unlike half of the genes in SSV1, including the minor capsid protein gene VP3, the VP1 gene does not tolerate deletion or transposon insertion. To determine the role of the VP1 gene and its proteolysis for virus function, we developed techniques for site-directed mutagenesis of the SSV1 genome and complemented deletion mutants with VP1 genes from other SSVs. By analyzing these mutants, we demonstrate that the N-terminus of the VP1 protein is required, but the N-terminus, or entire SSV1 VP1 protein, can be exchanged with VP1s from other SSVs. However, the conserved glutamate at the cleavage site is not essential for infectivity. Interestingly, viruses containing point mutations at this position generate mostly abnormal virions.