Natural Selection Plays an Important Role in Shaping the Codon Usage of Structural Genes of the Viruses Belonging to the Coronaviridae Family

Exploring Codon Usage Patterns and Influencing Factors in Ranavirus DNA Polymerase Genes

Ranaviruses, members of the genus Ranavirus within the family Iridoviridae, have become a significant concern for amphibian populations globally, along with other cold-blooded vertebrates, due to their emergence as a significant threat. We employed bioinformatics tools to examine the codon usage patterns in 61 DNA pol genes from Ranavirus, Lymphocystivirus, Megalocytivirus, and two unclassified ranaviruses, as no prior studies had been conducted on this topic. The results showed a slight or low level of codon usage bias (CUB) in the DNA pol genes of Ranavirus. Relative synonymous codon usage (RSCU) analysis indicated that the predominant codons favored in Ranavirus DNA pol genes terminate with C or G. Correlation analysis examining nucleotide content, third codon position, effective number of codons (ENC), correspondence analysis (COA), Aroma values, and GRAVY values indicated that the CUB across DNA pol genes could be influenced by both mutation pressure and natural selection. The neutrality plot indicated that natural selection is the primary factor driving codon usage. Furthermore, the analysis of the codon adaptation index (CAI) illustrated the robust adaptability of Ranavirus DNA pol genes to their hosts. Analysis of the relative codon deoptimization index (RCDI) suggested that Ranavirus DNA pol genes underwent greater selection pressure from their hosts. These findings will aid in comprehending the factors influencing the evolution and adaptation of Ranavirus to its hosts.

Differences between phytophagous and predatory species in Pentatomidae based on the mitochondrial genome

Pentatomidae includes many species of significant economic value as plant pests and biological control agents. The feeding habits of Pentatomidae are closely related to their energy metabolism and ecological adaptations. In this study, we sequenced the mitochondrial genomes of 12 Asopinae species using the next-generation sequencing to explore the effect of dietary changes on mitochondrial genome evolution. Notably, all sequences were double-stranded circular DNA molecules containing 37 genes and one control region. We then compared and analyzed the mitochondrial genome characteristics of phytophagous and predatory bugs. Notably, no significant difference was observed in the length of the mitochondrial genomes between the predatory and phytophagous bugs. However, the AT content was higher in the mitochondrial genomes of phytophagous bugs than that of predatory bugs. Moreover, phytophagous bugs prefer codon usage patterns ending in A/T compared with predatory bugs. The evolution rate of predatory bugs was lower than that of phytophagous bugs. The phylogenetic relationships across phytophagous bugs' lineages were largely consistent at depth nodes based on different datasets and tree-reconstructing methods, and strongly supported the monophyly of predatory bugs. Additionally, the estimated divergence times indicated that Pentatomidae explosively radiated in the Early Cretaceous. Subsequently, the subfamily Asopinae and the genus Menida diverged in the Late Cretaceous. Our research results provide data supporting for the evolutionary patterns and classification of Pentatomidae.

Comparative genome wise analysis of codon usage of Staphylococcus Genus

The genus Staphylococcus encompasses a diverse array of bacteria with significant implications for human health, including disreputable pathogens such as Staphylococcus aureus and Staphylococcus epidermidis. Understanding the genetic composition and codon usage patterns of Staphylococcus species is crucial for unraveling their evolutionary dynamics, adaptive strategies, and pathogenic potential. In this study, we conducted a comprehensive analysis of codon usage patterns across 48 species within the Staphylococcus genus. Our findings uncovered variations in genomic G-C content across Staphylococcus species, impacting codon usage preferences, with a notable preference for A/T-rich codons observed in pathogenic strains. This preference for A/T-rich codons suggests an energy-saving strategy in pathogenic organisms. Analysis of dinucleotide pair expression patterns unveiled insights into genomic dynamics, with overrepresented codon pairs reflecting trends in dinucleotide expression across genomes. Additionally, a significant correlation between CAI and genomic G-C content underscored the intricate relationship between codon usage patterns and gene expression strategies. Amino acid usage analysis highlighted preferences for energetically cheaper amino acids, suggesting adaptive strategies promoting energy efficiency. This comprehensive analysis sheds light on the evolutionary dynamics and adaptive mechanisms employed by Staphylococcus species, providing valuable insights into their pathogenic potential and clinical implications. Understanding these genomic features is crucial for devising strategies to combat staphylococcal infections and improve public health outcomes.

Analysis of codon usage bias of exonuclease genes in invertebrate iridescent viruses

Invertebrate iridescent viruses (IIVs) are double-stranded DNA viruses that belong to the Iridoviridae family. IIVs result diseases that vary in severity from subclinical to lethal in invertebrate hosts. Codon usage bias (CUB) analysis is a versatile method for comprehending the genetic and evolutionary aspects of species. In this study, we analyzed the CUB in 10 invertebrate iridescent viruses exonuclease genes by calculating and comparing the nucleotide contents, effective number of codons (ENC), codon adaptation index (CAI), relative synonymous codon usage (RSCU), and others. The results revealed that IIVs exonuclease genes are rich in A/T. The ENC analysis displayed a low codon usage bias in IIVs exonuclease genes. ENC-plot, neutrality plot, and parity rule 2 plot demonstrated that besides mutational pressure, other factors like natural selection, dinucleotide content, and aromaticity also contributed to CUB. The findings could enhance our understanding of the evolution of IIVs exonuclease genes.

Codon usage bias of secretory protein in Fusarium oxysporum f. sp. cubense tropical race 4

Banana Fusarium oxysporum f. sp. cubense tropical race 4 (Foc-TR4) is a highly destructive pathogen that infects nearly all major banana cultivars and has a tendency to spread further. Secreted proteins play a crucial role in the process of Fusarium wilt infection in bananas. In this study, we analyzed the codon usage bias (CUB) of the Foc-TR4 classical secretory protein genome for the first time and observed a strong bias toward codons ending with C. We found that 572 out of the 14,543 amino acid sequences in the Foc-TR4 genome exhibited characteristics of classical secretory proteins. The CUB was largely influenced by selection optimization pressure, as indicated by the ENC value and neutral plot analysis. Among the identified codons, such as UCC and CCC, 11 were found to be optimal for Foc-TR4 gene expression. Codons with higher GC content and a C base in the third position showed greater selectivity. The CUB in the secretory proteins encoded by Foc-TR4 provides insights into their evolutionary patterns, contributing to the development and screening of novel and effective antifungal drugs.

Reexamination of Aspergillus cristatus phylogeny in dark tea: Characteristics of the mitochondrial genome

To enhance our understanding of Aspergillus cristatus, an important functional microorganism, the characteristics of its mitochondrial genome were analyzed and compared with related species. The mitochondrial genome of A. cristatus was determined to be 77,649 bp in length, with 15 protein-coding regions. Notably, its length surpassed that of the other species, primarily attributable to the intron length. Gene order exhibited significant variations, with greater conservation observed in the genus Penicillium compared to Aspergillus. Phylogenetic tree analyses indicated that the genera Aspergillus and Penicillium are closely related but monophyletic. Furthermore, the phylogenetic tree constructed based on protein-coding genes effectively distinguished all strains with high branching confidence. This approach provides a robust reflection of the evolutionary relationship between A. cristatus and its related species, offering potential for the development of molecular markers suitable for Aspergillus and Penicillium.

Optimization and Deoptimization of Codons in SARS-CoV-2 and Related Implications for Vaccine Development

The spread of coronavirus disease 2019 (COVID-19), caused by severe respiratory syndrome coronavirus 2 (SARS-CoV-2), has progressed into a global pandemic. To date, thousands of genetic variants have been identified among SARS-CoV-2 isolates collected from patients. Sequence analysis reveals that the codon adaptation index (CAI) values of viral sequences have decreased over time but with occasional fluctuations. Through evolution modeling, it is found that this phenomenon may result from the virus's mutation preference during transmission. Using dual-luciferase assays, it is further discovered that the deoptimization of codons in the viral sequence may weaken protein expression during virus evolution, indicating that codon usage may play an important role in virus fitness. Finally, given the importance of codon usage in protein expression and particularly for mRNA vaccines, it is designed several codon-optimized Omicron BA.2.12.1, BA.4/5, and XBB.1.5 spike mRNA vaccine candidates and experimentally validated their high levels of expression. This study highlights the importance of codon usage in virus evolution and provides guidelines for codon optimization in mRNA and DNA vaccine development.

Natural selection plays a significant role in governing the codon usage bias in the novel SARS-CoV-2 variants of concern (VOC)

The ongoing prevailing COVID-19 pandemic caused by SARS-CoV-2 is becoming one of the major global health concerns worldwide. The SARS-CoV-2 genome encodes spike (S) glycoprotein that plays a very crucial role in viral entry into the host cell via binding of its receptor binding domain (RBD) to the host angiotensin converting enzyme 2 (ACE2) receptor. The continuously evolving SARS-CoV-2 genome results in more severe and transmissible variants characterized by the emergence of novel mutations called 'variants of concern' (VOC). The currently designated alpha, beta, gamma, delta and omicron VOC are the focus of this study due to their high transmissibility, increased virulence, and concerns for decreased effectiveness of the available vaccines. In VOC, the spike (S) gene and other non-structural protein mutations may affect the efficacies of the approved COVID-19 vaccines. To understand the diversity of SARS-CoV-2, several studies have been performed on a limited number of sequences. However, only a few studies have focused on codon usage bias (CUBs) pattern analysis of all the VOC strains. Therefore, to evaluate the evolutionary divergence of all VOC S-genes, we performed CUBs analysis on 300,354 sequences to understand the evolutionary relationship with its adaptation in different hosts, i.e., humans, bats, and pangolins. Base composition and RSCU analysis revealed the presence of 20 preferred AU-ended and 10 under-preferred GC-ended codons. In addition, CpG was found to be depleted, which may be attributable to the adaptive response by viruses to escape from the host defense process. Moreover, the ENC values revealed a higher bias in codon usage in the VOC S-gene. Further, the neutrality plot analysis demonstrated that S-genes analyzed in this study are under 83.93% influence of natural selection, suggesting its pivotal role in shaping the CUBs. The CUBs pattern of S-genes was found to be very similar among all the VOC strains. Interestingly, we observed that VOC strains followed a trend of antagonistic codon usage with respect to the human host. The identified CUBs divergence would help to understand the virus evolution and its host adaptation, thus help design novel vaccine strategies against the emerging VOC strains. To the best of our knowledge, this is the first report for identifying the evolution of CUBs pattern in all the currently identified VOC.

Insight into Codon Utilization Pattern of Tumor Suppressor Gene EPB41L3 from Different Mammalian Species Indicates Dominant Role of Selection Force

Simple Summary

The present study envisaged the codon usage pattern analysis of tumor suppressor gene EPB41L3 for the human, brown rat, domesticated cattle, and Sumatran orangutan. Most amino acids are coded by more than one synonymous codon, but they are used in a biased manner. The codon usage bias results from multiple factors like compositional properties, dinucleotide abundance, neutrality, parity, tRNA pool, etc. Understanding codon bias is central to fields as diverse as molecular evolution, gene expressivity, protein translation, and protein folding. This kind of studies is important to see the effects of various evolutionary forces on codon usage. The present study indicated that the selection force is dominant over other forces shaping codon usage in the envisaged organisms.

Abstract

Uneven codon usage within genes as well as among genomes is a usual phenomenon across organisms. It plays a significant role in the translational efficiency and evolution of a particular gene. EPB41L3 is a tumor suppressor protein-coding gene, and in the present study, the pattern of codon usage was envisaged. The full-length sequences of the EPB41L3 gene for the human, brown rat, domesticated cattle, and Sumatran orangutan available at the NCBI were retrieved and utilized to analyze CUB patterns across the selected mammalian species. Compositional properties, dinucleotide abundance, and parity analysis showed the dominance of A and G whilst RSCU analysis indicated the dominance of G/C-ending codons. The neutrality plot plotted between GC12 and GC3 to determine the variation between the mutation pressure and natural selection indicated the dominance of selection pressure (R = 0.926; p < 0.00001) over the three codon positions across the gene. The result is in concordance with the codon adaptation index analysis and the ENc-GC3 plot analysis, as well as the translational selection index (P2). Overall selection pressure is the dominant pressure acting during the evolution of the EPB41L3 gene.

Mutational Asymmetries in the SARS-CoV-2 Genome May Lead to Increased Hydrophobicity of Virus Proteins

The genomic diversity of SARS-CoV-2 has been a focus during the ongoing COVID-19 pandemic. Here, we analyzed the distribution and character of emerging mutations in a data set comprising more than 95,000 virus genomes covering eight major SARS-CoV-2 lineages in the GISAID database, including genotypes arising during COVID-19 therapy. Globally, the C>U transitions and G>U transversions were the most represented mutations, accounting for the majority of single-nucleotide variations. Mutational spectra were not influenced by the time the virus had been circulating in its host or medical treatment. At the amino acid level, we observed about a 2-fold excess of substitutions in favor of hydrophobic amino acids over the reverse. However, most mutations constituting variants of interests of the S-protein (spike) lead to hydrophilic amino acids, counteracting the global trend. The C>U and G>U substitutions altered codons towards increased amino acid hydrophobicity values in more than 80% of cases. The bias is explained by the existing differences in the codon composition for amino acids bearing contrasting biochemical properties. Mutation asymmetries apparently influence the biochemical features of SARS CoV-2 proteins, which may impact protein-protein interactions, fusion of viral and cellular membranes, and virion assembly.