Directional mutation pressure and neutral molecular evolution

Insights into the evolutionary forces that shape the codon usage in the viral genome segments encoding intrinsically disordered protein regions

Intrinsically disordered regions/proteins (IDRs) are abundant across all the domains of life, where they perform important regulatory roles and supplement the biological functions of structured proteins/regions (SRs). Despite the multifunctionality features of IDRs, several interrogations on the evolution of viral genomic regions encoding IDRs in diverse viral proteins remain unreciprocated. To fill this gap, we benchmarked the findings of two most widely used and reliable intrinsic disorder prediction algorithms (IUPred2A and ESpritz) to a dataset of 6108 reference viral proteomes to unravel the multifaceted evolutionary forces that shape the codon usage in the viral genomic regions encoding for IDRs and SRs. We found persuasive evidence that the natural selection predominantly governs the evolution of codon usage in regions encoding IDRs by most of the viruses. In addition, we confirm not only that codon usage in regions encoding IDRs is less optimized for the protein synthesis machinery (transfer RNAs pool) of their host than for those encoding SRs, but also that the selective constraints imposed by codon bias sustain this reduced optimization in IDRs. Our analysis also establishes that IDRs in viruses are likely to tolerate more translational errors than SRs. All these findings hold true, irrespective of the disorder prediction algorithms used to classify IDRs. In conclusion, our study offers a novel perspective on the evolution of viral IDRs and the evolutionary adaptability to multiple taxonomically divergent hosts.

Evolutionary and codon usage preference insights into spike glycoprotein of SARS-CoV-2

Interaction of SARS-CoV-2 spike glycoprotein with the ACE2 cell receptor is very crucial for virus attachment to human cells. Selected mutations in SARS-CoV-2 S-protein are reported to strengthen its binding affinity to mammalian ACE2. The N501T mutation in SARS-CoV-2-CTD furnishes better support to hotspot 353 in comparison with SARS-CoV and shows higher affinity for receptor binding. Recombination analysis exhibited higher recombination events in SARS-CoV-2 strains, irrespective of their geographical origin or hosts. Investigation further supports a common origin among SARS-CoV-2 and its predecessors, SARS-CoV and bat-SARS-like-CoV. The recombination events suggest a constant exchange of genetic material among the co-infecting viruses in possible reservoirs and human hosts before SARS-CoV-2 emerged. Furthermore, a comprehensive analysis of codon usage bias (CUB) in SARS-CoV-2 revealed significant CUB among the S-genes of different beta-coronaviruses governed majorly by natural selection and mutation pressure. Various indices of codon usage of S-genes helped in quantifying its adaptability in other animal hosts. These findings might help in identifying potential experimental animal models for investigating pathogenicity for drugs and vaccine development experiments.

Effect of genome composition and codon bias on infectious bronchitis virus evolution and adaptation to target tissues

Background

Infectious bronchitis virus (IBV) is one of the most relevant viruses affecting the poultry industry, and several studies have investigated the factors involved in its biological cycle and evolution. However, very few of those studies focused on the effect of genome composition and the codon bias of different IBV proteins, despite the remarkable increase in available complete genomes. In the present study, all IBV complete genomes were downloaded (n = 383), and several statistics representative of genome composition and codon bias were calculated for each protein-coding sequence, including but not limited to, the nucleotide odds ratio, relative synonymous codon usage and effective number of codons. Additionally, viral codon usage was compared to host codon usage based on a collection of highly expressed genes in IBV target and nontarget tissues.

Results

The results obtained demonstrated a significant difference among structural, non-structural and accessory proteins, especially regarding dinucleotide composition, which appears under strong selective forces. In particular, some dinucleotide pairs, such as CpG, a probable target of the host innate immune response, are underrepresented in genes coding for pp1a, pp1ab, S and N. Although genome composition and dinucleotide bias appear to affect codon usage, additional selective forces may act directly on codon bias. Variability in relative synonymous codon usage and effective number of codons was found for different proteins, with structural proteins and polyproteins being more adapted to the codon bias of host target tissues. In contrast, accessory proteins had a more biased codon usage (i.e., lower number of preferred codons), which might contribute to the regulation of their expression level and timing throughout the cell cycle.

Conclusions

The present study confirms the existence of selective forces acting directly on the genome and not only indirectly through phenotype selection. This evidence might help understanding IBV biology and in developing attenuated strains without affecting the protein phenotype and therefore immunogenicity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07559-5.

Comparative chloroplast genomes: insights into the evolution of the chloroplast genome of Camellia sinensis and the phylogeny of Camellia

BACKGROUND

Chloroplast genome resources can provide useful information for the evolution of plant species. Tea plant (Camellia sinensis) is among the most economically valuable member of Camellia. Here, we determined the chloroplast genome of the first natural triploid Chinary type tea ('Wuyi narcissus' cultivar of Camellia sinensis var. sinensis, CWN) and conducted the genome comparison with the diploid Chinary type tea (Camellia sinensis var. sinensis, CSS) and two types of diploid Assamica type teas (Camellia sinensis var. assamica: Chinese Assamica type tea, CSA and Indian Assamica type tea, CIA). Further, the evolutionary mechanism of the chloroplast genome of Camellia sinensis and the relationships of Camellia species based on chloroplast genome were discussed.

RESULTS

Comparative analysis showed the evolutionary dynamics of chloroplast genome of Camellia sinensis were the repeats and insertion-deletions (indels), and distribution of the repeats, indels and substitutions were significantly correlated. Chinese tea and Indian tea had significant differences in the structural characteristic and the codon usage of the chloroplast genome. Analysis of sequence characterized amplified region (SCAR) using sequences of the intergenic spacers (trnE/trnT) showed none of 292 different Camellia sinensis cultivars had similar sequence characteristic to triploid CWN, but the other four Camellia species did. Estimations of the divergence time showed that CIA diverged from the common ancestor of two Assamica type teas about 6.2 Mya (CI: 4.4-8.1 Mya). CSS and CSA diverged to each other about 0.8 Mya (CI: 0.4-1.5 Mya). Moreover, phylogenetic clustering was not exactly consistent with the current taxonomy of Camellia.

CONCLUSIONS

The repeat-induced and indel-induced mutations were two important dynamics contributed to the diversification of the chloroplast genome in Camellia sinensis, which were not mutually exclusive. Chinese tea and Indian tea might have undergone different selection pressures. Chloroplast transfer occurred during the polyploid evolution in Camellia sinensis. In addition, our results supported the three different domestication origins of Chinary type tea, Chinese Assamica type tea and Indian Assamica type tea. And, the current classification of some Camellia species might need to be further discussed.

Genetic and evolutionary characterization of avian paramyxovirus type 4 in China

As an economically important poultry pathogen, avian paramyxovirus serotype 4 (APMV-4) frequently reported and isolated from domestic and wild birds particularly waterfowls worldwide. However, evolutionary dynamics of APMV-4 based on genomic characteristics is lacking. In this study, APMV-4 strain designated JX-G13 was isolated from oropharyngeal and cloacal swab samples of wild birds in China. Phylogenetic analysis revealed APMV-4 strains were divided into four genetic genotypes and China isolates were mainly clustered into Genotype I. The MCMC tree indicated that APMV-4 diverged about 104 years ago with the evolutionary rate of 1.2927 × 10^-3 substitutions/site/year. BSP analysis suggested that the effective population size of APMV-4 exhibited a steady state and decreased slowly after 2013. The F gene of APMV-4 was considered relatively conserved among isolates based on nucleotide diversity analysis. Although the F gene was under purifying selection, two positions (5 and 21) located in 3'-UTR were subject to positive selection. Our study firstly presented the evolutionary assessments on the genetic diversity of circulating APMV-4 from wild birds and domestic poultry.

Molecular features similarities between SARS-CoV-2, SARS, MERS and key human genes could favour the viral infections and trigger collateral effects

In December 2019, rising pneumonia cases caused by a novel β-coronavirus (SARS-CoV-2) occurred in Wuhan, China, which has rapidly spread worldwide, causing thousands of deaths. The WHO declared the SARS-CoV-2 outbreak as a public health emergency of international concern, since then several scientists are dedicated to its study. It has been observed that many human viruses have codon usage biases that match highly expressed proteins in the tissues they infect and depend on the host cell machinery for the replication and co-evolution. In this work, we analysed 91 molecular features and codon usage patterns for 339 viral genes and 463 human genes that consisted of 677,873 codon positions. Hereby, we selected the highly expressed genes from human lung tissue to perform computational studies that permit to compare their molecular features with those of SARS, SARS-CoV-2 and MERS genes. The integrated analysis of all the features revealed that certain viral genes and overexpressed human genes have similar codon usage patterns. The main pattern was the A/T bias that together with other features could propitiate the viral infection, enhanced by a host dependant specialization of the translation machinery of only some of the overexpressed genes. The envelope protein E, the membrane glycoprotein M and ORF7 could be further benefited. This could be the key for a facilitated translation and viral replication conducting to different comorbidities depending on the genetic variability of population due to the host translation machinery. This is the first codon usage approach that reveals which human genes could be potentially deregulated due to the codon usage similarities between the host and the viral genes when the virus is already inside the human cells of the lung tissues. Our work leaded to the identification of additional highly expressed human genes which are not the usual suspects but might play a role in the viral infection and settle the basis for further research in the field of human genetics associated with new viral infections. To identify the genes that could be deregulated under a viral infection is important to predict the collateral effects and determine which individuals would be more susceptible based on their genetic features and comorbidities associated.

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.

Analysis of codon usage patterns and influencing factors in rice tungro bacilliform virus

Rice tungro bacilliform virus (RTBV) belongs to genus Tungrovirus within the family Caulimoviridae harbors circular double-stranded DNA (dsDNA). Rice tungro disease (RTD) caused by RTBV, responsible for severe rice yield losses in South and Southeast Asia. Here, we performed a systematic evolutionary and codon usage bias (CUB) analysis of RTBV genome sequences. We analysed different bioinformatics techniques to calculate the nucleotide compositions, the relative synonymous codon usage (RSCU), and other indices. The results indicated slightly or low codon usage bias in RTBV isolates. Mutation and natural selection pressures have equally contributed to this low codon usage bias. Additionally, multiple factors such as host, geographical distribution also affect codon usage patterns in RTBV genomes. RSCU analysis revealed that RTBV shows mutation bias and prefers A and U ended codons to code amino acids. Codon usage patterns of RTBV were also found to be influenced by its host. This indicates that RTBV have evolved codon usage patterns that are specific to its host. The findings from this study are expected to increase our understanding of factors leading to viral evolution and fitness with respect to hosts and the environment.

Analysis of synonymous codon usage of transcriptome database in Rheum palmatum

Background

Rheum palmatum is an endangered and important medicinal plant in Asian countries, especially in China. However, there is little knowledge about the codon usage bias for R. palmatum CDSs. In this project, codon usage bias was determined based on the R. palmatum 2,626 predicted CDSs from R. palmatum transcriptome.

Methods

In this study, all codon usage bias parameters and nucleotide compositions were calculated by Python script, Codon W, DNA Star, CUSP of EMBOSS.

Results

The average GC and GC3 content are 46.57% and 46.6%, respectively, the results suggested that there exists a little more AT than GC in the R. palmatum genes, and the codon bias of R. palmatum genes preferred to end with A/T. We concluded that the codon bias in R. palmatum was affect by nucleotide composition, mutation pressure, natural selection, gene expression levels, and the mutation pressure is the prominent factor. In addition, we figured out 28 optimal codons and most of them ended with A or U. The project here can offer important information for further studies on enhancing the gene expression using codon optimization in heterogeneous expression system, predicting the genetic and evolutionary mechanisms in R. palmatum.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

The nucleotide usages significantly impact synonymous codon usage in Mycoplasma hyorhinis

Five annotated genomes of Mycoplasma hyorhinis were analyzed for clarifying evolutionary dynamics driving the overall codon usage pattern. Information entropy used for estimating nucleotide usage pattern at the gene level indicates that multiple evolutionary dynamics participate in forcing nucleotide usage bias at every codon position. Moreover, nucleotide usage bias directly contributes to synonymous codon usage biases with two different extremes. The overrepresented codons tended to have A/T in the third codon position, and the underrepresented codons strongly used G/C in the third position. Furthermore, correspondence analysis and neutrality plot reflect an obvious interplay between mutation pressure and natural selection mediating codon usage in M. hyorhinis genome. Due to significant bias in usages between A/T and G/C at the gene level, different selective forces have been proposed to contribute to codon usage preference in M. hyorhinis genome, including nucleotide composition constraint derived from mutation pressure, translational selection involved in natural selection, and strand-specific mutational bias represented by different nucleotide skew index. The systemic analyses of codon usage for M. hyorhinis can enable us to better understand the mechanisms of evolution in this species.

Genome-wide analysis of the synonymous codon usage pattern of Streptococcus suis

Streptococcus suis (S. suis) is a gram-positive coccus that causes disease in humans and animals. The codon usage pattern of bacteria reveals a range of evolutionary changes that assist them to enhance tolerance to environments. To better understand the genetic features during the evolution of S. suis, we performed codon usage analysis. Nine pathogenic strains of different serotypes and different geographical distribution were analyzed to better understand the differences in their evolutionary process. Nucleotide compositions and relative synonymous codon usage (RSCU) analysis revealed that A/T-ending codons are dominant in S. suis. Neutrality analysis, correspondence analysis and ENC-plot results revealed that natural selection is the predominant element prompting codon usage. Cluster analysis based on RSCU was roughly consistent with the dendrogram rooted genomic BLAST analysis. Comparison of synonymous codon usage pattern between S. suis and susceptible hosts (H. sapiens and S. scrofa) revealed that the codon usage of S. suis is separated from the synonymous codon usage of susceptible hosts. The CAI values implied that S. suis includes a series of predicted highly expressed coding sequences contained in metabolism and transcriptional regulation, revealing the necessity of this pathogen to deal with various environmental conditions. The study of codon usage in S. suis may provide evidence involving the molecular evolution of bacteria and a better understanding of evolutionary relationships between S. suis and its corresponding hosts.

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

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

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

Viruses belonging to the Coronaviridae family have a single-stranded positive-sense RNA with a poly-A tail. The genome has a length of ~29.9 kbps, which encodes for genes that are essential for cell survival and replication. Different evolutionary constraints constantly influence the codon usage bias (CUB) of different genes. A virus optimizes its codon usage to fit the host environment on which it savors. This study is a comprehensive analysis of the CUB for the different genes encoded by viruses of the Coronaviridae family. Different methods including relative synonymous codon usage (RSCU), an Effective number of codons (ENc), parity plot 2, and Neutrality plot, were adopted to analyze the factors responsible for the genetic evolution of the Coronaviridae family. Base composition and RSCU analyses demonstrated the presence of A-ended and U-ended codons being preferred in the 3rd codon position and are suggestive of mutational selection. The lesser ENc value for the spike ‘S’ gene suggests a higher bias in the codon usage of this gene compared to the other structural genes. Parity plot 2 and neutrality plot analyses demonstrate the role and the extent of mutational and natural selection towards the codon usage pattern. It was observed that the structural genes of the Coronaviridae family analyzed in this study were at the least under 84% influence of natural selection, implying a major role of natural selection in shaping the codon usage.

Analysis of codon usage patterns in citrus based on coding sequence data

BACKGROUND

Codon usage is an important determinant of gene expression levels that can help us understand codon biology, evolution and mRNA translation of species. The majority of previous codon usage studies have focused on single species analysis, although few studies have focused on the species within the same genus. In this study, we proposed a multispecies codon usage analysis workflow to reveal the genetic features and correlation in citrus.

RESULTS

Our codon usage analysis workflow was based on the GC content, GC plot, and relative synonymous codon usage value of each codon in 8 citrus species. This approach allows for the comparison of codon usage bias of different citrus species. Next, we performed cluster analysis and obtained an overview of the relationship in citrus. However, traditional methods cannot conduct quantitative analysis of the correlation. To further estimate the correlation among the citrus species, we used the frequency profile to construct feature vectors of each species. The Pearson correlation coefficient was used to quantitatively analyze the distance among the citrus species. This result was consistent with the cluster analysis.

CONCLUSIONS

Our findings showed that the citrus species are conserved at the genetic level and demonstrated the existing genetic evolutionary relationship in citrus. This work provides new insights into codon biology and the evolution of citrus and other plant species.

Genome-wide analysis of codon usage pattern in herpesviruses and its relation to evolution

The preferential use of a specific codon, out of a group of synonymous codons encoding the same amino acid, in a gene transcript results from the bias in codon choice. Various evolutionary forces namely mutation pressure and natural selection influence the pattern of codon usage i.e. distinct for each gene/genome. We investigated the pattern of codon usage of eight human herpesvirus genomes and compared them with two other herpesvirus genomes namely murine herpesvirus 68 and bovine herpesvirus type 1.1 to elucidate its compositional features, pattern of codon usage across the genomes and report the differences of codon usage pattern of human herpesviruses from that of other two other viruses. We also identified the similarity of the codon usage of human herpesviruses with its host (human). The genes were found to be CG rich in HHV2, HHV3, HHV4, HHV6, HHV7 and BH genomes while TA rich in HHV1, HHV5, HHV8 and MH genomes. The codon usage bias (CUB) of genes was low. A highly significant correlation was found among compositional contents depicting the role of mutational pressure along with natural selection in framing CUB. Several more frequently used codons as well as less frequently used codons were identified to be similar between each human virus and its host (human), while murine herpesvirus 68 and bovine herpesvirus type 1.1 genomes did not possess similar adaptation strategy as human herpesviruses to human (host), thus we could conclude that viral CUB might have been shaped as per their host's nature for better surveillance. Neutrality plot revealed mutational pressure mostly influenced the CUB of HHV1, HHV8 and MH viruses, while natural selection had a major impact in the CUB of HHV2, HHV3, HHV4, HHV5, HHV6, HHV7 and BH genomes.

Host Adaptation and Evolutionary Analysis of Zaire ebolavirus: Insights From Codon Usage Based Investigations

Ebola virus (EBOV) has caused several outbreaks as the consequence of spillover events from zoonotic sources and has resulted in huge death tolls. In spite of considerable progress, a thorough know-how regarding EBOV adaptation in various host species and detailed information about the potential reservoirs of EBOV still remains obscure. The present study was executed to examine the patterns of codon usage and its associated influence in the adaptation of EBOV to potential hosts that dwell in Africa, the origin of the viral outbreaks. Correspondence analysis (CA) revealed that the codon usage signature in EBOV is a complex interplay of factors including compositional bias and natural selection, with the latter having a more pronounced impact. Low codon usage bias in EBOV indicates a flexibility of the viruses in adapting to diverse range of hosts with different codon usage architectures. EBOV adaptation in potential hosts, as estimated by codon adaptation index (CAI) and relative codon deoptimization index (RCDI), revealed that the viruses were relatively better adapted to African primates than other mammals examined, which might account for the high fatality rate of primates owing to EBOV infection. Bats have been speculated as natural reservoirs of EBOV. In the present analysis it was interesting to note that EBOV displayed lower degrees of adaptation, as estimated by CAI and RCDI, with bats in comparison to the primate hosts. Lower degrees of adaptation might contribute to long-term co-existence and circulation of the viral pathogens in bat populations. Codon usage patterns of EBOV isolates associated with different outbreaks varied significantly, with discrete patterns between the West and Central African isolates. Additional evolutionary analyses indicated that the West African Epidemic began with an initial spillover infection and there was more than one population of EBOV circulating in the natural reservoir in the Democratic Republic of the Congo. The present study yields valuable information regarding the possible circulation of EBOV in various African mammals.

Comprehensive codon usage analysis of rice black-streaked dwarf virus based on P8 and P10 protein coding sequences

Rice black-streaked dwarf virus (RBSDV) belongs to the genus Fijivirus of the family Reoviridae and is an important pathogen that damages rice, maize and wheat worldwide. Previously, several reports have described the genetic variation and population structure of RBSDV. However, the details of the evolutionary changes, synonymous codon usage patterns and host adaptation of the virus are largely unclear. Here, we performed a detailed analysis of the codon usage and host adaptability of RBSDV based on 130 full-length P8 and 234 full-length P10 sequences. Infrequent recombination and frequent segment reassortment influence the genomic evolution of RBSDV. Our phylogenetic analysis found three and four lineages based on the P8 and P10 non-recombinant sequences respectively. We found relatively stable and conserved genomic composition with lower codon usage choice in the RBSDV P8 and P10 protein coding sequences. Both ENC-plot and neutrality-plot analyses showed that natural selection is the key factor that shapes the codon usage pattern of RBSDV. Codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses indicated strong correlation between RBSDV and rice rather than maize, wheat or Laodelphax striatellus. Our study provides deep insight into the evaluation of the codon usage pattern and adaptive evolution of RBSDV based on P8 and P10 sequences and should be taken into consideration for the prevention and control of this virus.

Codon bias analyses on thyroid carcinoma genes

BACKGROUND

Thyroid carcinoma is one of the most common cancers in the world. Although the genetics of thyroid carcinoma was intensively studied, new mechanisms could be involved in its development as the codon bias. In this paper, we studied the codon bias of thyroid-cancer genes, considering not only the sequences but also the synonymous mutations.

METHODS

Different measures and statistical analyses were employed to characterize the thyroid-cancer genes. We considered classical measures as RSCU and ENC, the compositional and protein characteristics, but also the codon bias landscape via the %MinMax algorithm.

RESULTS

The compositional analyses highlighted two groups of thyroid cancer genes according to the GC% and GC3% content. The ENC did not show a clear codon bias in the genes. Differently, the RSCU analyses showed interesting codons that could play an important role in the development of thyroid cancer as the codon Ser-tcG. Furthermore, interesting synonymous mutations were detected that could affect the codon bias. The codon bias landscape detected genes enriched in rare codons as AKAP9 and KTN1. A cluster analysis based on %MinMax classified the thyroid cancer genes in four different groups according to the distribution of rare/frequent codons in the sequence.

CONCLUSIONS

This is the first study that analyzed the codon bias in thyroid cancer genes based also on synonymous mutations. This study provided different hints that should be further investigated by wet-lab validation and that it could open new scenarios in the understanding the molecular mechanisms involved in thyroid cancer development based on codon bias.

Evolution of Transmissible Gastroenteritis Virus (TGEV): A Codon Usage Perspective

Transmissible gastroenteritis virus (TGEV) is a coronavirus associated with diarrhea and high mortality in piglets. To gain insight into the evolution and adaptation of TGEV, a comprehensive analysis of phylogeny and codon usage bias was performed. The phylogenetic analyses of maximum likelihood and Bayesian inference displayed two distinct genotypes: genotypes I and II, and genotype I was classified into subtypes Ia and Ib. The compositional properties revealed that the coding sequence contained a higher number of A/U nucleotides than G/C nucleotides, and that the synonymous codon third position was A/U-enriched. The principal component analysis based on the values of relative synonymous codon usage (RSCU) showed the genotype-specific codon usage patterns. The effective number of codons (ENC) indicated moderate codon usage bias in the TGEV genome. Dinucleotide analysis showed that CpA and UpG were over-represented and CpG was under-represented in the coding sequence of the TGEV genome. The analyses of Parity Rule 2 plot, ENC-plot, and neutrality plot displayed that natural selection was the dominant evolutionary driving force in shaping codon usage preference in genotypes Ia and II. In addition, natural selection played a major role, while mutation pressure had a minor role in driving the codon usage bias in genotype Ib. The codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses suggested that genotype I might be more adaptive to pigs than genotype II. Current findings contribute to understanding the evolution and adaptation of TGEV.

Insights into Genomic Epidemiology, Evolution, and Transmission Dynamics of Genotype VII of Class II Newcastle Disease Virus in China

Newcastle disease virus (NDV) is distributed worldwide and has caused significant losses to the poultry industry. Almost all virulent NDV strains belong to class II, among which genotype VII is the predominant genotype in China. However, the molecular evolution and phylodynamics of class II genotype VII NDV strains in China remained largely unknown. In this study, we identified 13 virulent NDV including 11 genotype VII strains and 2 genotype IX strains, from clinical samples during 1997 to 2019. Combined NDV sequences submitted to GenBank, we investigate evolution, and transmission dynamics of class II NDVs in China, especially genotype VII strains. Our results revealed that East and South China have the most genotypic diversity of class II NDV, and East China might be the origin of genotype VII NDVs in China. In addition, genotype VII NDVs in China are presumably transmitted by chickens, as the virus was most prevalent in chickens. Furthermore, codon usage analysis revealed that the F genes of genotype VII NDVs have stronger adaptation in chickens, and six amino acids in this gene are found under positive selection via selection model analysis. Collectively, our results revealed the genetic diversity and evolutionary dynamics of genotype VII NDVs in China, providing important insights into the epidemiology of these viruses in China.

A Crosstalk on Codon Usage in Genes Associated with Leukemia

Leukemia is the outcome of aggregation of damaged white blood cells. Several genes were reported to be associated with the pathogenesis of leukemia. These genes were computationally analyzed to decipher their codon usage bias (CUB) and to identify the prime factors influencing the codon usage profile as no work was reported yet. The mean values of synonymous codon usage order (SCUO) parameter indicated low CUB of the genes. Significant positive association of SCUO with overall GC and positional GCs might signal the presence of mutational pressure. However, neutrality plot suggested the dominant role of natural selection across the genes. Along with natural selection, the role of mutation pressure was also prominent and that might be responsible for lower CUB (SCUO = 0.19) of genes. Low translational speed might permit accuracy in the process. A strong inverse relationship of translational rate was observed with CUB of genes and folding energy.

Analysis of the Codon Usage Pattern of HA and NA Genes of H7N9 Influenza A Virus

Novel H7N9 influenza virus transmitted from birds to human and, since March 2013, it has caused five epidemic waves in China. Although the evolution of H7N9 viruses has been investigated, the evolutionary changes associated with codon usage are still unclear. Herein, the codon usage pattern of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA), was studied to understand the evolutionary changes in relation to host, epidemic wave, and pathogenicity. Both genes displayed a low codon usage bias, with HA higher than NA. The codon usage was driven by mutation pressure and natural selection, although the main contributing factor was natural selection. Additionally, the codon adaptation index (CAI) and deoptimization (RCDI) illustrated the strong adaptability of H7N9 to Gallus gallus. Similarity index (SiD) analysis showed that Homo sapiens posed a stronger selection pressure than Gallus gallus. Thus, we assume that this may be related to the gradual adaptability of the virus to human. In addition, the host strong selection pressure was validated based on CpG dinucleotide content. In conclusion, this study analyzed the usage of codons of two genes of H7N9 and expanded our understanding of H7N9 host specificity. This aids into the development of control measures against H7N9 influenza virus.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Substitution scoring matrices for proteins - An overview

Sequence analysis is the primary and simplest approach to discover structural, functional and evolutionary details of related proteins. All the alignment based approaches of sequence analysis make use of amino acid substitution matrices, and the accuracy of the results largely depends on the type of scoring matrices used to perform alignment tasks. An amino acid substitution matrix is a 20 × 20 matrix in which the individual elements encapsulate the rates at which each of the 20 amino acid residues in proteins are substituted by other amino acid residues over time. In contrast to most globular/ordered proteins whose amino acids composition is considered as standard, there are several classes of proteins (e.g., transmembrane proteins) in which certain types of amino acid (e.g., hydrophobic residues) are enriched. These compositional differences among various classes of proteins are manifested in their underlying residue substitution frequencies. Therefore, each of the compositionally distinct class of proteins or protein segments should be studied using specific scoring matrices that reflect their distinct residue substitution pattern. In this review, we describe the development and application of various substitution scoring matrices peculiar to proteins with standard and biased compositions. Along with most commonly used standard matrices (PAM, BLOSUM, MD and VTML) that act as default parameters in various homologs search and alignment tools, different substitution scoring matrices specific to compositionally distinct class of proteins are discussed in detail.

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.

Comprehensive Analysis of Codon Usage on Porcine Astrovirus

Porcine astrovirus (PAstV), associated with mild diarrhea and neurological disease, is transmitted in pig farms worldwide. The purpose of this study is to elucidate the main factors affecting codon usage to PAstVs. Phylogenetic analysis showed that the subtype PAstV-5 sat at the bottom of phylogenetic tree, followed by PAstV-3, PAstV-1, PAstV-2, and PAstV-4, indicating that the five existing subtypes (PAstV1-PAstV5) may be formed by multiple differentiations of PAstV ancestors. A codon usage bias was found in the PAstVs-2,3,4,5 from the analyses of effective number of codons (ENC) and relative synonymous codon usage (RSCU). Nucleotides A/U are more frequently used than nucleotides C/G in the genome CDSs of the PAstVs-3,4,5. Codon usage patterns of PAstV-5 are dominated by mutation pressure and natural selection, while natural selection is the main evolutionary force that affects the codon usage pattern of PAstVs-2,3,4. The analyses of codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) showed the codon usage similarities between the PAstV and animals might contribute to the broad host range and the cross-species transmission of astrovirus. Our results provide insight into understanding the PAstV evolution and codon usage patterns.

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

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

Codon usage pattern and its influencing factors for mitochondrial CO genes among different classes of Arthropoda

Analysis of codon usage bias (CUB) is very much important in perceiving the knowledge of molecular biology, the discovery of a new gene, designing of transgenes and evolution of gene. In this study, we analyzed compositional features and codon usage of MT-CO (COI, COII and COIII) genes among the classes of Arthropoda to explore the pattern of CUB as no research work was reported yet. Nucleotide composition analysis in CO genes suggested that the genes were AT-rich in all the four classes of Arthropoda. CUB was low in all the classes of Arthropoda for MT-CO genes as revealed from a high effective number of codons (ENC). We also found that the evolutionary forces namely mutation pressure and natural selection were the key influencing factors in CUB among MT-CO genes as revealed by correlation analysis between overall nucleotide composition and nucleotide composition at the 3rd codon position. Correspondence analysis suggested that the pattern of CUB was different among the classes of Arthropoda. Further, it was revealed from the neutrality plot that natural selection had a dominant role while mutation pressure exhibited a minor role in structuring the pattern of codon usage in all the classes of Arthropoda across COI, COII and COIII genes.

Codon Usage Optimization in the Prokaryotic Tree of Life: How Synonymous Codons Are Differentially Selected in Sequence Domains with Different Expression Levels and Degrees of Conservation

The prokaryotic genomes—the current heritage of the most ancient life forms on earth—are comprised of diverse gene sets, all characterized by varied origins, ancestries, and spatial-temporal expression patterns. Such genetic diversity has for a long time raised the question of how cells shape their coding strategies to optimize protein demands (i.e., product abundance) and accuracy (i.e., translation fidelity) through the use of the same genetic code in genomes with GC contents that range from less than 20 to more than 80%. Here, we present evidence on how codon usage is adjusted in the prokaryotic tree of life and on how specific biases have operated to improve translation. Through the use of proteome data, we characterized conserved and variable sequence domains in genes of either high or low expression level and quantitated the relative weight of efficiency and accuracy—as well as their interaction—in shaping codon usage in prokaryotes.

Prokaryote genomes exhibit a wide range of GC contents and codon usages, both resulting from an interaction between mutational bias and natural selection. In order to investigate the basis underlying specific codon changes, we performed a comprehensive analysis of 29 different prokaryote families. The analysis of core gene sets with increasing ancestries in each family lineage revealed that the codon usages became progressively more adapted to the tRNA pools. While, as previously reported, highly expressed genes presented the most optimized codon usage, the singletons contained the less selectively favored codons. The results showed that usually codons with the highest translational adaptation were preferentially enriched. In agreement with previous reports, a C bias in 2- to 3-fold pyrimidine-ending codons, and a U bias in 4-fold codons occurred in all families, irrespective of the global genomic GC content. Furthermore, the U biases suggested that U₃-mRNA–U₃₄-tRNA interactions were responsible for a prominent codon optimization in both the most ancestral core and the highly expressed genes. A comparative analysis of sequences that encode conserved (cr) or variable (vr) translated products, with each one being under high (HEP) and low (LEP) expression levels, demonstrated that the efficiency was more relevant (by a factor of 2) than accuracy to modeling codon usage. Finally, analysis of the third position of codons (GC3) revealed that in genomes with global GC contents higher than 35 to 40%, selection favored a GC3 increase, whereas in genomes with very low GC contents, a decrease in GC3 occurred. A comprehensive final model is presented in which all patterns of codon usage variations are condensed in four distinct behavioral groups.

Generalizing rate heterogeneity across sites in statistical phylogenetics

Phylogenetics uses alignments of molecular sequence data to learn about evolutionary trees relating species. Along branches, sequence evolution is modelled using a continuous-time Markov process characterized by an instantaneous rate matrix. Early models assumed the same rate matrix governed substitutions at all sites of the alignment, ignoring variation in evolutionary pressures. Substantial improvements in phylogenetic inference and model fit were achieved by augmenting these models with multiplicative random effects that describe the result of variation in selective constraints and allow sites to evolve at different rates which linearly scale a baseline rate matrix. Motivated by this pioneering work, we consider an extension using a quadratic, rather than linear, transformation. The resulting models allow for variation in the selective coefficients of different types of point mutation at a site in addition to variation in selective constraints.

We derive properties of the extended models. For certain non-stationary processes, the extension gives a model that allows variation in sequence composition, both across sites and taxa. We adopt a Bayesian approach, describe an MCMC algorithm for posterior inference and provide software. Our quadratic models are applied to alignments spanning the tree of life and compared with site-homogeneous and linear models.

Analysis of codon usage pattern of mitochondrial ND genes in Platyhelminthes

Codon usage bias (CUB) is the nonrandom usage of synonymous codons in which some codons are more preferred to others.CUB can be determined by mutation pressure and selection. Various approaches have been used to understand the pattern of CUB in the mitochondrial ND (MT-ND or ND) genes involved in complex I of respiratory chain in five different classes of Platyhelminthes as no work was reported yet. The present study revealed that the CUB varies across MT-ND genes and the coding sequences showed the richness of A and T. Correspondence analysis implied the effect of mutational pressure and also the pattern of codon usage was different in different classes of platyhelminthes for MT-ND genes. Highly significant correlation was observed between overall nucleotide composition and its 3rd codon position in most of the homogeneous nucleotides such as A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% and also some significant correlations observed among heterogeneous nucleotides in all the five classes for MT-ND genes suggested the role of mutational pressure as well as natural selection in affecting the CUB. Neutrality plot suggested that the contributions of natural selection and mutational pressure varied across different classes of platyhelminthes and also differed in different MT-ND genes.

Mycobacterium lepromatosis genome exhibits unusually high CpG dinucleotide content and selection is key force in shaping codon usage

Mycobacterium lepromatosis was identified as a causative agent for leprosy in the year 2008 in the United States and later more cases were identified in Canada, Singapore, Brazil, and Myanmar. It is known to cause diffuse lepromatosis leprosy among humans. Since it is invasive, the mortality rates are higher in comparison to the M. leprae. At genomic level, there exists 90.9% similarity between M. lepromatosis and M. leprae. Codon usage analysis based on analyses of 228 coding sequences (CDSs) of M. lepromatosis, revealed that the genome is GC rich. Among the total 16 dinucleotides, CpG dinucleotide possesses the highest dinucleotide frequency in M. lepromatosis, that is strikingly an unobvious observation since higher CpG is associated with higher proinflammatory cytokine production and NF-κB activation that eventually leads to high pathogenicity. To evade immune response, CpG content is generally less in pathogens. The unusually high CpG content can be explained by the fact that the nucleotide composition of M. lepromatosis is CG rich. Various forces interplay to shape codon usage pattern of any organism including selection; mutation, nucleotide composition as well as GC biased gene conversion. To understand the interplay between various forces; neutrality, parity, Nc-GC3 (Effective number of codons-GC content at 3rd position of the codon), aromaticity (AROMO) and the general average hydropathicity score (GRAVY) analyses have been carried out. The analyses revealed that selection force is the major contributory force. Along with the selection; mutation, nucleotide composition as well as GC biased gene conversion also play role in shaping codon usage bias in M. lepromatosis. This is the first report on the codon usage in M. lepromatosis.

Codon usage trend in genes associated with obesity

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

Analysis of compositional properties and codon usage bias of mitochondrial CYB gene in anura, urodela and gymnophiona

We delineated the pattern of synonymous codon usage bias (CUB) and its determinants in mitochondrial CYB gene of respiratory chain across different amphibian groups namely orders anura, urodela and gymnophiona. We observed that CUB was low in CYB gene of amphibia. The gymnophionans had comparatively high bias followed by urodeles and anurans. The codons namely TCA, CCA, CAA, CGA, TGA, AAA and ACA were over-represented in all three orders. The codons such as GCC and TCC were over-represented in anura whereas in urodela, the over-represented codons were TTA, CTA, ATA, GTA, GAA, GGA and GCA. In gymnophiona, GCC, TTA, CTA, ATA, GTA, GAA and GGA codons were over-represented. The regression analysis between effective number of codons (ENC) and nucleobase at the 3rd position revealed that nucleobase A and C influenced CUB positively in order anura, while in urodela and gymnophiona, nucleobase A and T influenced the CUB positively. Mutation pressure and natural selection mutually illustrate the CUB of CYB gene (complex III gene) of amphibia as elucidated by correlation analysis between 3rd nucleotide in a codon and overall nucleotide content of the gene. However, neutrality plot showed that natural selection was the dominant evolutionary factor of CUB.

Genetic changes and host adaptability in sugarcane mosaic virus based on complete genome sequences

Sugarcane mosaic virus (SCMV), a member of the genus Potyvirus in the family Potyviridae, is an important pathogen that causes mosaic diseases in maize, sugarcane, canna and other graminaceous species worldwide. Previously, several reports have showed the genetic variation and population structure of SCMV. However, the evolutionary dynamics, synonymous codon usage pattern and adaptive evolution of the virus is unclear. In this study, we performed comprehensive analyses of phylodynamics, composition bias and codon usage of SCMV using 108 complete genomic sequences. Our phylogenetic analysis found six host- and geographically confined phylogenetic lineages within the SCMV non-recombinant isolates. We found a relatively stable and conserved genomic composition with a lower codon usage choice in the SCMV protein coding sequences. Mutation pressure and natural selection have shaped the codon usage patterns of the SCMV protein coding sequences with natural selection being the dominant factor. The codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses revealed a stronger correlation between SCMV and maize than between SCMV and sugarcane or canna. Our study is the first to evaluate the codon usage pattern of SCMV based on complete sequences and may provide a better understanding of the origin of SCMV and its evolutionary patterns for future research.

Analysis of Codon Usage and Nucleotide Bias in Middle East Respiratory Syndrome Coronavirus Genes

The Middle East Respiratory Syndrome (MERS) is an emerging disease caused by a recently identified human coronavirus (CoV). Over 2494 laboratory-confirmed cases and 858 MERS-related deaths have been reported from 27 countries. MERS-CoV has been associated with a high case fatality rate, especially in patients with pre-existing conditions. Despite the fatal nature of MERS-CoV infection, a comprehensive study to explore its evolution and adaptation in different hosts is lacking. We performed codon usage analyses on 4751 MERS-CoV genes and determined underlying forces that affect the codon usage bias in the MERS-CoV genome. The current analyses revealed a low but highly conserved, gene-specific codon usage bias in the MERS-CoV genome. The codon usage bias is mainly shaped by natural selection, while mutational pressure emerged as a minor factor affecting codon usage in some genes. Other contributory factors included CpG dinucleotide bias, physical and chemical properties of encoded proteins and gene length. Results reported in this study provide considerable insights into the molecular evaluation of MERS-CoV and could serve as a theoretical basis for optimizing MERS-CoV gene expression to study the functional relevance of various MERS-CoV proteins. Alternatively, an attenuated vaccine strain containing hundreds of silent mutations could be engineered. Codon de-optimization will not affect the amino acid sequence or antigenicity of a vaccine strain, but the sheer number of mutations would make viral reversion to a virulent phenotype extremely unlikely.

Evidence of multifaceted functions of codon usage in translation within the model beetle Tribolium castaneum

Synonymous codon use is non-random. Codons most used in highly transcribed genes, often called optimal codons, typically have high gene counts of matching tRNA genes (tRNA abundance) and promote accurate and/or efficient translation. Non-optimal codons, those least used in highly expressed genes, may also affect translation. In multicellular organisms, codon optimality may vary among tissues. At present, however, tissue specificity of codon use remains poorly understood. Here, we studied codon usage of genes highly transcribed in germ line (testis and ovary) and somatic tissues (gonadectomized males and females) of the beetle Tribolium castaneum. The results demonstrate that: (i) the majority of optimal codons were organism-wide, the same in all tissues, and had numerous matching tRNA gene copies (Opt-codon_↑tRNAs), consistent with translational selection; (ii) some optimal codons varied among tissues, suggesting tissue-specific tRNA populations; (iii) wobble tRNA were required for translation of certain optimal codons (Opt-codon_wobble), possibly allowing precise translation and/or protein folding; and (iv) remarkably, some non-optimal codons had abundant tRNA genes (Nonopt-codon_↑tRNAs), and genes using those codons were tightly linked to ribosomal and stress-response functions. Thus, Nonopt-codon_↑tRNAs codons may regulate translation of specific genes. Together, the evidence suggests that codon use and tRNA genes regulate multiple translational processes in T. castaneum.

Synonymous codon usage and context analysis of genes associated with pancreatic cancer

Pancreatic cancer is a fatal disorder which originates in pancreas. Its mortality rate is increasing with time. Some studies also reported that pancreatic cancer would be ranked 2nd by the year 2030. Codon usage bias (CUB) arises when synonymous codons for each amino acid are not used randomly in the coding sequences of genes. We used bioinformatic methods to analyze the compositional properties, codon context and codon usage trend of the genes associated with pancreatic cancer as no work was reported yet. From the base composition analysis, the pancreatic cancer genes were found to be GC-rich and at the 3rd codon position the G/C ending codons were more preferred to A/T ending ones. The CUB was low in genes associated with pancreatic cancer. Correspondence analysis proposed that other than base constraints, CUB might also be affected by some other factors such as natural selection. Moreover, results of correlation analysis indicated that CUB and various GC contents i.e. GC, GC1, GC2, GC3 played important role in the release of free energy by transcripts of the genes associated with pancreatic cancer. The low compAI values of coding sequences suggested a low translation rate of the genes.

Codon Usage and Phenotypic Divergences of SARS-CoV-2 Genes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which first occurred in Wuhan (China) in December of 2019, causes a severe acute respiratory illness with a high mortality rate, and has spread around the world. To gain an understanding of the evolution of the newly emerging SARS-CoV-2, we herein analyzed the codon usage pattern of SARS-CoV-2. For this purpose, we compared the codon usage of SARS-CoV-2 with that of other viruses belonging to the subfamily of Orthocoronavirinae. We found that SARS-CoV-2 has a high AU content that strongly influences its codon usage, which appears to be better adapted to the human host. We also studied the evolutionary pressures that influence the codon usage of five conserved coronavirus genes encoding the viral replicase, spike, envelope, membrane and nucleocapsid proteins. We found different patterns of both mutational bias and natural selection that affect the codon usage of these genes. Moreover, we show here that the two integral membrane proteins (matrix and envelope) tend to evolve slowly by accumulating nucleotide mutations on their corresponding genes. Conversely, genes encoding nucleocapsid (N), viral replicase and spike proteins (S), although they are regarded as are important targets for the development of vaccines and antiviral drugs, tend to evolve faster in comparison to the two genes mentioned above. Overall, our results suggest that the higher divergence observed for the latter three genes could represent a significant barrier in the development of antiviral therapeutics against SARS-CoV-2.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Distinctive features gleaned from the comparative genomes analysis of clinical and non-clinical isolates of Klebsiella pneumoniae

It is of interest to describe the distinctive features gleaned from the comparative genome analysis of clinical and non-clinical isolates of Klebsiella pneumoniae. The core genome of K. pneumoinae consisted of 3568 genes. Comparative genome analysis shows that mdtABCD, toxin-antitoxin systems are unique to clinical isolates and catB, benA, and transporter genes for citrate utilization are exclusive to non-clinical isolates. We further noted aromatic compound degrading genes in non-clinical isolates unlike in the later isolates. We grouped 88 core genes into 3 groups linked to infections, drug-resistance or xenobiotic metabolism using codon usage variation analysis. It is inferred using the neutrality plot analysis of GC12 with GC3 that codon usage variation is dominant over mutation pressure. Thus, we document data to distinguish clinical and non-clinical isolates of K. pneumoniae using comparative genomes analysis for understanding of genome diversity during speciation.

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.

Selection constraints determine preference for A/U-ending codons in Taxus contorta

Unequal utilization of synonymous codons is a well-known phenomenon among living organisms. This phenomenon plays a major role in the enhancement of the accuracy and efficiency of translation. Gymnosperms are rarely paid attention in this aspect. Understanding the degree of and determining the forces influencing codon usage bias (CUB) in Taxus contorta, an endangered Himalayan gymnosperm, will prove useful in interpreting the evolutionary characteristics of this species. Using RNAseq data, 93 790 assembled transcripts were clustered into 32 701 unigenes. Around 13 061 full-length sequences were utilized for the analysis of CUB. Compositional properties showed that GC-content ranged from 28.76% to 65.22%, with an average value of 44.28%, suggesting an AT-rich genome. The mean effective number of codons (ENC) value revealed that CUB is not strong in T. contorta. The preferred codons tended to be A/U ending, whereas the avoided codons tended to be G/C ending. A P2 index of 0.54 and a Mutation Responsive Index (MRI) value of -0.02 in addition to the results revealed by the neutrality, ENC, and parity plots showed that natural selection is a predominating factor governing CUB. Mutational pressure, gene length, hydropathiciy, aromaticity, and nucleotide composition influence CUB weakly.