Next Article in Journal
Genome-Wide Association Study Identifies 12 Loci Associated with Body Weight at Age 8 Weeks in Korean Native Chickens
Previous Article in Journal
Successful Genetic Screening and Creating Awareness of Familial Hypercholesterolemia and Other Heritable Dyslipidemias in the Netherlands
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 8
10.3390/genes12081169
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB

by
Xin Li
,
Xiaocen Wang
,
Pengtao Gong
,
Nan Zhang
,
Xichen Zhang
and
Jianhua Li
*
Key Laboratory of Zoonosis Research, Ministry of Education, College of Veterinary Medicine, Jilin University, Changchun 130062, China
*
Author to whom correspondence should be addressed.
Genes 2021, 12(8), 1169; https://doi.org/10.3390/genes12081169
Submission received: 29 June 2021 / Revised: 24 July 2021 / Accepted: 27 July 2021 / Published: 29 July 2021
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite-induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duodenalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G. duodenalis CUB formation. This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB. These results will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

1. Introduction

Codon usage bias is a widespread feature among both prokaryotes and eukaryotes, which describes synonymous codons that are not used at the same frequency in the process of gene translation [1,2,3]. CUB is widely distributed and affected by tRNA abundance [4], nucleotide composition [5], translational processes [6], gene function [7,8], protein structure [9], hydrophobicity [10], environment temperature [11], etc. Among the known factors affecting CUB in different organisms, constraints on composition and translation selection are believed to be the primary causes for the differences in CUB among genes in different species [12]. CUB can be used to optimize the expression of foreign genes in given host cells. In addition, CUB could also provide clues about the evolution and environmental adaptation of various species [13].
Analysis of synonymous codon usage bias has various important applications, such as degenerate primer design [14], heterologous gene expression [15], species origin determination [16], the prediction of the expression level of genes [17,18], and the prediction of gene functions [19]. CUB has been studied in various species, such as Taenia saginata [20], Taenia multiceps [21], Plasmodium falciparum [22], Entamoeba histolytica [23], Microsporidia [24], Caenorhabditis, Drosophila, Arabidopsis [25], Streptomyces [26], Borrelia burgdorferi [27], and Saccharomyces cerevisiae [28,29]. Giardiasis is a worldwide epidemic water source diarrhea disease, which can infect various mammals, including human beings [30,31,32]. Previous studies of G. duodenalis codon usage only analyzed with small number of genes [33,34].
This study investigates the CUB of G. duodenalis from transcriptome data, which provides useful information for elucidation of the mechanism of synonymous CUB, design of gene vaccine, and better control strategy of G. duodenalis.

2. Methods

2.1. Transcriptome Data

Nine thousand, seven hundred and forty-seven coding sequences (CDS) of G. duodenalis (Giardia Assemblage A isolate WB) obtained from GiardiaDB were investigated based on the transcriptome profiling of Oscar Franzén et al. [35]. To minimize outliers caused by small size, CDS with sizes below 300 bp were eliminated. Thus, a total of 5968 CDS were selected for the analysis.

2.2. Indices of Codon Usage

A codon usage table of G. duodenalis was obtained from Codon Usage Database via Optimizer software [36]. CDS sequences codon usage of G. duodenalis genes was analyzed by Codon W 1.4.4 software and EMBOSS online tools (https://www.bioinformatics.nl/emboss-explorer/, accessed on 26 July 2021). 17 CDS sequences codon usage of G. duodenalis highly expressed variant-specific surface proteins (VSPs) [37] were analyzed by EMBOSS online tools (https://www.bioinformatics.nl/emboss-explorer/, accessed on 26 July 2021). The analysis parameters included an effective number of codons (ENC), the codon adaptation index (CAI), relative synonymous codon usage (RSCU), GC content (guanine and cytosine content), A3s, T3s, G3s, C3s, and GC3s. Nucleotide composition and its frequency on the third synonymous codon can be used to reflect the codon usage preference of CDS sequences. The RSCU value is the observed codon frequency divided by the frequency expected under the same assumption for amino acid synonymous codons, and it is an index for studying differences in synonymous codon usage among genes. The RSCU values were measured according to the previous study [12]. ENC values reflect the number of codon types used in a gene, and its value generally ranged from 20 (when only one codon was used) to 61 (when all codons are used equally), and the ENC values were measured as previously described [38]. CAI refers to the consistency between the usage frequency of synonymous codons and optimal codons in the coding region, with a value of 0–1. CAI is used to estimate the degree of CUB, which is preferred in highly expressed genes. Higher CAI values mean that CUB may be stronger, and the potential expression level may be higher, and the CAI values were calculated as previously described [39].

2.3. Neutrality Plot

The effect of selection on CUB was usually measured by neutrality plot analysis. In this method, the average GC content of GC12 and GC3s was calculated. The scatter plot was drawn with GC3s as an independent variable and GC12 as a dependent variable. The points representing genes are distributed on or near the diagonal, indicating that the codon usage pattern is greatly affected by mutation; on the contrary, the smaller the slope of the scatter formation curve is, or even parallels to the horizontal axis, indicating that codon usage pattern is greatly influenced by environmental selection [40].

2.4. ENC Plot

The ENC plot of ENC values plotted against GC3s values is used to analyze the influence of base composition on the codon usage in a genome [41]. A standard curve is used to show the functional relationship between ENC and GC3s values under mutation pressure rather than selection pressure. In the ENC plot correlation analysis, GC3s were used as the independent variable and ENC as the dependent variable to construct the scatter diagram and to analyze the correlation between ENC and GC3s. In addition, according to the CUB, the standard curve was constructed under the condition of mutation pressure, but not selection pressure. If the predicted ENC value is on or near the standard curve, it represents that the CUB is mainly influenced by mutation rather than selection pressure; if the predicted ENC value is far below the standard curve, it represents that the codon composition is mainly affected by selection pressure [41].

2.5. PR2 Bias Plot Analysis

Parity rule 2 analysis, also known as parity rule analysis, is a method to study the base composition of codons. If the gene is not under the pressure of mutation or environmental selection, the internal composition of the base is A = T, C = G. However, due to the influence of gene mutation and environmental selection pressure, the usage of G and C in the genome coding sequence is often uneven, especially the third codon deviates from the rule of intrachain equivalence. In this method, amino acids encoded by four synonymous codons were analyzed, and the calculated results of G3/ (G3 + C3) and A3/ (A3 + T3) were plotted. The coordinates (0.5, 0.5) represent the PR2 principle (A = T, C = G). The distance and position of the scattered points from the center indicate the degree and direction of the gene deviation from the rules [42].

2.6. Determination of Optimal Codons

Based on the CAI values, 5% of the total genes with extremely high and low CAI values were regarded as high and low datasets, respectively. Codon usage was compared using a Chi-squared contingency test of the two groups, and codons whose frequency of usage was significantly higher (p < 0.01) in highly-expressed genes than those with low levels of expression were defined as optimal codons [43].

2.7. Correspondence Analysis (COA)

The connection between variables and samples was widely analyzed by Multivariate statistical analysis. COA has been widely used to study codon usage variation. CondonW was used to analyze RSCU values by COA analysis, COA was performed on RSCU values using to compared the intra-genomic variation of 59 informative codons partitioned along 59 orthogonal axes (excluding Met, Trp, and stop codons) [44].

2.8. Statistical Analysis

The indices of codon usage were analyzed by CondonW1.4.4 software. Microsoft Excel and SPSS 19.0 were used to analyze the correlation based on Spearman’s rank correlation.

3. Results

3.1. Nucleotide Contents of G. duodenalis Genes

The nucleotide contents of G. duodenalis CDSs (expressed as % GC) were shown in Figure 1. The results suggested that the GC content of 5968 G. duodenalis genes exhibited a distinctly unimodal distribution. The GC contents of G. duodenalis genes varied from 39.1% to 81.1% (SD = 0.0523), and the GC contents of the 5968 CDSs were mainly ranged from 45% to 55%. To understand the distribution of nucleotides, we studied the content of GC and GC3s. GC1, GC2 and GC3 were 53.38%, 41.79% and 51.67%, respectively. This result suggested that the GC2 was different from GC1, and GC3—GC2 was the lowest among the three codon positions. The mean value of GC contents of all codons was 48.95%.
To study the relationship between the three codon positions, the neutrality graph (GC12 versus GC3) of the G. duodenalis gene was constructed. Our result demonstrated that the GC3s of G. duodenalis genes was widely distributed (26.17% to 99.32%), and there was a significant correlation between GC12 and GC3 (r = 0.283, p < 0.0001) (Figure 2), indicating that mutation might affect the codon preference in G. duodenalis genome.

3.2. Codon Usage in G. duodenalis

The synonymous codon usage patterns of G. duodenalis were shown in Table 1. The G + C content of the G. duodenalis genome was 49.1%, which indicated that the G. duodenalis genome was a little AT-rich. The total codon usage was biased towards G-and C-terminal codons (31 codons were common codons, and 16/31 frequently used codons ended with G or C). These results indicated that gene composition restriction played a significant role in the formation of codon usage variation in the G. duodenalis genome. In addition, we analyzed 17 CDS sequences codon usage of G. duodenalis highly expressed VSPs genes based on the actual protein levels [37], and 12/15 most frequently used codons were biased towards G-and C-terminal codons (Table 2). At last, a codon usage table of G. duodenalis was obtained from Codon Usage Database via Optimizer software, and 9/10 most frequently used codons were biased towards G-and C-terminal codons (Table 3). These results further confirm the accuracy of our predicted codon usage bias.

3.3. Relation between ENC and GC3

ENC values were analyzed according to the fraction of GC3s (Figure 3) to clarify the relationship between nucleotide composition and G. duodenalis CUB. The ENC values ranged from 24 to 61, suggesting significant differences in codon bias among these genes. As shown in Figure 3, most of the points were clustered near the expected ENC curve, which indicated that the ENC of most genes was close to the expected ENC value based on their GC3. In addition, the ENC of some points was lower than the expected curve, indicating that the codon usage was also affected by other factors beyond mutation pressure.
(ENCexp-ENCobs)/ENCexp was calculated to estimate the observed and expected ENC values more accurately. We found that (ENCexp-ENCobs)/ENCexp was mainly located in 0–0.05, and the values of most genes were located in −0.05–0.15 (Figure 4). The results suggested that the ENCs of most genes were slightly different from the expected ENCs values in GC3s. The observed ENCs values of most genes were near expected ENCs in GC3s, although the ENCs values of some genes were much lower.

3.4. Correspondence Analysis

The differences in synonymous codon usage among G. duodenalis genes were further studied by RSCU correspondence analysis. We determined four major contributors as axis 1–4. Axis 1 and axis 2 accounted for 14.81% and 5.30% of the total variance, respectively, while axis 3 and axis 4 accounted for 4.05% and 3.46% of the total variance, respectively, suggesting that axis 1 and axis 2 were the main contributors to G. duodenalis CUB. Axis 1 and axis 2 showed each gene in Figure 5A. To clarify the effect of gene GC contents on CUB, the gene GC contents were color coded. The genes with GC content more than or equal to 60% were shown in green, while those with GC content less than 45% were shown in red. GC content ranged from 45% to 60% were shown in blue. The results showed that the high and low GC contents of G. duodenalis genes could be separated by the primary axis. In addition, as shown in Figure 5B, we found that the terminal codons of different bases could be separated along two axes. It seemed that the separation of the first axis codon was mainly due to the frequency difference of G/C and A/T terminal codons. Further calculation showed a significant correlation between the GC content of each gene and its position on the first axis (r = 0.1154, p < 0.0001). The location of axis 1 gene was positively correlated with GC3s (r = 0.1160, p < 0.0001) and negatively correlated with ENC (r = −0.0626, p < 0.0001). These results indicated that the genes with higher GC and GC3 contents and lower ENC content on the left side of the first axis showed stronger codon bias, which demonstrated that the main factor affecting the CUB of G. duodenalis was the nucleotide composition.
To investigate different gene codon usages, we selected hydrophobic genes, aromatic genes, ribosomal genes, and other genes from 5968 genes. Figure 5C showed the distribution of these four categories of genes. Multivariate analysis of variance indicated that there were statistically significant differences in codon usage among these different genes (p < 0.01).

3.5. PR2-Bias Plot Analysis

To study whether high bias genes restrict the selection of biased codons, the relationship between A/G purines and C/T pyrimidines in amino acid was analyzed by PR2 bias plot. Figure 6 showed that genes in the upper left quadrant had low expression and the genes in the lower right had high expression. We demonstrated that G. duodenalis preferred to use G and T rather than use C and A (Figure 6), which suggested that mutation bias, selection, and other factors were involved in G. duodenalis codon usage bias.

3.6. Role of Gene Expression Level and Encoded Protein Size Synonymous CUB

To investigate the relationship between CUB and gene expression level, the correlation coefficient between CAI and ENC was calculated and analyzed. The CAI values were used to evaluate G. duodenalis genes expression level, which ranged from 0.157 to 0.912 (mean = 0.318, SD = 0.07942). As shown in Figure 7, our results indicated that genes expression level were significantly negatively correlated with genes position along axis 1 (r = 0.6811, −0.1086, p < 0.0001), while CAI value showed positive relationship significantly with GC3s and GC content (r = 0.9022, 0.6486, respectively, p < 0.0001). These results demonstrated that highly expressed genes had a strong CUB and preferred to choose G or C codons in synonymous positions.
As shown in Figure 8, correlation analysis between protein size and axis 1 values indicated that the 3 correlation coefficients (r = −0.0817, 0.1822, −0.1254, respectively, p < 0.01) all significantly correlated, suggesting that genes with higher expression levels had a smaller size.

3.7. Optimal Codons

Translational optimal codons of G. duodenalis were represented by the average value of RSCU in high/low expressed genes (Table 4). Chi-square test showed that 26 codons were optimal codons, and the frequency of these codons was significantly higher in highly expressed genes (p < 0.01), which all end in G or C, indicating that G. duodenalis preferred to use G or C ending synonymous codons.

4. Discussion

CUB widely exists among both prokaryotes and eukaryotes. This is an interesting and complex phenomenon in the process of biological evolution. Previous studies have proposed hypotheses trying to explain the origin of CUB, among which the selection-mutation-drift balance model and the neutral theory are the most influential ones [12,45], which considers that CUB is determined by the balance between mutation pressure, genetic drift, and weak selection [12,46], while neutral theory believes that mutations of degenerate coding cites should be neutral selection, which leads to random synonymous codon usage selection [45]. However, studies have also suggested that many others factors could affect CUB, such as GC-content [47,48], gene size [25], gene expression level [25,49,50], and gene recombination rate [47,49,51,52]. Furthermore, RNA and protein structure [41,53,54,55], intron length [56], evolutionary age of the genes [57], population size [58], the aromaticity and the hydrophobicity of the coding proteins [10,59] have all been found to be influencing factors.
A previous analysis of G. duodenalis codon usage was restricted, which only considered eight genes, and yet it seems that the codon usage of G. duodenalis has been quite heterogeneous [33]. Another analysis of G. duodenalis codon usage investigated 65 genes, and 21 codons were the optimal codons, which were all end in C or G and almost exclusively used in the highly expressed genes [34], which was similar to our conclusion. However, the CUB of G. duodenalis has not been fully studied yet. In the present study, we performed a more comprehensive analysis based on the whole transcriptome data and found that multiple factors influence shaping G. duodenalis CUB, such as mutation pressure, selection, gene expression, and compositional constraints, and protein size.
Generally speaking, nucleotide composition is one of the most important effect factors in forming codon usage, while GC content reflects the overall trend of codon mutation [60]. Some GC-rich organisms, such as Triticum Aestivum, Oryzasativa, Bacteria, Archea, and Fungi [61,62], have been proved that they tend to use G or C ending codons. In contrast, some AT-rich organisms, such as Onchocerca volvulus, Mycoplasma capricolum, and P. falciparum [63,64,65], have been shown that they tend to A or T ending codons. In this study, we demonstrated that the average GC content among the 5968 G. duodenalis genes was 49.1%. Although the genome of G. duodenalis seems to be slightly AT-rich, the usage of all codons is biased towards G or C-terminal codons (Table 2), which is similar to those in T. saginata [20] and T. multiceps [21].
Previous studies have demonstrated that over-expressed genes are expressed more frequently, and produced more protein than other genes [22]. Moreover, preferred codons were more frequently used in highly expressed genes than other genes [66]. ENC represents the species independent synonymous bias in genes [38,67]. The genes expression level of protein coding genes can be classified according to ENC values. Highly expressed genes show less ENC value, while the lowly expressed genes show more ENC value. In the present study, the ENC values of G. duodenalis genes ranged from 24–61, suggesting that the codon bias among these genes was significantly different.
CAI is a method of recognizing differential gene expression by codon selection. Highly expressed genes exhibit a high tendency to use certain codons and tend to use them frequently [68,69]. In the present study, the low CAI of G. duodenalis indicated that highly expressed genes faced greater translation selection pressure in shaping CUB, which is similar to the codon usage in T. saginata [20] and P. falciparum [70].
It has been confirmed that various organisms, such as T. multiceps [21], T. saginata [20], S. cerevisiae [29], Silenelatifolia [71], Caenorhabditis elegans [25], and Arabidopsis thaliana [25], showed significant negative correlations between gene size and CUB. Interestingly, our research suggested that G. duodenalis also showed a negative correlation between gene size and CUB, suggesting that selection restriction made G. duodenalis tend to produce smaller proteins with similar functions to larger proteins, thus reducing the energy consumption of producing specific functional proteins [72].
Identifying the optimal codons could provide an effective means for rational codon usage rearrangement and evolutionary molecular genetics research [73,74,75]. The optimal codons tend to reflect the GC and AT content of the genomes [62,76]. In this study, 26 codons were identified as optimal codons, ending with G or C. This phenomenon is similar to previous study of G. duodenalis [34] and other eukaryotic genomes, such as T. multiceps [21] and T. saginata [20]. Identifying optimal codons in the G. duodenalis genome might provide valuable information for genetic engineering and evolutionary study.
Our study revealed the pattern of CUB in the G. duodenalis genome and its influencing factors. Our results indicated that the CUB of G. duodenalis seemed to be a complex equilibrium under different pressures: Natural selection, mutation, GC content, gene expression level, and protein size. Interestingly, all 26 optimal codons ended with G or C, which would be useful for cloning and expression of foreign genes in G. duodenalis. Together, our study elucidated the codon usage pattern of G. duodenalis and provided useful information for genetic engineering and evolutionary studies in this primitive eukaryote.

5. Conclusions

This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB, which will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

Author Contributions

X.L., X.W. and J.L. drafted the main manuscript and performed the data analysis; X.L., N.Z. and J.L. were responsible for experimental design; P.G., X.Z. and J.L. were responsible for guiding and manuscript revisions. All authors have read and agreed to the published version of the manuscript.

Funding

This Research was funded by China Postdoctoral Science Foundation (No.2019M651213) and National Natural Science Foundation of China (No.31772732 and 31672288).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Excluded this statement.

Acknowledgments

We thank Tingzhang Wang for analysis guidance.

Conflicts of Interest

The authors declare that they have no competing financial interests.

Abbreviations

GC1GC-content at the first codon positions
GC2GC-content at the second codon positions
GC12The average of GC1 and GC2
GC3GC-content at the third codon positions
GC3sFrequency of either a G or C at the third codon position of synonymous codons

References

  1. Akashi, H.; Eyrewalker, A. Translational selection and molecular evolution. Curr. Opin. Genet. Dev. 1998, 8, 688–693. [Google Scholar] [CrossRef]
  2. Akashi, H. Gene expression and molecular evolution. Curr. Opin. Genet. Dev. 2001, 11, 660–666. [Google Scholar] [CrossRef]
  3. Duret, L. Evolution of synonymous codon usage in metazoans. Curr. Opin. Genet. Dev. 2002, 12, 640–649. [Google Scholar] [CrossRef]
  4. Ikemura, T. Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes: A proposal for a synonymous codon choice that is optimal for the E. coli translational system. J. Mol. Biol. 1981, 146, 1–21. [Google Scholar] [CrossRef]
  5. Osawa, S.; Ohama, T.; Yamao, F.; Muto, A.; Jukes, T.H.; Ozeki, H.; Umesono, K. Directional Mutation Pressure and Transfer RNA in Choice of the Third Nucleotide of Synonymous Two-Codon Sets. Proc. Natl. Acad. Sci. USA 1988, 85, 1124–1128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sharp, P.M.; Li, W.H. Codon usage in regulatory genes in Escherichia coli does not reflect selection. Nucleic Acids Res. 1986, 14, 7737–7749. [Google Scholar] [CrossRef] [Green Version]
  7. Chiapello, H.; Al, E. Codon usage and gene function are related in sequences of Arabidopsis thaliana. Gene 1998, 209, GC1–GC38. [Google Scholar] [CrossRef]
  8. Moriyama, E.N.; Powell, J.R. Gene length and codon usage bias in Drosophila melanogaster, Saccharomyces cerevisiae and Escherichia coli. Nucleic Acids Res. 1998, 26, 3188–3193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Oresic, M.; Shalloway, D. Specific correlations between relative synonymous codon usage and protein secondary structure. J. Mol. Biol. 1998, 281, 31–48. [Google Scholar] [CrossRef] [PubMed]
  10. Romero, H.; Zavala, A.; Musto, H. Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces. Nucleic Acids Res. 2000, 28, 2084–2090. [Google Scholar] [CrossRef] [Green Version]
  11. Sau, K.; Deb, A. Temperature influences synonymous codon and amino acid usage biases in the phages infecting extremely thermophilic prokaryotes. Silico Biol. 2009, 9, 1–9. [Google Scholar] [CrossRef]
  12. Sharp, P.M.; Li, W.H. An evolutionary perspective on synonymous codon usage in unicellular organisms. J. Mol. Evol. 1986, 24, 28–38. [Google Scholar] [CrossRef]
  13. Angellotti, M.C.; Bhuiyan, S.B.; Chen, G.; Wan, X.F. CodonO: Codon usage bias analysis within and across genomes. Nucleic Acids Res. 2007, 35, 132–136. [Google Scholar] [CrossRef]
  14. Zheng, Y.; Zhao, W.M.; Wang, H.; Zhou, Y.B.; Luan, Y.; Qi, M.; Cheng, Y.Z.; Tang, W.; Liu, J.; Yu, H. Codon usage bias in Chlamydia trachomatis and the effect of codon modification in the MOMP gene on immune responses to vaccination. Biochem. Cell Biol.-Biochim. Biol. Cell 2007, 85, 218–226. [Google Scholar] [CrossRef] [PubMed]
  15. Kane, J.F. Effects of rare codon clusters on high-level expression of heterologous proteins in Escherichia coli. Curr. Opin. Biotechnol. 1995, 6, 494–500. [Google Scholar] [CrossRef]
  16. Ahn, I.; Jeong, B.J.; Bae, S.E.; Jin, J.; Son, H.S. Genomic Analysis of Influenza A Viruses, including Avian Flu (H5N1) Strains. Eur. J. Epidemiol. 2006, 21, 511–519. [Google Scholar] [CrossRef] [PubMed]
  17. Naya, H.; Romero, H.; Carels, N.; Zavala, A.; Musto, H. Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii. FEBS. Lett. 2001, 501, 127–130. [Google Scholar] [CrossRef] [Green Version]
  18. Gupta, S.K.; Bhattacharyya, T.K.; Ghosh, T.C. Synonymous Codon Usage in Lactococcus lactis: Mutational Bias Versus Translational Selection. J. Biomol. Struct. Dyn. 2004, 21, 527–535. [Google Scholar] [CrossRef]
  19. Lin, K.; Kuang, Y.; Joseph, J.S.; Kolatkar, P.R. Conserved codon composition of ribosomal protein coding genes in Escherichia coli, Mycobacterium tuberculosis and Saccharomyces cerevisiae: Lessons from supervised machine learning in functional genomics. Nucleic Acids Res. 2002, 30, 2599–2607. [Google Scholar] [CrossRef]
  20. Yang, X.; Luo, X.; Cai, X. Analysis of codon usage pattern in Taenia saginata based on a transcriptome dataset. Parasites Vectors 2014, 7, 527. [Google Scholar] [CrossRef]
  21. Huang, X.; Xu, J.; Chen, L.; Wang, Y.; Gu, X.; Peng, X.; Yang, G. Analysis of transcriptome data reveals multifactor constraint on codon usage in Taenia multiceps. BMC Genom. 2017, 18, 308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Peixoto, L.; Fernández, V.; Musto, H. The effect of expression levels on codon usage in Plasmodium falciparum. Parasitology 2004, 128, 245–251. [Google Scholar] [CrossRef] [PubMed]
  23. Ghosh, T.C.; Gupta, S.K.; Majumdar, S. Studies on codon usage in Entamoeba histolytica. Int. J. Parasitol. 2000, 30, 715–722. [Google Scholar] [CrossRef]
  24. Xiang, H.; Zhang, R.; Iii, R.R.B.; Liu, T.; Li, Z.; Pombert, J.F.; Zhou, Z. Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes. PLoS ONE 2015, 10, e0129223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Duret, L.; Mouchiroud, D. Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis. Proc. Natl. Acad. Sci. USA 1999, 96, 4482–4487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wright, F.; Bibb, M.J. Codon usage in the G+C-rich Streptomyces genome. Gene 1992, 113, 55–65. [Google Scholar] [CrossRef]
  27. Mcinerney, J.O. Replicational and Transcriptional Selection on Codon Usage in Borrelia burgdorferi. Proc. Natl. Acad. Sci. USA 1998, 95, 10698–10703. [Google Scholar] [CrossRef] [Green Version]
  28. Sharp, P.M.; Cowe, E. Synonymous codon usage in Saccharomyces cerevisiae. Yeast 1991, 7, 657–678. [Google Scholar] [CrossRef]
  29. Kliman, R.M.; Irving, N.; Santiago, M. Selection conflicts, gene expression, and codon usage trends in yeast. J. Mol. Evol. 2003, 57, 98–109. [Google Scholar] [CrossRef] [PubMed]
  30. Savioli, L.; Smith, H.; Thompson, A. Giardia and Cryptosporidium join the ‘Neglected Diseases Initiative. Trends Parasitol. 2006, 22, 203–208. [Google Scholar] [CrossRef]
  31. Geurden, T.; Vercruysse, J.; Claerebout, E. Is Giardia a significant pathogen in production animals? Exp. Parasitol. 2010, 124, 98–106. [Google Scholar] [CrossRef] [PubMed]
  32. Adam, R.D. Biology of Giardia lamblia. Clin. Microbiol. Rev. 2001, 14, 447–475. [Google Scholar] [CrossRef] [Green Version]
  33. Char, S.; Farthing, M.J. Codon usage in Giardia lamblia. J. Protozool. 1992, 39, 642–644. [Google Scholar] [CrossRef] [PubMed]
  34. Lafay, B.; Sharp, P.M. Synonymous codon usage variation among Giardia lamblia genes and isolates. Mol. Biol. Evol. 1999, 16, 1484–1495. [Google Scholar] [CrossRef] [Green Version]
  35. Franzén, O.; Jerlström-Hultqvist, J.; Einarsson, E.; Ankarklev, J.; Ferella, M.; Andersson, B.; Svärd, S.G. Transcriptome profiling of Giardia intestinalis using strand-specific RNA-seq. PLoS Comput. Biol. 2013, 9, e1003000. [Google Scholar] [CrossRef] [Green Version]
  36. Pere, P.; Guzmán, E.; Romeu, A.; Garcia-Vallvé, S. OPTIMIZER: A web server for optimizing the codon usage of DNA sequences. Nucleic Acids Res. 2007, 35, W126–W131. [Google Scholar]
  37. Müller, J.; Braga, S.; Uldry, A.C.; Heller, M.; Müller, N. Comparative proteomics of three Giardia lamblia strains: Investigation of antigenic variation in the post-genomic era. Parasitology 2020, 147, 1008–1018. [Google Scholar] [CrossRef]
  38. Wright, F. The ‘effective number of codons’ used in a gene. Gene 1990, 87, 23–29. [Google Scholar] [CrossRef]
  39. Sharp, P.M.; Li, W.H. The codon Adaptation Index—A measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15, 1281–1295. [Google Scholar] [CrossRef] [Green Version]
  40. Sueoka, N. Directional mutation pressure and neutral molecular evolution. Proc. Natl. Acad. Sci. USA 1988, 85, 2653–2657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Hartl, D.L.; Moriyama, E.N.; Sawyer, S.A. Selection intensity for codon bias. Genetics 1994, 138, 227–234. [Google Scholar] [CrossRef]
  42. Sueoka, N. Near Homogeneity of PR2-Bias Fingerprints in the Human Genome and Their Implications in Phylogenetic Analyses. J. Mol. Evol. 2001, 53, 469–476. [Google Scholar] [CrossRef]
  43. Liu, Q. Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans. Biosystems 2006, 85, 99–106. [Google Scholar] [CrossRef]
  44. Greenacre, M.J. Theory and applications of correspondence analysis. J. Am. Stat. Assoc. 1984, 80, 1067. [Google Scholar]
  45. Nakamura, Y.; Gojobori, T.T. Codon usage tabulated from the international DNA sequence databases. Nucleic Acids Res. 1998, 26, 334. [Google Scholar] [CrossRef] [Green Version]
  46. Bulmer, M. Are codon usage patterns in unicellular organisms determined by selection-mutation balance? J. Evol. Biol. 1988, 1, 15–26. [Google Scholar] [CrossRef]
  47. Comeron, J.M.; Kreitman, M.; Aguadé, M. Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila. Genetics 1999, 151, 239–249. [Google Scholar] [CrossRef]
  48. Marais, G.; Mouchiroud, D.; Duret, L. Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes. Proc. Natl. Acad. Sci. USA. 2001, 98, 5688–5692. [Google Scholar] [CrossRef] [Green Version]
  49. Hey, J.; Kliman, R.M. Interactions between natural selection, recombination and gene density in the genes of Drosophila. Genetics 2002, 160, 595–608. [Google Scholar] [CrossRef]
  50. Stenico, M.; Lloyd, A.T.; Sharp, P.M. Codon usage in Caenorhabditis elegans: Delineation of translational selection and mutational biases. Nucleic Acids Res. 1994, 22, 2437–2446. [Google Scholar] [CrossRef] [Green Version]
  51. Kliman, R.M.; Hey, J. Hill-Robertson interference in Drosophila melanogaster: Reply to Marais, Mouchiroud and Duret. Genet. Res. 2003, 81, 89–90. [Google Scholar] [CrossRef] [PubMed]
  52. Marais, G. Hill-Robertson Interference is a Minor Determinant of Variations in Codon Bias Across Drosophila melanogaster and Caenorhabditis elegans Genomes. Mol. Biol. Evol. 2002, 19, 1399–1406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Chen, Y.; Carlini, D.B.; Baines, J.F.; Parsch, J.; Braverman, J.M.; Tanda, S.; Stephan, W. RNA secondary structure and compensatory evolution. Genes Genet. Syst. 1999, 74, 271–286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Carlini, D.B.; Chen, Y.; Stephan, W. The relationship between third-codon position nucleotide content, codon bias, mRNA secondary structure and gene expression in the drosophilid alcohol dehydrogenase genes Adh and Adhr. Genetics 2001, 159, 623–633. [Google Scholar] [CrossRef]
  55. Orešič, M.; Dehn, M.H.H.; Korenblum, D.H.H.; Shalloway, D.H.H. Tracing Specific Synonymous Codon–Secondary Structure Correlations Through Evolution. J. Mol. Evol. 2003, 56, 473–484. [Google Scholar] [CrossRef]
  56. Vinogradov, A.E. Intron length and codon usage. J. Mol. Evol. 2001, 52, 310. [Google Scholar] [CrossRef] [Green Version]
  57. Prat, Y.; Fromer, M.; Linial, N.; Linial, M. Codon usage is associated with the evolutionary age of genes in metazoan genomes. Bmc Evol. Biol. 2009, 9, 285. [Google Scholar] [CrossRef] [Green Version]
  58. Berg, O.G. Selection intensity for codon bias and the effective population size of Escherichia coli. Genetics 1996, 142, 1379–1382. [Google Scholar] [CrossRef]
  59. Rispe, C.; Delmotte, F.; van Ham, R.C.; Moya, A. Mutational and selective pressures on codon and amino acid usage in Buchnera, endosymbiotic bacteria of aphids. Genome Res. 2004, 14, 44–53. [Google Scholar] [CrossRef] [Green Version]
  60. Shang, M.Z.; Liu, F.; Hua, J.P.; Wang, K.B. Analysis on Codon Usage of Chloroplast Genome of Gossypium hirsutum. Scientia Agricultura Sinica 2011, 44, 245–253. [Google Scholar]
  61. Kawabe, A.; Miyashita, N.T. Patterns of codon usage bias in three dicot and four monocot plant species. Genes Genet. Syst. 2003, 78, 343–352. [Google Scholar] [CrossRef] [Green Version]
  62. Hershberg, R.; Petrov, D. A: General Rules for Optimal Codon Choice. PLoS Genet. 2009, 5, e1000556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Saul, A.; Battistutta, D. Codon usage in Plasmodium falciparum. Mol. Biochem. Parasitol. 1988, 27, 35–42. [Google Scholar] [CrossRef]
  64. Milhon, J.L.; Tracy, J.W. Updated Codon Usage in Schistosoma. Exp. Parasitol. 1995, 80, 353–356. [Google Scholar] [CrossRef] [PubMed]
  65. Muto, A.; Yamao, F.; Osawa, S. The genome of Mycoplasma capricolum. Prog. Nucleic Acid Res. Mol. Biol. 1987, 34, 29–58. [Google Scholar]
  66. Quax, T.E.F.; Claassens, N.J.; Söll, D.; Oost, J.V.D. Codon Bias as a Means to Fine-Tune Gene Expression. Mol. Cell 2015, 59, 149–161. [Google Scholar] [CrossRef] [Green Version]
  67. Fuglsang, A. The ‘effective number of codons’ revisited. Biochem. Biophys. Res. Commun. 2004, 317, 957–964. [Google Scholar] [CrossRef]
  68. Supek, F. The Code of Silence: Widespread Associations Between Synonymous Codon Biases and Gene Function. J. Mol. Evol. 2016, 82, 65–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Carbone, A.; Zinovyev, A.; Képès, F. Codon adaptation index as a measure of dominating codon bias. Bioinformatics 2003, 19, 2005–2015. [Google Scholar] [CrossRef] [Green Version]
  70. Gajbhiye, S.; Patra, P.K.; Yadav, M.K. New insights into the factors affecting synonymous codon usage in human infecting Plasmodium species. Acta Trop. 2017, 176, 29–33. [Google Scholar] [CrossRef]
  71. Qiu, S.; Bergero, R.; Zeng, K.; Charlesworth, D. Patterns of Codon Usage Bias in Silene latifolia. Mol. Biol. Evol. 2011, 28, 771–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Moriyama, E.N.; Powell, J.R. Codon Usage Bias and tRNA Abundance in Drosophila. J. Mol. Evol. 1997, 45, 514–523. [Google Scholar] [CrossRef] [PubMed]
  73. Ko, H.J.; Ko, S.Y.; Kim, Y.J.; Lee, E.G.; Cho, S.N.; Kang, C.Y. Optimization of codon usage enhances the immunogenicity of a DNA vaccine encoding mycobacterial antigen Ag85B. Infect. Immun. 2005, 73, 5666–5674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Peng, R.; Yao, Q.; Xiong, A.; Cheng, Z.; Li, Y. Codon-modifications and an endoplasmic reticulum-targeting sequence additively enhance expression of an Aspergillus phytase gene in transgenic canola. Plant Cell Rep. 2006, 25, 124–132. [Google Scholar] [CrossRef] [PubMed]
  75. Rouwendal, G.J.A.; Mendes, O.; Wolbert, E.J.H.; Boer, A.D.D. Enhanced expression in tobacco of the gene encoding green fluorescent protein by modification of its codon usage. Plant Mol. Biol. 1997, 33, 989–999. [Google Scholar] [CrossRef]
  76. Rao, Y.; Wu, G.; Wang, Z.; Chai, X.; Nie, Q.; Zhang, X. Mutation Bias is the Driving Force of Codon Usage in the Gallus gallusgenome. DNA Res. Int. J. Rapid Publ. Rep. Genes Genomes 2011, 18, 499–512. [Google Scholar]
Genes 12 01169 g001 550
Figure 1. Distribution of GC contents in G. duodenalis genes.
Figure 1. Distribution of GC contents in G. duodenalis genes.
Genes 12 01169 g001
Genes 12 01169 g002 550
Figure 2. Neutrality plots analysis of GC12 and GC3 of G. duodenalis transcriptome.
Figure 2. Neutrality plots analysis of GC12 and GC3 of G. duodenalis transcriptome.
Genes 12 01169 g002
Genes 12 01169 g003 550
Figure 3. Distribution of ENC and GC3s of G. duodenalis genes. ENC is plotted against GC3s. The solid line (green) shows the expected ENC value if the CUB is caused by GC3s only.
Figure 3. Distribution of ENC and GC3s of G. duodenalis genes. ENC is plotted against GC3s. The solid line (green) shows the expected ENC value if the CUB is caused by GC3s only.
Genes 12 01169 g003
Genes 12 01169 g004 550
Figure 4. Frequency distribution of the ENC ratio.
Figure 4. Frequency distribution of the ENC ratio.
Genes 12 01169 g004
Genes 12 01169 g005 550
Figure 5. Correspondence analysis of RSCU for G. duodenalis genes. (A). The first two axes show the distribution of G. duodenalis genes. Green dots represent G + C content ≥ 60%; Blue dots represent G + C content ≥ 45%, but less than 60%; Red dots represent G + C content ≤ 45%. (B). Panel A shows the distribution of codons on the same two axes. Red dots indicate A and T ending codons, and green dots indicate G and C ending codons. (C). Red dots represent ribosomal genes, yellow dots represent genes with an Aromo value ≥0.15, green dots represent genes with a Gravy value higher than 5, and blue dots represent other genes.
Figure 5. Correspondence analysis of RSCU for G. duodenalis genes. (A). The first two axes show the distribution of G. duodenalis genes. Green dots represent G + C content ≥ 60%; Blue dots represent G + C content ≥ 45%, but less than 60%; Red dots represent G + C content ≤ 45%. (B). Panel A shows the distribution of codons on the same two axes. Red dots indicate A and T ending codons, and green dots indicate G and C ending codons. (C). Red dots represent ribosomal genes, yellow dots represent genes with an Aromo value ≥0.15, green dots represent genes with a Gravy value higher than 5, and blue dots represent other genes.
Genes 12 01169 g005
Genes 12 01169 g006 550
Figure 6. PR2-bias plot [A3/(A3 + T3) against G3/(G3 + C3)]. Red circle represents the average position for each plot.
Figure 6. PR2-bias plot [A3/(A3 + T3) against G3/(G3 + C3)]. Red circle represents the average position for each plot.
Genes 12 01169 g006
Genes 12 01169 g007 550
Figure 7. Relationship between ENC and gene expression level for.
Figure 7. Relationship between ENC and gene expression level for.
Genes 12 01169 g007
Genes 12 01169 g008 550
Figure 8. Relationship between ENC and encoded protein length for G. duodenalis.
Figure 8. Relationship between ENC and encoded protein length for G. duodenalis.
Genes 12 01169 g008
Table
Table 1. Codon usage in G. duodenalis.
Table 1. Codon usage in G. duodenalis.
AACodonNRSCUAACodonNRSCU
PheUUU59,5501.07SerUCU73,3181.49
UUC52,1730.93 UCC48,3160.98
LeuUUA31,1520.55 UCA44,6680.91
UUG46,1110.82 UCG32,4960.66
CUU82,2861.47ProCCU41,7551.11
CUC68,9321.23 CCC36,3590.97
CUA44,0310.78 CCA43,7751.17
CUG64,3051.15 CCG28,1580.75
IleAUU64,2661.06ThrACU51,1201.01
AUC62,8511.04 ACC48,6070.96
AUA54,8490.90 ACA64,1621.27
MetAUG71,8241.00 ACG38,7780.77
ValGUU53,7111.12AlaGCU66,0071.04
GUC53,6661.12 GCC64,5371.02
GUA33,6230.70 GCA82,0611.29
GUG50,2861.05 GCG41,0210.65
TyrUAU53,4401.01CysUGU32,7430.82
UAC52,1430.99 UGC46,7571.18
HisCAU35,2690.90ArgCGU26,1860.90
CAC43,2631.10 CGC35,5901.22
GlnCAA50,0960.79 CGA22,0410.76
CAG76,0571.21 CGG24,9100.85
AsnAAU61,6100.97SerAGU38,4080.78
AAC66,0241.03 AGC57,3891.17
LysAAA50,6230.63ArgAGA34,3341.18
AAG109,7171.37 AGG31,9721.10
AspGAU85,4571.00GlyGGU32,9100.77
GAC85,4761.00 GGC49,2001.15
GluGAA72,2470.78 GGA47,8741.12
GAG112,0611.22 GGG40,5860.95
N: The number of codons; the frequently used codons of G. duodenalis are displayed in bold.
Table
Table 2. 15 most frequently used codons of G. duodenalis highly expressed VSPs genes.
Table 2. 15 most frequently used codons of G. duodenalis highly expressed VSPs genes.
Codon Amino Acid Fraction Frequency Number
UGC Cys 0.730 88.230 1042
AAG Lys 0.760 60.203 711
GGC Gly 0.367 40.898 483
GAC Asp 0.575 35.309 417
AAC Asn 0.679 34.208 404
GCC Ala 0.299 32.769 387
ACG Thr 0.325 32.769 387
UGU Cys 0.270 32.684 386
GAG Glu 0.661 30.821 364
GGA Gly 0.268 29.890 353
ACC Thr 0.284 28.620 338
GCG Ala 0.257 28.196 333
GGG Gly 0.242 26.926 318
GAU Asp 0.425 26.080 308
AGC Ser 0.367 26.080 308
Bold, frequently used codons ended with G or C; Italics, frequently used codons ended with A or U.
Table
Table 3. The codon usage table of Giardia intestinalis was obtained from the Codon Usage Database via Optimizer software.
Table 3. The codon usage table of Giardia intestinalis was obtained from the Codon Usage Database via Optimizer software.
CAAFRA.FRE.NCAAFRA.FRE.NCAAFRA.FRE.NCAAFRA.FRE.N
UUUF0.4414.62723UCUS0.2317.73293UAUY0.4414.32658UGUC0.3512.62349
UUCF0.5618.93516UCCS0.1914.12632UACY0.5618.53448UGCC0.6523.84430
UUAL0.075.91095UCAS0.139.91849UAA*0.530.7137UGAW0.070.6108
UUGL0.1110.01861UCGS0.129.21721UAG*0.470.7122UGGW0.937.41378
CUUL0.2421.03091CCUP0.2511.22077CAUH0.367.61421CGUR0.167.91462
CUCL0.2824.64578CCCP0.2711.92215CACH0.6413.42493CGCR0.2813.72552
CUAL0.109.11702CCAP0.2410.51963CAAQ0.3312.02228CGAR0.104.8891
CUGL0.2017.53266CCGP0.2310.21898CAGQ0.6724.54551CGGR0.104.8896
AUUI0.3217.33217ACUT0.2314.92776AAUN0.4016.53074AGUS0.128.91658
AUCI0.4423.74418ACCT0.2718.13377AACN0.6024.44548AGCS0.2216.43061
AUAI0.2513.52505ACAT0.2717.93331AAAK0.2313.82574AGAR0.167.71433
AUGM1.0021.33961ACGT0.2315.32846AAGK0.7745.38433AGGR0.199.41749
GUUV0.2616.23009GCUA0.2319.33583GAUD0.4224.64575GGUG0.1711.62155
GUCV0.3622.74229GCCA0.3025.34707GACD0.5833.66253GGCG0.3423.14291
GUAV0.138.01483GCAA0.2823.24318GAAE0.3118.93516GGAG0.2516.83128
GUGV0.2515.62895GCGA0.1915.92965GAGE0.6941.57719GGGG0.2315.72917
Bold, 10 most frequently used codons. C, Codon; AA, amino acid; FRA., fraction; FRE., frequency; N, number.
Table
Table 4. Translational optimal codons of G. duodenalis.
Table 4. Translational optimal codons of G. duodenalis.
AACodonHigh RSCU (N)Low RSCU (N)AACodonHigh RSCU (N)Low RSCU (N)
PheUUU0.36 (632)1.24 (1846)SerUCU0.88 (972)1.58 (3140)
UUC *1.64 (2845)0.76 (1132) UCC *1.66 (1833)0.75 (1484)
LeuUUA0.07 (94)0.99 (1832) UCA0.27 (301)1.20 (2392)
UUG0.27 (377)0.96 (1777) UCG *1.13 (1249)0.53 (1044)
CUU0.96 (1350)1.31 (2421) AGU0.32 (353)0.98 (1943)
CUC *2.96 (4153)0.76 (1399) AGC *1.74 (1929)0.96 (1915)
CUA0.14 (190)1.06 (1962)ProCCU0.64 (807)1.21 (1684)
CUG *1.61 (2262)0.92 (1690) CCC *1.57 (1974)0.76 (1063)
IleAUU0.51 (801)1.16 (2187) CCA0.38 (480)1.43 (1992)
AUC *2.16 (3418)0.73 (1369) CCG *1.41 (1778)0.59 (824)
AUA0.33 (529)1.11 (2088)ThrACU0.50 (729)1.21 (2058)
MetAUG1.00 (2296)1.00 (2331) ACC *1.30 (1914)0.77 (1306)
ValGUU0.60 (1048)1.13 (1611) ACA0.62 (917)1.42 (2413)
GUC *2.27 (3977)0.74 (1066) ACG *1.58 (2319)0.60 (1013)
GUA0.17 (297)1.11 (1590)AlaGCU0.52 (1271)1.22 (2344)
GUG0.96 (1687)1.02 (1461) GCC *1.70 (4180)0.77 (1487)
TyrUAU0.40 (638)1.21 (1682) GCA0.63 (1551)1.47 (2826)
UAC *1.60 (2521)0.79 (1095) GCG *1.15 (2823)0.54 (1042)
HisCAU0.32 (338)1.13 (1582)CysUGU0.32 (702)1.05 (1125)
CAC *1.68 (1779)0.87 (1221) UGC *1.68 (3650)0.95 (1015)
GlnCAA0.21 (335)1.07 (2489)TrpUGG1.00 (823)1.00 (749)
CAG *1.79 (2878)0.93 (2171)ArgCGU0.57 (547)0.87 (1036)
AsnAAU0.39 (763)1.15 (2189) CGC *2.58 (2471)0.75 (890)
AAC *1.61 (3105)0.85 (1618) CGA0.25 (242)0.93 (1105)
LysAAA0.18 (548)0.92 (2284) CGG0.91 (868)0.89 (1056)
AAG *1.82 (5508)1.08 (2680) AGA0.37 (355)1.57 (1859)
AspGAU0.45 (1310)1.23 (3040) AGG *1.32 (1264)0.99 (1168)
GAC *1.55 (4509)0.77 (1913)GlyGGU0.41 (833)0.96 (1045)
GluGAA0.23 (721)1.00 (2804) GGC *1.87 (3831)0.91 (987)
GAG *1.77 (5474)1.00 (2795) GGA0.54 (1106)1.30 (1407)
GGG *1.18 (2417)0.83 (898)
Comparison of codon usage frequency between high and low expression genes of G. duodenalis. The optimal codons were determined by a Chi-square contingency test. * indicates that the frequency of the codons is much higher (p < 0.01). AA, amino acid; N, number of codons.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, X.; Wang, X.; Gong, P.; Zhang, N.; Zhang, X.; Li, J. Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB. Genes 2021, 12, 1169. https://doi.org/10.3390/genes12081169

AMA Style

Li X, Wang X, Gong P, Zhang N, Zhang X, Li J. Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB. Genes. 2021; 12(8):1169. https://doi.org/10.3390/genes12081169

Chicago/Turabian Style

Li, Xin, Xiaocen Wang, Pengtao Gong, Nan Zhang, Xichen Zhang, and Jianhua Li. 2021. "Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB" Genes 12, no. 8: 1169. https://doi.org/10.3390/genes12081169

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop