Codon usage of highly expressed genes affects proteome-wide translation efficiency

Prime mass amino acids: A new numbers based classification of significance to mass spectrometry and protein biology

The nominal mass of amino acid residues calculated from their elemental compositions are defined by prime numbers far more often than chance, and such residues appear to play an important role in the formation and biology of proteins. It is proposed therefore that consideration be given to classifying prime mass amino acids as such, beyond the more common, familiar definitions associated with the other physicochemical properties of amino acids including charged or non-charged, hydrophobic or hydrophilic, polar or non-polar, acidic or basic, aliphatic or aromatic. Greater focus could also be given to such residues during peptide and protein sequencing with mass spectrometry and the construction of structural maps, given their predominantly hydrophobic character and thus their role in protein folding and transmembrane domains. The use of prime numbers to define amino acids based on the sum of the atomic masses from their elemental compositions invokes other recent interest and observations whereby prime numbers were organised in a way that mirrors electrons arranged within the orbitals of an atom. The article links number theory with both the physical and biological sciences, and mass spectrometry, for the first time.

Lost in translation: conserved amino acid usage despite extreme codon bias in foraminifera

Analyses of codon usage in eukaryotes suggest that amino acid usage responds to GC pressure so AT-biased substitutions drive higher usage of amino acids with AT-ending codons. Here, we combine single-cell transcriptomics and phylogenomics to explore codon usage patterns in foraminifera, a diverse and ancient clade of predominantly uncultivable microeukaryotes. We curate data from 1,044 gene families in 49 individuals representing 28 genera, generating perhaps the largest existing dataset of data from a predominantly uncultivable clade of protists, to analyze compositional bias and codon usage. We find extreme variation in composition, with a median GC content at fourfold degenerate silent sites below 3% in some species and above 75% in others. The most AT-biased species are distributed among diverse non-monophyletic lineages. Surprisingly, despite the extreme variation in compositional bias, amino acid usage is highly conserved across all foraminifera. By analyzing nucleotide, codon, and amino acid composition within this diverse clade of amoeboid eukaryotes, we expand our knowledge of patterns of genome evolution across the eukaryotic tree of life.IMPORTANCEPatterns of molecular evolution in protein-coding genes reflect trade-offs between substitution biases and selection on both codon and amino acid usage. Most analyses of these factors in microbial eukaryotes focus on model species such as Acanthamoeba, Plasmodium, and yeast, where substitution bias is a primary contributor to patterns of amino acid usage. Foraminifera, an ancient clade of single-celled eukaryotes, present a conundrum, as we find highly conserved amino acid usage underlain by divergent nucleotide composition, including extreme AT-bias at silent sites among multiple non-sister lineages. We speculate that these paradoxical patterns are enabled by the dynamic genome structure of foraminifera, whose life cycles can include genome endoreplication and chromatin extrusion.

Codon usage patterns and genomic variation analysis of chloroplast genomes provides new insights into the evolution of Aroideae

Aroideae is an important subfamily of the Araceae family and contains many plants with medicinal and edible value. It is difficult to identify and classify Aroideae species accurately on the basis of morphology alone because of their polymorphic phenotypic traits. The chloroplast genome (CPG) is useful for studying on plant taxonomy and phylogeny, and the analysis of codon usage bias (CUB) in CPGs provides further insights into the intricate phylogenetic relationships among Aroideae. The results showed that the codon third position of the chloroplast genome coding sequence in Aroideae was rich in A and T, with a GC content of 37.91%. The ENC-plot and PR2-plot revealed that the codon usage bias of Aroideae was influenced by multiple factors, with natural selection as the dominant factor. Thirteen to twenty optimal codons ending in A/T were identified in 61 Aroideae species. Additionally, the comparative analysis of CPGs revealed that two single copy regions and non-coding regions were variable in Aroideae. Eight highly divergent regions (Pi > 0.064) were identified (ndhF, rpl32, ccsA, ndhE, ndhG, ndhF-rpl32, ccsA-ndhD, and ndhE-ndhG) , in which ndhE have the potential to serve as a reliable DNA marker to discriminate chloroplasts in Aroideae subfamily. Furthermore, the maximum likelihood-based phylogenetic trees constructed from complete chloroplast genomes and protein-coding sequences presented similar topologies. Principal component clustering analysis based on relative synonymous codon usage values (RSCUs) revealed that Calla was clearly deviated from Montrichardia and Anubias, and that Alocasia was closer to Colocasieae than to Arisaemateae. These findings suggest that the use of RSCU for clustering analysis could offer new theoretical support for species classification and evolution. Our research could provide a theoretical foundation for the chloroplast genetic engineering, taxonomy, and phylogenetic relationships of Aroideae chloroplasts.

Decoding the Chloroplast Genome of Korean endemic plant Acer okamotoanum: Comparative Genomics, Phylogenetic Insights, and Potential for Marker Development

Acer okamotoanum, a medicinally significant endemic plant of Korea, has seen limited genomic research. To address this gap, we conducted a comprehensive sequencing and analysis of its chloroplast genome. The assembled genome is 156,242 bp in length, with typical quadripartite structure, consisting of a large single-copy region, a small single-copy region, and two inverted repeat regions. It contains 130 genes, including 85 protein-coding, 37 tRNA, and 8 rRNA genes. Sixteen genes have a single intron, while clpP and ycf3 possess two introns each. Additionally, 17 genes are duplicated within the inverted repeat regions. The genome analysis revealed 92 Simple Sequence Repeats (SSRs), predominantly located in intergenic regions, with a bias toward A/T-rich codons. Comparative analysis with five closely related Acer species highlighted a highly conserved genomic structure, but also revealed differences in SSRs and repeat sequences. Hypervariable regions, such as rpl32-trnL and ycf1, were identified as potential molecular markers for phylogenetic and population studies. Phylogenetic analysis involving 37 chloroplast genomes confirmed the monophyly of the Acer genus and placed A. okamotoanum within the Platanoidea section, closely related to A. truncatum. This study improves the understanding of A. okamotoanum's genomic structure, offering insights for phylogenetic analysis, marker development, and conservation efforts.

Fitness and adaptive evolution of a Rhodococcus sp. harboring dioxin-catabolic plasmids

Catabolic plasmids are critical factors in the degradation of recalcitrant xenobiotics, such as dioxins. Understanding the persistence and evolution of native catabolic plasmids is pivotal for controlling their function in microbial remediation. Here, we track the fitness and evolution of Rhodococcus sp. strain p52 harboring dioxin-catabolic plasmids under nonselective conditions without contaminant. Growth curve analysis and competition experiments demonstrated that pDF01 imposed fitness costs, whereas pDF02 conferred fitness benefits. During stability tests, pDF01 tended to be lost from the population, while pDF02 maintained at least one copy in the cell until proliferation of the 400th generation. Genome-wide gene expression profiling combined with codon usage bias analysis revealed that the high expression of pDF01 genes involved in dibenzofuran catabolism and regulation caused metabolic burdens. In contrast, potential cooperation between the pDF02-encoded short-chain dehydrogenase/reductase family oxidoreductase and the redox cofactor mycofactocin, which synthetic genes are located on the chromosome, may explain the benefit of pDF02. The fitness cost imposed by pDF01 was alleviated during adaptive evolution and was associated with the transcriptional downregulation of the dibenzofuran degradation genes on pDF01, and the global regulation of genome-wide gene expression involving basic metabolism, transport, and signal transduction. This study broadens our understandings on the persistence and evolution of dioxin-catabolic mega-plasmids, thus paving the way for the bioremediation of recalcitrant xenobiotic pollution in the environment.

Comparative genomics of thermosynechococcaceae and thermostichaceae: insights into codon usage bias

Members of the families Thermosynechococcaceae and Thermostichaceae are well-known unicellular thermophilic cyanobacteria and a non-thermophilic genus Pseudocalidococcus was newly classified into the former. Analysis of the codon usage bias (CUB) of cyanobacterial species inhabiting different thermal and non-thermal niches will benefit the understanding of their genetic and evolutionary characteristics. Herein, the CUB and codon context patterns of protein-coding genes were systematically analyzed and compared between members of the two families. Overall, the nucleotide composition and CUB indices were found to differ between thermophiles and non-thermophiles. The thermophiles showed a higher G/C content in the codon base composition and tended to end with G/C compared to the non-thermophiles. Correlation analysis indicated significant associations between codon base composition and CUB indices. The results of the effective number of codons, parity-rule 2, neutral and correspondence analyses indicated that mutational pressure and natural selection primarily account for CUB in these cyanobacterial species, but the primary driving forces exhibit variation among genera. Moreover, the optimal codons identified based on relative synonymous codon usage values were found to differ among genera and even within genera. In addition, codon context pattern analysis revealed the specificity of the sequence context of start and stop codons among genera. Intriguingly, the clustering of codon context patterns appeared to be more related to thermotolerance than to phylogenomic relationships. In conclusion, this study facilitates the understanding of the characteristics and sources of variation of CUB and the evolution of the surveyed cyanobacterial clades with different thermotolerance and provides insights into their adaptation to different environments.

Predicting gene sequences with AI to study codon usage patterns

Selective pressure acts on the codon use, optimizing multiple, overlapping signals that are only partially understood. We trained AI models to predict codons given their amino acid sequence in the eukaryotes Saccharomyces cerevisiae and Schizosaccharomyces pombe and the bacteria Escherichia coli and Bacillus subtilis to study the extent to which we can learn patterns in naturally occurring codons to improve predictions. We trained our models on a subset of the proteins and evaluated their predictions on large, separate sets of proteins of varying lengths and expression levels. Our models significantly outperformed naïve frequency-based approaches, demonstrating that there are learnable dependencies in evolutionary-selected codon usage. The prediction accuracy advantage of our models is greater for highly expressed genes and is greater in bacteria than eukaryotes, supporting the hypothesis that there is a monotonic relationship between selective pressure for complex codon patterns and effective population size. In S. cerevisiae and bacteria, our models were more accurate for longer proteins, suggesting that the learned patterns may be related to cotranslational folding. Gene functionality and conservation were also important determinants that affect the performance of our models. Finally, we showed that using information encoded in homologous proteins has only a minor effect on prediction accuracy, perhaps due to complex codon-usage codes in genes undergoing rapid evolution. Our study employing contemporary AI methods offers a unique perspective and a deep-learning-based prediction tool for evolutionary-selected codons. We hope that these can be useful to optimize codon usage in endogenous and heterologous proteins.

Aminoglycoside tolerance in Vibrio cholerae engages translational reprogramming associated with queuosine tRNA modification

Tgt is the enzyme modifying the guanine (G) in tRNAs with GUN anticodon to queuosine (Q). tgt is required for optimal growth of Vibrio cholerae in the presence of sub-lethal aminoglycoside concentrations. We further explored here the role of the Q34 in the efficiency of codon decoding upon tobramycin exposure. We characterized its impact on the overall bacterial proteome, and elucidated the molecular mechanisms underlying the effects of Q34 modification in antibiotic translational stress response. Using molecular reporters, we showed that Q34 impacts the efficiency of decoding at tyrosine TAT and TAC codons. Proteomics analyses revealed that the anti-SoxR factor RsxA is better translated in the absence of tgt. RsxA displays a codon bias toward tyrosine TAT and overabundance of RsxA leads to decreased expression of genes belonging to SoxR oxidative stress regulon. We also identified conditions that regulate tgt expression. We propose that regulation of Q34 modification in response to environmental cues leads to translational reprogramming of transcripts bearing a biased tyrosine codon usage. In silico analysis further identified candidate genes which could be subject to such translational regulation, among which DNA repair factors. Such transcripts, fitting the definition of modification tunable transcripts, are central in the bacterial response to antibiotics.

Insight into the codon usage patterns and adaptation of Tembusu Virus

Since its emergence in 2010, Tembusu virus (TMUV) has rapidly spread across poultry populations in Asia, leading to substantial economic losses in these areas. Here, we investigate the codon usage patterns (CUPs) underlying TMUV's adaptation and evolutionary dynamics within host environments. Phylogenetic and compositional analyses consistently classify TMUV into four evolutionary lineages-Clusters 1, 2, 3, and ancestral TMUV-with Cluster 2 emerging as the dominant lineage. Codon adaptation index (CAI) analysis reveals that this lineage of TMUV show best adapted to the CUPs of duck than other lineages, underscoring the role of natural selection in shaping viral evolution, a finding in line with evidence that CUPs in the TMUV genome is predominantly shaped by natural selection. Furthermore, TMUV exhibits markedly higher adaptation to the CUPs of poultry hosts (duck, goose, and chicken) compared to potential host humans or vector mosquito. Thus, species-specific adaptability to the host environment may be a reason account for the distinct infectivity and clinic outcome of TMUV acted on hosts. Analysis of dinucleotide distribution reveals significant suppression of CpG and UpA dinucleotides in the TMUV genome, reflecting adaptive pressures to evade vertebrate immune responses. During transmission, TMUV shows increasing alignment with host CUPs and a continuous reduction in CpG dinucleotides, potentially enhancing its fitness within host microenvironments. This work advances our understanding of the basic biology underlying TMUV epidemiology, pathogenicity, and species-specific adaptation.

Compare Analysis of Codon Usage Bias of Nuclear Genome in Eight Sapindaceae Species

Codon usage bias (CUB) refers to the different frequencies with which various codons are utilized within a genome. Examining CUB is essential for understanding genome structure, function, and evolution. However, little was known about codon usage patterns and the factors influencing the nuclear genomes of eight ecologically significant Sapindaceae species widely utilized for food and medicine. In this study, an analysis of nucleotide composition revealed a higher A/T content and showed a preference for A/T at the third codon position in the eight species of Sapindaceae. A correspondence analysis of relative synonymous codon usage explained only part of the variation, suggesting that not only natural selection but also various other factors contribute to selective constraints on codon bias in the nuclear genomes of the eight Sapindaceae species. Additionally, ENC-GC3 plot, PR2-Bias, and neutrality plot analyses indicated that natural selection exerted a greater influence than mutation pressure across these eight species. Among the eight Sapindaceae species, 16 to 26 optimal codons were identified, with two common high-frequency codons: AGA (encoding Arg) and GCU (encoding Ala). The clustering heat map, which included the 8 Sapindaceae species and 13 other species, revealed two distinct clusters corresponding to monocots and dicots. This finding suggested that CUB analysis was particularly effective in elucidating evolutionary relationships at the family level. Collectively, our results emphasized the distinct codon usage characteristics and unique evolutionary traits of the eight Sapindaceae species.

Investigating the role of RNA-binding protein Ssd1 in aneuploidy tolerance through network analysis

RNA-binding proteins (RBPs) play critical cellular roles by mediating various stages of RNA life cycles. Ssd1, an RBP with pleiotropic effects, has been implicated in aneuploidy tolerance in Saccharomyces cerevisiae but its mechanistic role remains unclear. Here, we used a network-based approach to inform on Ssd1's role in aneuploidy tolerance, by identifying and experimentally perturbing a network of RBPs that share mRNA targets with Ssd1. We identified RBPs whose bound mRNA targets significantly overlap with Ssd1 targets. For 14 identified RBPs, we then used a genetic approach to generate all combinations of genotypes for euploid and aneuploid yeast with an extra copy of chromosome XII, with and without SSD1 and/or the RBP of interest. Deletion of 10 RBPs either exacerbated or alleviated the sensitivity of wild-type and/or ssd1Δ cells to chromosome XII duplication, in several cases indicating genetic interactions with SSD1 in the context of aneuploidy. We integrated these findings with results from a global overexpression screen that identified genes whose duplication complements ssd1Δ aneuploid sensitivity. The resulting network points to a subgroup of proteins with shared roles in translational repression and P-body formation, implicating these functions in aneuploidy tolerance. Our results reveal a role for new RBPs in aneuploidy tolerance and support a model in which Ssd1 mitigates translation-related stresses in aneuploid cells.

Modeling and measuring how codon usage modulates the relationship between burden and yield during protein overexpression in bacteria

Excess utilization of translational resources is a critical source of burden on cells engineered to over-express exogenous proteins. To improve protein yields and genetic stability, researchers often use codon optimization strategies that improve translational efficiency by matching an exogenous gene's codon usage with that of the host organism's highly expressed genes. Despite empirical data that shows the benefits of codon optimization, little is known quantitatively about the relationship between codon usage bias and the burden imposed by protein overexpression. Here, we develop and experimentally evaluate a stochastic gene expression model that considers the impact of codon usage bias on the availability of ribosomes and different tRNAs in a cell. In agreement with other studies, our model shows that increasing exogenous protein expression decreases production of native cellular proteins in a linear fashion. We also find that the slope of this relationship is modulated by how well the codon usage bias of the exogenous gene and the host's genes match. Strikingly, we predict that an overoptimization domain exists where further increasing usage of optimal codons worsens yield and burden. We test our model by expressing sfGFP and mCherry2 from constructs that have a wide range of codon optimization levels in Escherichia coli . The results agree with our model, including for an mCherry2 gene sequence that appears to lose expression and genetic stability from codon overoptimization. Our findings can be leveraged by researchers to predict and design more optimal cellular systems through the use of more nuanced codon optimization strategies.

Further Evidence for Strong Nonneutrality of Yeast Synonymous Mutations

Although synonymous mutations are commonly assumed neutral or nearly so, recent years have seen reports of fitness effects of synonymous mutations detected under laboratory conditions. In a previous study, we used genome editing to construct thousands of yeast mutants each carrying a synonymous or nonsynonymous mutation in one of 21 genes, and discovered that most synonymous and most nonsynonymous mutations are deleterious. A concern was raised that this observation could be caused by the fitness effects of potential CRISPR/Cas9 off-target edits and/or secondary mutations, and an experiment that would be refractory to such effects was proposed. Using genome sequencing, we here show that no CRISPR/Cas9 off-target editing occurred, although some mutants did carry secondary mutations. Analysis of mutants with negligible effects from secondary mutations and new data collected from the proposed experiment confirms the original conclusion. These findings, along with other reports of fitness effects of synonymous mutations from both case and systematic studies, necessitate a paradigm shift from assuming (near) neutrality of synonymous mutations.

SARS-CoV-2 Displays a Suboptimal Codon Usage Bias for Efficient Translation in Human Cells Diverted by Hijacking the tRNA Epitranscriptome

Codon bias analysis of SARS-CoV-2 reveals suboptimal adaptation for translation in human cells it infects. The detailed examination of the codons preferentially used by SARS-CoV-2 shows a strong preference for LysAAA, GlnCAA, GluGAA, and ArgAGA, which are infrequently used in human genes. In the absence of an adapted tRNA pool, efficient decoding of these codons requires a 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2) modification at the U34 wobble position of the corresponding tRNAs (tLysUUU; tGlnUUG; tGluUUC; tArgUCU). The optimal translation of SARS-CoV-2 open reading frames (ORFs) may therefore require several adjustments to the host's translation machinery, enabling the highly biased viral genome to achieve a more favorable "Ready-to-Translate" state in human cells. Experimental approaches based on LC-MS/MS quantification of tRNA modifications and on alteration of enzymatic tRNA modification pathways provide strong evidence to support the hypothesis that SARS-CoV-2 induces U34 tRNA modifications and relies on these modifications for its lifecycle. The conclusions emphasize the need for future studies on the evolution of SARS-CoV-2 codon bias and its ability to alter the host tRNA pool through the manipulation of RNA modifications.

Link Between Individual Codon Frequencies and Protein Expression: Going Beyond Codon Adaptation Index

An important role of a particular synonymous codon composition of a gene in its expression level is well known. There are a number of algorithms optimizing codon usage of recombinant genes to maximize their expression in host cells. Nevertheless, the underlying mechanism remains unsolved and is of significant relevance. In the realm of modern biotechnology, directing protein production to a specific level is crucial for metabolic engineering, genome rewriting and a growing number of other applications. In this study, we propose two new simple statistical and empirical methods for predicting the protein expression level from the nucleotide sequence of the corresponding gene: Codon Expression Index Score (CEIS) and Codon Productivity Score (CPS). Both of these methods are based on the influence of each individual codon in the gene on the overall expression level of the encoded protein and the frequencies of isoacceptors in the species. Our predictions achieve a correlation level of up to r = 0.7 with experimentally measured quantitative proteome data of Escherichia coli, which is superior to any previously proposed methods. Our work helps understand how codons determine protein abundances. Based on these methods, it is possible to design proteins optimized for expression in a particular organism.

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

Genome-wide impact of codon usage bias on translation optimization in Drosophila melanogaster

Accuracy and efficiency are fundamental to mRNA translation. Codon usage bias is widespread across species. Despite the long-standing association between optimized codon usage and improved translation, our understanding of its evolutionary basis and functional effects remains limited. Drosophila is widely used to study codon usage bias, but genome-scale experimental data are scarce. Using high-resolution mass spectrometry data from Drosophila melanogaster, we show that optimal codons have lower translation errors than nonoptimal codons after accounting for these biases. Genomic-scale analysis of ribosome profiling data shows that optimal codons are translated more rapidly than nonoptimal codons. Although we find no long-term selection favoring synonymous mutations in D. melanogaster after diverging from D. simulans, we identify signatures of positive selection driving codon optimization in the D. melanogaster population. These findings expand our understanding of the functional consequences of codon optimization and serve as a foundation for future investigations.

High-level expression of codon-optimized Taq DNA polymerase under the control of rhaBAD promoter

A DNA polymerase from Thermus aquaticus remains the most popular among DNA polymerases. It was widely applied in various fields involving the application of polymerase chain reaction (PCR), implying the high commercial value of this enzyme. For this reason, an attempt to obtain a high yield of Taq DNA polymerase is continuously conducted. In this study, the l-rhamnose-inducible promoter rhaBAD was utilized due to its ability to produce recombinant protein under tight control in E. coli expression system. Instead of full-length Taq polymerase, an N-terminal deletion of Taq polymerase was selected. To obtain a high-level expression, we attempted to optimize the codon by reducing the rare codon and GC content, and in a second attempt, we optimized the culture conditions for protein expression. The production of Taq polymerase using the optimum culture condition improved the level of expression by up to 3-fold. This approach further proved that a high level of recombinant protein expression could be achieved by yielding a purified Taq polymerase of about 8.5 mg/L of culture. This is the first research publication on the production of Taq polymerase with N-terminal deletion in E. coli with the control of the rhaBAD promoter system.

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria species.

Investigating the role of RNA-binding protein Ssd1 in aneuploidy tolerance through network analysis

RNA-binding proteins (RBPs) play critical cellular roles by mediating various stages of RNA life cycles. Ssd1, an RBP with pleiotropic effects, has been implicated in aneuploidy tolerance in Saccharomyces cerevisiae but its mechanistic role remains unclear. Here we used a network-based approach to inform on Ssd1's role in aneuploidy tolerance, by identifying and experimentally perturbing a network of RBPs that share mRNA targets with Ssd1. We identified RBPs whose bound mRNA targets significantly overlap with Ssd1 targets. For 14 identified RBPs, we then used a genetic approach to generate all combinations of genotypes for euploid and aneuploid yeast with an extra copy of chromosome XII, with and without SSD1 and/or the RBP of interest. Deletion of 10 RBPs either exacerbated or alleviated the sensitivity of wild-type and/or ssd1 Δ cells to chromosome XII duplication, in several cases indicating genetic interactions with SSD1 in the context of aneuploidy. We integrated these findings with results from a global over-expression screen that identified genes whose duplication complements ssd1 Δ aneuploid sensitivity. The resulting network points to a sub-group of proteins with shared roles in translational repression and p-body formation, implicating these functions in aneuploidy tolerance. Our results reveal a role for new RBPs in aneuploidy tolerance and support a model in which Ssd1 mitigates translation-related stresses in aneuploid cells.

Codon usage bias in yeasts and its correlation with gene expression, growth temperature, and protein structure

Codon usage bias (CUB) has been described in viruses, prokaryotes, and eukaryotes and has been linked to several cellular and environmental factors, such as the organism's growth temperature, gene expression levels, and regulation of protein synthesis and folding. Most of the studies in this area have been conducted in bacteria and higher eukaryotes, in some cases with different results. In this study, a comparative analysis of CUB in yeasts isolated from cold and template environments was performed in order to evaluate the correlation of CUB with yeast optimal temperature of growth (OTG), gene expression levels, cellular function, and structure of encoded proteins. Among the main findings, highly expressed ORFs tend to have a more similar CUB within and between yeasts, and a direct correlation between codons ending in C and expression level was generally found. A low correspondence between CUB and OTG was observed, with an inverse correlation for some codons ending in C. The clustering of yeasts based on their CUB partially aligns with their OTG, being more consistent for yeasts with lower OTG. In most yeasts, the abundance of preferred codons was generally lower at the 5' end of ORFs, higher in segments encoding beta strand, lower in segments encoding extracellular and transmembrane regions, and higher in "translation" and "energy metabolism" pathways, especially in highly expressed ORFs. Based on our findings, it is suggested that the abundance and distribution of preferred and non-preferred codons along mRNAs contribute to proper protein folding and functionality by regulating protein synthesis rates, becoming a more important factor under conditions that require faster protein synthesis in yeasts.

Modifications of the 5' region of the CASPONTM tag's mRNA further enhance soluble recombinant protein production in Escherichia coli

BACKGROUND

Escherichia coli is one of the most commonly used host organisms for the production of biopharmaceuticals, as it allows for cost-efficient and fast recombinant protein expression. However, challenging proteins are often produced with low titres or as inclusion bodies, and the manufacturing process needs to be developed individually for each protein. Recently, we developed the CASPONTM technology, a generic fusion tag-based platform process for high-titer soluble expression including a standardized downstream processing and highly specific enzymatic cleavage of the fusion tag. To assess potential strategies for further improvement of the N-terminally fused CASPONTM tag, we modified the 5'UTR and 5' region of the tag-coding mRNA to optimize the ribosome-mRNA interactions.

RESULTS

In the present work, we found that by modifying the 5'UTR sequence of a pET30acer plasmid-based system, expression of the fusion protein CASPONTM-tumour necrosis factor α was altered in laboratory-scale carbon-limited fed-batch cultivations, but no significant increase in expression titre was achieved. Translation efficiency was highest for a construct carrying an expression enhancer element and additionally possessing a very favourable interaction energy between ribosome and mRNA (∆Gtotal). However, a construct with comparatively low transcriptional efficiency, which lacked the expression enhancer sequence and carried the most favourable ∆Gtotal tested, led to the highest recombinant protein formation alongside the reference pET30a construct. Furthermore, we found, that by introducing synonymous mutations within the nucleotide sequence of the T7AC element of the CASPONTM tag, utilizing a combination of rare and non-rare codons, the free folding energy of the nucleotides at the 5' end (-4 to + 37) of the transcript encoding the CASPONTM tag increased by 6 kcal/mol. Surprisingly, this new T7ACrare variant led to improved recombinant protein titres by 1.3-fold up to 5.3-fold, shown with three industry-relevant proteins in lab-scale carbon limited fed-batch fermentations under industrially relevant conditions.

CONCLUSIONS

This study reveals some of the complex interdependencies between the ribosome and mRNA that govern recombinant protein expression. By modifying the 5'UTR to obtain an optimized interaction energy between the mRNA and the ribosome (ΔGtotal), transcript levels were changed, highlighting the different translation efficiencies of individual transcripts. It was shown that the highest recombinant titre was not obtained by the construct with the most efficient translation but by a construct with a generally high transcript amount coupled with a favourable ΔGtotal. Furthermore, an unexpectedly high potential to enhance expression by introducing silent mutations including multiple rare codons into the 5'end of the CAPONTM tag's mRNA was identified. Although the titres of the fusion proteins were dramatically increased, no formation of inclusion bodies or negative impact on cell growth was observed. We hypothesize that the drastic increase in titre is most likely caused by better ribosomal binding site accessibility. Our study, which demonstrates the influence of changes in ribosome-mRNA interactions on protein expression under industrially relevant production conditions, opens the door to the applicability of the new T7ACrare tag in biopharmaceutical industry using the CASPONTM platform process.

The RdRp genotyping of SARS-CoV-2 isolated from patients with different clinical spectrum of COVID-19

BACKGROUND

The evolution of SARS-CoV-2 has been observed from the very beginning of the fight against COVID-19, some mutations are indicators of potentially dangerous variants of the virus. However, there is no clear association between the genetic variants of SARS-CoV-2 and the severity of COVID-19. We aimed to analyze the genetic variability of RdRp in correlation with different courses of COVID-19.

RESULTS

The prospective study included 77 samples of SARS-CoV-2 isolated from outpatients (1st degree of severity) and hospitalized patients (2nd, 3rd and 4th degree of severity). The retrospective analyses included 15,898,266 cases of SARS-CoV-2 genome sequences deposited in the GISAID repository. Single-nucleotide variants were identified based on the four sequenced amplified fragments of SARS-CoV-2. The analysis of the results was performed using appropriate statistical methods, with p < 0.05, considered statistically significant. Additionally, logistic regression analysis was performed to predict the strongest determinants of the observed relationships. The number of mutations was positively correlated with the severity of the COVID-19, and older male patients. We detected four mutations that significantly increased the risk of hospitalization of COVID-19 patients (14676C > T, 14697C > T, 15096 T > C, and 15279C > T), while the 15240C > T mutation was common among strains isolated from outpatients. The selected mutations were searched worldwide in the GISAID database, their presence was correlated with the severity of COVID-19.

CONCLUSION

Identified mutations have the potential to be used to assess the increased risk of hospitalization in COVID-19 positive patients. Experimental studies and extensive epidemiological data are needed to investigate the association between individual mutations and the severity of COVID-19.

Codon usage bias and phylogenetic analysis of chloroplast genome in 36 gracilariaceae species

Gracilariaceae is a group of marine large red algae and main source of agar with important economic and ecological value. The codon usage patterns of chloroplast genomes in 36 species from Graciliaceae show that GC range from 0.284 to 0.335, the average GC3 range from 0.135 to 0.243 and the value of ENC range from 35.098 to 42.327, which indicates these genomes are rich in AT and prefer to use codons ending with AT in these species. Nc plot, PR2 plot, neutrality plot analyses and correlation analysis indicate that these biases may be caused by multiple factors, such as natural selection and mutation pressure, but prolonged natural selection is the main driving force influencing codon usage preference. The cluster analysis and phylogenetic analysis show that the differentiation relationship of them is different and indicate that codons with weak or unbiased preferences may also play an irreplaceable role in these species' evolution. In addition, we identified 26 common high-frequency codons and 8-18 optimal codons all ending in A/U in these 36 species. Our results will not only contribute to carrying out transgenic work in Gracilariaceae species to maximize the protein yield in the future, but also lay a theoretical foundation for further exploring systematic classification of them.

Specific codons control cellular resources and fitness

As cellular engineering progresses from simply overexpressing proteins to imparting complex phenotypes through multigene expression, judicious appropriation of cellular resources is essential. Since codon use is degenerate and biased, codons may control cellular resources at a translational level. We investigate how partitioning transfer RNA (tRNA) resources by incorporating dissimilar codon usage can drastically alter interdependence of expression level and burden on the host. By isolating the effect of individual codons' use during translation elongation while eliminating confounding factors, we show that codon choice can trans-regulate fitness of the host and expression of other heterologous or native genes. We correlate specific codon usage patterns with host fitness and derive a coding scheme for multigene expression called the Codon Health Index (CHI, χ). This empirically derived coding scheme (χ) enables the design of multigene expression systems that avoid catastrophic cellular burden and is robust across several proteins and conditions.

Characterization of the chloroplast genome of Gleditsia species and comparative analysis

The genus Gleditsia has significant medicinal and economic value, but information about the chloroplast genomic characteristics of Gleditsia species has been limited. Using the Illumina sequencing, we assembled and annotated the whole chloroplast genomes of seven Gleditsia species (Gleditsia sinensis, Gleditsia japonica var. delavayi (G. delavayi), G. fera, G. japonica, G. microphylla, Fructus Gleditsiae Abnormalis (Zhū Yá Zào), G. microphylla mutant). The assembled genomes revealed that Gleditsia species have a typical circular tetrad structure, with genome sizes ranging from 162,746 to 170,907 bp. Comparative genomic analysis showed that most (65.8-75.8%) of the abundant simple sequence repeats in Gleditsia and Gymnocladus species were located in the large single copy region. The Gleditsia chloroplast genome prefer T/A-ending codons and avoid C/G-ending codons, positive selection was acting on the rpoA, rpl20, atpB, ndhA and ycf4 genes, most of the chloroplast genes of Gleditsia species underwent purifying selection. Expansion and contraction of the inverted repeat (IR)/single copy (SC) region showed similar patterns within the Gleditsia genus. Polymorphism analysis revealed that coding regions were more conserved than non-coding regions, and the IR region was more conserved than the SC region. Mutational hotspots were mostly found in intergenic regions such as "rps16-trnQ", "trnT-trnL", "ndhG-ndhI", and "rpl32-trnL" in Gleditsia. Phylogenetic analysis showed that G. fera is most closely related to G. sinensis,G. japonica and G. delavayi are relatively closely related. Zhū Yá Zào can be considered a bud mutation of the G. sinensis. The albino phenotype of G. microphylla mutant is not caused by variations in the chloroplast genome, and that the occurrence of the albino phenotype may be due to mutations in chloroplast-related genes involved in splicing or localization functions. This study will help us enhance our exploration of the genetic evolution and geographical origins of the Gleditsia genus.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Codon usage characterization and phylogenetic analysis of the mitochondrial genome in Hemerocallis citrina

BACKGROUND

Hemerocallis citrina Baroni is a traditional vegetable crop widely cultivated in eastern Asia for its high edible, medicinal, and ornamental value. The phenomenon of codon usage bias (CUB) is prevalent in various genomes and provides excellent clues for gaining insight into organism evolution and phylogeny. Comprehensive analysis of the CUB of mitochondrial (mt) genes can provide rich genetic information for improving the expression efficiency of exogenous genes and optimizing molecular-assisted breeding programmes in H. citrina.

RESULTS

Here, the CUB patterns in the mt genome of H. citrina were systematically analyzed, and the possible factors shaping CUB were further evaluated. Composition analysis of codons revealed that the overall GC (GCall) and GC at the third codon position (GC3) contents of mt genes were lower than 50%, presenting a preference for A/T-rich nucleotides and A/T-ending codons in H. citrina. The high values of the effective number of codons (ENC) are indicative of fairly weak CUB. Significant correlations of ENC with the GC3 and codon counts were observed, suggesting that not only compositional constraints but also gene length contributed greatly to CUB. Combined ENC-plot, neutrality plot, and Parity rule 2 (PR2)-plot analyses augmented the inference that the CUB patterns of the H. citrina mitogenome can be attributed to multiple factors. Natural selection, mutation pressure, and other factors might play a major role in shaping the CUB of mt genes, although natural selection is the decisive factor. Moreover, we identified a total of 29 high-frequency codons and 22 optimal codons, which exhibited a consistent preference for ending in A/T. Subsequent relative synonymous codon usage (RSCU)-based cluster and mt protein coding gene (PCG)-based phylogenetic analyses suggested that H. citrina is close to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in evolutionary terms, reflecting a certain correlation between CUB and evolutionary relationships.

CONCLUSIONS

There is weak CUB in the H. citrina mitogenome that is subject to the combined effects of multiple factors, especially natural selection. H. citrina was found to be closely related to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in terms of their evolutionary relationships as well as the CUB patterns of their mitogenomes. Our findings provide a fundamental reference for further studies on genetic modification and phylogenetic evolution in H. citrina.

An Improved Method for Eliminating or Creating Intragenic Bacterial Promoters

Recent advances in genomic refactoring have been hindered by the ever-present complication of internal or cryptic transcriptional regulation. Typical approaches to these features have been to randomize or perform mass alterations to the gene sequences thought to contain the regulatory motifs; however, this approach can cause problems by altering translational speeds, introducing long distance DNA-DNA interaction effects, and inducing RNA toxicity. Previously, we developed a rational design approach named COdon Restrained Promoter SilEncing (CORPSE) which takes externally identified promoter sequences and uses position-specific scoring matrices as proxy promoter strengths to make minimal changes to promoter sequences to disable their activity. Additionally, through inverting our system we were also able to modify weak internal promoters to increase their activity. In this chapter, we augment our previous process with the biophysical model Promoter Calculator v1.0 developed by LaFleur et al. to combine promoter identification and activity prediction, with our algorithm to silently modify promoter sequences, to provide more robust promoter elimination and creation.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis

As an economically important tree, Gleditsia sinensis Lam. is widely planted. A lack of background genetic information on G. sinensis hinders molecular breeding. Based on PacBio single-molecule real-time (SMRT) sequencing and analysis of G. sinensis, a total of 95,183 non-redundant transcript sequences were obtained, of which 93,668 contained complete open reading frames (ORFs), 2,858 were long non-coding RNAs (LncRNAs) and 18,855 alternative splicing (AS) events were identified. Genes orthologous to different Gleditsia species pairs were identified, stress-related genes had been positively selected during the evolution. AGA, AGG, and CCA were identified as the universal optimal codon in the genus of Gleditsia. EIF5A was selected as a suitable fluorescent quantitative reference gene. 315 Cytochrome P450 monooxygenases (CYP450s) and 147 uridine diphosphate (UDP)-glycosyltransferases (UGTs) were recognized through the PacBio SMRT transcriptome. Randomized selection of GsIAA14 for cloning verified the reliability of the PacBio SMRT transcriptome assembly sequence. In conclusion, the research data lay the foundation for further analysis of the evolutionary mechanism and molecular breeding of Gleditsia.

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Bioinformatic Analysis Reveals the Role of Translation Elongation Efficiency Optimisation in the Evolution of Ralstonia Genus

Translation efficiency modulates gene expression in prokaryotes. The comparative analysis of translation elongation efficiency characteristics of Ralstonia genus bacteria genomes revealed that these characteristics diverge in accordance with the phylogeny of Ralstonia. The first branch of this genus is a group of bacteria commonly found in moist environments such as soil and water that includes the species R. mannitolilytica, R. insidiosa, and R. pickettii, which are also described as nosocomial infection pathogens. In contrast, the second branch is plant pathogenic bacteria consisting of R. solanacearum, R. pseudosolanacearum, and R. syzygii. We found that the soil Ralstonia have a significantly lower number and energy of potential secondary structures in mRNA and an increased role of codon usage bias in the optimization of highly expressed genes' translation elongation efficiency, not only compared to phytopathogenic Ralstonia but also to Cupriavidus necator, which is closely related to the Ralstonia genus. The observed alterations in translation elongation efficiency of orthologous genes are also reflected in the difference of potentially highly expressed gene' sets' content among Ralstonia branches with different lifestyles. Analysis of translation elongation efficiency characteristics can be considered a promising approach for studying complex mechanisms that determine the evolution and adaptation of bacteria in various environments.

Codon usage pattern of the ancestor of green plants revealed through Rhodophyta

Rhodophyta are among the closest known relatives of green plants. Studying the codons of their genomes can help us understand the codon usage pattern and characteristics of the ancestor of green plants. By studying the codon usage pattern of all available red algae, it was found that although there are some differences among species, high-bias genes in most red algae prefer codons ending with GC. Correlation analysis, Nc-GC3s plots, parity rule 2 plots, neutrality plot analysis, differential protein region analysis and comparison of the nucleotide content of introns and flanking sequences showed that the bias phenomenon is likely to be influenced by local mutation pressure and natural selection, the latter of which is the dominant factor in terms of translation accuracy and efficiency. It is worth noting that selection on translation accuracy could even be detected in the low-bias genes of individual species. In addition, we identified 15 common optimal codons in seven red algae except for G. sulphuraria for the first time, most of which were found to be complementary and bound to the tRNA genes with the highest copy number. Interestingly, tRNA modification was found for the highly degenerate amino acids of all multicellular red algae and individual unicellular red algae, which indicates that highly biased genes tend to use modified tRNA in translation. Our research not only lays a foundation for exploring the characteristics of codon usage of the red algae as green plant ancestors, but will also facilitate the design and performance of transgenic work in some economic red algae in the future.

Growth-dependent Gene Expression Variation Influences the Strength of Codon Usage Biases

The most highly expressed genes in microbial genomes tend to use a limited set of synonymous codons, often referred to as "preferred codons." The existence of preferred codons is commonly attributed to selection pressures on various aspects of protein translation including accuracy and/or speed. However, gene expression is condition-dependent and even within single-celled organisms transcript and protein abundances can vary depending on a variety of environmental and other factors. Here, we show that growth rate-dependent expression variation is an important constraint that significantly influences the evolution of gene sequences. Using large-scale transcriptomic and proteomic data sets in Escherichia coli and Saccharomyces cerevisiae, we confirm that codon usage biases are strongly associated with gene expression but highlight that this relationship is most pronounced when gene expression measurements are taken during rapid growth conditions. Specifically, genes whose relative expression increases during periods of rapid growth have stronger codon usage biases than comparably expressed genes whose expression decreases during rapid growth conditions. These findings highlight that gene expression measured in any particular condition tells only part of the story regarding the forces shaping the evolution of microbial gene sequences. More generally, our results imply that microbial physiology during rapid growth is critical for explaining long-term translational constraints.

Growth rate is a dominant factor predicting the rhizosphere effect

The root microbiome is shaped by plant root activity, which selects specific microbial taxa from the surrounding soil. This influence on the microorganisms and soil chemistry in the immediate vicinity of the roots has been referred to as the rhizosphere effect. Understanding the traits that make bacteria successful in the rhizosphere is critical for developing sustainable agriculture solutions. In this study, we compared the growth rate potential, a complex trait that can be predicted from bacterial genome sequences, to functional traits encoded by proteins. We analyzed 84 paired rhizosphere- and soil-derived 16S rRNA gene amplicon datasets from 18 different plants and soil types, performed differential abundance analysis, and estimated growth rates for each bacterial genus. We found that bacteria with higher growth rate potential consistently dominated the rhizosphere, and this trend was confirmed in different bacterial phyla using genome sequences of 3270 bacterial isolates and 6707 metagenome-assembled genomes (MAGs) from 1121 plant- and soil-associated metagenomes. We then identified which functional traits were enriched in MAGs according to their niche or growth rate status. We found that predicted growth rate potential was the main feature for differentiating rhizosphere and soil bacteria in machine learning models, and we then analyzed the features that were important for achieving faster growth rates, which makes bacteria more competitive in the rhizosphere. As growth rate potential can be predicted from genomic data, this work has implications for understanding bacterial community assembly in the rhizosphere, where many uncultivated bacteria reside.

Developmental Selection and the Perception of Mutation Bias

The notion that mutations are random relative to their fitness effects is central to the Neo-Darwinian view of evolution. However, a recent interpretation of the patterns of mutation accumulation in the genome of Arabidopsis thaliana has challenged this notion, arguing for the presence of a targeted DNA repair mechanism that causes a nonrandom association of mutation rates and fitness effects. Specifically, this mechanism was suggested to cause a reduction in the rates of mutations on essential genes, thus lowering the rates of deleterious mutations. Central to this argument were attempts to rule out selection at the population level. Here, we offer an alternative and parsimonious interpretation of the patterns of mutation accumulation previously attributed to mutation bias, showing how they can instead or additionally be caused by developmental selection, that is selection occurring at the cellular level during the development of a multicellular organism. Thus, the depletion of deleterious mutations in A. thaliana may indeed be the result of a selective process, rather than a bias in mutation. More broadly, our work highlights the importance of considering development in the interpretation of population-genetic analyses of multicellular organisms, and it emphasizes that efforts to identify mechanisms involved in mutational biases should explicitly account for developmental selection.

Codon Usage Bias Analysis in Macronuclear Genomes of Ciliated Protozoa

Ciliated protozoa (ciliates) are unicellular eukaryotes, several of which are important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome of ciliates can promote a better understanding of the genetic mode and evolutionary history of these organisms and help optimize codons to improve gene editing efficiency in model ciliates. In this study, the following indices were calculated: the guanine-cytosine (GC) content, the frequency of the nucleotides at the third position of codons (T3, C3, A3, G3), the effective number of codons (ENc), GC content at the 3rd position of synonymous codons (GC3s), and the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, Neutrality plot analysis, ENc plot analysis, and correlation analysis were employed to explore the main influencing factors of CUB. The results showed that the GC content in the MAC genomes of each of 21 ciliate species, the genomes of which were relatively complete, was lower than 50%, and the base compositions of GC and GC3s were markedly distinct. Synonymous codon analysis revealed that the codons in most of the 21 ciliates ended with A or T and four codons were the general putative optimal codons. Collectively, our results indicated that most of the ciliates investigated preferred using the codons with anof AT-ending and that codon usage bias was affected by gene mutation and natural selection.

Growth-dependent gene expression variation influences the strength of codon usage biases

The most highly expressed genes in microbial genomes tend to use a limited set of synonymous codons, often referred to as "preferred codons." The existence of preferred codons is commonly attributed to selection pressures on various aspects of protein translation including accuracy and/or speed. However, gene expression is condition-dependent and even within single-celled organisms transcript and protein abundances can vary depending on a variety of environmental and other factors. Here, we show that growth rate-dependent expression variation is an important constraint that significantly influences the evolution of gene sequences. Using large-scale transcriptomic and proteomic data sets in Escherichia coli and Saccharomyces cerevisiae, we confirm that codon usage biases are strongly associated with gene expression but highlight that this relationship is most pronounced when gene expression measurements are taken during rapid growth conditions. Specifically, genes whose relative expression increases during periods of rapid growth have stronger codon usage biases than comparably expressed genes whose expression decreases during rapid growth conditions. These findings highlight that gene expression measured in any particular condition tells only part of the story regarding the forces shaping the evolution of microbial gene sequences. More generally, our results imply that microbial physiology during rapid growth is critical for explaining long-term translational constraints.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Application of codon usage and context analysis in genes up- or down-regulated in neurodegeneration and cancer to combat comorbidities

Introduction

Neurodegeneration and cancer present in comorbidities with inverse effects due to the expression of genes and pathways acting in opposition. Identifying and studying the genes simultaneously up or downregulated during morbidities helps curb both ailments together.

Methods

This study examines four genes. Three of these (Amyloid Beta Precursor Protein (APP), Cyclin D1 (CCND1), and Cyclin E2 (CCNE2) are upregulated, and one protein phosphatase 2 phosphatase activator (PTPA) is simultaneously downregulated in both disorders. We investigated molecular patterns, codon usage, codon usage bias, nucleotide bias in the third codon position, preferred codons, preferred codon pairs, rare codons, and codon context.

Results

Parity analysis revealed that T is preferred over A, and G is preferred over C in the third codon position, suggesting composition plays no role in nucleotide bias in both the upregulated and downregulated gene sets and that mutational forces are stronger in upregulated gene sets than in downregulated ones. Transcript length influenced the overall %A composition and codon bias, and the codon AGG exerted the strongest influence on codon usage in both the upregulated and downregulated gene sets. Codons ending in G/C were preferred for 16 amino acids, and glutamic acid-, aspartic acid-, leucine-, valine-, and phenylalanine-initiated codon pairs were preferred in all genes. Codons CTA (Leu), GTA (Val), CAA (Gln), and CGT (Arg) were underrepresented in all examined genes.

Discussion

Using advanced gene editing tools such as CRISPR/Cas or any other gene augmentation technique, these recoded genes may be introduced into the human body to optimize gene expression levels to augment neurodegeneration and cancer therapeutic regimens simultaneously.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

Plasmid Costs Explain Plasmid Maintenance, Irrespective of the Nature of Compensatory Mutations

Conjugative plasmids often carry virulence and antibiotic-resistant genes. Therefore, understanding the behavior of these extra-chromosomal DNA elements gives insights into their spread. Bacteria frequently replicate slower after plasmids' entry, an observation inconsistent with the plasmids' ubiquity in nature. Several hypotheses explain the maintenance of plasmids among bacterial communities. However, the numerous combinations of bacterial species and strains, plasmids, and environments claim a robust elucidatory mechanism of plasmid maintenance. Previous works have shown that donor cells already adapted to the plasmid may use the plasmid as a 'weapon' to compete with non-adapted plasmid-free cells. Computer simulations corroborated this hypothesis with a wide range of parameters. Here we show that donor cells benefit from harboring conjugative plasmids even if compensatory mutations in transconjugant cells occur in the plasmid, not on chromosomes. The advantage's leading causes are as follows: mutations take time to appear, many plasmids remain costly, and re-transfer of mutated plasmids usually occurs in sites distant to the original donors, implying little competition between these cells. Research in previous decades cautioned against uncritical acceptance of the hypothesis that resistance cost helps to preserve antibiotics' effectiveness. This work gives a new twist to this conclusion by showing that costs help antibiotic-resistant bacteria to compete with plasmid-free cells even if compensatory mutations appear in plasmids.

G. de Lomana A, Gomes-Filho JV, Vêncio RZN, Moritz RL, Koide T, Baliga NS. A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

Variation in synonymous evolutionary rates in the SARS-CoV-2 genome

Introduction

Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Influential variants and mutants of this virus continue to emerge, and more effective virus-related information is urgently required for identifying and predicting new mutants. According to earlier reports, synonymous substitutions were considered phenotypically silent; thus, such mutations were frequently ignored in studies of viral mutations because they did not directly cause amino acid changes. However, recent studies have shown that synonymous substitutions are not completely silent, and their patterns and potential functional correlations should thus be delineated for better control of the pandemic.

Methods

In this study, we estimated the synonymous evolutionary rate (SER) across the SARS-CoV-2 genome and used it to infer the relationship between the viral RNA and host protein. We also assessed the patterns of characteristic mutations found in different viral lineages.

Results

We found that the SER varies across the genome and that the variation is primarily influenced by codon-related factors. Moreover, the conserved motifs identified based on the SER were found to be related to host RNA transport and regulation. Importantly, the majority of the existing fixed-characteristic mutations for five important virus lineages (Alpha, Beta, Gamma, Delta, and Omicron) were significantly enriched in partially constrained regions.

Discussion

Taken together, our results provide unique information on the evolutionary and functional dynamics of SARS-CoV-2 based on synonymous mutations and offer potentially useful information for better control of the SARS-CoV-2 pandemic.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Using protein-per-mRNA differences among human tissues in codon optimization

BACKGROUND

Codon usage and nucleotide composition of coding sequences have profound effects on protein expression. However, while it is recognized that different tissues have distinct tRNA profiles and codon usages in their transcriptomes, the effect of tissue-specific codon optimality on protein synthesis remains elusive.

RESULTS

We leverage existing state-of-the-art transcriptomics and proteomics datasets from the GTEx project and the Human Protein Atlas to compute the protein-to-mRNA ratios of 36 human tissues. Using this as a proxy of translational efficiency, we build a machine learning model that identifies codons enriched or depleted in specific tissues. We detect two clusters of tissues with an opposite pattern of codon preferences. We then use these identified patterns for the development of CUSTOM, a codon optimizer algorithm which suggests a synonymous codon design in order to optimize protein production in a tissue-specific manner. In human cell-line models, we provide evidence that codon optimization should take into account particularities of the translational machinery of the tissues in which the target proteins are expressed and that our approach can design genes with tissue-optimized expression profiles.

CONCLUSIONS

We provide proof-of-concept evidence that codon preferences exist in tissue-specific protein synthesis and demonstrate its application to synthetic gene design. We show that CUSTOM can be of benefit in biological and biotechnological applications, such as in the design of tissue-targeted therapies and vaccines.

Comparative Genome Analysis of Enterococcus cecorum Reveals Intercontinental Spread of a Lineage of Clinical Poultry Isolates

Enterococcus cecorum is an emerging pathogen responsible for osteomyelitis, spondylitis, and femoral head necrosis causing animal suffering and mortality and requiring antimicrobial use in poultry. Paradoxically, E. cecorum is a common inhabitant of the intestinal microbiota of adult chickens. Despite evidence suggesting the existence of clones with pathogenic potential, the genetic and phenotypic relatedness of disease-associated isolates remains little investigated. Here, we sequenced and analyzed the genomes and characterized the phenotypes of more than 100 isolates, the majority of which were collected over the last 10 years from 16 French broiler farms. Comparative genomics, genome-wide association studies, and the measured susceptibility to serum, biofilm-forming capacity, and adhesion to chicken type II collagen were used to identify features associated with clinical isolates. We found that none of the tested phenotypes could discriminate the origin of the isolates or the phylogenetic group. Instead, we found that most clinical isolates are grouped phylogenetically, and our analyses selected six genes that discriminate 94% of isolates associated with disease from those that are not. Analysis of the resistome and the mobilome revealed that multidrug-resistant clones of E. cecorum cluster into a few clades and that integrative conjugative elements and genomic islands are the main carriers of antimicrobial resistance. This comprehensive genomic analysis shows that disease-associated clones of E. cecorum belong mainly to one phylogenetic clade. IMPORTANCE Enterococcus cecorum is an important pathogen of poultry worldwide. It causes a number of locomotor disorders and septicemia, particularly in fast-growing broilers. Animal suffering, antimicrobial use, and associated economic losses require a better understanding of disease-associated E. cecorum isolates. To address this need, we performed whole-genome sequencing and analysis of a large collection of isolates responsible for outbreaks in France. By providing the first data set on the genetic diversity and resistome of E. cecorum strains circulating in France, we pinpoint an epidemic lineage that is probably also circulating elsewhere that should be targeted preferentially by preventive strategies in order to reduce the burden of E. cecorum-related diseases.

Nucleotide and codon usage biases involved in the evolution of African swine fever virus: A comparative genomics analysis

Since African swine fever virus (ASFV) replication is closely related to its host's machinery, codon usage of viral genome can be subject to selection pressures. A better understanding of codon usage can give new insights into viral evolution. We implemented information entropy and revealed that the nucleotide usage pattern of ASFV is significantly associated with viral isolation factors (region and time), especially the usages of thymine and cytosine. Despite the domination of adenine and thymine in the viral genome, we found that mutation pressure alters the overall codon usage pattern of ASFV, followed by selective forces from natural selection. Moreover, the nucleotide skew index at the gene level indicates that nucleotide usages influencing synonymous codon bias of ASFV are significantly correlated with viral protein hydropathy. Finally, evolutionary plasticity is proved to contribute to the weakness in synonymous codons with A- or T-end serving as optimal codons of ASFV, suggesting that fine-tuning translation selection plays a role in synonymous codon usages of ASFV for adapting host. Taken together, ASFV is subject to evolutionary dynamics on nucleotide selections and synonymous codon usage, and our detailed analysis offers deeper insights into the genetic characteristics of this newly emerging virus around the world.