Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic.

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

Strategies and Patterns of Codon Bias in Molluscum Contagiosum Virus

Trends associated with codon usage in molluscum contagiosum virus (MCV) and factors governing the evolution of codon usage have not been investigated so far. In this study, attempts were made to decipher the codon usage trends and discover the major evolutionary forces that influence the patterns of codon usage in MCV with special reference to sub-types 1 and 2, MCV-1 and MCV-2, respectively. Three hypotheses were tested: (1) codon usage patterns of MCV-1 and MCV-2 are identical; (2) SCUB (synonymous codon usage bias) patterns of MCV-1 and MCV-2 slightly deviate from that of human host to avoid affecting the fitness of host; and (3) translational selection predominantly shapes the SCUB of MCV-1 and MCV-2. Various codon usage indices viz. relative codon usage value, effective number of codons and codon adaptation index were calculated to infer the nature of codon usage. Correspondence analysis and correlation analysis were performed to assess the relative contribution of silent base contents and significance of codon usage indices in defining bias in codon usage. Among the tested hypotheses, only the second and third hypotheses were accepted.

Molecular Mimicry of Pathogenicity of Neisseria

Neisseria, a genus from beta-proteobacteria class, is of potent clinical importance. This genus contains both pathogenic and commensal strains. Gonorrhea and meningitis are two major diseases caused by pathogens belonging to this genus. With increased use of antimicrobial agents against these pathogens they have evolved the antimicrobial resistance (AMR) capacity making these diseases nearly untreatable. The set of anti-bacterial resistance genes (resistome) and genes associated with signal processing (secretomes) are crucial for the host-microbial interaction. With the virtue of whole genome sequences and computational biology it is now possible to study the genomic and proteomic riddles of Neisseria along with their comprehensive evolutionary and metabolic profiling. We have studied relative synonymous codon usage, amino acid usage, reverse ecology, comparative genomics, evolutionary analysis and pathogen-host (Neisseria-human) interaction through bioinformatics analysis. Our analysis revealed the co-evolution of Neisseria genomes with the human host. Moreover, co-occurrence of Neisseria and humans has been supported through reverse ecology analysis. A differential pattern of evolutionary rate of resistomes and secretomes was evident among the pathogenic and commensal strains. Comparative genomics supported the presence of virulent genes in both pathogenic and commensal strains of select genus. Our analysis also indicated a transition from commensal to pathogenic Neisseria strains through the long run of evolution.

Host Adaptation of Codon Usage in SARS-CoV-2 From Mammals Indicate Natural Selection

The outbreak of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections, spread across hosts from humans to animals, transmitting particularly effectively in mink. How SARS-CoV-2 selects and evolves in the host, and the differences in the evolution of different animals are still unclear. To analysis the mutation and codon usage bias of SARS-CoV-2 in infected humans and animals. The SARS-CoV-2 sequence in mink (Mink-SARS2) and binding energy with receptor were calculated compared with human. The relative synonymous codon usage of viral encoded gene was analyzed to characterize the differences and the evolutionary characteristics. A synonymous codon usage analysis showed that SARS-CoV-2 is optimized to adapt in the animals in which it is currently reported, and all of the animals showed decreased adaptability relative to that of humans, except for mink. The neutrality plot showed that the effect of natural selection on different SARS-CoV-2 sequences is stronger than mutation pressure. A binding affinity analysis indicated that the spike protein of the SARS-CoV-2 variant in mink showed a greater preference for binding with the mink receptor ACE2 than with the human receptor, especially as the mutation Y453F and N501T in Mink-SARS2 lead to improvement of binding affinity for mink receptor. In summary, mutations Y453F and N501T in Mink-SARS2 lead to improvement of binding affinity with mink receptor, indicating possible natural selection and current host adaptation. Monitoring the variation and codon bias of SARS-CoV-2 provides a theoretical basis for tracing the epidemic, evolution and cross-species spread of SARS-CoV-2.

Complete chloroplast genome of Hordeum brevisubulatum: Genome organization, synonymous codon usage, phylogenetic relationships, and comparative structure analysis

Background Hordeum brevisubulatum, known as fine perennial forage, is used for soil salinity improvement in northern China. Chloroplast (cp) genome is an ideal model for assessing its genome evolution and the phylogenetic relationships. We de novo sequenced and analyzed the cp genome of H. brevisubulatum, providing a fundamental reference for further studies in genetics and molecular breeding. Results The cp genome of H. brevisubulatum was 137,155 bp in length with a typical quadripartite structure. A total of 130 functional genes were annotated and the gene of accD was lost in the process of evolution. Among all the annotated genes, 16 different genes harbored introns and the genes of ycf3 and rps12 contained two introns. Parity rule 2 (PR2) plot analysis showed that majority of genes had a bias toward T over A in the coding strand in all five Hordeum species, and a slight G over C in the other four Hordeum species except for H. bogdanil. Additionally, 52 dispersed repeat sequences and 182 simple sequence repeats were identified. Moreover, some unique SSRs of each species could be used as molecular markers for further study. Compared to the other four Hordeum species, H. brevisubulatum was most closely related to H. bogdanii and its cp genome was relatively conserved. Moreover, inverted repeat regions (IRa and IRb) were less divergent than other parts and coding regions were relatively conserved compared to non-coding regions. Main divergence was presented at the SSC/IR border. Conclusions This research comprehensively describes the architecture of the H. brevisubulatum cp genome and improves our understanding of its cp biology and genetic diversity, which will facilitate biological discoveries and cp genome engineering.

Comparative analysis of codon usage patterns of FUT2 from different species

Enterotoxigenic E. coli is an important zoonotic pathogen causing diarrhea in human and newborn animals. α - (1,2) fucosyltransferase 2 (FUT2) is closely associated with the formation of pathogenic receptors of Enterotoxigenic E. coli. Codon usage bias analysis can help to better understand the molecular mechanisms and evolutionary relationships of a particular gene. In order to understand the codon usage pattern of FUT2 gene, FUT2 gene coding sequences of nine species were selected from GenBank database for calculating the nucleotide composition (GC content) and genetic indices including effective number of codons, relative synonymous codon usage and relative codon usage bias using R software, in order to analyze codon usage bias and base composition in FUT2 gene from different species. The results showed that the codon usage of FUT2 gene in different species was affected by GC bias, especially GC frequency at the third position of codon (GC3). Most of the optimal codons were biased towards the G/C-ending types. GCC, CUG, UCC, GUG and AUC showed the highest relative synonymous codon usage value among different species, belonging to the most dominant codons. The usage characteristic of the codens for FUT2 gene in Sus scrofa was similar to that of Bos taurus; Homo sapiens was similar to Pan troglodytes. Effective number of codons was significantly, negatively correlated with GC3, and the relative higher frequency of optimal codon implied that FUT2 genes from different species had a strong bias in codon usage.

The First Eight Mitogenomes of Leaf-Mining Dactylispa Beetles (Coleoptera: Chrysomelidae: Cassidinae) Shed New Light on Subgenus Relationships

The taxonomic classification of Dactylispa, a large genus of leaf-mining beetles, is problematic because it is currently based on morphology alone. Here, the first eight mitochondrial genomes of Dactylispa species, which were used to construct the first molecular phylogenies of this genus, are reported. The lengths of the eight mitogenomes range from 17,189 bp to 20,363 bp. All of the mitochondrial genomes include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 A + T-rich region. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively, which is a common phenomenon among animals. According to relative synonymous codon usage, UUA(L) has the highest frequency. With two Gonophorini species as the outgroup, mitogenome-based phylogenetic trees of the eight Dactylispa species were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the PCGs, tRNAs, and rRNAs. Two DNA-based phylogenomic inferences and one protein-based phylogenomic inference support the delimitation of the subgenera Dactylispa s. str. and Platypriella as proposed in the system of Chen et al. (1986). However, the subgenus Triplispa is not recovered as monophyletic. The placement of Triplispa species requires further verification and testing with more species. We also found that both adult body shape and host plant relationship might explain the subgeneric relationships among Dactylispa beetles to a certain degree.

Codon Usage Bias in Autophagy-Related Gene 13 in Eukaryotes: Uncovering the Genetic Divergence by the Interplay Between Nucleotides and Codon Usages

Synonymous codon usage bias is a universal characteristic of genomes across various organisms. Autophagy-related gene 13 (atg13) is one essential gene for autophagy initiation, yet the evolutionary trends of the atg13 gene at the usages of nucleotide and synonymous codon remains unexplored. According to phylogenetic analyses for the atg13 gene of 226 eukaryotic organisms at the nucleotide and amino acid levels, it is clear that their nucleotide usages exhibit more genetic information than their amino acid usages. Specifically, the overall nucleotide usage bias quantified by information entropy reflected that the usage biases at the first and second codon positions were stronger than those at the third position of the atg13 genes. Furthermore, the bias level of nucleotide ‘G’ usage is highest, while that of nucleotide ‘C’ usage is lowest in the atg13 genes. On top of that, genetic features represented by synonymous codon usage exhibits a species-specific pattern on the evolution of the atg13 genes to some extent. Interestingly, the codon usages of atg13 genes in the ancestor animals (Latimeria chalumnae, Petromyzon marinus, and Rhinatrema bivittatum) are strongly influenced by mutation pressure from nucleotide composition constraint. However, the distributions of nucleotide composition at different codon positions in the atg13 gene display that natural selection still dominates atg13 codon usages during organisms’ evolution.

Intragenomic variation in mutation biases causes underestimation of selection on synonymous codon usage

Patterns of non-uniform usage of synonymous codons (codon bias) varies across genes in an organism and across species from all domains of life. The bias in codon usage is due to a combination of both non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most population genetics models quantify the effects of mutation bias and selection on shaping codon usage patterns assuming a uniform mutation bias across the genome. However, mutation biases can vary both along and across chromosomes due to processes such as biased gene conversion, potentially obfuscating signals of translational selection. Moreover, estimates of variation in genomic mutation biases are often lacking for non-model organisms. Here, we combine an unsupervised learning method with a population genetics model of synonymous codon bias evolution to assess the impact of intragenomic variation in mutation bias on the strength and direction of natural selection on synonymous codon usage across 49 Saccharomycotina budding yeasts. We find that in the absence of a priori information, unsupervised learning approaches can be used to identify regions evolving under different mutation biases. We find that the impact of intragenomic variation in mutation bias varies widely, even among closely-related species. We show that the overall strength and direction of selection on codon usage can be underestimated by failing to account for intragenomic variation in mutation biases. Interestingly, genes falling into clusters identified by machine learning are also often physically clustered across chromosomes, consistent with processes such as biased gene conversion. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable mutation biases on codon frequencies.