Characterization of Synonymous Codon Usage Bias in the Pseudorabies Virus US1 Gene

In the present study, we examined the codon usage bias between pseudorabies virus （PRV） US1 gene and the USl-like genes of 20 reference alphaherpesviruses. Comparative analysis showed noticeable disparities of the synonymous codon usage bias in the 21 alphaherpesviruses, indicated by codon adaptation index, effective number of codons （ENc） and GC3s value. The codon usage pattern of PRV US1 gene was phylogenetically conserved and similar to that of the USl-like genes of the genus Varicellovirus of alphaherpesvirus, with a strong bias towards the codons with C and G at the third codon position. Cluster analysis of codon usage pattern of PRV US1 gene with its reference alphaherpesviruses demonstrated that the codon usage bias of USl-like genes of 21 alphaherpesviruses had a very close relation with their gene functions. ENc-plot revealed that the genetic heterogeneity in PRV US1 gene and the 20 reference alphaherpesviruses was constrained by G＋C content, as well as the gene length. In addition, comparison of codon preferences in the US1 gene of PRV with those ofE. coli, yeast and human revealed that there were 50 codons showing distinct usage differences between PRV and yeast, 49 between PRV and human, but 48 between PRV and E. coil Although there were slightly fewer differences in codon usages between E.coli and PRV, the difference is unlikely to be statistically significant, and experimental studies are necessary to establish the most suitable expression system for PRV US1. In conclusion, these results may improve our understanding of the evolution, pathogenesis and functional studies of PRV, as well as contributing to the area of herpesvirus research or even studies with other viruses.

Analysis of Synonymous Codon Usage Bias in 09H1N1

A novel subtype of influenza A virus 09H1N1 has rapidly spread across the world. Evolutionary analyses of this virus have revealed that 09H1N1 is a triple reassortant of segments from swine, avian and human influenza viruses. In this study, we investigated factors shaping the codon usage bias of 09H1N1 and carried out cluster analysis of 60 strains of influenza A virus from different subtypes based on their codon usage bias. We discovered that more preferentially used codons of 09H1N1 are A-ended or U-ended, and the intra-genomic codon usage bias of 09H1N1 is quite low. Base composition constraint, dinucleotide biases and translational selection are the main factors influencing the codon usage bias of 09H1N1. At the genome level, we find that the codon usage bias of 09H1N1 is similar to H1N1 (A/swine/Kansas/77778/2007H1N1), H9N2 from Asia, H1N2 from Asia and North America and H3N2 from North America. Our results provide insight for understanding the processes governing evolution, regulation of gene expression, and revealing the evolution of 09H1N1.

Epidemiological Surveillance: Genetic Diversity of Rotavirus Group A in the Pearl River Delta, Guangdong, China in 2019

Objective This study aimed to understand the epidemic status and phylogenetic relationships of rotavirus group A(RVA)in the Pearl River Delta region of Guangdong Province,China.Methods This study included individuals aged 28 days–85 years.A total of 706 stool samples from patients with acute gastroenteritis collected between January 2019 and January 2020 were analyzed for 17 causative pathogens,including RVA,using a Gastrointestinal Pathogen Panel,followed by genotyping,virus isolation,and complete sequencing to assess the genetic diversity of RVA.Results The overall RVA infection rate was 14.59%(103/706),with an irregular epidemiological pattern.The proportion of co-infection with RVA and other pathogens was 39.81%(41/103).Acute gastroenteritis is highly prevalent in young children aged 0–1 year,and RVA is the key pathogen circulating in patients 6–10 months of age with diarrhea.G9P[8](58.25%,60/103)was found to be the predominant genotype in the RVA strains,and the 41 RVA-positive strains that were successfully sequenced belonged to three different RVA genotypes in the phylogenetic analysis.Recombination analysis showed that gene reassortment events,selection pressure,codon usage bias,gene polymorphism,and post-translational modifications(PTMs)occurred in the G9P[8]and G3P[8]strains.Conclusion This study provides molecular evidence of RVA prevalence in the Pearl River Delta region of China,further enriching the existing information on its genetics and evolutionary characteristics and suggesting the emergence of genetic diversity.Strengthening the surveillance of genotypic changes and gene reassortment in RVA strains is essential for further research and a better understanding of strain variations for further vaccine development.

Analysis of Synonymous Codon Usage in Aeropyrum pernix K1 and Other Crenarchaeota Microorganisms

In this study, a comparative analysis of the codon usage bias was performed in Aeropyrum pernix K1 and two other phylogenetically related Crenarchaeota microorganisms （i.e., Pyrobaculum aerophilum str. IM2 and Sulfolobus acidocaldarius DSM 639）. The results indicated that the synonymous codon usage in A. pernix K1 was less biased, which was highly correlated with the GC3s value. The codon usage patterns were phylogenetically conserved among these Crenarchaeota microorganisms. Comparatively, it is the species function rather than the gene function that determines their gene codon usage patterns. A. pernix K1, P. aerophilum str. IM2, and S. acidocaldarius DSM 639 live in differently extreme conditions. It is presumed that the hving environment played an important role in determining the codon usage pattern of these microorganisms. Besides, there was no strain-specific codon usage among these microorganisms. The extent of codon bias in A. pernix K1 and S. acidocaldarius DSM 639 were highly correlated with the gene expression level, but no such association was detected in P. aerophilum str. IM2 genomes.

Synonymous codon usage in Methanosarcina mazei str.Goe1 and other Euryarchaeota microorganisms

A comparative analysis of the codon usage bias was conducted in Methanosarcina mazei str. Goel and two related Euryarchaeota microorganisms （Picrophilus torridus str. DSM 9790 and Natronomonas pharaonis str. DSM 2160）. Results revealed that synonymous codon usage in Methanosarcina mazei str. Goel was less biased, which was highly correlated with the GC3S value. And the codon usage patterns were phylogenetically conserved among those Euryarchaeota microorganisms. By employing a hierarchical clustering analysis, it can be seen that it is more the species than the gene function that determines their gene codon usage pattems. Considering that those microorganisms live in different environments where the pH conditions vary quite a lot, it can be presumed that their living environments, especially the pH conditions, play an important role in determining those microorganisms＇ codon usage pattems.

Molecular Evolution of Protein Sequences and Codon Usage in Monkeypox Viruses

The monkeypox virus(mpox virus,MPXV)epidemic in 2022 has posed a significant public health risk.Yet,the evolutionary principles of MPXV remain largely unknown.Here,we examined the evolutionary patterns of protein sequences and codon usage in MPXV.We first demonstrated the signal of positive selection in OPG027,specifically in the CladeⅠlineage of MPXV.Subsequently,we discovered accelerated protein sequence evolution over time in the variants responsible for the 2022 outbreak.Furthermore,we showed strong epistasis between amino acid substitutions located in different genes.The codon adaptation index(CAI)analysis revealed that MPXV genes tended to use more non-preferred codons compared to human genes,and the CAI decreased over time and diverged between clades,with CladeⅠ>Ⅱa andⅡb-A>Ⅱb-B.While the decrease in fatality rate among the three groups aligned with the CAI pattern,it remains unclear whether this correlation was coincidental or if the deoptimization of codon usage in MPXV led to a reduction in fatality rates.This study sheds new light on the mechanisms that govern the evolution of MPXV in human populations.

Comparative analysis of flower-meristem-identity gene APETALA2 (AP2) codon in different plant species

The flower-meristem-identity gene APETALA2(AP2), one of class-A genes, is involved in the establishment of the floral meristem and the forming of sepals and petals. Codon usage bias(CUB) identifies differences among species, meanwhile dynamic analysis of base composition can identify the molecular mechanisms and evolutionary relationships of a specific gene. In this study, eight coding sequences(CDS) of AP2 gene were selected from different plant species using the Gen Bank database. Their nucleotide composition(GC content), genetic index, relative synonymous codon usage(RSCU) and relative codon usage bias(RCUB) were calculated with R Software to compare codon bias and base composition dynamics of AP2 gene codon usage patterns in different plant species. The results showed that the usage of AP2 gene codons from different plant species were influened by GC bias, especially GC3 s. Overall, base composition analysis indicated that the usage frequency of codon AT in the gene coding sequence was higher than GC among AP2 gene CDS from different plant species. Furthermore, most AP2 gene CDSs ended with AT; AGA, GCU and UGU had relatively high RSCU values as the most dominant codon; the usage characteristic of the AP2 gene codon in Malus domestica was similar to that of Vitis vinifera; Paeonia lactiflora was similar to Paeonia suffruticosa and Solanum lycopersicum was similar to Petunia×hybrida. There was a moderate preference in the usage of AP2 gene codon among different plant species from relatively low frequency of optimal codon(Fop) values and high effective number of codons(ENC) value. This study has revealed the usage characteristics of the AP2 gene codon from the comparision of AP2 gene codon preference and base dynamics in different plant species and provide a platform for further study towards transgenic engineering and codon optimization.

Genome-wide Identification and Analysis of DNA Methyltransferases in Grape

DNA methylation is an important epigenetic regulation mechanism, which is catalyzed by DNA methyltransferases. In this study, eight DNA methyltransferase genes were identified in grape genome to analyze the selective pressure, gene expression and codon usage bias. The results showed grape DNA methyltransferase MET subfamily underwent relatively strong purifying selection during evolution, while chromomethylase CMT subfamily underwent positive selection during evolution. Under different abiotic(heat, drought or cold) stresses, the expression level of many grape DNA methyltransferase genes changed significantly. The expression level of these genes might be related with cis-regulatory elements of their promoters. The results of codon usage bias analysis showed that synonymous codon bias existed in grape DNA methyltransferase gene family, which might be affected by mutation pressure. These results laid a solid foundation for in-depth study of DNA methyltransferases in grape.

Lack of transcriptional coordination between mitochondrial and nuclear oxidative phosphorylation genes in the presence of two divergent mitochondrial genomes

In most eukaryotes,oxidative phosphorylation(OXPHOS)is the main energy production process and it involves both mitochondrial and nuclear genomes.The close interaction between the two genomes is critical for the coordinated function of the OXPHOS process.Some bivalves show doubly uniparental inheritance(DUI)of mitochondria,where two highly divergent mitochondrial genomes,one inherited through eggs(F-type)and the other through sperm(M-type),coexist in the same individual.However,it remains a puzzle how nuclear OXPHOS genes coordinate with two divergent mitochondrial genomes in DUI species.In this study,we compared transcription,polymorphism,and synonymous codon usage in the mitochondrial and nuclear OXPHOS genes of the DUI species Ruditapes philippinarum using sex-and tissue-specific transcriptomes.Mitochondrial and nuclear OXPHOS genes showed different transcription profiles.Strong co-transcription signal was observed within mitochondrial(separate for F-and M-type)and within nuclear OXPHOS genes but the signal was weak or absent between mitochondrial and nuclear OXPHOS genes,suggesting that the coordination between mitochondrial and nuclear OXPHOS subunits is not achieved transcriptionally.McDonald-Kreitman and frequency-spectrum based tests indicated that M-type OXPHOS genes deviated significantly from neutrality,and that F-type and M-type OXPHOS genes undergo different selection patterns.Codon usage analysis revealed that mutation bias and translational selection were the major factors affecting the codon usage bias in different OXPHOS genes,nevertheless,translational selection in mitochondrial OXPHOS genes appears to be less efficient than nuclear OXPHOS genes.Therefore,we speculate that the coordination between OXPHOS genes may involve post-transcriptional/translational regulation.

Preliminary analysis of the mitochondrial genome evolutionary pattern in primates

Since the birth of molecular evolutionary analysis,primates have been a central focus of study and mitochondrial DNA is well suited to these endeavors because of its unique features.Surprisingly,to date no comprehensive evaluation of the nucleotide substitution patterns has been conducted on the mitochondrial genome of primates.Here,we analyzed the evolutionary patterns and evaluated selection and recombination in the mitochondrial genomes of 44 Primates species downloaded from GenBank.The results revealed that a strong rate heterogeneity occurred among sites and genes in all comparisons.Likewise,an obvious decline in primate nucleotide diversity was noted in the subunit rRNAs and tRNAs as compared to the protein-coding genes.Within 13 protein-coding genes,the pattern of nonsynonymous divergence was similar to that of overall nucleotide divergence,while synonymous changes differed only for individual genes,indicating that the rate heterogeneity may result from the rate of change at nonsynonymous sites.Codon usage analysis revealed that there was intermediate codon usage bias in primate protein-coding genes,and supported the idea that GC mutation pressure might determine codon usage and that positive selection is not the driving force for the codon usage bias.Neutrality tests using site-specific positive selection from a Bayesian framework indicated no sites were under positive selection for any gene,consistent with near neutrality.Recombination tests based on the pairwise homoplasy test statistic supported complete linkage even for much older divergent primate species.Thus,with the exception of rate heterogeneity among mitochondrial genes,evaluating the validity assumed complete linkage and selective neutrality in primates prior to phylogenetic or phylogeographic analysis seems unnecessary.

Preliminary analysis of the mitochondrial genome evolutionary pattern in primates

Since the birth of molecular evolutionary analysis, primates have been a central focus of study and mitochondrial DNA is well suited to these endeavors because of its unique features. Surprisingly, to date no comprehensive evaluation of the nucleotide substitution patterns has been conducted on the mitochondrial genome of primates. Here, we analyzed the evolutionary patterns and evaluated selection and recombination in the mitochondrial genomes of 44 Primates species downloaded from Genl3ank. The results revealed that a strong rate heterogeneity occurred among sites and genes in all comparisons. Likewise, an obvious decline in primate nucleotide diversity was noted in the subunit rRNAs and tRNAs as compared to the protein-coding genes. Within 13 protein-coding genes, the pattern of nonsynonymous divergence was similar to that of overall nucleotide divergence, while synonymous changes differed only for individual genes, indicating that the rate heterogeneity may result from the rate of change at nonsynonymous sites. Codon usage analysis revealed that there was intermediate codon usage bias in primate protein-coding genes, and supported the idea that GC mutation pressure might determine codon usage and that positive selection is not the driving force for the codon usage bias. Neutrality tests using site-specific positive selection from a Bayesian framework indicated no sites were under positive selection for any gene, consistent with near neutrality. Recombination tests based on the pairwise homoplasy test statistic supported complete linkage even for much older divergent primate species. Thus, with the exception of rate heterogeneity among mitochondrial genes, evaluating the validity assumed complete linkage and selective neutrality in primates prior to phylogenetic or phylogeographic analysis seems unnecessary.

Phylogenetic Relationships Among Species Subgroups in the Drosophila saltans Group (Diptera: Drosophilidae): Can Morphology Solve a Molecular Conflict?

Proper phylogenetic reconstruction is crucial for understanding many evolutionary phenomena. In spite of the great success of molecular phylogenetics, DNA signal still may be limited by some intrinsic constraints such as codon usage bias. The phylogenetic relationships between the five species subgroups of the Drosophila saltans group are a good example of conflicting molecular phylogenies drawn from different genes due to an ancestral substitutional shift. Here, forty morphological characters were analyzed using the same set of species used in previous molecular studies, with at least a single representative of each subgroup. The cladistic analysis was in disagreement with most of the previous hypotheses in placing the sturtevanti subgroup as an early branch, whereas the four remaining subgroups form a well supported clade that can be further subdivided into two sister clades： one containing the cordata and the elliptica subgroups, whereas the second includes the parasaltans and the saltans subgroups. The molecular evolution （codon usage bias） of the saltans group were revised in light of the present finding. The analysis highlights the important role of morphology in phylogeny reconstruction and in understanding molecular evolutionary phenomena.

Effect of the Flanking Sequence Architecture of Translation Initiation AUG Codon on Gene Expression Level in Rice

The relationship between the codon usage bias, gene expression level and the AUG context(from -20 to +6 positions relative to the initiator AUG codon) was examined in 541unigene sequences of rice. A significant correlation for CAI values (codon adaptationindex) was observed at five nucleotide positions (-19, -18, -9, -4, +5), eight (-19, -18,-14, -9, -6, -4, -1, +5) for CPP (codon preference parameter), and seven (-18, -16, -15,-9, -7, -1, +6) for mRNA abundance in the flanking sequence of the initiator AUG codonrespectively, but a significantly positive correlation for both CAI and CPP at twopositions (-4 and +5), indicating that both those positions are evolutionally under thenatural selection constraint at the translational level. By site-directed mutagenesis atseven specific positions (-18, -16, -15, -9, -7, -1 and +6) for allergenic protein thathad the highest mRNA abundance in this study, its expression level decreased dramatically63.3 and 72.5% respectively, indicating the importance of those 7 positions for geneexpression. A highly positive correlation (r=0.625, P<0.01) between AUGCAI and GCcontent in the flanking sequence of the initiator AUG codon showed a more effectivehigher GC content on translation initiation efficiency. The strong preference for G orC at those 8 positions (-6, -5, -3, -2, -1, +4, +5 and +6) in the AUG context suggestedthat an important factor in modulation of the translation efficiency, as well assynonymous codon usage bias, particularly in highly expressed genes.

Dynamic Landscapes of tRNA Transcriptomes and Translatomes in Diverse Mouse Tissues

Although the function of tRNAs in the translational process is well established,it remains controversial whether tRNA abundance is tightly associated with translational efficiency(TE)in mammals.Moreover,how critically the expression of tRNAs contributes to the establishment of tissue-specific proteomes in mammals has not been well addressed.Here,we measured both tRNA expression using demethylase-tRNA sequencing(DM-tRNA-seq)and TE of mRNAs using ribosome-tagging sequencing(RiboTag-seq)in the brain,heart,and testis of mice.Remarkable variation in the expression of tRNA isodecoders was observed among different tissues.When the statistical effect of isodecoder-grouping on reducing variations is considered through permutating the anticodons,we observed an expected reduction in the variation of anticodon expression across all samples,an unexpected smaller variation of anticodon usage bias,and an unexpected larger variation of tRNA isotype expression at amino acid level.Regardless of whether or not they share the same anticodons,the isodecoders encoding the same amino acids are co-expressed across different tissues.Based on the expression of tRNAs and the TE of mRNAs,we find that the tRNA adaptation index(tAI)and TE are significantly correlated in the same tissues but not between tissues;and tRNA expression and the amino acid composition of translating peptides are positively correlated in the same tissues but not between tissues.We therefore hypothesize that the tissue-specific expression of tRNAs might be due to post-transcriptional mechanisms.This study provides a resource for tRNA and translation studies,as well as novel insights into the dynamics of tRNAs and their roles in translational regulation.

Functional genomic analysis of Hawaii marine metagenomes

Using high-throughput sequencing on metagenome to analyze marine microbial community, it is one of current main issues in the field of environmental microbe research. In this paper, we conducted the functional analysis on seven samples of metagenomic data from different depth seawater in Hawaii. The results of gene prediction and function annotation indicate that there are large amounts of potential novel genes of which functions remain unknown at present. Based on the gene annotation, codon usage bias is studied on ribosomal protein-related genes and shows an evident influence by the marine extreme environment. Furthermore, focusing on the marine environmental differences such as light intensity, dissolved oxygen, temperature and pressure among various depths, comparative analysis is carried out on related genes and metabolic pathways. Thus, the understanding as well as new insights into the correlation between marine environment and microbes are proposed at molecular level. Therefore, the studies herein afford a clue to reveal the special living strategies of microbial community from sea surface to deep sea.