Factors influencing synonymous codon and amino acid usage biases in Mimivirus

Multi-Chromosomal mitochondrial genome of medicinal plant Acorus tatarinowii (Acoraceae): Firstly reported from Acorales Order

Acorus tatarinowii Schott (A. tatarinowii), a well-known traditional Chinese medicinal plant renowned for its high medicinal value, but its mitochondrial genome (mitogenome) is still unexplored. In this study, we meticulously assembled the complete mitochondrial genome of A. tatarinowii using a combination of Illumina short reads and Oxford Nanopore long reads. Our findings revealed that A. tatarinowii possesses a complex chromosomal structural mitogenome, comprising two linear chromosomes and seven circular chromosomes. This mitogenome spans 1.81 Mb in length with a GC content of 38.29 %. Notably, it contained 24 unique mitochondrial core genes, seven unique variable genes, 17 tRNA genes, and three rRNA genes. Analyses of codon usage, most protein-coding genes (PCGs) exhibited a common codon usage preference, with RSCU values greater than 1, and the codon with the highest RSCU value was UAA(End, 1.90). We conducted a thorough analysis of repeat sequences, the distribution of repetitive sequences in nine mitochondrial chromosomes showed distinct patterns. Moreover, we identified 82 and 12 homologous fragments by comparing the sequences of chloroplast and nuclear genomes to the A. tatarinowii mitogenome, respectively. Lastly, We predicted a total of 234 potential RNA editing sites in 28 unique PCGs and discovered that the nad4 gene has been edited the most often, at 26 times. Our results contribute to the enrichment of mitochondrial genome resources for Acoraceae, and the mitogenome also can be used as a reference for other species.

Synonymous sites for accessibility around microRNA binding sites in bacterial spot and speck disease resistance genes of tomato

The major causes of mass tomato infections in both covered and open ground are agents of bacterial spot and bacterial speck diseases. MicroRNAs (miRNAs) are 16-21 nucleotides in length, non-coding RNAs that inhibit translation and trigger mRNA degradation. MiRNAs play a significant part in plant resistance to abiotic and biotic stresses by mediating gene regulation via post-transcriptional RNA silencing. In this study, we analyzed a collection of bacterial resistance genes of tomato and their binding sites for tomato miRNAs and Pseudomonas syringe pv. tomato miRNAs. Our study found that two genes, bacterial spot disease resistance gene (Bs4) and bacterial speck disease resistance gene (Prf), have a 7mer-m8 perfect seed match with miRNAs. Bs4 was targeted by one tomato miRNA (sly-miR9470-3p) and three Pseudomonas syringe pv. tomato miRNAs (PSTJ4_3p_27246, PSTJ4_3p_27246, and PSTJ4_3p_27246). Again, Prf gene was found to be targeted by two tomato miRNAs namely, sly-miR9469-5p and sly-miR9474-3p. The accessibility of the miRNA-target site and its flanking regions and the relationship between relative synonymous codon usage and tRNAs were compared. Strong access to miRNA targeting regions and decreased rate of translations suggested that miRNAs might be efficient in binding to their particular targets. We also found the existence of rare codons, which suggests that it could enhance miRNA targeting even more. The codon usage pattern analysis of the two genes revealed that both were AT-rich (Bs4 = 63.2%; Prf = 60.8%). We found a low codon usage bias in both genes, suggesting that selective restriction might regulate them. The silencing property of miRNAs would allow researchers to discover the involvement of plant miRNAs in pathogen invasion. However, the efficient validation of direct targets of miRNAs is an urgent need that might be highly beneficial in enhancing plant resistance to multiple pathogenic diseases.

Assembly of the Complete Mitochondrial Genome of Pereskia aculeata Revealed That Two Pairs of Repetitive Elements Mediated the Recombination of the Genome

Pereskia aculeata is a potential new crop species that has both food and medicinal (antinociceptive activity) properties. However, comprehensive genomic research on P. aculeata is still lacking, particularly concerning its organelle genome. In this study, P. aculeata was studied to sequence the mitochondrial genome (mitogenome) and to ascertain the assembly, informational content, and developmental expression of the mitogenome. The findings revealed that the mitogenome of P. aculeata is circular and measures 515,187 bp in length with a GC content of 44.05%. It contains 52 unique genes, including 33 protein-coding genes, 19 tRNA genes, and three rRNA genes. Additionally, the mitogenome analysis identified 165 SSRs, primarily consisting of tetra-nucleotides, and 421 pairs of dispersed repeats with lengths greater than or equal to 30, which were mainly forward repeats. Based on long reads and PCR experiments, we confirmed that two pairs of long-fragment repetitive elements were highly involved with the mitogenome recombination process. Furthermore, there were 38 homologous fragments detected between the mitogenome and chloroplast genome, and the longest fragment was 3962 bp. This is the first report on the mitogenome in the family Cactaceae. The decoding of the mitogenome of P. aculeata will provide important genetic materials for phylogenetic studies of Cactaceae and promote the utilization of species germplasm resources.

Genes Common in Primary Immunodeficiencies and Cancer Display Overrepresentation of Codon CTG and Dominant Role of Selection Pressure in Shaping Codon Usage

Primary immunodeficiencies (PIDs) are disorders of the immune system that involve faulty cellular, humoral, or both cellular and humoral functions. PIDs are present at the crossroad between infections, immune dysregulation, and cancers. A panel encompassing 42 genes involved in both PIDs and cancer has been investigated for the genes’ compositional properties, codon usage patterns, various forces affecting codon choice, protein properties, and gene expression profiles. In the present study, the codon choice of genes was found to be dependent upon the richness of the nucleotide; the viz AT nucleotide rich genome preferred AT ending codons. The dinucleotide TpA adversely affected protein expression, while CpG did not. The CTG codon was the most overrepresented codon in 80.95% of genes. Analysis of various protein properties, including GRAVY, AROMA, isoelectric point, aliphatic index, hydrophobicity, instability index, and numbers of acidic, basic, and neutral amino acid residues revealed that the hydrophobicity index, instability index, and numbers of acidic and basic amino acid residues are the factors affecting gene expression. Based on neutrality analysis, parity analysis, ENc-GC3 analysis, and regression analysis of nucleotides present at the first and third positions of the codon, it was determined that selection pressure, mutation pressure, and compositional constraints all participated in shaping codon usage. The study will help determine the various evolutionary forces acting on genes common to both PIDs and cancer. Codon usage analysis might be helpful in the future to augment both diseases simultaneously. The research also indicates a peculiar pattern adapted by a set of genes involved in any disease.

Insights into molecular structure, genome evolution and phylogenetic implication through mitochondrial genome sequence of Gleditsia sinensis

Gleditsia sinensis is an endemic species widely distributed in China with high economic and medicinal value. To explore the genomic evolution and phylogenetic relationships of G. sinensis, the complete mitochondrial (mt) genome of G. sinensis was sequenced and assembled, which was firstly reported in Gleditsia. The mt genome was circular and 594,121 bp in length, including 37 protein-coding genes (PCGs), 19 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. The overall base composition of the G. sinensis mt genome was 27.4% for A, 27.4% for T, 22.6% for G, 22.7% for C. The comparative analysis of PCGs in Fabaceae species showed that most of the ribosomal protein genes and succinate dehydrogenase genes were lost. In addition, we found that the rps4 gene was only lost in G. sinensis, whereas it was retained in other Fabaceae species. The phylogenetic analysis based on shared PCGs of 24 species (22 Fabaceae and 2 Solanaceae) showed that G. sinensis is evolutionarily closer to Senna species. In general, this research will provide valuable information for the evolution of G. sinensis and provide insight into the phylogenetic relationships within the family Fabaceae.

Altering Compositional Properties of Viral Genomes to Design Live-Attenuated Vaccines

Live-attenuated vaccines have been historically used to successfully prevent numerous diseases caused by a broad variety of RNA viruses due to their ability to elicit strong and perdurable immune-protective responses. In recent years, various strategies have been explored to achieve viral attenuation by rational genetic design rather than using classic and empirical approaches, based on successive passages in cell culture. A deeper understanding of evolutionary implications of distinct viral genomic compositional aspects, as well as substantial advances in synthetic biology technologies, have provided a framework to achieve new viral attenuation strategies. Herein, we will discuss different approaches that are currently applied to modify compositional features of viruses in order to develop novel live-attenuated vaccines.

Analysis of codon usage patterns of the chloroplast genome in Delphinium grandiflorum L. reveals a preference for AT-ending codons as a result of major selection constraints

Background

Codon usage bias analysis is a suitable strategy for identifying the principal evolutionary driving forces in different organisms. Delphinium grandiflorum L. is a perennial herb with high economic value and typical biological characteristics. Evolutionary analysis of D. grandiflorum can provide a rich resource of genetic information for developing hybridization resources of the genus Delphinium.

Methods

Synonymous codon usage (SCU) and related indices of 51 coding sequences from the D. grandiflorum chloroplast (cp) genome were calculated using Codon W, Cups of EMBOSS, SPSS and Microsoft Excel. Multivariate statistical analysis combined by principal component analysis (PCA), correspondence analysis (COA), PR2-plot mapping analysis and ENC plot analysis was then conducted to explore the factors affecting the usage of synonymous codons.

Results

The SCU bias of D. grandiflorum was weak and codons preferred A/T ending. A SCU imbalance between A/T and G/C at the third base position was revealed by PR2-plot mapping analysis. A total of eight codons were identified as the optimal codons. The PCA and COA results indicated that base composition (GC content, GC₃ content) and gene expression were important for SCU bias. A majority of genes were distributed below the expected curve from the ENC plot analysis and up the standard curve by neutrality plot analysis. Our results showed that with the exception of notable mutation pressure effects, the majority of genetic evolution in the D. grandiflorum cp genome might be driven by natural selection.

Discussions

Our results provide a theoretical foundation for elucidating the genetic architecture and mechanisms of D. grandiflorum, and contribute to enriching D. grandiflorum genetic resources.

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.

Codon usage pattern and its influencing factors in different genomes of hepadnaviruses

Codon usage bias (CUB) arises from the preference for a codon over codons for the same amino acid. The major factors contributing to CUB are evolutionary forces, compositional properties, gene expression, and protein properties. The present analysis was performed to investigate the compositional properties and the extent of CUB across the genomes of members of the family Hepadnaviridae, as previously no work using bioinformatic tools has been reported. The viral genes were found to be AT rich with low CUB. Analysis of relative synonymous codon usage (RSCU) was used to identify overrepresented and underrepresented codons for each amino acid. Correlation analysis of overall nucleotide composition and its composition at the third codon position suggested that mutation pressure might influence the CUB. A highly significant correlation was observed between GC12 and GC3 (r = 0.910, p < 0.01), indicating that directional mutation affected all three codon positions across the genome. Translational selection (P2) and mutational responsive index (MRI) values of genes suggested that mutation plays a more important role than translational selection in members of the family Hepadnaviridae.

Close evolutionary relationship between rice black-streaked dwarf virus and southern rice black-streaked dwarf virus based on analysis of their bicistronic RNAs

Background

Rice black-streaked dwarf virus (RBSDV) and Southern rice black-streaked dwarf virus (SRBSDV) seriously interfered in the production of rice and maize in China. These two viruses are members of the genus Fijivirus in the family Reoviridae and can cause similar dwarf symptoms in rice. Although some studies have reported the phylogenetic analysis on RBSDV or SRBSDV, the evolutionary relationship between these viruses is scarce.

Methods

In this study, we analyzed the evolutionary relationships between RBSDV and SRBSDV based on the data from the analysis of codon usage, RNA recombination and phylogenetic relationship, selection pressure and genetic characteristics of the bicistronic RNAs (S5, S7 and S9).

Results

RBSDV and SRBSDV showed similar patterns of codon preference: open reading frames (ORFs) in S7 and S5 had with higher and lower codon usage bias, respectively. Some isolates from RBSDV and SRBSDV formed a clade in the phylogenetic tree of S7 and S9. In addition, some recombination events in S9 occurred between RBSDV and SRBSDV.

Conclusions

Our results suggest close evolutionary relationships between RBSDV and SRBSDV. Selection pressure, gene flow, and neutrality tests also supported the evolutionary relationships.

Electronic supplementary material

The online version of this article (10.1186/s12985-019-1163-3) contains supplementary material, which is available to authorized users.

Differences in Molecular Characteristics of Segment 8 in Rice black-streaked dwarf virus and Southern rice black-streaked dwarf virus

Rice black-streaked dwarf virus (RBSDV) and Southern rice black-streaked dwarf virus (SRBSDV) cause maize rough dwarf disease (MRDD) and rice black-streaked dwarf disease (RBSDD) in China. RBSDV segment 8 (S8) contains the only deletion mutation in the genomes of these viruses, which are both members of the genus Fijivirus. To illuminate the molecular differences between the RBSDV and SRBSDV genomes and better understand the evolution of these viruses, and to determine which virus is specifically associated with MRDD and RBSDD in each region, S8 was analyzed in 66 virus isolates collected from 10 geographic locations in China and 14 S8 sequences obtained from the National Center for Biotechnology Information GenBank. Phylogenetic analysis showed that the pathogen associated with MRDD and RBSDD in the Yellow and Huai River valleys was RBSDV, whereas the pathogen associated with these diseases in Sanya was SRBSDV. Codon usage bias in S8 differed significantly between RBSDV and SRBSDV, as indicated by effective number of codons used by a gene (Nc) and GC values, Nc plots, and variation explained by the first axis in correspondence analysis. The nucleotide identities among these 66 RBSDV and SRBSDV isolates ranged from 66.2 to 68.2%, and were considerably lower than the nucleotide identities within RBSDV (from 94.1 to 99.9%) or SRBSDV (from 93.9 to 100%) isolates. Most S8 polymorphisms were identified in the region from 1,000 to 1,200 bp in RBSDV and in the region from 500 to 700 bp in SRBSDV. The difference in the lengths of RBSDV (1,936 bp) and SRBSDV (1,928 bp) was due to an 8-bp deletion in the 3'-untranslated region of SRBSDV. Six recombination events were detected in S8 in RBSDV and two recombination events were detected in S8 in SRBSDV. Recombination breakpoints were found within the region containing the deletion mutation in nine isolates. However, no recombination events were detected between RBSDV and SRBSDV. Both of these viruses were under negative and purifying selection, although the ratio of nonsynonymous mutations to synonymous mutations (Ka/Ks) for RBSDV S8 (0.0530) was not significantly lower than that of SRBSDV S8 (0.0823, P = 0.1550). We found that SRBSDV was more highly genetically differentiated (product of effective population size and the migration rate among populations < 1; values for the among-populations component of genetic variation or normalized variation > 0.33; and P values of the sequence statistic, the rank statistic, and the nearest-neighbor statistic < 0.01) than RBSDV. However, gene flow between RBSDV and SRBSDV was not frequent.

Evolutionary and genetic analysis of the VP2 gene of canine parvovirus

BACKGROUND

Canine parvovirus (CPV) type 2 emerged in 1978 in the USA and quickly spread among dog populations all over the world with high morbidity. Although CPV is a DNA virus, its genomic substitution rate is similar to some RNA viruses. Therefore, it is important to trace the evolution of CPV to monitor the appearance of mutations that might affect vaccine effectiveness.

RESULTS

Our analysis shows that the VP2 genes of CPV isolated from 1979 to 2016 are divided into six groups: GI, GII, GIII, GIV, GV, and GVI. Amino acid mutation analysis revealed several undiscovered important mutation sites: F267Y, Y324I, and T440A. Of note, the evolutionary rate of the CPV VP2 gene from Asia and Europe decreased. Codon usage analysis showed that the VP2 gene of CPV exhibits high bias with an ENC ranging from 34.93 to 36.7. Furthermore, we demonstrate that natural selection plays a major role compared to mutation pressure driving CPV evolution.

CONCLUSIONS

There are few studies on the codon usage of CPV. Here, we comprehensively studied the genetic evolution, codon usage pattern, and evolutionary characterization of the VP2 gene of CPV. The novel findings revealing the evolutionary process of CPV will greatly serve future CPV research.

Nucleotide composition and codon usage bias of SRY gene

The SRY gene is present within the sex-determining region of the Y chromosome which is responsible for maleness in mammals. The nonuniform usage of synonymous codons in the mRNA transcript for encoding a particular amino acid is the codon usage bias (CUB). Analysis of codon usage pattern is important to understand the genetic and molecular organisation of a gene. It also helps in heterologous gene expression, design of primer and synthetic gene. However, the analysis of codon usage bias of SRY gene was not yet studied. We have used bioinformatic tools to analyse codon usage bias of SRY gene across mammals. Codon bias index (CBI) indicated that the overall extent of codon usage bias was weak. The relative synonymous codon usage (RSCU) analysis suggested that most frequently used codons had an A or C at the third codon position. Compositional constraint played an important role in codon usage pattern as evident from correspondence analysis (CA). Significant correlation among nucleotides constraints indicated that both mutation pressure and natural selection affect the codon usage pattern. Neutrality plot suggested that natural selection might play a major role, while mutation pressure might play a minor role in codon usage pattern in SRY gene in different species of mammals.

Analysis of codon usage patterns in Ginkgo biloba reveals codon usage tendency from A/U-ending to G/C-ending

As one of the most ancient tree species, the codon usage pattern analysis of Ginkgo biloba is a useful way to understand its evolutionary and genetic mechanisms. Several studies have been conducted on angiosperms, but seldom on gymnosperms. Based on RNA-Seq data of the G. biloba transcriptome, amount to 17,579 unigenes longer than 300 bp were selected and analyzed from 68,547 candidates. The codon usage pattern tended towards more frequently use of A/U-ending codons, which showed an obvious gradient progressing from gymnosperms to dicots to monocots. Meanwhile, analysis of high/low-expression unigenes revealed that high-expression unigenes tended to use G/C-ending codons together with more codon usage bias. Variation of unigenes with different functions suggested that unigenes involving in environment adaptation use G/C-ending codons more frequently with more usage bias, and these results were consistent with the conclusion that the formation of G. biloba codon usage bias was dominated by natural selection.

Molecular Genetic Analysis and Evolution of Segment 7 in Rice Black-Streaked Dwarf Virus in China

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease or rice black-streaked dwarf disease and can lead to severe yield losses in maize and rice. To analyse RBSDV evolution, codon usage bias and genetic structure were investigated in 111 maize and rice RBSDV isolates from eight geographic locations in 2013 and 2014. The linear dsRNA S7 is A+U rich, with overall codon usage biased toward codons ending with A (A3s, S7-1: 32.64%, S7-2: 29.95%) or U (U3s, S7-1: 44.18%, S7-2: 46.06%). Effective number of codons (Nc) values of 45.63 in S7-1 (the first open reading frame of S7) and 39.96 in S7-2 (the second open reading frame of S7) indicate low degrees of RBSDV-S7 codon usage bias, likely driven by mutational bias regardless of year, host, or geographical origin. Twelve optimal codons were detected in S7. The nucleotide diversity (π) of S7 sequences in 2013 isolates (0.0307) was significantly higher than in 2014 isolates (0.0244, P = 0.0226). The nucleotide diversity (π) of S7 sequences in isolates from Jinan (0.0391) was higher than that from the other seven locations (P < 0.01). Only one S7 recombinant was detected in Baoding. RBSDV isolates could be phylogenetically classified into two groups according to S7 sequences, and further classified into two subgroups. S7-1 and S7-2 were under negative and purifying selection, with respective Ka/Ks ratios of 0.0179 and 0.0537. These RBSDV populations were expanding (P < 0.01) as indicated by negative values for Tajima's D, Fu and Li's D, and Fu and Li's F. Genetic differentiation was detected in six RBSDV subpopulations (P < 0.05). Absolute Fst (0.0790) and Nm (65.12) between 2013 and 2014, absolute Fst (0.1720) and Nm (38.49) between maize and rice, and absolute Fst values of 0.0085-0.3069 and Nm values of 0.56-29.61 among these eight geographic locations revealed frequent gene flow between subpopulations. Gene flow between 2013 and 2014 was the most frequent.

System analysis of synonymous codon usage biases in archaeal virus genomes

Recent studies of geothermally heated aquatic ecosystems have found widely divergent viruses with unusual morphotypes. Archaeal viruses isolated from these hot habitats usually have double-stranded DNA genomes, linear or circular, and can infect members of the Archaea domain. In this study, the synonymous codon usage bias (SCUB) and dinucleotide composition in the available complete archaeal virus genome sequences have been investigated. It was found that there is a significant variation in SCUB among different Archaeal virus species, which is mainly determined by the base composition. The outcome of correspondence analysis (COA) and Spearman׳s rank correlation analysis shows that codon usage of selected archaeal virus genes depends mainly on GC richness of genome, and the gene׳s function, albeit with smaller effects, also contributes to codon usage in this virus. Furthermore, this investigation reveals that aromaticity of each protein is also critical in affecting SCUB of these viral genes although it was less important than that of the mutational bias. Especially, mutational pressure may influence SCUB in SIRV1, SIRV2, ARV1, AFV1, and PhiCh1 viruses, whereas translational selection could play a leading role in HRPV1׳s SCUB. These conclusions not only can offer an insight into the codon usage biases of archaeal virus and subsequently the possible relationship between archaeal viruses and their host, but also may help in understanding the evolution of archaeal viruses and their gene classification, and more helpful to explore the origin of life and the evolution of biology.

•

The SCUB of archaeal virus genes depends mainly on GC richness of genome.

•

The mutational pressure is the main factor that influences SCUB.

•

The aromaticity of each protein is also critical in affecting SCUB.

•

The translational selection could play a leading role in HRPV1׳s SCUB.

•

The mode is helpful to explore the origin of life and the evolution of biology.

Codon usage, amino acid usage, transfer RNA and amino-acyl-tRNA synthetases in Mimiviruses

Mimiviruses are giant viruses that infect phagocytic protists, including Acanthamoebae spp., which were discovered during the past decade. They are the current record holder among viruses for their large particle and genome sizes. One group is composed of three lineages, referred to as A, B and C, which include the vast majority of the Mimiviridae members. Cafeteria roenbergensis virus represents a second group, though the Mimiviridae family is still expanding. We analyzed the codon and amino acid usages in mimiviruses, as well as both the transfer RNA (tRNA) and amino acyl-tRNA synthetases. We confirmed that the codon and amino acid usages of these giant viruses are highly dissimilar to those in their amoebal host Acanthamoeba castellanii and are instead correlated with the high adenine and thymine (AT) content of Mimivirus genomes. We further describe that the set of tRNAs and amino acyl-tRNA synthetases in mimiviruses is globally not adapted to the codon and amino acid usages of these viruses. Notwithstanding, Leu(TAA)tRNA, present in several Mimivirus genomes and in multiple copies in some viral genomes, may complement the amoebal tRNA pool and may contribute to accommodate the viral AT-rich codons. In addition, we found that the genes most highly expressed at the beginning of the Mimivirus replicative cycle have a nucleotide content more adapted to the codon usage in A.castellanii.

Analysis of codon use features of stearoyl-acyl carrier protein desaturase gene in Camellia sinensis

The stearoyl-acyl carrier protein desaturase (SAD) gene widely exists in all kinds of plants. In this paper, the Camellia sinensis SAD gene (CsSAD) sequence was firstly analyzed by Codon W, CHIPS, and CUSP programs online, and then compared with genomes of the tea plant, other species and SAD genes from 11 plant species. The results show that the CsSAD gene and the selected 73 of C. sinensis genes have similar codon usage bias. The CsSAD gene has a bias toward the synonymous codons with A and T at the third codon position, the same as the 73 of C. sinensis genes. Compared with monocotyledons such as Triticum aestivum and Zea mays, the differences in codon usage frequency between the CsSAD gene and dicotyledons such as Arabidopsis thaliana and Nicotiana tobacum are less. Therefore, A. thaliana and N. tobacum expression systems may be more suitable for the expression of the CsSAD gene. The analysis result of SAD genes from 12 plant species also shows that most of the SAD genes are biased toward the synonymous codons with G and C at the third codon position. We believe that the codon usage bias analysis presented in this study will be essential for providing a theoretical basis for discussing the structure and function of the CsSAD gene.

Analysis of synonymous codon usage patterns in torque teno sus virus 1 (TTSuV1)

Torque teno sus virus 1 (TTSuV1) is a novel virus that has been found widely distributed in the swine population in recent years. Analysis of codon usage can reveal much about the molecular evolution of TTSuV1. In this study, synonymous codon usage patterns and the key determinants in the coding region of 29 available complete TTSuV1 genome sequences were examined. By calculating the nucleotide content and relative synonymous codon usage (RSCU) of TTSuV1 coding sequences, we found that the preferentially used codons were mostly those ending with A or C nucleotides; less-used codons were mostly codons ending with U or G nucleotides, and these were mainly affected by composition constraints. Although there was a variation in codon usage bias among different TTSuV1 genomes, the codon usage bias and GC content in the TTSuV1 coding region was lower, which was mainly determined by the base composition in the third codon position and the effective number of codons (ENC) value. Moreover, the results of correspondence analysis (COA) indicated that the codon usage patterns of TTSuV1 isolated from different countries varied greatly and had significant differences. In addition, Spearman’s rank correlation analysis and an ENC plot revealed that apart from mutation pressure, which was critical in determining the codon usage pattern, other factors were involved in shaping the evolution of codon usage bias in TTSuV1, such as natural selection. Those results suggested that synonymous codon usage patterns of TTSuV1 genomes were the result of interaction between mutation pressure and natural selection. The information from this study not only provides important insights into the synonymous codon usage pattern of TTSuV1, but also helps to identify the main factors affecting codon usage by this virus.

Patterns and influencing factor of synonymous codon usage in porcine circovirus

Background

Analysis of codon usage can reveal much about the molecular evolution of the viruses. Nevertheless, little information about synonymous codon usage pattern of porcine circovirus (PCV) genome in the process of its evolution is available. In this study, to give a new understanding on the evolutionary characteristics of PCV and the effects of natural selection from its host on the codon usage pattern of the virus, Patterns and the key determinants of codon usage in PCV were examined.

Methods

We carried out comprehensive analysis on codon usage pattern in the PCV genome, by calculating relative synonymous codon usage (RSCU), effective number of codons (ENC), dinucleotides and nucleic acid content of the PCV genome.

Results

PCV genomes have relatively much lower content of GC and codon preference, this result shows that nucleotide constraints have a major impact on its synonymous codon usage. The results of the correspondence analysis indicate codon usage patterns of PCV of various genotypes, various subgenotypes changed greatly, and significant differences in codon usage patterns of Each virus of Circoviridae.There is much comparability between PCV and its host in their synonymous codon usage, suggesting that the natural selection pressure from the host factor also affect the codon usage patterns of PCV. In particular, PCV genotype II is in synonymous codon usage more similar to pig than to PCV genotype I, which may be one of the most important molecular mechanisms of PCV genotype II to cause disease. The calculations results of the relative abundance of dinucleotides indicate that the composition of dinucleotides also plays a key role in the variation found in synonymous codon usage in PCV. Furthermore, geographic factors, the general average hydrophobicity and the aromaticity may be related to the formation of codon usage patterns of PCV.

Conclusion

The results of these studies suggest that synonymous codon usage pattern of PCV genome are the result of interaction between mutation pressure and natural selection from its host. The information from this study may not only have theoretical value in understanding the characteristics of synonymous codon usage in PCV genomes, but also have significant value for the molecular evolution of PCV.

Categorizing host-dependent RNA viruses by principal component analysis of their codon usage preferences

Viruses have to exploit host transcription and translation mechanisms to replicate in a hostile host cellular environment, and therefore, it is likely that the infected host may impose pressure on viral evolution. In this study, we investigated differences in codon usage preferences among the highly mutable single strain RNA viruses which infect vertebrate or invertebrate hosts, respectively. We incorporate principal component analysis (PCA) and k-mean methods to clustering viruses infected with different type of hosts. The relative synonymous codon usage (RSCU) indices of all genes in 32 RNA viruses were calculated, and the correlation of the RSCU indices among different viruses was analyzed by the PCA. Our results show a positive correlation in codon usage preferences among viruses that target the same host category. Results of k-means clustering analysis further confirmed the statistical significance of this study, demonstrating that viruses infecting vertebrate hosts have different codon usage preferences to those of invertebrate viruses. Based on the analysis of the effective number of codons (ENC) in relation to the GC-content at the synonymous third codon position (GC3s), we further identified that mutational pressure was the dominant evolution driving force in making the different codon usage preferences. This study suggests a new and effective way to characterize host-dependent RNA viruses based on the codon usage pattern.

Codon Usage Patterns in Corynebacterium glutamicum: Mutational Bias, Natural Selection and Amino Acid Conservation

The alternative synonymous codons in Corynebacterium glutamicum, a well-known bacterium used in industry for the production of amino acid, have been investigated by multivariate analysis. As C. glutamicum is a GC-rich organism, G and C are expected to predominate at the third position of codons. Indeed, overall codon usage analyses have indicated that C and/or G ending codons are predominant in this organism. Through multivariate statistical analysis, apart from mutational selection, we identified three other trends of codon usage variation among the genes. Firstly, the majority of highly expressed genes are scattered towards the positive end of the first axis, whereas the majority of lowly expressed genes are clustered towards the other end of the first axis. Furthermore, the distinct difference in the two sets of genes was that the C ending codons are predominate in putatively highly expressed genes, suggesting that the C ending codons are translationally optimal in this organism. Secondly, the majority of the putatively highly expressed genes have a tendency to locate on the leading strand, which indicates that replicational and transciptional selection might be invoked. Thirdly, highly expressed genes are more conserved than lowly expressed genes by synonymous and nonsynonymous substitutions among orthologous genes fromthe genomes of C. glutamicum and C. diphtheriae. We also analyzed other factors such as the length of genes and hydrophobicity that might influence codon usage and found their contributions to be weak.

Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences

Viruses differ markedly in their specificity toward host organisms. Here, we test the level of general sequence adaptation that viruses display toward their hosts. We compiled a representative data set of viruses that infect hosts ranging from bacteria to humans. We consider their respective amino acid and codon usages and compare them among the viruses and their hosts. We show that bacteria-infecting viruses are strongly adapted to their specific hosts, but that they differ from other unrelated bacterial hosts. Viruses that infect humans, but not those that infect other mammals or aves, show a strong resemblance to most mammalian and avian hosts, in terms of both amino acid and codon preferences. In groups of viruses that infect humans or other mammals, the highest observed level of adaptation of viral proteins to host codon usages is for those proteins that appear abundantly in the virion. In contrast, proteins that are known to participate in host-specific recognition do not necessarily adapt to their respective hosts. The implication for the potential of viral infectivity is discussed.

Analysis of synonymous codon usage in 11 human bocavirus isolates

Human Bocavirus (HBoV) is a novel virus which can cause respiratory tract disease in infants or children. In this study, the codon usage bias and the base composition variations in the available 11 complete HBoV genome sequences have been investigated. Although, there is a significant variation in codon usage bias among different HBoV genes, codon usage bias in HBoV is a little slight, which is mainly determined by the base compositions on the third codon position and the effective number of codons (ENC) value. The results of correspondence analysis (COA) and Spearman's rank correlation analysis reveals that the G + C compositional constraint is the main factor that determines the codon usage bias in HBoV and the gene's function also contributes to the codon usage in this virus. Moreover, it was found that the hydrophobicity of each protein and the gene length are also critical in affecting these viruses’ codon usage, although they were less important than that of the mutational bias and the genes’ function. At last, the relative synonymous codon usage (RSCU) of 44 genes from these 11 HBoV isolates is analyzed using a hierarchical cluster method. The result suggests that genes with same function yet from different isolates are classified into the same lineage and it does not depend on geographical location. These conclusions not only can offer an insight into the codon usage patterns and gene classification of HBoV, but also may help in increasing the efficiency of gene delivery/expression systems.