Non-uniqueness of factors constraint on the codon usage in Bombyx mori

Analysis of Codon Usage Bias in the Plastid Genome of Diplandrorchis sinica (Orchidaceae)

In order to understand the bias and main affecting factors of codon usage in the plastid genome of Diplandrorchis sinica, which is a rare and endangered plant species in the Orchidaceae family, the complete plastid genome sequence of D. sinica was downloaded from the GenBank database and 20 protein-coding sequences that met the analysis requirements were finally selected. The GC content, length of the amino acid (Laa), relative synonymous codon usage (RSCU), and effective number of codon (ENC) of each gene and codon were calculated using the CodonW and EMBOSS online programs. Neutral plot analysis, ENC-plot analysis, PR2-plot analysis, and correspondence analysis were performed using Origin Pro 2024 software, and correlation analysis between various indicators was performed using SPSS 23.0 software. The results showed that the third base of the codon in the plastid genome of D. sinica was rich in A and T, with a GC3 content of 27%, which was lower than that of GC1 (45%) and GC2 (39%). The ENC value ranged from 35 to 57, with an average of 47. The codon usage bias was relatively low, and there was a significant positive correlation between ENC and GC3. There were a total of 32 codons with RSCU values greater than 1, of which 30 ended with either A or U. There were a total of nine optimal codons identified, namely, UCU, UCC, UCA, GCA, UUG, AUA, CGU, CGA, and GGU. This study indicated that the dominant factor affecting codon usage bias in the plastid genome of D. sinica was natural selection pressure, while the impact of base mutations was limited. The codon usage patterns were not closely related to gene types, and the distribution of photosynthetic system genes and ribosomal protein-coding gene loci was relatively scattered, indicating significant differences in the usage patterns of these gene codons. In addition, the codon usage patterns may not be related to whether the plant is a photosynthetic autotrophic or heterotrophic nutritional type. The results of this study could provide scientific references for the genomic evolution and phylogenetic research of plant species in the family Orchidaceae.

Understanding the nucleotide composition and patterns of codon usage in the expression of human oral cancer genes

Oral squamous cell carcinoma (OSCC) is primarily known as oral cancer (OC) that mostly occurs in mouth, lips and tongue. Mutations in some of the genes cause OC and some genes are risk factors for progression of OC. In this study, we analyzed the compositional features and pattern of codon usage in genes involved in OC using computational method as no work was reported yet. Compositional features suggested that the overall GC content was higher i.e. genes were GC rich. Effective number of codons (ENC) values ranged from 34.6 to 55.9 with a mean value of 49.03±4.22 representing low codon usage bias (CUB). Correspondence analysis (COA) suggested that the codon usage pattern was different in different genes. In genes associated with OC, highly significant correlation was observed between GC12 and GC3 (r=0.454, p<0.01) suggesting that directional mutation affected all the three codon positions. This is the first report on pattern of codon usage pattern on genes involved in OC, which not only alludes a new perspective for elucidating the mechanisms of biased usage of synonymous codons but also provide valuable clues for molecular genetic engineering.

Comparative analysis of the organelle genomes of Aconitum carmichaelii revealed structural and sequence differences and phylogenetic relationships

In this study, we conducted an assembly and analysis of the organelle genomes of Aconitum carmichaelii. Our investigation encompassed the examination of organelle genome structures, gene transfer events, and the environmental selection pressures affecting A. carmichaelii. The results revealed distinct evolutionary patterns in the organelle genomes of A. carmichaelii. Especially, the plastome exhibited a more conserved structure but a higher nucleotide substitution rate (NSR), while the mitogenome displayed a more complex structure with a slower NSR. Through homology analysis, we identified several instances of unidirectional protein-coding genes (PCGs) transferring from the plastome to the mitogenome. However, we did not observe any events which genes moved from the mitogenome to the plastome. Additionally, we observed multiple transposable element (TE) fragments in the organelle genomes, with both organelles showing different preferences for the type of nuclear TE insertion. Divergence time estimation suggested that rapid differentiation occurred in Aconitum species approximately 7.96 million years ago (Mya). This divergence might be associated with the reduction in CO2 levels and the significant uplift of the Qinghai-Tibet Plateau (QTP) during the late Miocene. Selection pressure analysis indicated that the dN/dS values of both organelles were less than 1, suggested that organelle PCGs were subject to purification selection. However, we did not detect any positively selected genes (PSGs) in Subg. Aconitum and Subg. Lycoctonum. This observation further supports the idea that stronger negative selection pressure on organelle genes in Aconitum results in a more conserved amino acid sequence. In conclusion, this study contributes to a deeper understanding of organelle evolution in Aconitum species and provides a foundation for future research on the genetic mechanisms underlying the structure and function of the Aconitum plastome and mitogenome.

Molecular Dissection of Herpes Simplex Virus Type 1 to Elucidate Molecular Mechanisms Behind Latency and Comparison of Its Codon Usage Patterns with Genes Modulated During Alzheimer's Disease as a Part of Host-Pathogen Interaction

BACKGROUND

Herpes simplex virus type 1 (HSV-1) is associated with Alzheimer's disease, which goes into a cycle of latency and reactivation. The present study was envisaged to understand the reasons for latency and specific molecular patterns present in the HSV-1.

OBJECTIVE

The objective is the molecular dissection of Herpes simplex virus type 1 to elucidate molecular mechanisms behind latency and compare its codon usage patterns with genes modulated during Alzheimer's disease as a part of host-pathogen interaction.

METHODS

In the present study, we tried to investigate the potential reasons for the latency of HSV-1 virus bioinformatically by determining the CpG patterns. Also, we investigated the codon usage pattern, the presence of rare codons, codon context, and protein properties.

RESULTS

The top 222 codon pairs graded based on their frequency in the HSV-1 genome revealed that with only one exception (CUG-UUU), all other codon pairs have codons ending with G/C. Considering it an extension of host-pathogen interaction, we compared HSV-1 codon usage with that of codon usage of genes modulated during Alzheimer's disease, and we found that CGT and TTT are only two codons that exhibited similar codon usage patterns and other codons showed statistically highly significant different codon preferences. Dinucleotide CpG tends to mutate to TpG, suggesting the presence of mutational forces and the imperative role of CpG methylation in HSV-1 latency.

CONCLUSIONS

Upon comparison of codon usage between HSV-1 and Alzheimer's disease genes, no similarities in codon usage were found as a part of host-pathogen interaction. CpG methylation plays an imperative role in latency HSV-1.

The m5 C methyltransferase NSUN2 promotes codon-dependent oncogenic translation by stabilising tRNA in anaplastic thyroid cancer

BACKGROUND

Translation dysregulation plays a crucial role in tumourigenesis and cancer progression. Oncogenic translation relies on the stability and availability of tRNAs for protein synthesis, making them potential targets for cancer therapy.

METHODS

This study performed immunohistochemistry analysis to assess NSUN2 levels in thyroid cancer. Furthermore, to elucidate the impact of NSUN2 on anaplastic thyroid cancer (ATC) malignancy, phenotypic assays were conducted. Drug inhibition and time-dependent plots were employed to analyse drug resistance. Liquid chromatography-mass spectrometry and bisulphite sequencing were used to investigate the m5 C methylation of tRNA at both global and single-base levels. Puromycin intake and high-frequency codon reporter assays verified the protein translation level. By combining mRNA and ribosome profiling, a series of downstream proteins and codon usage bias were identified. The acquired data were further validated by tRNA sequencing.

RESULTS

This study observed that the tRNA m5 C methyltransferase NSUN2 was up-regulated in ATC and is associated with dedifferentiation. Furthermore, NSUN2 knockdown repressed ATC formation, proliferation, invasion and migration both in vivo and in vitro. Moreover, NSUN2 repression enhanced the sensitivity of ATC to genotoxic drugs. Mechanically, NSUN2 catalyses tRNA structure-related m5 C modification, stabilising tRNA that maintains homeostasis and rapidly transports amino acids, particularly leucine. This stable tRNA has a substantially increased efficiency necessary to support a pro-cancer translation program including c-Myc, BCL2, RAB31, JUNB and TRAF2. Additionally, the NSUN2-mediated variations in m5C levels and different tRNA Leu iso-decoder families, partially contribute to a codon-dependent translation bias. Surprisingly, targeting NSUN2 disrupted the c-Myc to NSUN2 cycle in ATC.

CONCLUSIONS

This research revealed that a pro-tumour m5C methyltransferase, dynamic tRNA stability regulation and downstream oncogenes, c-Myc, elicits a codon-dependent oncogenic translation network that enhances ATC growth and formation. Furthermore, it provides new opportunities for targeting translation reprogramming in cancer cells.

Mutation and codon bias analysis of the spike protein of Omicron, the recent variant of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant is a heavily mutated virus and designated as a variant of concern. To investigate the codon usage pattern of this new variant, we performed mutation and codon bias analysis for Omicron as well as for its sub-lineages BA.1 and BA.2 and compared them with the original SARS-CoV-2 and the Delta variant sequences obtained in this study. Our results indicate that the sub-lineage BA.1 and BA.2 have up to 23 sites of difference on the spike protein, which have minimal impact on function. The Omicron variant and its sub-lineages have similar codon usage patterns and A/U ending codons appear to be preferred over G/C ending codons. The Omicron has a lower degree of codon usage bias in spite of evidence that natural selection, mutation pressure and dinucleotide abundance shape the codon usage bias of Omicron, with natural selection being more significant on BA.2 than the other sub-lineages of Omicron. The codon usage pattern of Omicron variant that we explored provides valid information for a clearer understanding of Omicron and its sub-lineages, which could find application in vaccine development and optimization.

Identification and characterization of codon usage pattern and influencing factors in HFRS-causing hantaviruses

Hemorrhagic fever with renal syndrome (HFRS) is an acute viral zoonosis carried and transmitted by infected rodents through urine, droppings, or saliva. The etiology of HFRS is complex due to the involvement of viral factors and host immune and genetic factors which hinder the development of potential therapeutic solutions for HFRS. Hantaan virus (HTNV), Dobrava-Belgrade virus (DOBV), Seoul virus (SEOV), and Puumala virus (PUUV) are predominantly found in hantaviral species that cause HFRS in patients. Despite ongoing prevention and control efforts, HFRS remains a serious economic burden worldwide. Furthermore, recent studies reported that the hantavirus nucleocapsid protein is a multi-functional protein and plays a major role in the replication cycle of the hantavirus. However, the precise mechanism of the nucleoproteins in viral pathogenesis is not completely understood. In the framework of the current study, various in silico approaches were employed to identify the factors influencing the codon usage pattern of hantaviral nucleoproteins. Based on the relative synonymous codon usage (RSCU) values, a comparative analysis was performed between HFRS-causing hantavirus and their hosts, suggesting that HTNV, DOBV, SEOV, and PUUV, were inclined to evolve their codon usage patterns that were comparable to those of their hosts. The results indicated that most of the overrepresented codons had AU-endings, which revealed that mutational pressure is the major force shaping codon usage patterns. However, the influence of natural selection and geographical factors cannot be ignored on viral codon usage bias. Further analysis also demonstrated that HFRS causing hantaviruses adapted host-specific codon usage patterns to sustain successful replication and transmission chains within hosts. To our knowledge, no study to date reported the factors influencing the codon usage pattern within hantaviral nucleoproteins. Thus, the proposed computational scheme can help in understanding the underlying mechanism of codon usage patterns in HFRS-causing hantaviruses which lend a helping hand in designing effective anti-HFRS treatments in future. This study, although comprehensive, relies on in silico methods and thus necessitates experimental validation for more solid outcomes. Beyond the identified factors influencing viral behavior, there could be other yet undiscovered influences. These potential factors should be targets for further research to improve HFRS therapeutic strategies.

Structural Characteristics and Phylogenetic Analysis of the Mitochondrial Genomes of Four Krisna Species (Hemiptera: Cicadellidae: Iassinae)

Krisna species are insects that have piercing-sucking mouthparts and belong to the Krisnini tribe in the Iassinae subfamily of leafhoppers in the Cicadellidae family. In this study, we sequenced and compared the mitochondrial genomes (mitogenomes) of four Krisna species. The results showed that all four mitogenomes were composed of cyclic double-stranded molecules and contained 13 protein-coding genes (PCGs) and 22 and 2 genes coding for tRNAs and rRNAs, respectively. Those mitogenomes exhibited similar base composition, gene size, and codon usage patterns for the protein-coding genes. The analysis of the nonsynonymous substitution rate (Ka)/synonymous substitution rate (Ks) showed that evolution occurred the fastest in ND4 and the slowest in COI. 13 PCGs that underwent purification selection were suitable for studying phylogenetic relationships within Krisna. ND2, ND6, and ATP6 had highly variable nucleotide diversity, whereas COI and ND1 exhibited the lowest diversity. Genes or gene regions with high nucleotide diversity can provide potential marker candidates for population genetics and species delimitation in Krisna. Analyses of parity and neutral plots showed that both natural selection and mutation pressure affected the codon usage bias. In the phylogenetic analysis, all subfamilies were restored to a monophyletic group; the Krisnini tribe is monophyletic, and the Krisna genus is paraphyletic. Our study provides novel insights into the significance of the background nucleotide composition and codon usage patterns in the CDSs of the 13 mitochondrial PCGs of the Krisna genome, which could enable the identification of a different gene organization and may be used for accurate phylogenetic analysis of Krisna species.

Deep decoding of codon usage strategies and host adaption preferences of soybean mosaic virus

Soybean mosaic virus (SMV) has threatened the global yield of Leguminosae crops, but the mechanism of its infection, spread, and evolution remains unknown. A systemic analysis of 107 SMV strains was performed to explore the genome-wide codon usage profile and the various factors influencing the codon usage patterns of SMV, which provides insight into its molecular evolution and elucidates its unknown host adaptation pattern. The overall nucleotide composition and correlation analysis revealed that the preferred synonymous codons mostly end with A/U. Clustering by RSCU value of each strain and phylogenetic tree analysis showed that the SMV isolates studied were divided into four clades, with a low overall extent of codon usage bias (CUB) in SMV. According to the ENC, PR2, neutrality plot, and correspondence analysis, natural selection of geographical diversity may play a critical role in the CUB. Higher adaptability was shown in Glycine with SMV and more pressure was received by clade III. These findings could not only provide valuable information about the overall codon usage pattern of the SMV genome, but could also aid in the clarification of the involved mechanisms that dominate the codon usage patterns and genetic evolution of the SMV genome.

Analysis of codon usage patterns in Haloxylon ammodendron based on genomic and transcriptomic data

Haloxylon ammodendron, a xero-halophytic shrub of Chenopodiaceae, is a dominant species in deserts, which has a strong drought and salt tolerance and plays an important role in sand fixation. However, the codon usage bias (CUB) in H. ammodendron is still unclear at present. In this study, the codon usage patterns of 38,657 coding sequences (CDSs) in the newly released whole-genome sequence data of H. ammodendron and 3,948 CDSs in the previously obtained transcriptome sequencing data were compared and analyzed. The results showed that the CDSs with the total guanineandcytosine(GC)content in the range of 40% ∼ 45% was the most in the genome and transcriptome. Among which, the GC1, GC2, and GC3 contents of genomic CDSs were 50.83%, 40.56%, and 40.23%, respectively, and those of CDSs in the transcriptome were 47.16%, 39.02%, and 39.59%, respectively. Therefore, the bases in H. ammodendron were rich in adenine and thymine, and the overallcodonusage was biasedtoward A- and U-ending codons. The analysis of neutrality plot, effective number of codon (ENC) plot, and parity rule 2 (PR2) bias plot showed that both natural selection and mutation pressure had great influences on the CUB of H. ammodendron, but natural selection was the most important determinant. Besides, gene expression level and the function and protein length of some specific genes also had influences on the codon usage pattern. Finally, a total of 25 common optimal codons were found in the genomic and transcriptomic data, and AU/GC-ending codons ratio was 24:1. It should be noted that the salt-tolerant unigenes had similar codon usage, and the highly expressed genes had higher usage frequency of optimal codons and lower GC content than the lowly expressed genes. In addition, there was no difference in the ENC values of salt-tolerant unigenes in H. ammodendron, and the expression level of the genes had no correlation with CAI. This study will help to elucidate the formation mechanism of H. ammodendron codon usage bias, and make contributions to the identification of new genes and the genetic engineering study on H. ammodendron.

Analysis of evolutionary and genetic patterns in structural genes of primate lentiviruses

Background

Primate lentiviruses (HIV1, HIV2, and Simian immunodeficiency virus [SIV]) cause immune deficiency, encephalitis, and infectious anemia in mammals such as cattle, cat, goat, sheep, horse, and puma.

Objective

This study was designed and conducted with the main purpose of confirming the overall codon usage pattern of primate lentiviruses and exploring the evolutionary and genetic characteristics commonly or specifically expressed in HIV1, HIV2, and SIV.

Methods

The gag, pol, and env gene sequences of HIV1, HIV2, and SIV were analyzed to determine their evolutionary relationships, nucleotide compositions, codon usage patterns, neutrality, selection pressure (influence of mutational pressure and natural selection), and viral adaptation to human codon usage.

Results

A strong ‘A’ bias was confirmed in all three structural genes, consistent with previous findings regarding HIV. Notably, the ENC-GC3s plot and neutral evolution analysis showed that all primate lentiviruses were more affected by selection pressure than by mutation caused by the GC composition of the gene, consistent with prior reports regarding HIV1. The overall codon usage bias of pol was highest among the structural genes, while the codon usage bias of env was lowest. The virus groups showing high codon bias in all three genes were HIV1 and SIVcolobus. The codon adaptation index (CAI) and similarity D(A, B) values indicated that although there was a high degree of similarity to human codon usage in all three structural genes of HIV, this similarity was not caused by translation pressure. In addition, compared with HIV1, the codon usage of HIV2 is more similar to the human codon usage, but the overall codon usage bias is lower.

Conclusion

The origin viruses of HIV (SIVcpz_gor and SIVsmm) exhibit greater similarity to human codon usage in the gag gene, confirming their robust adaptability to human codon usage. Therefore, HIV1 and HIV2 may have evolved to avoid human codon usage by selection pressure in the gag gene after interspecies transmission from SIV hosts to humans. By overcoming safety and stability issues, information from codon usage analysis will be useful for attenuated HIV1 vaccine development. A recoded HIV1 variant can be used as a vaccine vector or in immunotherapy to induce specific innate immune responses. Further research regarding HIV1 dinucleotide usage and codon pair usage will facilitate new approaches to the treatment of AIDS.

Strong Selectional Forces Fine-Tune CpG Content in Genes Involved in Neurological Disorders as Revealed by Codon Usage Patterns

Neurodegenerative disorders cause irreversible damage to the neurons and adversely affect the quality of life. Protein misfolding and their aggregation in specific parts of the brain, mitochondrial dysfunction, calcium load, proteolytic stress, and oxidative stress are among the causes of neurodegenerative disorders. In addition, altered metabolism has been associated with neurodegeneration as evidenced by reductions in glutamine and alanine in transient global amnesia patients, higher homocysteine-cysteine disulfide, and lower methionine decline in serum urea have been observed in Alzheimer's disease patients. Neurodegeneration thus appears to be a culmination of altered metabolism. The study's objective is to analyze various attributes like composition, physical properties of the protein, and factors like selectional and mutational forces, influencing codon usage preferences in a panel of genes involved directly or indirectly in metabolism and contributing to neurodegeneration. Various parameters, including gene composition, dinucleotide analysis, Relative synonymous codon usage (RSCU), Codon adaptation index (CAI), neutrality and parity plots, and different protein indices, were computed and analyzed to determine the codon usage pattern and factors affecting it. The correlation of intrinsic protein properties such as the grand average of hydropathicity index (GRAVY), isoelectric point, hydrophobicity, and acidic, basic, and neutral amino acid content has been found to influence codon usage. In genes up to 800 amino acids long, the GC3 content was highly variable, while GC12 content was relatively constant. An optimum CpG content is present in genes to maintain a high expression level as required for genes involved in metabolism. Also observed was a low codon usage bias with a higher protein expression level. Compositional parameters and nucleotides at the second position of codons played essential roles in explaining the extent of bias. Overall analysis indicated that the dominance of selection pressure and compositional constraints and mutational forces shape codon usage.

Human genes with codon usage bias similar to that of the nonstructural protein 1 gene of influenza A viruses are conjointly involved in the infectious pathogenesis of influenza A viruses

Molecular mechanisms of the non-structural protein 1 (NS1) in influenza A-induced pathological changes remain ambiguous. This study explored the pathogenesis of human infection by influenza A viruses (IAVs) through identifying human genes with codon usage bias (CUB) similar to NS1 gene of these viruses based on the relative synonymous codon usage (RSCU). CUB of the IAV subtypes H1N1, H3N2, H3N8, H5N1, H5N2, H5N8, H7N9 and H9N2 was analyzed and the correlation of RSCU values of NS1 sequences with those of the human genes was calculated. The CUB of NS1 was uneven and codons ending with A/U were preferred. The ENC-GC3 and neutrality plots suggested natural selection as the main determinant for CUB. The RCDI, CAI and SiD values showed that the viruses had a high degree of adaptability to human. A total of 2155 human genes showed significant RSCU-based correlation (p < 0.05 and r > 0.5) with NS1 coding sequences and was considered as human genes with CUB similar to NS1 gene of IAV subtypes. Differences and similarities in the subtype-specific human protein–protein interaction (PPI) networks and their functions were recorded among IAVs subtypes, indicating that NS1 of each IAV subtype has a specific pathogenic mechanism. Processes and pathways involved in influenza, transcription, immune response and cell cycle were enriched in human gene sets retrieved based on the CUB of NS1 gene of IAV subtypes. The present work may advance our understanding on the mechanism of NS1 in human infections of IAV subtypes and shed light on the therapeutic options.

Codon usage patterns across seven Rosales species

BACKGROUND

Codon usage bias (CUB) analysis is an effective method for studying specificity, evolutionary relationships, and mRNA translation and discovering new genes among various species. In general, CUB analysis is mainly performed within one species or between closely related species and no such study has been applied among species with distant genetic relationships. Here, seven Rosales species with high economic value were selected to conduct CUB analysis.

RESULTS

The results showed that the average GC1, GC2 and GC3 contents were 51.08, 40.52 and 43.12%, respectively, indicating that the A/T content is more abundant and the Rosales species prefer A/T as the last codon. Neutrality plot and ENc plot analysis revealed that natural selection was the main factor leading to CUB during the evolution of Rosales species. All 7 Rosales species contained three high-frequency codons, AGA, GTT and TTG, encoding Arg, Val and Leu, respectively. The 7 Rosales species differed in high-frequency codon pairs and the distribution of GC3, though the usage patterns of closely related species were more consistent. The results of the biclustering heat map among 7 Rosales species and 20 other species were basically consistent with the results of genome data, suggesting that CUB analysis is an effective method for revealing evolutionary relationships among species at the family or order level. In addition, chlorophytes prefer using G/C as ending codon, while monocotyledonous and dicotyledonous plants prefer using A/T as ending codon.

CONCLUSIONS

The CUB pattern among Rosales species was mainly affected by natural selection. This work is the first to highlight the CUB patterns and characteristics of Rosales species and provides a new perspective for studying genetic relationships across a wide range of species.

Codon usage signatures in Sabia and Chapare for host adaptation

Sabia and Chapare viruses in the Arenavirus family cause viral hemorrhagic fever among humans with a fatality rate of 30% with no treatment models. Therefore, it is of interest to document the codon usage, amino acid patterns and associated factors influencing the observed variations in Sabia and Chapare viruses for host adaptation. Multivariate statistical analysis revealed compositional constraint and host selection pressure influencing the viral codon usage patterns. These data suggests the codon usage signatures in Sabia and Chapare viruses for host adaptation in the human host implying its role in the rapid progression of the infection. Dinucleotides UpG and CpA were noted to be over-represented among the Sabia, Chapare viruses and human genomes. Strong restraint from the usage of CpG dinucleotides among viruses is linked with the molecular mimicry of the human immune system. Thus, the data reported from this study help in understanding the mechanism of viral adaptation inside the host genome for further consideration in drug discovery.

Comparative analysis of transcriptomic data shows the effects of multiple evolutionary selection processes on codon usage in Marsupenaeus japonicus and Marsupenaeus pulchricaudatus

BACKGROUND

Kuruma 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.

RESULTS

Using 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.

CONCLUSIONS

This 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.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Analysis of codon usage bias in mitochondrial CO gene among platyhelminthes

The phenomenon of non-uniform usage of the synonymous codons, where some codons are given more preference to others, is known as codon usage bias (CUB). CUB is known to be determined by two major evolutionary forces i.e. mutation pressure and selection. We used various approaches to understand the codon usage pattern in mitochondrial CO (MT-CO) genes involved in complex IV of the respiratory chain (RC) as no work was reported yet. Our present study revealed that CUB was relatively high and the coding sequences were rich in A and T. Correspondence analysis further indicated that A/T compositional properties under mutational pressure might be affecting the codon usage pattern and was different in different classes for MT-CO gene. A highly significant correlation between A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% in all the classes indicated that compositional constraints under mutational pressure and natural selection might affect the CUB. Neutrality plot indicated that both natural selection and mutational bias affected the CUB, where, natural selection played the major role as compared to mutational pressure.

Free energy of mRNA positively correlates with GC content in chloroplast transcriptomes of edible legumes

In the present study, the results of nucleotide composition analysis showed that the legume chloroplast (cp) transcriptomes were AT rich. From the neutrality plot, we observed that natural selection might have played a major role, while mutation pressure played a minor role in the CUB of cp transcriptomes. Highly significant (p < 0.05) negative correlation was found between mRNA free energy (mFE) and scaled chi-square for entire mRNA in Cicer arietinum and Lens culinaris suggesting that the release of higher energy by entire mRNA molecule might be associated with higher degree of codon usage bias in these two crop plants. Further, highly significant (p < 0.01, p < 0.05) positive correlation of mFE for entire mRNA was found with GC3 and that of mFE for 39 bases with GC, GC1, GC2 and GC3 contents among all the legumes. This indicated that higher GC content might induce the release of more free energy by cp transcriptomes.

Codon usage patterns and evolution of HSP60 in birds

Heat shock protein 60 (HSP60) is highly conserved from prokaryotic to eukaryotic organisms, acting as molecular chaperone and other vital biological functions. In spite of increasing knowledge of HSP60, its evolutionary mechanism on functional adaption is still far from completely understood. Moreover, analysis of codon usage bias (CUB) is a powerful tool to understand evolutionary association studies. However, so far, as we know, no scientific work on CUB of HSP60 in birds has been reported. In this study, we provide a comprehensive analysis on the codon usage and molecular evolution of HSP60 across birds. The results indicated that HSP60 had a weak codon usage bias with high ENC values (range from 52.66 to 61), low RSCU, and A/T-ending codons were mostly preferred. Meanwhile, it was considered that mutation, natural selection, and genetic drift combined to shape codon usage patterns with different strength proportions among various birds for HSP60. Then, the LRT tests suggested that different lineages of birds might be under similar selective pressures. Besides, the two positive selection sites (151 and 131) were detected and might undergo radical changes. These findings would contribute to understand function diversity and molecular evolution of HSP60 in birds.

An interplay between compositional constraint and natural selection dictates the codon usage pattern among select Galliformes

Galliformes are believed to be the first avian order that started living in human association and became domesticated. Members of this order ranged from common to rare species. Next-generation sequencing has availed researchers with the whole genome sequences of five Galliformes; chicken, helmeted Guinea fowl, turkey, Japanese quail, and peafowl. Bioinformatic analysis based on codon usage, evolution, and species-specific functional enrichment can provide some crucial information aiding proper understanding of their genomic strategies. In this study, we investigated the genomic features of chicken, helmeted guinea fowl, turkey, and Japanese quail. Their genomes were AT biased although the potentially highly expressed genes contained more GC than AT. Cytosine dominated the third position of frequently used optimal codons. Mutational pressures in the analyzed Galliformes were in the range of 0.2-0.6%. Neutrality plot, translational selection index, and mutational responsive index indicated the dominance of selection pressure over mutational pressure among Galliformes. A pair of di-nucleotides, TpA and CpG, was found to be used less frequently than others in protein-coding genes since both of them are associated with the conversion of euchromatin to heterochromatin. Functional enrichment analysis revealed the dominance of proteins associated with fundamental biological processes. In turkey, chicken and helmeted Guinea fowl proteins with immunity-boosting capacity prevailed along with proteins needed for signal transduction and maintenance of central dogma. Evolutionary analysis indicated a bias towards synonymous substitution than non-synonymous mutation.

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.

Analysis of codon usage patterns in citrus based on coding sequence data

BACKGROUND

Codon usage is an important determinant of gene expression levels that can help us understand codon biology, evolution and mRNA translation of species. The majority of previous codon usage studies have focused on single species analysis, although few studies have focused on the species within the same genus. In this study, we proposed a multispecies codon usage analysis workflow to reveal the genetic features and correlation in citrus.

RESULTS

Our codon usage analysis workflow was based on the GC content, GC plot, and relative synonymous codon usage value of each codon in 8 citrus species. This approach allows for the comparison of codon usage bias of different citrus species. Next, we performed cluster analysis and obtained an overview of the relationship in citrus. However, traditional methods cannot conduct quantitative analysis of the correlation. To further estimate the correlation among the citrus species, we used the frequency profile to construct feature vectors of each species. The Pearson correlation coefficient was used to quantitatively analyze the distance among the citrus species. This result was consistent with the cluster analysis.

CONCLUSIONS

Our findings showed that the citrus species are conserved at the genetic level and demonstrated the existing genetic evolutionary relationship in citrus. This work provides new insights into codon biology and the evolution of citrus and other plant species.

Codon usage signatures in the genus Cryptococcus: A complex interplay of gene expression, translational selection and compositional bias

The fungal genus Cryptococcus comprises of several diverse species. The pathogens forming Cryptococcus neoformans/ Cryptococcus gatti species complex are of immense clinical significance owing to the high frequency of infections and deaths globally. Three closely related non-pathogenic species namely, Cryptococcus amylolentus, Cryptococcus wingfieldii and Cryptococcus depauperatus are the non-pathogenic ancestral species from which pathogenic lineages have diverged. In the current study, a comprehensive analysis of factors influencing the codon and amino acid usage bias in six pathogenic and three non-pathogenic species was performed. Our results revealed that though compositional bias played a crucial role, translational selection and gene expression were the key determinants of codon usage variations. Analysis of relative dinucleotide abundance and codon context signatures revealed strict avoidance of TpA dinucleotide across genomes. Multivariate statistical analysis based on codon usage data resulted in discrete clustering of pathogens and non-pathogens which correlated with previous reports on their phylogenetic distribution.

Analysis of codon usage bias of chloroplast genes in Oryza species : Codon usage of chloroplast genes in Oryza species

The codon usage bias in chloroplast genes of Oryza species was low and AT rich. The pattern of codon usage was different among Oryza species and mainly influenced by mutation pressure and natural selection. Codon usage bias (CUB) is the unequal usage of synonymous codons in which some codons are more preferred to others in the coding sequences of genes. It shows a species-specific property. We studied the patterns of codon usage and the factors that influenced the CUB of protein-coding chloroplast (cp) genes in 18 Oryza species as no work was yet reported. The nucleotide composition analysis revealed that the overall GC content of cp genes in different species of Oryza was lower than 50%, i.e., Oryza cp genes were AT rich. Synonymous codon usage order (SCUO) suggested that CUB was weak in the cp genes of different Oryza species. A highly significant correlation was observed between overall nucleotides and its constituents at the third codon position suggesting that both, mutation pressure and natural selection, might influence the CUB. Correspondence analysis (COA) revealed that codon usage pattern differed across Oryza species. In the neutrality plot, a narrow range of GC3 distribution was recorded and some points were diagonally distributed in all the plots, suggesting that natural selection and mutation pressure might have influenced the CUB. The slope of the regression line was < 0.5, augmenting our inference that natural selection might have played a major role, while mutation pressure had a minor role in shaping the CUB of cp genes. The magnitudes of mutation pressure and natural selection on cp genes varied across Oryza species.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Comprehensive Analysis of Codon Usage on Porcine Astrovirus

Porcine astrovirus (PAstV), associated with mild diarrhea and neurological disease, is transmitted in pig farms worldwide. The purpose of this study is to elucidate the main factors affecting codon usage to PAstVs. Phylogenetic analysis showed that the subtype PAstV-5 sat at the bottom of phylogenetic tree, followed by PAstV-3, PAstV-1, PAstV-2, and PAstV-4, indicating that the five existing subtypes (PAstV1-PAstV5) may be formed by multiple differentiations of PAstV ancestors. A codon usage bias was found in the PAstVs-2,3,4,5 from the analyses of effective number of codons (ENC) and relative synonymous codon usage (RSCU). Nucleotides A/U are more frequently used than nucleotides C/G in the genome CDSs of the PAstVs-3,4,5. Codon usage patterns of PAstV-5 are dominated by mutation pressure and natural selection, while natural selection is the main evolutionary force that affects the codon usage pattern of PAstVs-2,3,4. The analyses of codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) showed the codon usage similarities between the PAstV and animals might contribute to the broad host range and the cross-species transmission of astrovirus. Our results provide insight into understanding the PAstV evolution and codon usage patterns.

Analysis of codon usage pattern of mitochondrial ND genes in Platyhelminthes

Codon usage bias (CUB) is the nonrandom usage of synonymous codons in which some codons are more preferred to others.CUB can be determined by mutation pressure and selection. Various approaches have been used to understand the pattern of CUB in the mitochondrial ND (MT-ND or ND) genes involved in complex I of respiratory chain in five different classes of Platyhelminthes as no work was reported yet. The present study revealed that the CUB varies across MT-ND genes and the coding sequences showed the richness of A and T. Correspondence analysis implied the effect of mutational pressure and also the pattern of codon usage was different in different classes of platyhelminthes for MT-ND genes. Highly significant correlation was observed between overall nucleotide composition and its 3rd codon position in most of the homogeneous nucleotides such as A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% and also some significant correlations observed among heterogeneous nucleotides in all the five classes for MT-ND genes suggested the role of mutational pressure as well as natural selection in affecting the CUB. Neutrality plot suggested that the contributions of natural selection and mutational pressure varied across different classes of platyhelminthes and also differed in different MT-ND genes.

Analysis of compositional properties and codon usage bias of mitochondrial CYB gene in anura, urodela and gymnophiona

We delineated the pattern of synonymous codon usage bias (CUB) and its determinants in mitochondrial CYB gene of respiratory chain across different amphibian groups namely orders anura, urodela and gymnophiona. We observed that CUB was low in CYB gene of amphibia. The gymnophionans had comparatively high bias followed by urodeles and anurans. The codons namely TCA, CCA, CAA, CGA, TGA, AAA and ACA were over-represented in all three orders. The codons such as GCC and TCC were over-represented in anura whereas in urodela, the over-represented codons were TTA, CTA, ATA, GTA, GAA, GGA and GCA. In gymnophiona, GCC, TTA, CTA, ATA, GTA, GAA and GGA codons were over-represented. The regression analysis between effective number of codons (ENC) and nucleobase at the 3rd position revealed that nucleobase A and C influenced CUB positively in order anura, while in urodela and gymnophiona, nucleobase A and T influenced the CUB positively. Mutation pressure and natural selection mutually illustrate the CUB of CYB gene (complex III gene) of amphibia as elucidated by correlation analysis between 3rd nucleotide in a codon and overall nucleotide content of the gene. However, neutrality plot showed that natural selection was the dominant evolutionary factor of CUB.

Evidence of multifaceted functions of codon usage in translation within the model beetle Tribolium castaneum

Synonymous codon use is non-random. Codons most used in highly transcribed genes, often called optimal codons, typically have high gene counts of matching tRNA genes (tRNA abundance) and promote accurate and/or efficient translation. Non-optimal codons, those least used in highly expressed genes, may also affect translation. In multicellular organisms, codon optimality may vary among tissues. At present, however, tissue specificity of codon use remains poorly understood. Here, we studied codon usage of genes highly transcribed in germ line (testis and ovary) and somatic tissues (gonadectomized males and females) of the beetle Tribolium castaneum. The results demonstrate that: (i) the majority of optimal codons were organism-wide, the same in all tissues, and had numerous matching tRNA gene copies (Opt-codon_↑tRNAs), consistent with translational selection; (ii) some optimal codons varied among tissues, suggesting tissue-specific tRNA populations; (iii) wobble tRNA were required for translation of certain optimal codons (Opt-codon_wobble), possibly allowing precise translation and/or protein folding; and (iv) remarkably, some non-optimal codons had abundant tRNA genes (Nonopt-codon_↑tRNAs), and genes using those codons were tightly linked to ribosomal and stress-response functions. Thus, Nonopt-codon_↑tRNAs codons may regulate translation of specific genes. Together, the evidence suggests that codon use and tRNA genes regulate multiple translational processes in T. castaneum.

Genome-wide codon usage pattern analysis reveals the correlation between codon usage bias and gene expression in Cuscuta australis

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis.

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

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

Selection constraints determine preference for A/U-ending codons in Taxus contorta

Unequal utilization of synonymous codons is a well-known phenomenon among living organisms. This phenomenon plays a major role in the enhancement of the accuracy and efficiency of translation. Gymnosperms are rarely paid attention in this aspect. Understanding the degree of and determining the forces influencing codon usage bias (CUB) in Taxus contorta, an endangered Himalayan gymnosperm, will prove useful in interpreting the evolutionary characteristics of this species. Using RNAseq data, 93 790 assembled transcripts were clustered into 32 701 unigenes. Around 13 061 full-length sequences were utilized for the analysis of CUB. Compositional properties showed that GC-content ranged from 28.76% to 65.22%, with an average value of 44.28%, suggesting an AT-rich genome. The mean effective number of codons (ENC) value revealed that CUB is not strong in T. contorta. The preferred codons tended to be A/U ending, whereas the avoided codons tended to be G/C ending. A P2 index of 0.54 and a Mutation Responsive Index (MRI) value of -0.02 in addition to the results revealed by the neutrality, ENC, and parity plots showed that natural selection is a predominating factor governing CUB. Mutational pressure, gene length, hydropathiciy, aromaticity, and nucleotide composition influence CUB weakly.

Insights into the nucleotide composition and codon usage pattern of human tumor suppressor genes

Tumor suppressor genes encode different proteins that inhibit the uncontrolled proliferation of cell growth and tumor development. To acquire clues for predicting gene expression level, it is essential to understand the codon usage bias (CUB) of genes to characterize genome which possesses its own compositional characteristics and unique coding sequences. We used bioinformatic tools to analyze the codon usage patterns of 637 human tumor suppressor genes as no work was reported earlier. The mean effective number of codons of these genes was 48, indicating low CUB. Our results exhibited a significant positive correlation among different nucleotide compositions and the codons ending with C base was most frequently used along with the most over-represented codon CTG and GTG codifying leucine and valine amino acid, respectively, in human tumor suppressor genes. The neutrality plot showed a significant positive correlation (Pearson, r = 0. 646; P < .01) suggesting that mutation on GC bias might affect the CUB. However, the linear regression coefficient of GC₁₂ on GC₃ in human tumor suppressor genes suggested that natural selection played a major role while mutation pressure played a minor role in the codon usage patterns of tumor suppressor genes in human. Our study would throw light into the factors that affect CUB and the codon usage patterns in the human tumor suppressor genes.

Compositional properties and codon usage pattern of mitochondrial ATP gene in different classes of Arthropoda

Codon usage bias (CUB) is defined as the usage of synonymous codons unequally for an amino acid in a gene transcript. It is influenced by both mutation pressure and natural selection and is a species-specific property. In our current study, we used bioinformatic methods to investigate the coding sequences of mitochondrial adenosine triphosphate gene (MT-ATP) in different classes of arthropoda to know the codon usage pattern of the gene as no work was described earlier. The analysis of compositional properties suggested that the gene is AT rich. The effective number of codons revealed the CUB of both ATP6 and ATP8 gene was moderate. Heat map showed that the codons ending with AT were negatively associated with GC3 while the codons ending with GC were positively associated with GC3 in all the classes of arthropoda. Correspondence study revealed that the pattern of codon usage of ATP6 and ATP8 genes differed across classes. Neutrality plot suggested the codon usage bias of these two genes in phylum arthropoda was influenced by both mutation pressure and natural selection.

The codon usage pattern of genes involved in ovarian cancer

In this study, we analyzed the compositional dynamics and codon usage pattern of genes involved in ovarian cancer (OC) using a computational method. Mutations in specific genes are associated with OC, and some genes are risk factors for progression of OC, but no work has been reported yet on the codon usage pattern of genes involved in OC. Nucleotide composition analysis of OC-related genes suggested that the overall GC content was higher than AT content; that is, the genes were GC rich. The improved effective number of codons indicated that the overall extent of codon usage bias of genes involved in OC was low. The codons AGC, CTG, ATC, ACC, GTG, and GCC were overrepresented, while the codons TCG, TTA, CTA, CCG, CAA, CGT, ATA, ACG, GTA, GTT, GCG, and GGT were underrepresented in the genes. Correspondence analysis suggested that the codon usage pattern was different in different genes. A highly significant correlation was observed between GC12 and GC3 (r = 0.587, P < 0.01) of genes, suggesting that directional mutation affected the three codon positions. Our report on the codon usage pattern of genes involved in OC includes a new perspective for elucidating the mechanisms of biased usage of synonymous codons, as well as providing useful clues for molecular genetic engineering.

Codon Usage Pattern of Genes Involved in Central Nervous System

Codon usage bias (CUB) is the non-uniform usage of synonymous codons in which some codons are more preferred to others in the transcript. Analysis of codon usage bias has applications in understanding the basics of molecular biology, genetics, gene expression, and molecular evolution. To understand the patterns of codon usage in genes involved in the central nervous system (CNS), we used bioinformatic approaches to analyze the protein-coding sequences of genes involved in the CNS. The improved effective number of codons (ENC) suggested that the overall codon usage bias was low. The relative synonymous codon usage (RSCU) revealed that the most frequently occurring codons had a G or C at the third codon position. The codons namely TCC, AGC, CTG, CAG, CGC, ATC, ACC, GTG, GCC, GGC, and CGG (average RSCU > 1.6) were over-represented. Both mutation pressure and natural selection might affect the codon usage pattern as evident from correspondence and parity plot analyses. The overall GC content (59.93) was higher than AT content, i.e., genes were GC-rich. The correlation of GC12 with GC3 suggested that mutation pressure might affect the codon usage pattern.

Dissimilar substitution rates between two strands of DNA influence codon usage pattern in some human genes

We illustrated the descriptive aspects of codon usage of some important human genes and their expression potential in E. coli. By comparing the results of various codon usage parameters, effects that are due to selection and mutational pressures have been deciphered. The variation in GC3s explains a significant proportion of the variation in codon usage patterns. The codons CGC, CGG, CTG and GCG showed strong positive correlation with GC3, which suggested that codon usage had been influenced by GC bias. We also found that ACC (Thr, RSCU-1.77), GCC (Ala, RSCU-1.67), CCC (Pro, RSCU-1.54), TCC (Ser, RSCU-1.47) were frequently used which signified that C was common at 2nd and 3rd codon positions. Correspondence analysis revealed that F1 axis had significant correlation with various GC contents suggesting that compositional properties under mutation pressure might affect codon usage bias. Nc-GC3 plot analysis suggested that both mutation pressure and natural selection might affect the codon usage bias which is also supported by neutrality plot analysis. The dinucleotide CT, TG and AG were significantly over-represented and CG, TA, AT, TT, and GT were underrepresented due to high rate of spontaneous mutation resulting from cytosine deamination.

Codon usage vis-a-vis start and stop codon context analysis of three dicot species

To understand the variation in genomic composition and its effect on codon usage, we performed the comparative analysis of codon usage and nucleotide usage in the genes of three dicots, Glycine max, Arabidopsis thaliana and Medicago truncatula. The dicot genes were found to be A/T rich and have predominantly A-ending and/or T-ending codons. GC3s directly mimic theusage pattern of global GC content. Relative synonymous codon usage analysis suggests that the high usage frequency of A/T over G/C mononucleotide containing codons in AT-rich dicot genome is due to compositional constraint as a factor of codon usage bias. Odds ratio analysis identified the dinucleotides TpG, TpC, GpA, CpA and CpT as over-represented, where, CpG and TpA as under-represented dinucleotides. The results of (NcExp-NcObs)/NcExp plot suggests that selection pressure other than mutation played a significant role in influencing the pattern of codon usage in these dicots. PR2 analysis revealed the significant role of selection pressure on codon usage. Analysis of varience on codon usage at start and stop site showed variation in codon selection in these sites. This study provides evidence that the dicot genes were subjected to compositional selection pressure.