Analysis for codon usage bias in membrane anchor of nonstructural protein 5A from BVDV

The nonstructural protein 5A (NS5A) of the bovine viral diarrhea virus (BVDV) is a monotopic membrane protein. This protein can anchor to the cell membrane by an in-plane amphipathic ⍺-helix, which participates in the viral replication complex. In this study, the effects of synonymous codon usage pattern of NS5A and the overall transfer RNA (tRNA) abundance in cells on the formation of the in-plane membrane anchor of NS5A were analyzed, based on NS5A coding sequences of different BVDV genotypes. BVDV NS5A coding sequences represent the most potential for BVDV genotyping. Moreover, the nucleotide usage of BVDV NS5A dominates the genotype-specific pattern of synonymous codon usage. There is an obvious relationship between synonymous codon usage bias and the spatial conformation of the in-plane membrane anchor. Furthermore, the overall tRNA abundance profiling displays that codon positions with a high level of tRNA abundance are more than ones with a low level of tRNA abundance in the in-plane membrane anchor, implying that high translation speed probably acts on the spatial conformation of in-plane membrane anchor of BVDV NS5A. These results give a new opinion on the effect of codon usage bias in the formation of the in-plane membrane anchor of BVDV NS5A.

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.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Codon Usage Bias in Autophagy-Related Gene 13 in Eukaryotes: Uncovering the Genetic Divergence by the Interplay Between Nucleotides and Codon Usages

Synonymous codon usage bias is a universal characteristic of genomes across various organisms. Autophagy-related gene 13 (atg13) is one essential gene for autophagy initiation, yet the evolutionary trends of the atg13 gene at the usages of nucleotide and synonymous codon remains unexplored. According to phylogenetic analyses for the atg13 gene of 226 eukaryotic organisms at the nucleotide and amino acid levels, it is clear that their nucleotide usages exhibit more genetic information than their amino acid usages. Specifically, the overall nucleotide usage bias quantified by information entropy reflected that the usage biases at the first and second codon positions were stronger than those at the third position of the atg13 genes. Furthermore, the bias level of nucleotide ‘G’ usage is highest, while that of nucleotide ‘C’ usage is lowest in the atg13 genes. On top of that, genetic features represented by synonymous codon usage exhibits a species-specific pattern on the evolution of the atg13 genes to some extent. Interestingly, the codon usages of atg13 genes in the ancestor animals (Latimeria chalumnae, Petromyzon marinus, and Rhinatrema bivittatum) are strongly influenced by mutation pressure from nucleotide composition constraint. However, the distributions of nucleotide composition at different codon positions in the atg13 gene display that natural selection still dominates atg13 codon usages during organisms’ evolution.

Tracing Brucella evolutionary dynamics in expanding host ranges through nucleotide, codon and amino acid usages in genomes

The host range of Brucella organisms has expanded from terrestrial and marine mammals to fish and amphibians. The high homology genomes of different Brucella organisms promote us to investigate evolutionary patterns for nucleotide, codon and amino acid usage patterns at gene levels among Brucella species. Although the similar patterns for nucleotide and synonymous codon usages exist in gene population, GC composition at the first codon position has significant correlations to that of the second and third codon positions, respectively, suggesting that nucleotide usages surrounding one codon influence synonymous codon usage patterns. Evolutionary patterns represented by synonymous codon and amino acid usages reflect host factor impacting Brucella speciation. As for genetic variations of important virulent factors involved with different biological functions, genes encoding lipoplysaccharides (LPSs) display more distinctive codon adaptation to Brucella than those of the BvrR/BvrS system and type IV secretion system. By Bayesian analysis, the polygenetic constructions for these genes of virulent factors shared by Brucella species display the purifying/positive selections and partially host factor in mediating genetic variations of these genes. The systemic analyses for nucleotide, synonymous codon and amino acid usages at gene level and genetic variations of important virulent factor genes display that host limitation influences either genetic characterizations at gene level or a particular gene involved in virulent factors of Brucella.Communicated by Ramaswamy H. Sarma.

Synonymous Dinucleotide Usage: A Codon-Aware Metric for Quantifying Dinucleotide Representation in Viruses

Distinct patterns of dinucleotide representation, such as CpG and UpA suppression, are characteristic of certain viral genomes. Recent research has uncovered vertebrate immune mechanisms that select against specific dinucleotides in targeted viruses. This evidence highlights the importance of systematically examining the dinucleotide composition of viral genomes. We have developed a novel metric, called synonymous dinucleotide usage (SDU), for quantifying dinucleotide representation in coding sequences. Our method compares the abundance of a given dinucleotide to the null hypothesis of equal synonymous codon usage in the sequence. We present a Python3 package, DinuQ, for calculating SDU and other relevant metrics. We have applied this method on two sets of invertebrate- and vertebrate-specific flaviviruses and rhabdoviruses. The SDU shows that the vertebrate viruses exhibit consistently greater under-representation of CpG dinucleotides in all three codon positions in both datasets. In comparison to existing metrics for dinucleotide quantification, the SDU allows for a statistical interpretation of its values by comparing it to a null expectation based on the codon table. Here we apply the method to viruses, but coding sequences of other living organisms can be analysed in the same way.

Significance and roles of synonymous codon usage in the evolutionary process of Proteus

Proteus spp. bacteria frequently serve as opportunistic pathogens that can infect many animals and show positive survival and existence in various natural environments. The evolutionary pattern of Proteus spp. is an unknown topic, which benefits understanding the different evolutionary dynamics for excellent bacterial adaptation to various environments. Here, the eight whole genomes of different Proteus species were analyzed for the interplay between nucleotide usage and synonymous codon usage. Although the orthologous average nucleotide identity and average nucleotide identity display the genetic diversity of these Proteus species at the genome level, the principal component analysis further shows that these species sustain the specific genetic niche at the aspect of synonymous codon usage patterns. Interestingly, although these Proteus species have A/T rich genes with underrepresented G (guanine) or C (cytosine) at the third codon positions and overrepresented A or T at these positions, some synonymous codons with A or T end are obviously suppressed in usage. The overall codon usage pattern reflected by the effective number of codons (ENC) has a significantly positive correlation with GC3 content (GC content at the third codon position), and ENC has a significantly negative correlation with the adaptation index for these species. These results suggest that the mutation pressure caused by nucleotide composition constraint serves as a dominant evolutionary dynamic driving evolutionary trend of Proteus spp., along with other selections related to natural selection, replication and fine-tune translation, and so on. Taken together, the analyses help to understand the evolutionary interplay between nucleotide and codon usage at the gene level of Proteus.

Evidence of evolutionary selection for cotranslational folding

Recent experiments and simulations have demonstrated that proteins can fold on the ribosome. However, the extent and generality of fitness effects resulting from cotranslational folding remain open questions. Here we report a genome-wide analysis that uncovers evidence of evolutionary selection for cotranslational folding. We describe a robust statistical approach to identify loci within genes that are both significantly enriched in slowly translated codons and evolutionarily conserved. Surprisingly, we find that domain boundaries can explain only a small fraction of these conserved loci. Instead, we propose that regions enriched in slowly translated codons are associated with cotranslational folding intermediates, which may be smaller than a single domain. We show that the intermediates predicted by a native-centric model of cotranslational folding account for the majority of these loci across more than 500 Escherichia coli proteins. By making a direct connection to protein folding, this analysis provides strong evidence that many synonymous substitutions have been selected to optimize translation rates at specific locations within genes. More generally, our results indicate that kinetics, and not just thermodynamics, can significantly alter the efficiency of self-assembly in a biological context.

Evolution of Synonymous Codon Usage in the Mitogenomes of Certain Species of Bilaterian Lineage with Special Reference to Chaetognatha

Chaetognatha is a minor phylum, comprising transparent marine invertebrates varying in size from 0.5 to 12 cm. The exact phylogenetic position of Chaetognatha in Metazoa has not been deciphered as some embryological characteristics place chaetognaths among deuterostomes and some morphological characteristics place these among protostomes. In this study, the major factors that drive synonymous codon usage bias (SCUB) in the mitogenomes of representative species of Chaetognatha and chosen species of other closely related phyla were analyzed. Spearman’s rank correlation analyses of nucleotide contents suggested that mutational pressure and selection were acting in all examined mitogenomes but with varying intensities. The quantification of SCUB using effective number of codons vs. GC composition at the third codon position (GC₃) plot suggested that mutational pressure due to GC compositional constraints might be one of the major influencing forces driving the SCUB in all chaetognaths except Sagitta enflata. However, neutrality plots revealed no significant correlation between GC₃ and cumulative GC content at first and second codon positions (GC₁₂) in all other species, except in Daphnia pulex. The parity rule 2 bias plot showed that significant compositional differences existed between C and G, as well as between A and T, contents in most of the protein-coding genes (PCGs) and, comparatively, A and T contents were used more proportionally than C and G contents in all chosen mitogenomes. Chi-square analysis revealed the presence of putative optimal codons in all species, except in S. enflata. The correspondence analysis identified that mutational pressure and selection act on the mitogenomes of the selected chaetognaths and other phyla with varying intensities. The cluster analysis based on relative synonymous codon usage (RSCU) values revealed that RSCU variations in the PCGs of mitogenomes of chaetognaths are more comparable with those of protostomes. Apart from mutational pressure and selection, certain unknown selective forces might be acting on the PCGs in the analyzed mitogenomes as the phenomenon of SCUB could not be explained by mutational pressure, by selection, or by both.

Estimating Gene Expression and Codon-Specific Translational Efficiencies, Mutation Biases, and Selection Coefficients from Genomic Data Alone

Extracting biologically meaningful information from the continuing flood of genomic data is a major challenge in the life sciences. Codon usage bias (CUB) is a general feature of most genomes and is thought to reflect the effects of both natural selection for efficient translation and mutation bias. Here we present a mechanistically interpretable, Bayesian model (ribosome overhead costs Stochastic Evolutionary Model of Protein Production Rate [ROC SEMPPR]) to extract meaningful information from patterns of CUB within a genome. ROC SEMPPR is grounded in population genetics and allows us to separate the contributions of mutational biases and natural selection against translational inefficiency on a gene-by-gene and codon-by-codon basis. Until now, the primary disadvantage of similar approaches was the need for genome scale measurements of gene expression. Here, we demonstrate that it is possible to both extract accurate estimates of codon-specific mutation biases and translational efficiencies while simultaneously generating accurate estimates of gene expression, rather than requiring such information. We demonstrate the utility of ROC SEMPPR using the Saccharomyces cerevisiae S288c genome. When we compare our model fits with previous approaches we observe an exceptionally high agreement between estimates of both codon-specific parameters and gene expression levels ([Formula: see text] in all cases). We also observe strong agreement between our parameter estimates and those derived from alternative data sets. For example, our estimates of mutation bias and those from mutational accumulation experiments are highly correlated ([Formula: see text]). Our estimates of codon-specific translational inefficiencies and tRNA copy number-based estimates of ribosome pausing time ([Formula: see text]), and mRNA and ribosome profiling footprint-based estimates of gene expression ([Formula: see text]) are also highly correlated, thus supporting the hypothesis that selection against translational inefficiency is an important force driving the evolution of CUB. Surprisingly, we find that for particular amino acids, codon usage in highly expressed genes can still be largely driven by mutation bias and that failing to take mutation bias into account can lead to the misidentification of an amino acid's "optimal" codon. In conclusion, our method demonstrates that an enormous amount of biologically important information is encoded within genome scale patterns of codon usage, accessing this information does not require gene expression measurements, but instead carefully formulated biologically interpretable models.

Synonymous codon usage affects the expression of wild type and F508del CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel composed of 1480 amino acids. The major mutation responsible for cystic fibrosis results in loss of amino acid residue, F508 (F508del). Loss of F508 in CFTR alters the folding pathway resulting in endoplasmic-reticulum-associated degradation. This study investigates the role of synonymous codon in the expression of CFTR and CFTR F508del in human HEK293 cells. DNA encoding the open reading frame (ORF) for CFTR containing synonymous codon replacements was expressed using a heterologous vector integrated into the genome. The results indicate that the codon usage greatly affects the expression of CFTR. While the promoter strength driving expression of the ORFs was largely unchanged and the mRNA half-lives were unchanged, the steady-state levels of the mRNA varied by as much as 30-fold. Experiments support that this apparent inconsistency is attributed to nonsense mediated decay independent of exon junction complex. The ratio of CFTR/mRNA indicates that mRNA containing native codons was more efficient in expressing mature CFTR as compared to mRNA containing synonymous high-expression codons. However, when F508del CFTR was expressed after codon optimization, a greater percentage of the protein escaped endoplasmic-reticulum-associated degradation resulting in considerable levels of mature F508del CFTR on the plasma membrane, which showed channel activity. These results indicate that codon usage has an effect on mRNA levels and protein expression, for CFTR, and likely on chaperone-assisted folding pathway, for F508del CFTR.

Analysis of codon usage preference in hemagglutinin genes of the swine-origin influenza A (H1N1) virus

BACKGROUND

The swine-origin influenza A (H1N1) virus (S-OIV) has come to the forefront since 2009 and was identified as a new reassortant strain. The hemagglutinin (HA) glycoprotein mediates virus binding, contains antigenic regions recognized by neutralizing antibodies, and is associated with viral cross-species infection and adaption. The comparison study of codon usage preferences in influenza viral genomes was less extensive. In this study, we used codon usage pattern analyses to validate the adaption and origins of S-OIV.

METHODS

Codon usage pattern was used to estimate the host adaption of S-OIVs. Phylogenetic analysis of the HA gene was conducted to understand the phylogeny of H1N1 viruses isolated from different hosts. Amino acid signature pattern on antigenic sites of HA was analyzed to understand the antigenic characteristics.

RESULTS

Results of phylogenetic analyses of HA gene indicate that S-OIVs group in identical clusters. The synonymous codon usage pattern analyses indicate that the effective number of codons versus GC content at the third codon position in the HA1 gene slightly differ from those in swine H1N1 and gradually adapted to human. Our data indicate that S-OIV evolution occurred according to positive selection within these antigenic regions. A comparison of antigenic site amino acids reveals similar signature patterns between S-OIV and 1918 human influenza strains.

CONCLUSION

This study proposes a new and effective way to gain a better understanding of the features of the S-OIV genome and evolutionary processes based on the codon usage pattern. It is useful to trace influenza viral origins and cross-species virus transmission.

Translational selection frequently overcomes genetic drift in shaping synonymous codon usage patterns in vertebrates

Synonymous codon usage patterns are shaped by a balance between mutation, drift, and natural selection. To date, detection of translational selection in vertebrates has proven to be a challenging task, obscured by small long-term effective population sizes in larger animals and the existence of isochores in some species. The consensus is that, in such species, natural selection is either completely ineffective at overcoming mutational pressures and genetic drift or perhaps is effective but so weak that it is not detectable. The aim of this research is to understand the interplay between mutation, selection, and genetic drift in vertebrates. We observe that although variation in mutational bias is undoubtedly the dominant force influencing codon usage, translational selection acts as a weak additional factor influencing synonymous codon usage. These observations indicate that translational selection is a widespread phenomenon in vertebrates and is not limited to a few species.

Mutation pressure shapes codon usage in the GC-Rich genome of foot-and-mouth disease virus

Foot-and-mouth disease (FMD) is economically the most important viral-induced livestock disease worldwide. In this study, we report the results of a survey of codon usage bias of FMD virus (FMDV) representing all seven serotypes (A, O, C, Asia 1, SAT 1, SAT 2, and SAT 3). Correspondence analysis, a commonly used multivariate statistical approach, was carried out to analyze synonymous codon usage bias. The analysis showed that the overall extent of codon usage bias in FMDV is low. Furthermore, the good correlation between the frequency of G + C at the synonymous third position of sense codons (GC3S) content at silent sites of each sequence and codon usage bias suggested that mutation pressure rather than natural (translational) selection is the most important determinant of the codon bias observed. In addition, other factors, such as the lengths of open reading frame (ORF) and the hydrophobicity of genes also influence the codon usage variation among the genomes of FMDV in a minor way. The result of phylogenetic analyses based on the relative synonymous codon usage (RSCU) values indicated a few obvious phylogenetic incongruities, which suggest that more FMDV genome diversity may exist in nature than is currently indicated. Our work might give some clues to the features of FMDV genome and some evolutionary information of this virus.