Analysis of Codon Usage Bias in the Streptococcus pneumoniae Pneumolysin Gene
Streptococcus pneumoniae pneumolysin is a key virulence factor belonging to the cholesterol‐dependent cytolysin family, enabling host cell lysis and immune evasion. While synonymous codon usage bias is known to fine‐tune virulence gene expression in pathogens, its role in pneumolysin remains uncharacterized. This study presents a comprehensive analysis of codon usage patterns in the pneumolysin gene across 420 curated coding sequences. We found a pronounced preference for A/U‐ending codons, significant underrepresentation of CpG dinucleotides, and moderate overall bias (effective number of codons, ENC = 50.28). Neutrality plot, parity rule 2 (PR2) bias, and ENC‐plot analyses collectively indicated that natural selection—not mutational pressure—is the dominant evolutionary force shaping this bias. Strikingly, pneumolysin’s codon usage showed a significant correlation with the abundant tRNA gene pool of its human host, suggesting an adaptive strategy that may minimize immunostimulation caused by bacterial mRNA release during infection. These findings reveal a balance between translational efficiency, proper protein folding, and immune evasion, thereby providing a functional understanding of pneumolysin evolution and a foundation for practical applications. These include guiding codon‐optimized heterologous expression for biochemical studies and enabling codon deoptimization for the design of safer live‐attenuated vaccines.
- Research Article
11
- 10.1186/s12985-019-1245-2
- Nov 14, 2019
- Virology Journal
BackgroundPolyomaviruses (PyVs) have a wide range of hosts, from humans to fish, and their effects on hosts vary. The differences in the infection characteristics of PyV with respect to the host are assumed to be influenced by the biochemical function of the LT-Ag protein, which is related to the cytopathic effect and tumorigenesis mechanism via interaction with the host protein.MethodsWe carried out a comparative analysis of codon usage patterns of large T-antigens (LT-Ags) of PyVs isolated from various host species and their functional domains and sequence motifs. Parity rule 2 (PR2) and neutrality analysis were applied to evaluate the effects of mutation and selection pressure on codon usage bias. To investigate evolutionary relationships among PyVs, we carried out a phylogenetic analysis, and a correspondence analysis of relative synonymous codon usage (RSCU) values was performed.ResultsNucleotide composition analysis using LT-Ag gene sequences showed that the GC and GC3 values of avian PyVs were higher than those of mammalian PyVs. The effective number of codon (ENC) analysis showed host-specific ENC distribution characteristics in both the LT-Ag gene and the coding sequences of its domain regions. In the avian and fish PyVs, the codon diversity was significant, whereas the mammalian PyVs tended to exhibit conservative and host-specific evolution of codon usage bias. The results of our PR2 and neutrality analysis revealed mutation bias or highly variable GC contents by showing a narrow GC12 distribution and wide GC3 distribution in all sequences. Furthermore, the calculated RSCU values revealed differences in the codon usage preference of the LT-AG gene according to the host group. A similar tendency was observed in the two functional domains used in the analysis.ConclusionsOur study showed that specific domains or sequence motifs of various PyV LT-Ags have evolved so that each virus protein interacts with host cell targets. They have also adapted to thrive in specific host species and cell types. Functional domains of LT-Ag, which are known to interact with host proteins involved in cell proliferation and gene expression regulation, may provide important information, as they are significantly related to the host specificity of PyVs.
- Research Article
11
- 10.1186/s12864-024-10825-x
- Oct 8, 2024
- BMC Genomics
The Halimedineae are marine green macroalgae that play crucial roles as primary producers in various habitats, including coral reefs, rocky shores, embayments, lagoons, and seagrass beds. Several tropical species have calcified thalli, which contribute significantly to the formation of coral reefs. In this study, we investigated the codon usage patterns and the main factors influencing codon usage bias in 16 chloroplast genomes of the suborder Halimedineae. Nucleotide composition analysis revealed that the codons of these species were enriched in A/U bases and preferred to end in A/U bases, and the distribution of GC content followed a trend of GC1 > GC2 > GC3. 30 optimal codons encoding 17 amino acids were identified, and most of the optimal codons and all of the over-expressed codons preferentially ended with A/U. The neutrality plot, effective number of codons (ENc) plot, and parity rule 2 (PR2) plot analysis indicated that natural selection played a major role in shaping codon usage bias of the most Halimedineae species. The genetic relationships based on their RSCU values and chloroplast protein-coding genes showed the closely related species have similar codon usage patterns. This study describes, for the first time, the codon usage patterns and characterization of Halimedineae chloroplast genomes, and provides new insights into the evolution of this suborder.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12864-024-10825-x.
- Research Article
11
- 10.1016/j.gene.2022.146842
- Aug 28, 2022
- Gene
Analysis of codon usage patterns in Haloxylon ammodendron based on genomic and transcriptomic data
- Research Article
88
- 10.1007/s11105-013-0691-z
- Dec 20, 2013
- Plant Molecular Biology Reporter
Codon usage bias (CUB) is an important evolutionary feature in a genome and has been widely documented from prokaryotes to eukaryotes. However, the significance of CUB in the Asteraceae family has not been well understood, with no Asteraceae species having been analyzed for this characteristic. Here, we use bioinformatics approaches to comparatively analyze the general patterns and influencing factors of CUB in five Asteraceae chloroplast (cp) genomes. The results indicated that the five genomes had similar codon usage patterns, showing a strong bias towards a high representation of NNA and NNT codons. Neutrality analysis showed that these cp genomes had a narrow GC distribution and no significant correlation was observed between GC12 and GC3. Parity Rule 2 (PR2) plot analysis revealed that purines were used more frequently than pyrimidines. Effective number of codons (ENc)-plot analysis showed that most genes followed the parabolic line of trajectory, but several genes with low ENc values lying below the expected curve were also observed. Furthermore, correspondence analysis of relative synonymous codon usage (RSCU) yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggested that both natural selection and mutational bias contributed to codon bias, while selection was the major force to shape the codon usage in these Asteraceae cp genomes. Our study, which is the first to investigate codon usage patterns in Asteraceae plastomes, will provide helpful information about codon distribution and variation in these species, and also shed light on the genetic and evolutionary mechanisms of codon biology within this family.
- Research Article
- 10.1007/s10528-025-11292-z
- Dec 8, 2025
- Biochemical genetics
Maple (Acer genus), deciduous perennial trees belong to the Sapindaceae family, thrive naturally across Asia, Europe and North America. In addition to their wide use in woodwork and pharmaceuticals, Acer species are favored choices for contemporary landscape design because of their large stature, intricate foliage patterns, vibrant colors, distinctive fruit shapes (paired winged samaras), and extensive genetic diversity. Understandably, they play a pivotal role in enhancing urban ecosystems and contributing to landscape construction. The aim of this study was to investigate the genomic characteristics and extent of codon usage bias (CUB) variation in the chloroplast (cp.) genome of ten sections of the genus Acer. Synonymous codon usage order (SCUO) indicated a weak codon usage bias (CUB) in the cp. genes of various Acer species. Notably, the third position of the codons exhibited a significant correlation with the overall nucleotide composition, suggesting that both natural selection and mutation pressure may have influenced the CUB. The Effective Number of Codon (ENC)-plot, Parity Rule 2 (PR2), and neutrality analyses showed that the codon bias of Acer cp. genetic material was influenced by several factors. While our analyses suggest a significant role for natural selection in shaping codon usage bias, as indicated by the ENC-plot analysis deviating from neutrality expectations, the contribution of mutation pressure and potentially genetic drift cannot be entirely excluded. These findings provide insights into codon usage patterns, which may be useful for future studies aiming to optimize codon usage for exogenous gene expression in the Acer cp. genome.
- Research Article
- 10.22438/jeb/47/2/mrn-5735
- Mar 15, 2026
- Journal of Environmental Biology
Aim: This study explores the codon usage bias (CUB) in the mitochondrial genes of two economically significant catfish species, Clarias batrachus and Heteropneustes fossilis, providing valuable insights for genetic and evolutionary research. Methodology: Effective Number of Codons (ENC), RSCU values, Compositional analysis, Correspondence and PR2 analysis, Neutrality plot, RNA editing etc., were conducted. Results: A low CUB was observed in both species, with average Effective Number of Codons (ENC) values 43.46 for C. batrachus and 44.08 for H. fossilis. In C. batrachus, 12 codons were over-represented, while 19 were under-represented; in H. fossilis, 11 codons were over-represented and 17 under-represented. Compositional analysis showed an AT-rich bias at third codon position in both the species. Correspondence and parity rule 2 (PR2) analyses suggested that mutational pressure and natural selection contribute to CUB. Neutrality plots indicated that natural selection plays a significant role than mutation in shaping CUB in the mitochondrial genes. Additionally, high C to T transitions was observed at RNA editing sites. Positive GRAVY scores pointed to the hydrophobic nature of mitochondrial proteins. Phylogenetic analysis revealed the evolutionary relationships between these species. Interpretation: This study provides valuable insights into the genetic evolution of these catfish species, with implications for aquaculture and evolutionary biology. Key words: Codon usage bias, Natural selection, Phylogenetic, RNA editing sites
- Research Article
29
- 10.1186/s12863-024-01191-4
- Jan 13, 2024
- BMC genomic data
BackgroundHemerocallis 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.ResultsHere, 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.ConclusionsThere 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.
- Research Article
7
- 10.15171/ijb.1156
- Dec 29, 2014
- Iranian Journal of Biotechnology
Background: Most of the amino acids are encoded by more than one codon, termed as synonymous codons. Synonymous codon usage is not random as it is unique to species. In each amino acid family, some synonymous codons are preferred and this is referred to as synonymous codon usage bias (SCUB). Trends associated with evolution of SCUB and factors influencing its diversification in plastomes of genetically distinct Oenothera plastomes have not been investigated so far. Objectives: In the present study, major forces that shape SCUB in Oenothera plastomes and putative preferred codons in the protein coding genes (PCG) of plastomes were identified. Materials and Methods: To unravel various features of SCUB across selected Oenothera plastomes, commonly used codon usage indices such as relative synonymous codon usage (RSCU), synonymous codon usage order (SCUO), effective number of codons (ENC) and codon adaptation index (CAI) were calculated. Correspondence analysis (COA) on RSCU was performed to identify various characteristics of SCUB across different PCG in Oenothera plastomes. Spearman’s rank correlation analysis was adopted to correlate nucleotide contents, codon usage indices and major axes of COA to find out critical parameters in shaping SCUB. Results: Mutational bias due to compositional constraints played crucial role in shaping SCUB as T3 and GC3 contents were in strong negative correlation with all axes of COA. Nevertheless, significant negative correlations between axis 1 and 3 with ENC and CAI respectively, in all species, and narrow distribution of GC contents in neutrality plot, indicate the role of natural selection. Hydropathy score of proteins was found to be influencing SCUB in O. glazioviana as it showed strong negative correlation with axis 2. Conclusion: We concluded that mutational pressure coupled with weak selection influenced SCUB in the examined plastomes of Oenothera. In addition, all examined species of Oenothera exist as disjunct populations in different parts of North America and these populations might have experienced genetic drift as random mutations in small populations that have been fixed over a period of time.
- Research Article
6
- 10.3390/su141711035
- Sep 4, 2022
- Sustainability
Codon usage bias (CUB) phenomenon varies with the species and even within the genes of the same species, where few codons are preferred more frequently than their other synonymous codons. It also categorizes the differences between species. Nucleotide compositional analysis reveals the molecular mechanisms of genes and the evolutionary relationship of a gene in dissimilar plant species. In the present study, three orthologous sequences of each FAE1 (FAE1.1, FAE1.2, and FAE1.3) and FAD2 (FAD2.1, FAD2.2, and FAD2.3) genes, from six Brassica species were accessed using the GenBank database. Further, CUB-related parameters such as nucleotide composition (AT and GC content), relative synonymous codon usage (RSCU), the effective number of codons (ENC), frequency of optimal codons (Fop), relative codon usage bias (RCBS), neutrality plot (GC12 vs. GC3), parity rule-2 [(A3/(A3 + T3) vs. (G3/(G3 + C3)], and correspondence analysis (COA) were analyzed to compare codon bias in U’s triangle Brassica species. The FAE1 genes were AT-biased and FAD2 genes were GC-biased across the studied Brassica species. RSCU values indicated that both the genes had moderate codon usage frequency for selected amino acids. The evolutionary study confirmed that codon usage preference is similar within the species grouped into the same cluster for FAE1; however, B. nigra performed differently for FAD2.2 orthologue. The high ENC value, low Fop, and RSCU value highlighted that FAE1 and FAD2 genes had a low level of gene expression and moderate preference for codon usage across the Brassicas. In addition, neutrality plot, parity rule, and correspondence analysis revealed that natural selection pressure had significantly contributed to CUB for FAE1 genes, whereas mutation and selection pressure occurred for FAD2 genes. This study would help to decode codon optimization, improve the level of expression of exogenous genes, and transgenic engineering to increase fatty acid profiling for the betterment of seed oil in Brassica species.
- Research Article
- 10.48129/kjs.11289
- Nov 11, 2021
- Kuwait Journal of Science
Comparative analysis of codon usage patterns of FUT2 from different species
- Research Article
30
- 10.3390/plants9111605
- Nov 18, 2020
- Plants
Biebersteiniaceae and Nitrariaceae, two small families, were classified in Sapindales recently. Taxonomic and phylogenetic relationships within Sapindales are still poorly resolved and controversial. In current study, we compared the chloroplast genomes of five species (Biebersteinia heterostemon, Peganum harmala, Nitraria roborowskii, Nitraria sibirica, and Nitraria tangutorum) from Biebersteiniaceae and Nitrariaceae. High similarity was detected in the gene order, content and orientation of the five chloroplast genomes; 13 highly variable regions were identified among the five species. An accelerated substitution rate was found in the protein-coding genes, especially clpP. The effective number of codons (ENC), parity rule 2 (PR2), and neutrality plots together revealed that the codon usage bias is affected by mutation and selection. The phylogenetic analysis strongly supported (Nitrariaceae (Biebersteiniaceae + The Rest)) relationships in Sapindales. Our findings can provide useful information for analyzing phylogeny and molecular evolution within Biebersteiniaceae and Nitrariaceae.
- Research Article
2
- 10.1016/j.bbrc.2023.03.030
- Mar 22, 2023
- Biochemical and Biophysical Research Communications
Ipomoea plants possess important commercial, medicinal, and ornamental value. Molecular and morphological studies have confirmed that most species of this genus exhibit similar phenotypes but complex phylogenetic relationships. To date, limited information is available on these evolutionary relationships. In this study, systematic analysis of diverse species from Ipomoea was used to elucidate the relationships in this genus. To this end, we employed the concept of codon usage bias (CUB) to analyze the codon usage bias of five Ipomoea species such as effective number of codons (ENC) and GC content at the third synonym codon position (GC3s). Three types of plots including ENC-GC3s, parity rule 2 (PR2) and neutrality plots were employed to discover the factors determining CUB, and the frequency of hydrogen bonds and nucleotide were calculated to dissect changes in GC content at the 5′-end of the coding sequence. Our results showed little distinctness in CUB among the five species, with a reduction of hydrogen bonds content at the 5′-end (with similar changes in cytosines). In addition, optimal codons of Ipomoea aquatica ended with G or C, different from those of the other four species, which ended in A or T. These results may be useful for exploring the evolutionary relationships among this group, and for understanding the reasons for the variation among Ipomoea species.
- Research Article
35
- 10.3389/fmicb.2022.1037598
- Dec 14, 2022
- Frontiers in Microbiology
Monkeypox virus (MPXV) causing multi-country outbreak-2022 is related to viruses caused outbreak-2017–2018 in West Africa. Still not fully understood which proteins of the MPXV discovered in Nigeria in 2017 have mutated through different lineages to the extent that it could cause a multi-country outbreak in 2022; similarly, codon usage bias, host adaptation indices, and the role of selection or mutation pressure in the mutated genes are also not fully studied. Here we report that according to the available sequence data this monkeypox virus acquires point mutations in multiple proteins in each period, and these point mutations accumulate and become a virus that can root outbreak-2022. Viruses exported from Nigeria to Singapore, Israel, and the United Kingdom in 2018–2019 were developed as evolutionary ancestors to B.1 viruses (MPXVs causing multi-country outbreak-2022) through MPXV/United States/2021/MD virus. Although these exported viruses have different amino acid mutations in different proteins, amino acid mutations in 10 proteins are common among them. The MPXV-United Kingdom-P2 virus evolved with only mutations in these 10 proteins and further evolved into MPXV/United States/2021/MD with amino acid mutations in 26 (including amino acid mutations in 10 proteins of the MPXV-United States-P2) proteins. It is noteworthy that specific amino acid mutations in these 22/26 (presence in MPXV/United States/2021/MD) proteins are present in B.1 viruses. Further, analysis of Relative Synonymous Codon Usage (RSCU), Synonymous Codon Usage Fraction (SCUF), and Effective Number of Codons (ENc) revealed codon usage bias in genes that exhibited nucleotide mutations in lineage B.1. Also, host adaptation indices analyzes such as Codon Adaptation Index (CAI), Expected-CAI (eCAI), Relative Codon Deoptimization Index (RCDI) and Expected value for the RCDI (eRCDI) analyzes reveal that the genes that demonstrated nucleotide mutations in lineage B.1 are favorable for human adaptation. Similarly, ENc-GC3s plot, Neutrality plot, and Parity Rule 2 (PR2)-bias plot analyzes suggest a major role of selection pressure than mutation pressure in the evolution of genes displaying nucleotide mutations in lineage B.1. Overall, from 2017 to 2022, MPXV’s mutation and spread suggests that this virus continues to evolve through point mutation in the genes according to the available sequence data.
- Research Article
10
- 10.1515/ap-2017-0090
- Jan 20, 2017
- Acta Parasitologica
Schistosoma japonicum and Ascaris suum are considered as the major parasites of human which cause various life threatening diseases such as schistomiasis and ascariasis. The codon usage bias (CUB) is known as the phenomenon of more usage of a specific codon than the other synonymous codons for an amino acid. The factors that influence the codon usage bias are mutation pressure, natural selection, gene expression, gene length, GC content, RNA stability, recombination rates, codon position etc. Here we had used various bioinformatic tools and statistical analyses to understand the compositional features, expression level and codon usage bias in the genes of these two species.After estimating the effective number of codon (ENC) in both the species, codon usage bias was found to be low and gene expression was high. The nucleobase A and T were used most often than C and G. From neutrality plot and correspondence analysis it was found that both natural selection and mutation pressure played an important role in shaping the codon usage pattern of both species. Moreover, natural selection played a major role while mutation pressure played a minor role in shaping the codon usage bias in S. japonicum and A.suum. This is the first report on the codon usage biology in S. japonicum and A.suum, and the factors influencing their codon usage bias. These results are expected to be useful for genetic engineering and evolutionary studies.
- Research Article
13
- 10.1186/s12864-021-08106-y
- Oct 30, 2021
- BMC Genomics
BackgroundKuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data.ResultsUsing CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω > 1) had a higher correlation coefficient than that of experienced purifying selection (ω < 1). Parity Rule 2 (PR2) and effective number of codons (ENc) plot analysis showed that the codon usage patterns of both species were influenced by both mutations and selection. Moreover, the average observed ENc value was lower than the expected value for both species, suggesting that factors other than GC may play roles in these phenomena. The results of multispecies clustering based on codon preference were consistent with traditional classification.ConclusionsThis study provides a relatively comprehensive understanding of the correlations among codon usage bias, gene expression, and selection pressures of CDSs for M. japonicus and M. pulchricaudatus. The genetic evolution was driven by mutations and selection pressure. Moreover, the results point out new insights into the specificities and evolutionary characteristics of the two Marsupenaeus species.