BACKGROUND/OBJECTIVES: Codon usage bias affects gene expression and translation efficiency across species. The effective number of codons (ENC) and GC content influence codon preference, often displaying unimodal or bimodal distributions. This study investigates the correlation between ENC and GC rankings across species and how their relationship affects codon usage distributions.

METHODS

I analyzed nuclear-encoded genes from 17 species representing six kingdoms: one bacteria (), three fungi (, , and ), one archaea (), three protists (, , and ),), three plants (, , and ), and six animals (, , , , , and ). Genes in all 17 species were ranked by GC content and ENC, and correlations were assessed. I examined how adding or subtracting these rankings influenced their overall distribution in a new method that I call Two-Rank Order Normalization or TRON. The equation, TRON = SUM(ABS((GC rank:GC rank) - (ENC rank:ENC rank))/(N/3), where (GC rank:GC rank) is a rank-order series of GC rank, (ENC rank:ENC rank) is a rank-order series ENC rank, sorted by the rank-order series GC rank. The denominator of TRON, N/3, is the normalization factor because it is the expected value of the sum of the absolute value of GC rank-ENC rank for all genes if GC rank and ENC rank are not correlated.

RESULTS

ENC and GC rankings are positively correlated (i.e., ENC increases as GC increases) in AT-rich species such as honeybees (R = 0.60, slope = 0.78) and wasps (R = 0.52, slope = 0.72) and negatively correlated (i.e., ENC decreases as GC increases) in GC-rich species such as humans (R = 0.38, slope = -0.61) and rice (R = 0.59, slope = -0.77). Second, the GC rank-ENC rank distributions change from unimodal to bimodal as GC content increases in the 17 species. Third, the GC rank+ENC rank distributions change from bimodal to unimodal as GC content increases in the 17 species. Fourth, the slopes of the correlations (GC versus ENC) in all 17 species are negatively correlated with TRON (R = 0.98) (see Graphic Abstract).

CONCLUSIONS

The correlation between ENC rank and GC rank differs among species, shaping codon usage distributions in opposite ways depending on whether a species' nuclear-encoded genes are AT-rich or GC-rich. Understanding these patterns might provide insights into translation efficiency, epigenetics mediated by CpG DNA methylation, epitranscriptomics of RNA modifications, RNA secondary structures, evolutionary pressures, and potential applications in genetic engineering and biotechnology.