Directional Mutation Pressure Research Articles

De Novo Assembly and Characterization of the Transcriptome of an Omnivorous Camel Cricket (Tachycines meditationis).

Tachycines meditationis (Orthoptera: Rhaphidophoridae: Tachycines) is a widely distributed insect in eastern Asia. This species is common in urban environments, and its unique omnivorous diet may contribute to its success in various habitats. However, molecular studies on the species are scarce. Here, we obtained the first transcriptome sequence of T. meditationis and performed preliminary analyses to test whether the evolution of coding sequences fits the expectations based on the species' ecology. We retrieved 476,495 effective transcripts and annotated 46,593 coding sequences (CDS). We analysed the codon usage and found that directional mutation pressure was the leading cause of codon usage bias in this species. This genome-wide relaxed codon usage pattern in T. meditationis is surprising, given the potentially large population size of this species. Moreover, despite the omnivorous diet, the chemosensory genes of this species do not exhibit codon usage deviating significantly from the genome-level pattern. They also do not seem to experience more gene family expansion than other cave cricket species do. A thorough search for rapidly evolved genes using the dN/dS value showed that genes associated with substance synthesis and metabolic pathways, such as retinol metabolism, aminoacyl-tRNA biosynthesis, and fatty acid metabolism, underwent species-specific positive selection. While some results seem to contradict the species ecology, our transcriptome assembly provides a valuable molecular resource for future studies on camel cricket evolution and molecular genetics for feeding ecology in insects, in general.

A detailed comparative analysis of codon usage bias in Alongshan virus

Alongshan virus (ALSV) is an emerging tick-borne pathogen that infects humans, causing febrile disease. ALSV uses Ixodes Persulcatus ticks to infect humans with a wide range of signs, from asymptomatic to encephalitis-like syndrome. There is an increasing public health concern about the ALSV infection. To get insight into the impacts of viral relations with their hosts on viral ability, survival, and evasion from hosts immune systems remain unknown. The codon usage is a driving force in viral genome evolution; therefore, we enrolled 41 ALSV strains in codon usage analysis to elucidate the molecular evolutionary dynamics of ALSV. The results indicate that the overall codon usage among ALSV isolates is relatively similar and slightly biased. Base compositions for the cds were in order of G >A >C >U and in the third position of codons G3 >A3 >C3 >T3. The RSCU values revealed that the more frequently used codons were mostly GC ended. Different codon preferences in ALSV genes in relation to codon usage of H. sapiens and Ixodes Persulcatus genes were found. Neutrality plot was determined to reveal the superiority of natural selection over directional mutation pressure in causing CUB based on GC12 versus GC3 contents. The results of these studies suggest that the emergence of ALSV in China, Russia and Finland may also be reflected in ALSV codon usage. Altogether, the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of ALSV.

Composition, codon usage pattern, protein properties, and influencing factors in the genomes of members of the family Anelloviridae.

The present study was carried out on 62 genome sequences of members of the family Anelloviridae, as there have been no reports of genome analysis of these DNA viruses using a bioinformatics approach. The genes were found to be rich in AC content with low codon usage bias (CUB). Relative synonymous codon usage (RSCU) values identified the preferred codons for each amino acid in the family. The codon AGA was overrepresented, while the codons TCG, TTG, CGG, CGT, ACG, GCG and GAT were underrepresented in all of the genomes. A significant correlation was found between the effective number of codons (ENC) and base constraints, indicating that compositional properties might have influenced the CUB. A highly significant correlation was observed between the overall base content and the base content at the third codon position, indicating that mutations might have affected the CUB. A highly significant positive correlation was observed between GC12 and GC3 (r = 0.904, p < 0.01), which indicated that directional mutation pressure influenced all three codon positions. A neutrality plot revealed that the contribution of mutation and natural selection in determining the CUB was 58.6% and 41.4%, respectively.Electronic supplementary materialThe online version of this article (10.1007/s00705-020-04890-2) contains supplementary material, which is available to authorized users.

Host-directed editing of the SARS-CoV-2 genome

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.

Codon usage bias in the H gene of canine distemper virus

Canine distemper virus (CDV), a non-segmented single negative-stranded RNA (ssRNA), is the etiological agent of canine distemper. Canine distemper is a highly contagious and lethal viral disease in domestic dogs and wild carnivores. Study of the evolution of CDV presents an essential key to improve the vaccine efficacy. In this study, a total of 328 full-length CDV hemagglutinin (H) gene sequences were subjected to phylogenetic, amino acid mutations, and codon usage analysis. In accordance with previous study, CDV genotypes consisted of fifteen lineages. The unique amino acid substitution sites in each CDV lineages have been identified for the first time, including America-1 (Q330H), America-2 (I585S), Asia-1 (A359V), Asia-2 (H61R), Asia-3 (P108Q), Asia-4 (K213T), India-1/Asia-5(S497P), Arctic (S20L), Africa-1(N489S), Colombian (V41I), EWL (I44V), Europe (D560E), Europe-1/South America-1(K161Q), South America-2 (R580Q), and East African (S214A). Codon usage analysis indicated that H gene exhibited low codon usage bias and further neutrality plot analysis demonstrated that natural selection played a dominated role in driving CPV evolution. The effective number of codons (ENC) plots show that all the different sequences are below the standard curve, indicating that mutational pressure is not the only factor affecting CUB but other forces, including natural selection. The neutrality analysis showed that the slope of the regression line was 0.1501, indicating natural selection dominates directional mutation pressure in driving the codon usage pattern. In addition, nucleotide composition, relative synonymous codon usage value, dinucleotide content, and geographical distribution have been proven to influence the codon usage bias of the CDV H gene. The novel findings enhanced the understanding of CDV evolution.

Compositional Features and Codon Usage Pattern of Genes Associated with Anxiety in Human.

Codon usage bias (CUB) is the unequal usage of synonymous codon; some codons are more preferred than others. CUB analysis has applications in understanding the molecular organization of genome, genetics, gene expression, and molecular evolution. Bioinformatic approach was used to analyze the protein-coding sequences of genes involved in the anxiety to understand the patterns of codon usage as no work was reported yet. The improved effective number of codons (Nc) values ranged from 43.55 to 55.06, with a mean of 44.57, suggested that the overall CUB was low for genes associated with anxiety. The overall GC and AT content was 54.76 and 45.24, respectively. Relative synonymous codon usage (RSCU) analysis revealed that most frequently used codon ended mostly with C or G. The over-represented codons in genes associated with anxiety were CTG, ATC, GTG, AGC, ACC, and GCC, while under-represented codons were TTA, CTT, CTA, ATA, GTT, GTA, TCG, CCG, GCG, CAA, and CGT. Correlation analysis was performed between overall nucleotide composition and its 3rd codon positions, and observed highly significant (p < 0.01) correlation between them suggested that both mutation pressure and natural selection might affect the pattern of CUB. The highly significant correlation (0.598**, p < 0.01) was also observed between GC12 with GC3 suggested that directional mutation pressure might acted on all codon positions for genes associated with anxiety.

Codon usage trend in genes associated with obesity.

Obesity is not only a social menace but also an economic burden as it reduces productivity and increases health care cost. We used bioinformatic tools to analyze the CUB of obesity associated genes and compared with housekeeping genes (control) to explore the similarities and differences between two data sets as no work was reported yet. The mean effective number of codons (ENC) in genes associated with obesity and housekeeping gene was 50.45 and 52.03 respectively, indicating low CUB. The relative synonymous codon usage (RSCU) suggested that codons namely CTG and GTG were over-represented in both obesity and housekeeping genes while under-represented codons were TCG, TTA, CTA, CCG, CAA, CGT, ATA, ACG, GTA and GCG in obesity genes and TCG, TTA, CCG, ATA, ACG, GTA, and GCG in housekeeping genes. t test analysis suggested that 11 codons namely TTA (Leu), TTG (Leu), CCG (Pro), CAC (His), CAA (Gln), CAG (Gln), CGT (Arg), AGA (Arg), ATA (Ile), ATT (Ile) and GCG (Ala) were significantly differed (p < 0.05 or p < 0.01) between obesity and housekeeping genes. Highly significant correlation was observed between GC12 and GC3 in obesity and housekeeping genes i.e. r = 0.580** and r = 0.498** (p < 0.01) respectively indicating the effect of directional mutation pressure present in all codon positions.

Analysis of codon usage in the mitochondrion genome of Spirodela polyrhiza

In this study, investigation into the characteristics of codon bias and factors in shaping the codon usage pattern of the mitochondrial genome of Spirodela polyrhiza L. Schleid (strain 7498) was performed. S. polyrhiza had a very high GC content of 45.7% in the mitochondrial genome. Compared with other monocots whose mitochondrial genomes have been sequenced, the GC contents of mitochondrial genome of monocots seem to have a decreasing trend as the molecular evolution, and then finally trend to stable at about 44%. The mitochondrion of S. polyrhiz had extremely narrow GC3s distributions from 0.212 to 0.433, if taking no account of SpMatR (0.559) which was the only protein-coding genes with G/C-biased codon usages. The neutrality plot (GC12 vs. GC3s) indicated that there was no correlation between GC12 and GC3s, suggesting there was little directional mutation pressure on the mitochondrial genome of S. polyrhiza. The ENC-plot (ENC vs. GC3s) showed that there were only a small number of genes lying on the expected curve, but a majority of the points with low ENC values lied well below the expected curve. According to the frequency distributions of the ENC ratio, the dataset displayed a single peak, and the most of genes located into a narrow region of ENC ratios between 0–0.1. Context-dependent mutation was analyzed by Chi square test for contingency tables at the four fold degenerate (FFD) sites, suggesting that the frequencies of N3-position bases were definitely not independent of the N2-position bases. Correspondence analysis (COA) of relative synonymous codon usage (RSCU) yielded a first axis which accounted for only13.23% of the total variation. These results suggested that although other factors might play a role in the mitochondrion evolution of S. polyrhiza, mutation press still was the major factor in shaping the codon usage. Meantime, six codons including UUA(Leu), CAU(His), CGU(Arg), AUU(Ile), GCU(Ala) and GAA(Glu) were determined as the “optimal codons” (p < 0.01). Among them, the codon UUA(Leu), its usage was significantly more often in highly biased genes than in lowly.

Chronic obstructive pulmonary disease: A crosstalk on nucleotide compositional dynamics and codon usage patterns of the genes involved in disease.

Chronic obstructive pulmonary disease (COPD), a lung disease, affects a large number of people worldwide, leading to death. Here, we analyzed the compositional features and trends of codon usage of the genes influencing COPD to understand molecular biology, genetics, and evolutionary relationships of these genes as no work was reported yet. Coding sequences of COPD genes were found to be rich in guanine-cytosine (GC) content. A high value (34-60) of the effective number of codons of the genes indicated low codon usage bias (CUB). Correspondence analysis suggested that the COPD genes were distinct in their codon usage patterns. Relative synonymous codon usage values of codons differed between the more preferred codons and the less-preferred ones. Correlation analysis between overall nucleotides and those at third codon position revealed that mutation pressure might influence the CUB of the genes. The high correlation between GC12 and GC3 signified that directional mutation pressure might have operated at all the three codon positions in COPD genes.

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.

Analysis of base and codon usage by rubella virus

Rubella virus (RUBV), a small, plus-strand RNA virus that is an important human pathogen, has the unique feature that the GC content of its genome (70%) is the highest (by 20%) among RNA viruses. To determine the effect of this GC content on genomic evolution, base and codon usage were analyzed across viruses from eight diverse genotypes of RUBV. Despite differences in frequency of codon use, the favored codons in the RUBV genome matched those in the human genome for 18 of the 20 amino acids, indicating adaptation to the host. Although usage patterns were conserved in corresponding genes in the diverse genotypes, within-genome comparison revealed that both base and codon usages varied regionally, particularly in the hypervariable region (HVR) of the P150 replicase gene. While directional mutation pressure was predominant in determining base and codon usage within most of the genome (with the strongest tendency being towards C's at third codon positions), natural selection was predominant in the HVR region. The GC content of this region was the highest in the genome (>80%), and it was not clear if selection at the nucleotide level accompanied selection at the amino acid level. Dinucleotide frequency analysis of the RUBV genome revealed that TpA usage was lower than expected, similar to mammalian genes; however, CpG usage was not suppressed, and TpG usage was not enhanced, as is the case in mammalian genes.

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.

Rossiter Henry Crozier 1943 - 2009

Ross Crozier, population geneticist and leader in the study of the evolutionary genetics of social insects, was born on 4 January 1943 in Jodhpur, India. He died of a heart attack in his office at James Cook University in Townsville on 12 November 2009. He is survived by his wife Yuen Ching Kok, who was his inseparable companion and collaborator in life as in the laboratory. Crozier was a pioneer in the application of molecular genetic markers to the analysis of social insect populations, and generated much of the theory that made these analyses possible. Ross and Ching Crozier produced the first sequence of the honey bee mitochondrial genome?the second insect mitochondria to be fully sequenced. From the sequence Crozier produced fundamental insights into the nature of DNA evolution, particularly directional mutation pressure towards particular nucleotides. Crozier contributed massively to the development of kin selection theory, which remains the most potent explanatory theory for the evolution of social behaviour in insects.

Adaptive Threonine Increase in Transmembrane Regions of Mitochondrial Proteins in Higher Primates

BackgroundThe mitochondrial (mt) gene tree of placental mammals reveals a very strong acceleration of the amino acid (AA) replacement rate and a change in AA compositional bias in the lineage leading to the higher primates (simians), in contrast to the nuclear gene tree. Whether this acceleration and compositional bias were caused by adaptive evolution at the AA level or directional mutation pressure at the DNA level has been vigorously debated.Methodology/Principal FindingsOur phylogenetic analysis indicates that the rate acceleration in the simian lineage is accompanied by a marked increase in threonine (Thr) residues in the transmembrane helix regions of mt DNA-encoded proteins. This Thr increase involved the replacement of hydrophobic AAs in the membrane interior. Even after accounting for lack of independence due to phylogeny, a regression analysis reveals a statistical significant positive correlation between Thr composition and longevity in primates.Conclusion/SignificanceBecause crucial roles of Thr and Ser in membrane proteins have been proposed to be the formation of hydrogen bonds enhancing helix-helix interactions, the Thr increase detected in the higher primates might be adaptive by serving to reinforce stability of mt proteins in the inner membrane. The correlation between Thr composition in the membrane interior and the longevity of animals is striking, especially because some mt functions are thought to be involved in aging.

Revisiting the directional mutation pressure theory: The analysis of a particular genomic structure in Leishmania major

This paper reports the existence of a significant negative correlation between GC12 and GC3 in the recently sequenced genome of Leishmania major. This result contradicts the previous evidence that the compositional correlations between codon positions are universal. Moreover, it challenges the interpretation of the GC12 vs. GC3 linear regression slope as the relative neutrality of GC12, within the framework of the directional mutation pressure theory [Sueoka, N., 1988. Directional mutational pressure and neutral molecular evolution. P Natl Acad Sci USA 85, 2653–2657.]. The analysis of the codon usage pattern for L. major shows that codon choice is most likely influenced by both mutation pressure and translational selection. Dinucleotide frequencies were also analysed; our results do not support the existence of an unusual neighbour-dependent mutation bias in this genome. We developed two evolutionary models that could explain the origin of the negative GC12/GC3 correlation. The first model is based on the effect of translational selection on the GC3 the second one is based on a potential mutation bias combined with purifying selection at the amino-acid level. Both models predict a negative GC12/GC3 correlation at the equilibrium. The potential implications of these results for this aspect of the directional mutation pressure theory are discussed. We conclude that the particular case of L. major should lead to a careful reevaluation of several hypotheses of this theory. The origin of the negative GC12/GC3 correlation remains for now an open question.

Insect mitochondrial genomics 2: the complete mitochondrial genome sequence of a giant stonefly,Pteronarcys princeps, asymmetric directional mutation bias, and conserved plecopteran A+T-region elements

Mitochondrial (mt) genome sequences of insects are receiving renewed attention in molecular phylogentic studies, studies of mt-genome rearrangement, and other unusual molecular phenomena, such as translational frameshifting. At present, the basal neopteran lineages are poorly represented by mt-genome sequences. Complete mt-genome sequences are available in the databases for only the Orthoptera and Blatteria; 9 orders are unrepresented. Here, we present the complete mt-genome sequence of a giant stonefly, Pteronarcys princeps (Plecoptera; Pteronarcyidae). The 16,004 bp genome is typical in its genome content, gene organisation, and nucleotide composition. The genome shows evidence of strand-specific mutational biases, correlated with the time between the initiation of leading and the initiation of lagging strand replication. Comparisons with other insects reveal that this trend is seen in other insect groups, but is not universally consistent among sampled mt-genomes. The A+T region is compared with that of 2 stoneflies in the family Peltoperlidae. Conserved stem-loop structures and sequence blocks are noted between these distantly related families.

Retroposed Elements as Archives for the Evolutionary History of Placental Mammals

Reconstruction of the placental mammalian (eutherian) evolutionary tree has undergone diverse revisions, and numerous aspects remain hotly debated. Initial hierarchical divisions based on morphology contained many misgroupings due to features that evolved independently by similar selection processes. Molecular analyses corrected many of these misgroupings and the superordinal hierarchy of placental mammals was recently assembled into four clades. However, long or rapid evolutionary periods, as well as directional mutation pressure, can produce molecular homoplasies, similar characteristics lacking common ancestors. Retroposed elements, by contrast, integrate randomly into genomes with negligible probabilities of the same element integrating independently into orthologous positions in different species. Thus, presence/absence analyses of these elements are a superior strategy for molecular systematics. By computationally scanning more than 160,000 chromosomal loci and judiciously selecting from only phylogenetically informative retroposons for experimental high-throughput PCR applications, we recovered 28 clear, independent monophyly markers that conclusively verify the earliest divergences in placental mammalian evolution. Using tests that take into account ancestral polymorphisms, multiple long interspersed elements and long terminal repeat element insertions provide highly significant evidence for the monophyletic clades Boreotheria (synonymous with Boreoeutheria), Supraprimates (synonymous with Euarchontoglires), and Laurasiatheria. More importantly, two retropositions provide new support for a prior scenario of early mammalian evolution that places the basal placental divergence between Xenarthra and Epitheria, the latter comprising all remaining placentals. Due to its virtually homoplasy-free nature, the analysis of retroposon presence/absence patterns avoids the pitfalls of other molecular methodologies and provides a rapid, unequivocal means for revealing the evolutionary history of organisms.

The Response of Amino Acid Frequencies to Directional Mutation Pressure in Mitochondrial Genome Sequences Is Related to the Physical Properties of the Amino Acids and to the Structure of the Genetic Code

The frequencies of A, C, G, and T in mitochondrial DNA vary among species due to unequal rates of mutation between the bases. The frequencies of bases at fourfold degenerate sites respond directly to mutation pressure. At first and second positions, selection reduces the degree of frequency variation. Using a simple evolutionary model, we show that first position sites are less constrained by selection than second position sites and, therefore, that the frequencies of bases at first position are more responsive to mutation pressure than those at second position. We define a measure of distance between amino acids that is dependent on eight measured physical properties and a similarity measure that is the inverse of this distance. Columns 1, 2, 3, and 4 of the genetic code correspond to codons with U, C, A, and G in their second position, respectively. The similarity of amino acids in the four columns decreases systematically from column 1 to column 2 to column 3 to column 4. We then show that the responsiveness of first position bases to mutation pressure is dependent on the second position base and follows the same decreasing trend through the four columns. Again, this shows the correlation between physical properties and responsiveness. We determine a proximity measure for each amino acid, which is the average similarity between an amino acid and all others that are accessible via single point mutations in the mitochondrial genetic code structure. We also define a responsiveness for each amino acid, which measures how rapidly an amino acid frequency changes as a result of mutation pressure acting on the base frequencies. We show that there is a strong correlation between responsiveness and proximity, and that both these quantities are also correlated with the mutability of amino acids estimated from the mtREV substitution rate matrix. We also consider the variation of base frequencies between strands and between genes on a strand. These trends are consistent with the patterns expected from analysis of the variation among genomes.

Codon usage patterns in Nematoda: analysis based on over 25 million codons in thirty-two species

A codon usage table for 32 nematode species is presented and suggests that total genomic GC content drives codon usage.

Comparative Analysis of the Base Composition and Codon Usages in Fourteen Mycobacteriophage Genomes

To study the possible codon usage and base composition variation in the bacteriophages, fourteen mycobacteriophages were used as a model system here and both the parameters in all these phages and their plating bacteria, M. smegmatis had been determined and compared. As all the organisms are GC-rich, the GC contents at third codon positions were found in fact higher than the second codon positions as well as the first + second codon positions in all the organisms indicating that directional mutational pressure is strongly operative at the synonymous third codon positions. Nc plot indicates that codon usage variation in all these organisms are governed by the forces other than compositional constraints. Correspondence analysis suggests that: (i) there are codon usage variation among the genes and genomes of the fourteen mycobacteriophages and M. smegmatis, i.e, codon usage patterns in the mycobacteriophages is phage-specific but not the M. smegmatis-specific; (ii) synonymous codon usage patterns of Barnyard, Che8, Che9d, and Omega are more similar than the rest mycobacteriophages and M. smegmatis; (iii) codon usage bias in the mycobacteriophages are mainly determined by mutational pressure; and (iv) the genes of comparatively GC rich genomes are more biased than the GC poor genomes. Translational selection in determining the codon usage variation in highly expressed genes can be invoked from the predominant occurrences of C ending codons in the highly expressed genes. Cluster analysis based on codon usage data also shows that there are two distinct branches for the fourteen mycobacteriophages and there is codon usage variation even among the phages of each branch.