The Influence of the Selection at the Amino Acid Level on Synonymous Codon Usage from the Viewpoint of Alternative Genetic Codes.

Synonymous codon usage can be influenced by mutations and/or selection, e.g., for speed of protein translation and correct folding. However, this codon bias can also be affected by a general selection at the amino acid level due to differences in the acceptance of the loss and generation of these codons. To assess the importance of this effect, we constructed a mutation-selection model model, in which we generated almost 90,000 stationary nucleotide distributions produced by mutational processes and applied a selection based on differences in physicochemical properties of amino acids. Under these conditions, we calculated the usage of fourfold degenerated (4FD) codons and compared it with the usage characteristic of the pure mutations. We considered both the standard genetic code (SGC) and alternative genetic codes (AGCs). The analyses showed that a majority of AGCs produced a greater 4FD codon bias than the SGC. The mutations producing more thymine or adenine than guanine and cytosine increased the differences in usage. On the other hand, the mutational pressures generating a lot of cytosine or guanine with a low content of adenine and thymine decreased this bias because the nucleotide content of most 4FD codons stayed in the compositional equilibrium with these pressures. The comparison of the theoretical results with those for real protein coding sequences showed that the influence of selection at the amino acid level on the synonymous codon usage cannot be neglected. The analyses indicate that the effect of amino acid selection cannot be disregarded and that it can interfere with other selection factors influencing codon usage, especially in AT-rich genomes, in which AGCs are usually used.

The structure of the genetic code as an optimal graph clustering problem.

The standard genetic code (SGC) is the set of rules by which genetic information is translated into proteins, from codons, i.e. triplets of nucleotides, to amino acids. The questions about the origin and the main factor responsible for the present structure of the code are still under a hot debate. Various methodologies have been used to study the features of the code and assess the level of its potential optimality. Here, we introduced a new general approach to evaluate the quality of the genetic code structure. This methodology comes from graph theory and allows us to describe new properties of the genetic code in terms of conductance. This parameter measures the robustness of codon groups against the potential changes in translation of the protein-coding sequences generated by single nucleotide substitutions. We described the genetic code as a partition of an undirected and unweighted graph, which makes the model general and universal. Using this approach, we showed that the structure of the genetic code is a solution to the graph clustering problem. We presented and discussed the structure of the codes that are optimal according to the conductance. Despite the fact that the standard genetic code is far from being optimal according to the conductance, its structure is characterised by many codon groups reaching the minimum conductance for their size. The SGC represents most likely a local minimum in terms of errors occurring in protein-coding sequences and their translation.

The influence of local environmental factors in southwestern Poland on the abundance of Ixodes ricinus and prevalence of infection with Borrelia burgdorferi s.l. and B. miyamotoi.

Ticks are important ectoparasites and vectors of pathogens that cause disease in humans and animals. The natural habitat of Ixodes ricinus ticks is forests, which are convenient habitats to search for hosts, including reservoir hosts, and therefore can be an important habitat source of tick-borne pathogens. The aim of the study was to assess the usefulness of detailed forest habitat-type maps to estimate the tick-borne risk at a local scale (Lower Silesia, SW Poland). For the purposes of estimating tick abundance, we used the land cover maps available from the Forest Data Bank. For I. ricinus collection, nine sites located in three forest habitat types were chosen: broadleaf forest, mixed broadleaf and coniferous forest and coniferous forest. Ticks were collected once a month from April to June 2018 and 2019 using the standard flagging method. At each of the nine sites, ticks were collected in four plots, of 100 m2 each. Tick abundance was analysed using general linear mixed models (GLMM). A total of 2196 (10.1/100 m2) ticks were collected, including 2093 Ixodes ricinus (95.3%; 9.6/100 m2), 46 Dermacentor reticulatus (2.1%; 0.2/100 m2) and 57 Haemaphysalis concinna (2.6%; 0.3/100 m2). Among the collected I. ricinus were 589 larvae (28.1%; 2.7/100 m2), 1261 nymphs (60.3%; 5.8/100 m2), 128 females (6.1%; 0.6/100 m2) and 115 males (5.5%; 0.5/100 m2). We found a highly significant effect of forest habitat type on the density of ticks for broadleaf forest (coefficient = 1.87267, p-value = 2.79e - 07). Additionally, a significant influence of air temperature and relative humidity on the abundance of ticks was observed. During spring, the peak activity of I. ricinus was recorded in May and June. For DNA amplification of Borrelia burgdorferi s.l., a nested PCR method was used. Out of 494 I. ricinus, 83 (16.8%) were positive for Borrelia spp. The RFLP method showed the occurrence of five species including four belonging to the B. burgdorferi s.l. complex: B. afzelii (30.1%), B. garinii (38.6%), B. valaisiana (2.4%) and B. lusitaniae (18.1%). Furthermore, B. miyamotoi (9.6%), a species belonging to bacteria that cause relapsing fever as well as co-infection of B. miyamotoi/B.lusitaniae (1.2%) were found. The differences in the infection level of Borrelia spp. between broadleaf forest and mixed broadleaf and coniferous forest were statistically significant.

Model of Genetic Code Structure Evolution under Various Types of Codon Reading.

The standard genetic code (SGC) is a set of rules according to which 64 codons are assigned to 20 canonical amino acids and stop coding signal. As a consequence, the SGC is redundant because there is a greater number of codons than the number of encoded labels. This redundancy implies the existence of codons that encode the same genetic information. The size and organization of such synonymous codon blocks are important characteristics of the SGC structure whose evolution is still unclear. Therefore, we studied possible evolutionary mechanisms of the codon block structure. We conducted computer simulations assuming that coding systems at early stages of the SGC evolution were sets of ambiguous codon assignments with high entropy. We included three types of reading systems characterized by different inaccuracy and pattern of codon recognition. In contrast to the previous study, we allowed for evolution of the reading systems and their competition. The simulations performed under minimization of translational errors and reduction of coding ambiguity produced the coding system resistant to these errors. The reading system similar to that present in the SGC dominated the others very quickly. The survived system was also characterized by low entropy and possessed properties similar to that in the SGC. Our simulation show that the unambiguous SGC could emerged from a code with a lower level of ambiguity and the number of tRNAs increased during the evolution.

The influence of different types of translational inaccuracies on the genetic code structure.

BACKGROUND: The standard genetic code is a recipe for assigning unambiguously 21 labels, i.e. amino acids and stop translation signal, to 64 codons. However, at early stages of the translational machinery development, the codons did not have to be read unambiguously and the early genetic codes could have contained some ambiguous assignments of codons to amino acids. Therefore, the goal of this work was to obtain the genetic code structures which could have evolved assuming different types of inaccuracy of the translational machinery starting from unambiguous assignments of codons to amino acids. RESULTS: We developed a theoretical model assuming that the level of uncertainty of codon assignments can gradually decrease during the simulations. Since it is postulated that the standard code has evolved to be robust against point mutations and mistranslations, we developed three simulation scenarios assuming that such errors can influence one, two or three codon positions. The simulated codes were selected using the evolutionary algorithm methodology to decrease coding ambiguity and increase their robustness against mistranslation. CONCLUSIONS: The results indicate that the typical codon block structure of the genetic code could have evolved to decrease the ambiguity of amino acid to codon assignments and to increase the fidelity of reading the genetic information. However, the robustness to errors was not the decisive factor that influenced the genetic code evolution because it is possible to find theoretical codes that minimize the reading errors better than the standard genetic code.

Many alternative and theoretical genetic codes are more robust to amino acid replacements than the standard genetic code.

We evaluated the differences between the standard genetic code (SGC) and its known alternative variants in terms of the consequences of amino acids replacements. Furthermore, the properties of all the possible theoretical genetic codes, which differ from the SGC by one, two or three changes in codon assignments were also tested. Although the SGC is closer to the best theoretical codes than to the worst ones due to the minimization of amino acid replacements, from 10% to 27% of the all possible theoretical codes minimize the effect of these replacements better than the SGC. Interestingly, many types of codon reassignments observed in the alternative codes are also responsible for the substantial robustness to amino acid replacements. As many as 18 out of 21 alternatives perform better than the SGC under the assumed optimization criteria. These findings suggest that not all reassignments in the alternative codes are neutral and some of them could be selected to reduce harmful effects of mutations or translation of protein-coding sequences. The results also imply that the standard genetic code can be improved in this respect by a quite small number of changes, which are in fact realized in its variants. It would mean that the tendency to minimize mutational errors was not the main force that drove the evolution of the SGC.

The optimality of the standard genetic code assessed by an eight-objective evolutionary algorithm.

BACKGROUND: The standard genetic code (SGC) is a unique set of rules which assign amino acids to codons. Similar amino acids tend to have similar codons indicating that the code evolved to minimize the costs of amino acid replacements in proteins, caused by mutations or translational errors. However, if such optimization in fact occurred, many different properties of amino acids must have been taken into account during the code evolution. Therefore, this problem can be reformulated as a multi-objective optimization task, in which the selection constraints are represented by measures based on various amino acid properties. RESULTS: To study the optimality of the SGC we applied a multi-objective evolutionary algorithm and we used the representatives of eight clusters, which grouped over 500 indices describing various physicochemical properties of amino acids. Thanks to that we avoided an arbitrary choice of amino acid features as optimization criteria. As a consequence, we were able to conduct a more general study on the properties of the SGC than the ones presented so far in other papers on this topic. We considered two models of the genetic code, one preserving the characteristic codon blocks structure of the SGC and the other without this restriction. The results revealed that the SGC could be significantly improved in terms of error minimization, hereby it is not fully optimized. Its structure differs significantly from the structure of the codes optimized to minimize the costs of amino acid replacements. On the other hand, using newly defined quality measures that placed the SGC in the global space of theoretical genetic codes, we showed that the SGC is definitely closer to the codes that minimize the costs of amino acids replacements than those maximizing them. CONCLUSIONS: The standard genetic code represents most likely only partially optimized systems, which emerged under the influence of many different factors. Our findings can be useful to researchers involved in modifying the genetic code of the living organisms and designing artificial ones.

New insights into the criteria of functional heterozygosity of the Apis mellifera complementary sex determining gene-Discovery of a functional allele pair differing by a single amino acid.

The complementary sex determiner (csd) gene is responsible for controlling the sex-determination molecular switch in western honey bees (Apis mellifera): bees that are heterozygous for csd develop into females, whereas bees that are hemizygous or homozygous develop into males. The homozygous diploid males are destroyed at an early stage of their development. It has been proposed that the minimal number of amino acid differences between two csd alleles needed to fully determine femaleness is five and it has also been shown that smaller differences may result in forming an evolutionary intermediate that is not fully capable of female determination, but has increased fitness compared to the homozygous genotype. In this study, we have implemented a terminal restriction length polymorphism-based method of identifying and distinguishing paternal alleles in a given bee colony and assigning them to a particular maternal allele in order to gather information on large number of functional csd pairs and also to identify, to some extent, genotypes that are underrepresented or absent in bee colonies. The main finding of this study is the identification of a fully functional genotype consisting of csd alleles that differed from each other by a one amino acid position. The individuals carrying this genotype expressed only female-specific transcripts of feminizer and double-sex genes. By comparing the sequences differences between the csd pair identified in our study with those described earlier, we conclude that functional heterozygosity of the csd gene is dependent not only on the number of the amino acid differences but also on the sequence context and position of the change. The discovery of a functional allele pair differing by a single amino acid also implies that the generation of a new csd specificity may also occur during a single mutation step with no need for evolutionary intermediates accumulating further mutations.

Models of genetic code structure evolution with variable number of coded labels.

It is assumed that at the early stage of cell evolution its translation machinery was characterized by high noise, i.e. ambiguous assignment of codons to amino acids in the genetic code, which initially encoded only few amino acids. Next, during its evolution new amino acids were added to this code. Taking into account this facts, we investigated theoretical models of genetic code's structure, which evolved from a set of ambiguous codons assignments into a coding system with a low level of uncertainty. We considered three types of translational inaccuracies assuming a different number of fixed codon positions. We applied a modified version of evolutionary algorithm for finding the genetic codes that the most effectively reduced the initial uncertainty in the assignment of codons to encoded labels, i.e. amino acids and a stop translation signal. We examined codes with the number of labels from four to 22. Our results indicated that the quality of genetic code structure is strongly dependent on the number of encoded labels as well as the type of translational mechanism. The more strict assignments of codon to the labels was preferred by the codes encoding more number of labels. The results showed that a smaller degeneracy of codes evolved from a more tolerant coding with the stepwise addition of coded amino acids to the genetic code. The distribution of codon groups in the standard genetic code corresponds well to the translation model assuming two fixed codon positions, whereas the six-codon groups can be relics form previous stages of evolution when the code characterized by a greater uncertainty.

Some theoretical aspects of reprogramming the standard genetic code.

Reprogramming of the standard genetic code to include non-canonical amino acids (ncAAs) opens new prospects for medicine, industry, and biotechnology. There are several methods of code engineering, which allow us for storing new genetic information in DNA sequences and producing proteins with new properties. Here, we provided a theoretical background for the optimal genetic code expansion, which may find application in the experimental design of the genetic code. We assumed that the expanded genetic code includes both canonical and non-canonical information stored in 64 classical codons. What is more, the new coding system is robust to point mutations and minimizes the possibility of reversion from the new to old information. In order to find such codes, we applied graph theory to analyze the properties of optimal codon sets. We presented the formal procedure in finding the optimal codes with various number of vacant codons that could be assigned to new amino acids. Finally, we discussed the optimal number of the newly incorporated ncAAs and also the optimal size of codon groups that can be assigned to ncAAs.

Basic principles of the genetic code extension.

Compounds including non-canonical amino acids (ncAAs) or other artificially designed molecules can find a lot of applications in medicine, industry and biotechnology. They can be produced thanks to the modification or extension of the standard genetic code (SGC). Such peptides or proteins including the ncAAs can be constantly delivered in a stable way by organisms with the customized genetic code. Among several methods of engineering the code, using non-canonical base pairs is especially promising, because it enables generating many new codons, which can be used to encode any new amino acid. Since even one pair of new bases can extend the SGC up to 216 codons generated by a six-letter nucleotide alphabet, the extension of the SGC can be achieved in many ways. Here, we proposed a stepwise procedure of the SGC extension with one pair of non-canonical bases to minimize the consequences of point mutations. We reported relationships between codons in the framework of graph theory. All 216 codons were represented as nodes of the graph, whereas its edges were induced by all possible single nucleotide mutations occurring between codons. Therefore, every set of canonical and newly added codons induces a specific subgraph. We characterized the properties of the induced subgraphs generated by selected sets of codons. Thanks to that, we were able to describe a procedure for incremental addition of the set of meaningful codons up to the full coding system consisting of three pairs of bases. The procedure of gradual extension of the SGC makes the whole system robust to changing genetic information due to mutations and is compatible with the views assuming that codons and amino acids were added successively to the primordial SGC, which evolved minimizing harmful consequences of mutations or mistranslations of encoded proteins.

Optimization of the standard genetic code in terms of two mutation types: Point mutations and frameshifts.

The distinct structure and universality of the standard genetic code (SGC) have fascinated the scientists ever since the first amino acid assignments were discovered. There are several hypotheses trying to explain the origin and evolution of this code. One of them postulates that the SGC evolved to minimize harmful effects of amino acid replacements in proteins, caused by mutations and translational errors. Many investigations concerning this hypothesis have already been carried out, but they were focused mainly on the consequences of single-point mutations. Therefore, we decided to check the influence of other types of mutations, i.e. insertions and deletions, on the robustness to amino acid replacements of the SGC. Such mutations cause shifts in the reading frame during the translation process which result in more harmful consequences in coded proteins than in the case of single-point mutations. We applied a multi-objective optimization algorithm to find the best and the worst genetic codes, regarding their robustness to both single-point and frameshift mutations, for various amino acid properties. Then we compared the features of the found codes with the properties of the standard genetic code. The results show that the SGC is not fully optimized for minimizing the effects of frameshift mutations but it is, nevertheless, much closer to the best solutions than to the worst ones. It implies that a certain tendency to minimize the costs of amino acids replacements resulting from various kinds of mutations is present in the standard genetic code.

Correction: Optimization of the standard genetic code according to three codon positions using an evolutionary algorithm.

[This corrects the article DOI: 10.1371/journal.pone.0201715.].

Optimization of the standard genetic code according to three codon positions using an evolutionary algorithm.

Many biological systems are typically examined from the point of view of adaptation to certain conditions or requirements. One such system is the standard genetic code (SGC), which generally minimizes the cost of amino acid replacements resulting from mutations or mistranslations. However, no full consensus has been reached on the factors that caused the evolution of this feature. One of the hypotheses suggests that code optimality was directly selected as an advantage to preserve information about encoded proteins. An important feature that should be considered when studying the SGC is the different roles of the three codon positions. Therefore, we investigated the robustness of this code regarding the cost of amino acid replacements resulting from substitutions in these positions separately and the sum of these costs. We applied a modified evolutionary algorithm and included four models of the genetic code assuming various restrictions on its structure. The SGC was compared both with the codes that minimize the objective function and those that maximize it. This approach allowed us to place the SGC in the global space of possible codes, which is a more appropriate and unbiased comparison than that with randomly generated codes because they are characterized by relatively uniform amino acid assignments to codons. The SGC appeared to be well optimized at the global scale, but its individual positions were not fully optimized because there were codes that were optimized for only one codon position and simultaneously outperformed the SGC at the other positions. We also found that different code structures may lead to the same optimality and that random codes can show a tendency to minimize costs under some of the genetic code models. Our results suggest that the optimality of SGC could be a by-product of other processes.

Role of recombination and faithfulness to partner in sex chromosome degeneration.

Sex determination in mammals is strongly linked to sex chromosomes. In most cases, females possess two copies of X chromosome while males have one X and one Y chromosome. It is assumed that these chromosomes originated from a pair of homologous autosomes, which diverged when recombination between them was suppressed. However, it is still debated why the sex chromosomes stopped recombining and how this process spread out over most part of the chromosomes. To study this problem, we developed a simulation model, in which the recombination rate between the sex chromosomes can freely evolve. We found that the suppression of recombination between the X and Y is spontaneous and proceeds very quickly during the evolution of population, which leads to the degeneration of the Y in males. Interestingly, the degeneration happens only when mating pairs are unfaithful. This evolutionary strategy purifies the X chromosome from defective alleles and leads to the larger number of females than males in the population. In consequence, the reproductive potential of the whole population increases. Our results imply that both the suppression of recombination and the degeneration of Y chromosome may be associated with reproductive strategy and favoured in polygamous populations with faithless mating partners.

The influence of habitat preferences on shell morphology in ecophenotypes of Trochulus hispidus complex.

Trochulus hispidus and T. sericeus are hairy snails widely distributed in Europe. They differ in shell morphology and are usually found in various land habitats. However, their morphology does not match genetic distance as they do not form distinct clades. Therefore, it is interesting to determine to what extent environmental factors can control their phenotypes. We analysed the morphological traits and many environmental features of their habitats to find relationships between these parameters and explain ecological reasons for this plasticity. We found many statistically significant correlations between morphological traits and environmental variables. Illumination, forestation, precipitation and temperature occurred the most important features discriminating habitats of these snails. It turned out that T. sericeus prefers forests and moist shaded places, while T. hispidus chooses more dry habitats and open areas exposed to the sun. T. sericeus is also probably more tolerant to low and variable temperatures. The hair durability is also correlated with their habitats: the shell of T. hispidus is mostly hairless but hairs almost always cover the shell of T. sericeus. These results support the hypothesis that the lack of hairs is associated with the loss of a potential adaptive function due to the change from wet to dry habitats. The hairs facilitate the adherence of snails to herbaceous plants during feeding when the humidity levels are high. The morphological divergence of T. hispidus and T. sericeus is the result of phenotypic plasticity and selection associated with the habitat, which affect both the shell shape and the hair durability. Since T. hispidus and T. sericeus do to not represent separate biological species and their variability has no genetic basis, they should be considered as ecophenotypes. This and our previous studies suggest that phenotypic plasticity in widely distributed Trochulus species is quite common and may have been of ancestral origin.

Representations of Search Spaces in the Problem of Mutational Pressure Optimization According to Protein-Coding Sequences.

The proper representation of the search space is the fundamental step in every optimization task, because it has a decisive impact on the quality of potential solutions. In particular, this problem appears when the search spaces are nonstandard and complex, with the large number of candidate solutions that differ from classical forms usually investigated. One of such spaces is the set of continuous-time, homogenous, and stationary Markov processes. They are commonly used to describe biological phenomena, for example, mutations in DNA sequences and their evolution. Because of the complexity of these processes, the representation of their search space is not an easy task but it is important for effective solving of the biological problems. One of them is optimality of mutational pressure acting on protein-coding sequences. Therefore, we described three representations of the search spaces and proposed several specific evolutionary operators that are used in evolutionary-based optimization algorithms to solve the biological problem of mutational pressure optimality. In addition, we gave a general formula for the fitness function, which can be used to measure the quality of potential solutions. The structures of these solutions are based on two models of DNA evolution described by substitution-rate matrices, which are commonly used in phylogenetic analyzes. The proposed representations have been successfully utilized in various issues, and the obtained results are very interesting from a biological point of view. For example, they show that mutational pressures are, to some extent, optimized to minimize cost of amino acid substitutions in proteins.

The Impact of Selection at the Amino Acid Level on the Usage of Synonymous Codons.

There are two main forces that affect usage of synonymous codons: directional mutational pressure and selection. The effectiveness of protein translation is usually considered as the main selectional factor. However, biased codon usage can also be a byproduct of a general selection at the amino acid level interacting with nucleotide replacements. To evaluate the validity and strength of such an effect, we superimposed >3.5 billion unrestricted mutational processes on the selection of nonsynonymous substitutions based on the differences in physicochemical properties of the coded amino acids. Using a modified evolutionary optimization algorithm, we determined the conditions in which the effect on the relative codon usage is maximized. We found that the effect is enhanced by mutational processes generating more adenine and thymine than guanine and cytosine, as well as more purines than pyrimidines. Interestingly, this effect is observed only under an unrestricted model of nucleotide substitution, and disappears when the mutational process is time-reversible. Comparison of the simulation results with data for real protein coding sequences indicates that the impact of selection at the amino acid level on synonymous codon usage cannot be neglected. Furthermore, it can considerably interfere, especially in AT-rich genomes, with other selections on codon usage, e.g., translational efficiency. It may also lead to difficulties in the recognition of other effects influencing codon bias, and an overestimation of protein coding sequences whose codon usage is subjected to adaptational selection.

Optimization of amino acid replacement costs by mutational pressure in bacterial genomes.

Mutations are considered a spontaneous and random process, which is important component of evolution because it generates genetic variation. On the other hand, mutations are deleterious leading to non-functional genes and energetically costly repairs. Therefore, one can expect that the mutational pressure is optimized to simultaneously generate genetic diversity and preserve genetic information. To check if empirical mutational pressures are optimized in these ways, we compared matrices of nucleotide mutation rates derived from bacterial genomes with their best possible alternatives that minimized or maximized costs of amino acid replacements associated with differences in their physicochemical properties (e.g. hydropathy and polarity). It should be noted that the studied empirical nucleotide substitution matrices and the costs of amino acid replacements are independent because these matrices were derived from sites free of selection on amino acid properties and the amino acid costs assumed only amino acid physicochemical properties without any information about mutation at the nucleotide level. Obtained results indicate that the empirical mutational matrices show a tendency to minimize costs of amino acid replacements. It implies that bacterial mutational pressures can evolve to decrease consequences of amino acid substitutions. However, the optimization is not full, which enables generation of some genetic variability.

Uneven distribution of complementary sex determiner (csd) alleles in Apis mellifera population.

The complementary sex determiner (csd) gene determines the sex of the western honey bee (Apis mellifera L.). Bees that are heterozygous at the csd locus develop into females; whereas hemizygous bees develop into males. The co-occurrence of two identical csd alleles in a single diploid genome leads to the genetic death of the bee. Thus, the maintenance of csd diversity in the population is favoured. The number and distribution of csd alleles is particularly interesting in light of the recent decline in the honey bee population. In this study, we analysed the distribution of csd alleles in two Polish populations separated by about 100 km. We analysed the maternal alleles of 193 colonies and found 121 different alleles. We also analysed the distribution and frequency of the alleles, and found that they are distributed unevenly. We show that the methods that have been used so far to estimate the total worldwide number of csd alleles have significantly underestimated their diversity. We also show that the uneven distribution of csd alleles is caused by a large number of infrequent alleles, which most likely results from the fact that these alleles are generated very frequently.