Next Article in Journal
Long-Term EEG Component Analysis Method Based on Lasso Regression
Next Article in Special Issue
Machine Learning-Based Prediction of the Seismic Bearing Capacity of a Shallow Strip Footing over a Void in Heterogeneous Soils
Previous Article in Journal
Algorithms for Bidding Strategies in Local Energy Markets: Exhaustive Search through Parallel Computing and Metaheuristic Optimization
Previous Article in Special Issue
A New Constructive Heuristic Driven by Machine Learning for the Traveling Salesman Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of the Codon Table to Quantify the Evolutionary Role of Random Mutations

Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
Submission received: 18 August 2021 / Revised: 14 September 2021 / Accepted: 15 September 2021 / Published: 17 September 2021
(This article belongs to the Special Issue 2021 Selected Papers from Algorithms Editorial Board Members)

Abstract

:
The various biases affecting RNA mutations during evolution is the subject of intense research, leaving the extent of the role of random mutations undefined. To remedy this lacuna, using the codon table, the number of codons representing each amino acid was correlated with the amino acid frequencies in different branches of the evolutionary tree. The correlations were seen to increase as evolution progressed. Furthermore, the number of RNA mutations that resulted in a given amino acid mutation were found to be correlated with several widely used amino acid similarity tables (used in sequence alignments). These correlations were seen to increase when the observed codon usage was factored in.

1. Introduction

It is a fundamental tenet of evolutionary biology that evolution is driven by mutations in the carrier of the genetic code, the DNA, and, consequently, the RNA. While a simplistic view of Darwinian evolution would assume all mutations to be equally probable, the current work in this field is dominated by characterizing various biases in the mutation rates. A recent paper reviewed the research on these biases indicated that biases are categorized as mutation bias, i.e., the difference in the rate of changes at the RNA level (e.g., C–G bias [1]), and composition bias, i.e., the difference in the fitness of the mutated structure, showing bias at the protein level, and showed that they act synergistically [2]. Another recent work [3] emphasized the non-randomness of mutations.
The aim of this paper is to divert the attention from the various biases, and instead attempts to quantify the randomness of the mutations observed during evolution. For this purpose, the results of the mutations will be examined at the protein level. In support of the protein-centric view of mutation, it should be noted that while mutations at the RNA level are interesting biochemical problems, the biological significance of the mutations manifests itself at the protein level.
It has been recognized that if evolution were driven by unbiased random mutations, then the amino acid (AA) frequencies in proteins should fully correlate with the number of codons that code each AA; a correlation was indeed found, but it was not very strong [4]. The number of codons were also found negatively correlated with the AA size and were correlated with the AAs hydrostatic pressure asymmetry index [5]. This paper will examine the codon multiplicity–AA frequency correlation according to families of organisms.
One important use of sequence alignment algorithms is for the detection of the evolutionary relations of proteins. Sequence alignments are based on so-called substitution matrices that measure the similarity of various AAs; in the evolutionary context they should measure the likelihood of one AA to be mutated into another one. The more the RNA mutation are unbiasedly random, the more these substitution matrices should be correlated with the likelihood of AA mutations resulting from unbiased random RNA mutations. To get at these correlations, the codons representing each AA will be used to derive the unbiased likelihood of AA mutations resulting from single or double nucleotide changes followed by the calculation of the correlation between the mutation likelihoods derived from the codon table and various AA substitution matrices. The larger the correlation, the larger the role of randomness in the AA mutations.

2. Materials and Methods

The AA frequencies were obtained for fungi, protozoa, invertebrate, vertebrate (mammals), vertebrate (non-mammals), and plants from the work of Gaur [6]; Tables 1 and 2 of [6] show the AA frequencies for membrane and non-membrane proteins, respectively. The codon table used is from Wikipedia, using the AA coding established for the majority of organisms, with the understanding that there are some organisms where different coding is used. The codon usage frequencies cfreq(i,j,k) were obtained from the GenScript (GenScript Inc., Piscataway, NJ, USA) website (https://www.genscript.com/tools/codon-frequency-table, accessed on 16 September 2021), which provides usage frequencies for a number of organisms, including homo sapiens. The list of AA substitution matrices is described in [7] and was downloaded from the AAindex database (http://www.genome.jp/dbget-bin/www_bfind?aaindex, accessed on 16 September 2021).
The number of single (M = 2) or double (M = 1) RNA mutations, possibly weighted, that change AA a to AA b, m(a,b), is obtained as follows
m ( a , b ) = i , j , k = 1 4 δ ( a , a a ( i , j , k ) )       i , j , k = 1 4 δ ( δ ( i , i ) + δ ( j , j ) + δ ( k , k ) , M ) δ ( b , a a ( i , j , k ) )       W ( i , j , k , i , j , k )
where δ(i,j) is the Kroenecker delta; aa(i,j,k) is the amino acid number (≤20) coded by the nucleotides i, j, and k; and W(i,j,k,i′,j′,k’) is a weight depending on the codons (i,j,k) and (i′,j′,k′) involved.
All of the calculations described in this paper were performed by the program codon.f, available at the URL https://mezeim01.u.hpc.mssm.edu/codon (accessed on 16 September 2021).

3. Results

3.1. Correlation between Codon Multiplicities and AA Propensities

The Pearson correlations (r) between the codon multiplicities and propensities of the AAs (shown in Table 1) were calculated for protozoa, invertebrates, fungi, plants, vertebrates (non-mammals), and vertebrates (mammals); the calculations were performed separately for the membrane and non-membrane proteins. In addition, all correlations were also calculated using only the 11 AAs that were found to be proximal to each other in [4]; for this subset, the correlations were stronger. With the observation in mind that in the evolutionary tree both fungi and plants are on a separate branch of the branch formed by protozoa, invertebrates, and vertebrates, we can conclude that the correlations increased as evolution progressed.

3.2. The Number of RNA Mutations per AA Mutation

Table 2 and Table 3 show for each pair of AAs, the number of single and double RNA mutations, respectively, that result in the AA mutation. First, these data reflect the observation that single mutations are likely to result in small chemical changes [8,9]. While it is not surprising that not all of the 20 × 19/2 = 180 AA mutations can be achieved by a single RNA mutation, it is important to notice that a significant number of AA mutations cannot be achieved even by double nucleotide changes.

3.3. Correlation between the Number of RNA Mutations and AA Substitution Matrices

The larger the role of random mutations in the genetic evolution, the more the number of mutations per AA substitution should correlate with the AA similarity measures. The observation suggests the examination of the correlation between the number of RNA mutations that can change a given AA into another one and the corresponding elements of the various substitution matrices used in the AA sequence alignments.
For the calculation of these correlations, 98 substitution matrices were downloaded from the AAindex database [7]. The list of matrices involves several variants of the two widely used matrix types: PAM [10,11,12,13] and BLOSUM [14]. PAM stands for “point accepted mutation”, the matrices are based on scoring all amino acid positions in the related sequences; the different versions are based on different timeframes for the mutations. BLOSUM stands for “BLOcks SUbstitution Matrix”, and the matrices are based on substitutions and conserved positions in blocks; the closer the compared sequences are, the higher the number version of the matrix that is supposed to be used. It is important to note that these substitution matrices were developed from data at the AA level without considering the requisite RNA mutations. Several of the rest were designed for specific situations. Most, but not all, were symmetric.
Table 4 shows the Pearson correlations between the matrix of the number of mutations, both unweighted and weighted (with the codon usage frequency from the human genome), as shown in Table 2 and Table 3, respectively, and the PAM and BLOSUM matrices; the correlations with the full set of 98 substitution matrices are provided in the Supplementary Materials. For the matrices that were not symmetric, the correlations were calculated for the lower triangular matrix. In addition, the correlation calculations were also performed using only AA mutations that could be achieved with only one or two RNA mutations.

4. Discussion

The variation in the number of codons coding for a given amino acid suggested a number of ways to quantify the role of random mutations in genetic evolution. In particular, the work of Mittal and Jayaram was refined to study the correlation between the number of codons and the propensity of the corresponding AAs to calculate it for different branches of the evolutionary tree and the correlations between the number of mutations changing each AA to other AAs with the corresponding substitution matrices used in sequence alignments. Significant correlations were indeed found, confirming that random mutations indeed form one contribution to genetic evolution, but are not enough to account for it completely. At the same time, the results confirmed the existence of significant biases, as referred to in the Introduction.
The interesting observation from the present calculations is the increase of correlation moving up in the evolutionary tree; this observation also holds for the subset of 11 AAs defined by Mittal and Jayaram. If one accepts the role of randomness of mutations, then this observation provides experimental support for the idea of Darwinian evolution (not that it lacks other, more convincing support). If one accepts Darwinian evolution, then this observation provides experimental support and a measure of the role of random mutations in the evolution.
The correlations found also raise an intriguing question. It has been shown recently that a significant contribution to the foldability of an AA sequence is its adherence to the non-uniform distribution of AAs [15]. This leads to the following question: was there, at the early stages of genetic evolution, a mechanism that biased the establishment of the codon table to use more codons for the more frequent AAs? It is not impossible, as it has already been suggested, that the establishment of the genetic code was influenced by the chemistry of the AAs [16], as it has been found to be likely that the genetic code was established at a later stage of evolution [17].
The fact that the codon table as finalized during evolution ended up with a number of AA mutations that can only be achieved with three RNA mutation can be considered an imperfection. It could have been the result of either too fast evolution (like an optimization process using simulated annealing that was not run slow enough) or of competing requirements (like minimizing the number of RNA mutations, which would result in a drastic AA switch).
In the study examining the correlations between the number of RNA mutations per AA mutations and the substitution matrices, the results were different for AA mutations that required single and those that required double RNA mutations. For single RNA mutations, the results showed significant correlations, both when using all matrix elements and when excluding the matrix elements that correspond to AA mutations that cannot be achieved with a single RNA mutation, although the correlations were weaker for the latter. Furthermore, it is gratifying that the correlations increased when the experimental codon usage probabilities were taken into account. The increased correlations provide additional confirmation for the role of random mutations.
For the double RNA mutations, however, the results were different. When using all matrix elements and without using frequency weighting, the correlations were small, but at least positive, but they become negative when the usage frequencies were factored in. Omitting the “forbidden” AA mutations, contrary to the single RNA mutation case, the correlations increased, but remained small. This suggests that when an alignment indicates AA mutations requiring more than one RNA mutation, it is unlikely to be the result of just a pair of mutations.
One explanation for the low or even negative correlations, using the assumption of random double mutations, could be that such AA changes are the result of insertions and/or frame shifts; this could also explain the occurrence of AA mutations requiring three RNA mutations (which is even more unlikely). The incorporation of this idea could be an alignment algorithm that is done simultaneously at the RNA and the AA level, which could verify the presence of a frame shift.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/a14090270/s1: the correlations of the mutation numbers with the full set of 98 substitution matrices.

Funding

This research received no external funding. It was supported in part through the computational resources and staff expertise provided by the Department of Scientific Computing at the Icahn School of Medicine at Mount Sinai.

Acknowledgments

George Rose, B. Jayaram, and Aditya Mittal are thanked for their helpful suggestions.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. International Human Genome Sequencing Consortium. Initial sequencing and analysis of the human genome. Nature 2001, 409, 860–921. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Cano, A.V.; Payne, J.L.P. Mutation bias interacts with composition bias to influence adaptive evolution. PLoS Comput. Biol. 2020, 16, e1008296. [Google Scholar] [CrossRef] [PubMed]
  3. Caporale, L.H.; Doyle, J. In darwinian evolution, feedback from natural selection leads to biased mutations. Ann. N. Y. Acad. Sci. 2013, 1305, 18–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mittal, A.; Jayaram, B. A possible molecular metric for biological evolvability. J. Biosci. 2012, 37, 573–577. [Google Scholar] [CrossRef] [PubMed]
  5. Giulio, M.D. The origin of the genetic code: Theories and their relationships, a review. BioSystems 2005, 80, 175–184. [Google Scholar] [CrossRef] [PubMed]
  6. Gaur, R.K. Amino acid frequency distribution among eukaryotic proteins. IIOAB J. 2014, 5, 6–11. [Google Scholar]
  7. Kawashima, S.; Pokrowski, P.; Pokarowska, M.; Kolinski, A.; Katayama, T.; Kanehisa, M. AAindex: Amino acid index database, progress report 2008. Nucleic Acids Res. 2008, 36, D202–D205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Giulio, M.D. The extension reached by the minimization of the polarity distances during the evolution of the genetic code. J. Mol. Evol. 1989, 29, 288–293. [Google Scholar] [CrossRef] [PubMed]
  9. Wong, J.T.-F. Role of minimization of chemical distances between amino acids in the evolution of the genetic code. Proc. Natl. Acad. Sci. USA 1980, 77, 1083–1086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Dayhoff, M.O.; Schwartz, R.M.; Orcutt, B.C. A model of evolutionary change in proteins. Atlas Protein Seq. Struct. 1972, 5, 89–99. [Google Scholar]
  11. Altschul, S.F. Amino acid substitution matrices from an information theoretic perspective. J. Mol. Biol. 1991, 219, 555–565. [Google Scholar] [CrossRef]
  12. Benner, S.A.; Cohen, M.A.; Gonnet, G.H. Amino acid substitution during functionally constrained divergent evolution of protein sequences. Protein Eng. 1994, 7, 1323–1332. [Google Scholar] [CrossRef] [Green Version]
  13. Jones, D.T.; Taylor, W.R.; Thornton, J.M. The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 1992, 8, 275–282. [Google Scholar] [CrossRef]
  14. Henikoff, S.; Henikoff, J.G. Amino acid substitution matrices from protein blocks. Proc. Natl. Acad. Sci. USA 1992, 89, 10915–10919. [Google Scholar] [CrossRef] [Green Version]
  15. Mezei, M. On predicting foldability of a protein from its sequence. Proteins 2020, 88, 355–356. [Google Scholar] [CrossRef]
  16. Wong, J.T.-F. A co-evolution theory of the genetic code. Proc. Nat. Acad. Sci. USA 1975, 72, 1909–1912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Giulio, M.D. Reflections on the origin of the genetic code: A hypothesis. J. Theor. Biol. 1998, 191, 191–196. [Google Scholar] [CrossRef]
Table 1. Pearson correlation r between amino acid frequencies and codon multiplicities.
Table 1. Pearson correlation r between amino acid frequencies and codon multiplicities.
All 20Top 11
MembraneNon-MembraneMembraneNon-Membrane
Vertebrates (mammals)0.870.810.970.98
Vertebrates (non-mammals)0.760.730.930.97
Plants0.720.710.880.81
Fungi0.670.750.860.94
Invertebrates0.520.520.820.88
Protozoa0.500.230.780.70
Legend: membrane and non-membrane: correlations calculated for all membrane proteins and all non-membrane proteins, respectively; All 20 and Top 11: correlations calculated for all 20 AAs and for the top 11 AAs in the ranking of [4], respectively.
Table 2. Number of single mutations that result in each AA mutation.
Table 2. Number of single mutations that result in each AA mutation.
Amino AcidAmino Acid CodeNumber of Single Mutations
GLYG12
ALAA412
VALV4412
LEUL00618
ILEI00346
SERS2402214
THRT04003612
ASPD22200002
GLUE222000042
ASNN0000222202
GLNQ00020000202
LYSK000010202422
HISH0002000202402
ARGR600416200022218
PHEF002622000000002
TYRY0000020202002022
TRPW10010100000002000
CYSC200004000000022222
METM0012301000010100000
PROP040404400020240000012
STP 100303002022020422004
GAVLISTDENQKHRFYWCMPSTP
Table 3. Number of double mutations that result in each AA mutation.
Table 3. Number of double mutations that result in each AA mutation.
Amino AcidAmino Acid CodeNumber of Double Mutations
GLYG0
ALAA120
VALV12120
LEUL661812
ILEI339140
SERS101461274
THRT412469180
ASPD66622420
GLUE666412200
ASNN2222446240
GLNQ22281224240
LYSK222456642020
HISH2226242242040
ARGR1866188121024684612
PHEF226648220200220
TYRY2226282242442220
TRPW31120510101104220
CYSC6226282202002102200
METM1134033012110220100
PROP4124143141222222614221210
STP 5336211343435486214132
GAVLISTDENQKHRFYWCMPSTP
Table 4. Correlation between PAM or BLOSUM matrices and the number of mutation paths.
Table 4. Correlation between PAM or BLOSUM matrices and the number of mutation paths.
All AA Mutations IncludedAA Mutations Requiring Two/Three RNA Mutations Excluded
Single MutationsDouble MutationsSingle MutationsDouble Mutations
Code 1Corr_n 2Corr_w 2Corr_nCorr_wCorr_nCorr_wCorr_nCorr_w
ALTS910101 30.5080.5400.080−0.0910.3890.3640.3600.172
BENS940101 40.4880.6260.073−0.0750.3060.5070.3780.250
BENS940102 50.4690.5770.051−0.0900.3330.4860.3460.230
BENS940103 60.4470.518−.007−0.1140.3630.4620.2700.221
DAYM780301 70.4490.4800.0550.0850.3380.3260.3070.167
JOND920103 80.4760.5960.057−0.0960.3180.4860.3760.237
JOND940101 90.4040.591−0.043−0.1110.1690.4440.2160.243
DAYM780302 100.5620.5650.071−0.1040.4810.3500.3390.247
HENS920101 110.4870.491−0.093−0.2240.4200.4070.2220.125
HENS920102 120.5060.518−0.072−0.2060.4400.4500.2350.139
HENS920103 130.5210.532−0.061−0.2050.4410.4470.2520.139
HENS920104 140.4830.498−0.103−0.2140.4130.4010.1140.187
1 Code: the AAindex database identifier; 2 Corr_n, Corr_w: number of mutation paths calculated with W = 1 and W = cfreq(i,j,k) × cfreq(i,j′,k′), respectively; 3 PAM-120 [11]; 4 PAM-log-odds_6.4–8.7 [12]; 5 PAM-log-odds_22–29 [12]; 6 PAM-log-odds_74–100 [12]; 7 PAM-250-log-odds [10]; 8 PAM-250_PET91 [13]; 9 PAM-250_TM [13]; 10 PAM-40–log-odds [10]; 11 BLOSUM45 [14]; 12 BLOSUM62 [14]; 13 BLOSUM80 [14]; 14 BLOSUM50 [14].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mezei, M. Use of the Codon Table to Quantify the Evolutionary Role of Random Mutations. Algorithms 2021, 14, 270. https://0-doi-org.brum.beds.ac.uk/10.3390/a14090270

AMA Style

Mezei M. Use of the Codon Table to Quantify the Evolutionary Role of Random Mutations. Algorithms. 2021; 14(9):270. https://0-doi-org.brum.beds.ac.uk/10.3390/a14090270

Chicago/Turabian Style

Mezei, Mihaly. 2021. "Use of the Codon Table to Quantify the Evolutionary Role of Random Mutations" Algorithms 14, no. 9: 270. https://0-doi-org.brum.beds.ac.uk/10.3390/a14090270

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