Codon Adaptation Index Values Research Articles

The evolutionary and genetic patterns of African swine fever virus

African swine fever (ASF) is a serious animal disease, and has spread to Africa, Europe and Asia, causing massive economic losses. African swine fever virus (ASFV) is transmitted from a reservoir host (warthog) to domestic pigs via a sylvatic cycle (transmission between warthogs and soft ticks) and a domestic cycle (transmission between domestic pigs) and survives by expressing a variety of genes related to virus–host interactions. We evaluated differences in codon usage patterns among ASFV genotypes and clades and explored the common and specific evolutionary and genetic characteristics of ASFV sequences. We analysed the evolutionary relationships, nucleotide compositions, codon usage patterns, selection pressures (mutational pressure and natural selection) and viral adaptation to host codon usage based on the coding sequences (CDS) of key functional genes of ASFV. AT bias was detected in the six genes analysed, irrespective of clade. The AT bias of genes (A224L, A179L, EP153R) encoding proteins involved in interaction with host cells after infection was high; among them, the AT bias of EP153R was the greatest at 78.3%. A large number of overrepresented codons were identified in EP153R, whereas there were no overrepresented codons with a relative synonymous codon usage (RSCU) value of ≥3 in B646L. In most genes, the pattern of selection pressure was similar for each clade, but in EP153R, diverse patterns of selection pressure were captured within the same clade and genotype. As a result of evaluating host adaptation based on the codon adaptation index (CAI), for B646L, E183L, CP204L and A179L, the codon usage patterns in all sequences were more similar to tick than domestic pig or wild boar. However, EP153R showed the lowest average CAI value of 0.52 when selecting tick as a reference set. The genes analysed in this study showed different magnitudes of selection pressure at the clade and genotype levels, which is likely to be related to the function of the encoded proteins and may determine key evolutionary traits of viruses, such as the level of genetic variation and host range. The diversity of codon adaptations at the genetic level in ASFV may account for differences in translational selection in ASFV hosts and provides insight into viral host adaptation and co-evolution.

Investigation of CaMV-host co-evolution through synonymous codon pattern.

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

DERNA Enables Pareto Optimal RNA Design.

The design of an RNA sequence that encodes an input target protein sequence is a crucial aspect of messenger RNA (mRNA) vaccine development. There are an exponential number of possible RNA sequences for a single target protein due to codon degeneracy. These potential RNA sequences can assume various secondary structure conformations, each with distinct minimum free energy (MFE), impacting thermodynamic stability and mRNA half-life. Furthermore, the presence of species-specific codon usage bias, quantified by the codon adaptation index (CAI), plays a vital role in translation efficiency. While earlier studies focused on optimizing either MFE or CAI, recent research has underscored the advantages of simultaneously optimizing both objectives. However, optimizing one objective comes at the expense of the other. In this work, we present the Pareto Optimal RNA Design problem, aiming to identify the set of Pareto optimal solutions for which no alternative solutions exist that exhibit better MFE and CAI values. Our algorithm DEsign RNA (DERNA) uses the weighted sum method to enumerate the Pareto front by optimizing convex combinations of both objectives. We use dynamic programming to solve each convex combination in time and space. Compared with a CDSfold, previous approach that only optimizes MFE, we show on a benchmark data set that DERNA obtains solutions with identical MFE but superior CAI. Moreover, we show that DERNA matches the performance in terms of solution quality of LinearDesign, a recent approach that similarly seeks to balance MFE and CAI. We conclude by demonstrating our method's potential for mRNA vaccine design for the SARS-CoV-2 spike protein.

A detailed analysis of codon usages bias and influencing factors in the nucleocapsid gene of Nipah Virus

Several outbreaks of Nipah Virus (NiV) have recently been reported in various parts of the world including India. The 18.2-kb size single-stranded negative-sense RNA genome of NiV has six genes that encode nine proteins. Among nine, the nucleoprotein (N) protein is the most abundant viral protein. The relative availability of N protein is a determining factor for the activation of genome encapsidation, replicase activity, and regulating viral RNA synthesis. In this study, we have conducted a codon usage analysis of N protein-encoding gene (N gene) of NiV. The relative synonymous codon usage (RSCU) values, in combination with an ENC value of 51.048, represented low codon usage bias in the N gene. The effect of mutational pressure on codon usage bias was confirmed by significant correlations of GC3s, G3s, C3s, A3s, U3s, and ENC values with whole nucleotide contents (GC%, G%, C%, A%, and U%). Correlation study of GC3s, G3s, C3s, A3s, and U3s with axis values of correspondence analysis (CA) also supported the role of mutational pressure. The correlation study of Gravy values with GC3s, G3s, C3s, A3s, and U3s revealed the presence of natural selection in addition to mutational pressure on codon usage bias. Moreover, NiV codon adaptation index (CAI) value higher than their corresponding expected CAI (eCAI) values against human (CAI, 0.726; eCAI, 0.715), pig (CAI, 0.837; eCAI, 0.819), and bat (CAI, 0.762; eCAI, 0.750) also indicated natural selection play role on codon usage bias. Additionally, geographical distribution, and evolutionary processes also influenced the codon usage bias to some extent.

Base composition, adaptation, and evolution of goose astroviruses: codon-based investigation

Goose astroviruses (GoAstVs) are causative agents that account for fatal infection of goslings characterized by visceral urate deposition, resulting in severe economic losses in major goose-producing regions in China since 2017. In this study, we sought to unravel the intrinsic properties associated with adaptation and evolution in the host environment of GoAstVs. Consistent results from phylogenetic analysis and correspondence analysis performed on the codon usage patterns (CUPs) reveal 2 clusters of GoAstVs, namely, GoAstV-1 and GoAstV-2. However, multiple similar compositional characteristics were found, despite the high divergence between GoAstV-1 and GoAstV-2. Studies on the base composition of GoAstVs reveal an A/U bias, indicating a compositional constraint, while natural selection prevailed in determining the CUPs in the virus genome based on our neutrality plot analysis, reflecting high adaptive pressure to fit the host environment. Codon adaptation index (CAI) analysis revealed a higher degree of fitness to the CUPs of the corresponding host for GoAstVs than avian influenza virus and betacoronaviruses, which may be a favorable factor contributing to the high pathogenicity and wide distribution of GoAstVs in goslings. In addition, GoAstVs were less adapted to ducks and chickens, with significantly lower CAI values than to geese, which may be a reason for the different prevalence of GoAstVs among these species. Extensive investigations on dinucleotide distribution revealed a significant suppression of the CpG and UpA motifs in the virus genome, which may facilitate adaptation to the host's innate immune system by evading surveillance. In addition, our study reported the trends of increasing fitness to the host's microenvironment for GoAstVs through increasing adaptation to host CUPs and ongoing reduction of CpG motifs in the virus genome. The present analysis deepens our understanding of the basic biology, pathogenesis, adaptation and evolutionary pattern of GoAstVs, and contributes to the development of novel antiviral strategies.

Intracellular gene transfer and codon usage of cytoplasmic genomes in green plants

Abstract Intracellular gene transfer is widely recognized as one of the most important driving forces for species evolution. Here we investigated transferred cytoplasmic motifs in green plants including spore-bearing plants and seed-bearing plants (hereafter termed spore plants and seed plants) . Our analyses revealed that gene transfer in spore plants was characterized by shorter motifs than that of seed plants. Several spore species did not exhibit intracellular gene transfer. Meanwhile, high frequency transferred tRNA genes were identified with average values of minimum free energy at moderate level. From the chloroplast to the mitochondrial genome, trnP was found to have transferred with high frequency in green plants. In gene transfer from the mitochondrial to the chloroplast genome, trnN was found to be a highly transferred gene. We observed that several tRNA genes including trnF, trnW, and trnN were involved in bidirectional transfer, which may be related to application strategy of functional protein-coding genes in a plant’s adaptive evolution. Codon Adaptation Index (CAI) analysis showed that codon usage was unbalanced in spore and seed plants. CAI values for seed plants were higher than those for spore plants in general, which may reveal rapid divergence adaptability of codon usage in the former. These results provide novel insights into gene transfer and codon usage within cytoplasmic genomes.

Designing a multi-epitopes subunit vaccine against human herpes virus 6A based on molecular dynamics and immune stimulation

Human Herpesvirus 6A (HHV-6A) is a prevalent virus associated with various clinical manifestations, including neurological disorders, autoimmune diseases, and promotes tumor cell growth. HHV-6A is an enveloped, double-stranded DNA virus with a genome of approximately 160–170 kb containing a hundred open-reading frames. An immunoinformatics approach was applied to predict high immunogenic and non-allergenic CTL, HTL, and B cell epitopes and design a multi-epitope subunit vaccine based on HHV-6A glycoprotein B (gB), glycoprotein H (gH), and glycoprotein Q (gQ). The stability and correct folding of the modeled vaccines were confirmed through molecular dynamics simulation. Molecular docking found that the designed vaccines have a strong binding network with human TLR3, with Kd values of 1.5E-11 mol/L, 2.6E-12 mol/L, 6.5E-13 mol/L, and 7.1E-11 mol/L for gB-TLR3, gH-TLR3, gQ-TLR3, and the combined vaccine-TLR3, respectively. The codon adaptation index values of the vaccines were above 0.8, and their GC content was around 67 % (normal range 30–70 %), indicating their potential for high expression. Immune simulation analysis demonstrated robust immune responses against the vaccine, with approximately 650,000/ml combined IgG and IgM antibody titer. This study lays a strong foundation for developing a safe and effective vaccine against HHV-6A, with significant implications for treating associated conditions.

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content, protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.

Designing a novel mRNA vaccine against Vibrio harveyi infection in fish: an immunoinformatics approach.

Vibrio harveyi belongs to the family Vibrionaceae of class Gammaproteobacteria. Around 12 Vibrio species can cause gastroenteritis (gastrointestinal illness) in humans. A large number of bacterial particles can be found in the infected cells, which may cause death. Despite these devastating complications, there is still no cure or vaccine for the bacteria. As a result, we used an immunoinformatics approach to develop a multi-epitope vaccine against the most pathogenic hemolysin gene of V. harveyi. The immunodominant T- and B-cell epitopes were identified using the hemolysin protein. We developed a vaccine employing three possible epitopes: cytotoxic T-lymphocytes, helper T-lymphocytes, and linear B-lymphocyte epitopes, after thorough testing. The vaccine was developed to be antigenic, immunogenic, and non-allergenic, as well as have a better solubility. Molecular dynamics simulation revealed significant structural stiffness and binding stability. In addition, the immunological simulation generated by computers revealed that the vaccination might elicit immune reactions Escherichia coli K12 as a model, codon optimization yielded ideal GC content and a higher codon adaptation index value, which was then included in the cloning vector pET2+ (a). Altogether, our experiment implies that the proposed peptide vaccine might be a good option for vibriosis prophylaxis.

Analysis of codon usage bias of classical swine fever virus.

Background and Aim:Classical swine fever (CSF), caused by CSF virus (CSFV), is a highly contagious disease in pigs causing 100% mortality in susceptible adult pigs and piglets. High mortality rate in pigs causes huge economic loss to pig farmers. CSFV has a positive-sense RNA genome of 12.3 kb in length flanked by untranslated regions at 5’ and 3’ end. The genome codes for a large polyprotein of 3900 amino acids coding for 11 viral proteins. The 1300 codons in the polyprotein are coded by different combinations of three nucleotides which help the infectious agent to evolve itself and adapt to the host environment. This study performed and employed various methods/techniques to estimate the changes occurring in the process of CSFV evolution by analyzing the codon usage pattern.Materials and Methods:The evolution of viruses is widely studied by analyzing their nucleotides and coding regions/codons using various methods. A total of 115 complete coding regions of CSFVs including one complete genome from our laboratory (MH734359) were included in this study and analysis was carried out using various methods in estimating codon usage bias and evolution. This study elaborates on the factors that influence the codon usage pattern.Results:The effective number of codons (ENC) and relative synonymous codon usage showed the presence of codon usage bias. The mononucleotide (A) has a higher frequency compared to the other mononucleotides (G, C, and T). The dinucleotides CG and CC are underrepresented and overrepresented. The codons CGT was underrepresented and AGG was overrepresented. The codon adaptation index value of 0.71 was obtained indicating that there is a similarity in the codon usage bias. The principal component analysis, ENC-plot, Neutrality plot, and Parity Rule 2 plot produced in this article indicate that the CSFV is influenced by the codon usage bias. The mutational pressure and natural selection are the important factors that influence the codon usage bias.Conclusion:The study provides useful information on the codon usage analysis of CSFV and may be utilized to understand the host adaptation to virus environment and its evolution. Further, such findings help in new gene discovery, design of primers/probes, design of transgenes, determination of the origin of species, prediction of gene expression level, and gene function of CSFV. To the best of our knowledge, this is the first study on codon usage bias involving such a large number of complete CSFVs including one sequence of CSFV from India.

Codon usage analysis of zoonotic coronaviruses reveals lower adaptation to humans by SARS-CoV-2.

Since 2002, the world has witnessed major outbreaks of acute respiratory illness by three zoonotic coronaviruses (CoVs), which differ from each other in pathogenicity. Reasons for the lower pathogenicity of SARS-CoV-2 than the other two zoonotic coronaviruses, SARS-CoV and MERS-CoV, are not well understood. We herein compared the codon usage patterns of the three zoonotic CoVs causing severe acute respiratory syndromes and four human-specific CoVs (NL63, 229E, OC43, and HKU1) causing mild diseases. We found that the seven viruses have different codon usages, with SARS-CoV-2 having the lowest effective number of codons (ENC) among the zoonotic CoVs. Human codon adaptation index (CAI) analysis revealed that the CAI value of SARS-CoV-2 is the lowest among the zoonotic CoVs. The ENC and CAI values of SARS-CoV-2 were more similar to those of the less-pathogenic human-specific CoVs. To further investigate adaptive evolution within SARS-CoV-2, we examined codon usage patterns in 3573 genomes of SARS-CoV-2 collected over the initial 4 months of the pandemic. We showed that the ENC values and the CAI values of SARS-CoV-2 were decreasing over the period. The low ENC and CAI values could be responsible for the lower pathogenicity of SARS-CoV-2. While mutational pressure appears to shape codon adaptation in the overall genomes of SARS-CoV-2 and other zoonotic CoVs, the E gene of SARS-CoV-2, which has the highest codon usage bias, appears to be under strong natural selection. Data from the study contribute to our understanding of the pathogenicity and evolution of SARS-CoV-2 in humans.

The characteristic of the synonymous codon usage and phylogenetic analysis of hepatitis B virus

BackgroundHepatitis B virus (HBV) infection is a crucial medical issue worldwide. The dependence of HBV replication on host cell machineries and their co-evolutionary interactions prompt the codon usage pattern of viral genes to translation selection and mutation pressure.ObjectiveThe evolutionary characteristics of HBV and the natural selection effects of the human genome on the codon usage characteristics were analyzed to provide a basis for medication development for HBV infection.MethodsThe codon usage pattern of sequences from different HBV genotypes of our isolates and reference HBV genome sequences downloaded from the National Center for Biotechnology Information (NCBI) database were analyzed by computing the relative synonymous codon usage (RSCU), nucleotide content, codon adaptation index (CAI) and the effective number of codons (ENC).ResultsThe highest ENC values were observed in the C genotypes, followed by the B genotypes. The ENC values indicated a weak codon usage bias (CUB) in HBV genome. The number of codons differentially used between the three genotypes was markedly higher than that of similarly used codons. High CAI values indicated a good adaptability of HBV to its host. The ENC plot indicated the occurrence of mutational pressure in the three genotypes. The mean Ka/Ks ratios in the three genotypes were lower than 1, which indicated a negative selection pressure. The CAI and GC3% plot indicated the existence of CUB in the HBV genome.ConclusionsNucleotide composition, mutation bias, negative selection and mutational pressure are key factors influencing the CUB and phylogenetic diversity in HBV genotypes. The data provided here could be useful for developing drugs for HBV infection.

Gene transfer and nucleotide sequence evolution by Gossypium cytoplasmic genomes indicates novel evolutionary characteristics.

The DNA fragments transferred among cotton cytoplasmic genomes are highly differentiated. The wild D group cotton species have undergone much greater evolution compared with cultivated AD group. Cotton (Gossypium spp.) is one of the most economically important fiber crops worldwide. Gene transfer, nucleotide evolution, and the codon usage preferences in cytoplasmic genomes are important evolutionary characteristics of high plants. In this study, we analyzed the nucleotide sequence evolution, codon usage, and transfer of cytoplasmic DNA fragments in Gossypium chloroplast (cp) and mitochondrial (mt) genomes, including the A genome group, wild D group, and cultivated AD group of cotton species. Our analyses indicated that the differences in the length of transferred cytoplasmic DNA fragments were not significant in mitochondrial and chloroplast sequences. Analysis of the transfer of tRNAs found that trnQ and nine other tRNA genes were commonly transferred between two different cytoplasmic genomes. The Codon Adaptation Index values showed that Gossypium cp genomes prefer A/T-ending codons. Codon preference selection was higher in the D group than the other two groups. Nucleotide sequence evolution analysis showed that intergenic spacer sequences were more variable than coding regions and nonsynonymous mutations were clearly more common in cp genomes than mt genomes. Evolutionary analysis showed that the substitution rate was much higher in cp genomes than mt genomes. Interestingly, the D group cotton species have undergone much faster evolution compared with cultivated AD groups, possibly due to the selection and domestication of diverse cotton species. Our results demonstrate that gene transfer and differential nucleotide sequence evolution have occurred frequently in cotton cytoplasmic genomes.

Expression of Lon-like Protease Gene from Lactobacillus plantarum IIA-1A5 in Escherichia coli BL21(DE3)

Proteases are one of most important and abundant enzymes produced by the biotechnology industry, for scientific, physiological and industrial application and dominates of the whole enzyme market. Lactobacillus plantarum IIA-1A5 is an Indonesian lactic acid bacteria (LAB) isolated from beef Peranakan Ongole cattle. Preliminary analysis on its whole genome sequence indicated that this strain harbours some genes involved in protein degradation and might be promising to be further applied. This study aims to optimize the gene sequence of a lon-like protease of L. plantarum IIA-1A5 for heterologous expression system. The Lon-like gene expression system is made using genes that have been optimized first in silico. pET-28a(+), E. coli BL21(DE3), Nde1 and BamH1 were used in this study as a expression vector, a host and retriction enzyme, respectively. Molecular weight was validated using SDS-PAGE and expasy.org software. The results showed that optimization increased codon adaptation index value (CAI) and GC content to 0.97 and 56.57%, respectively, which were suitable for the E. coli expression system. The Lon-like IIA gene was successfully expressed in the cell cytoplasm by induction of 1 mM Isopropyl β-D-1-thiogalactopyranoside (IPTG) at 37 °C. As many as 88% of Lon-like IIA codons were distributed in the 91-100 quality group. Lon-like IIA was successfully expressed in a host cell induced with 1 mM IPTG at 37oC . IPTG induction was performed at the 3rd hour of incubation with OD600 0.59. In addition, Lon-like IIA molecular weight was detected approximately 43 kDa.

Exploring the codon patterns between CCD and NCED genes among different plant species

Plant carotenoid cleavage oxygenase (CCO) is an enzyme which catalyzes carotenoids to apocarotenoid products that are involved in several vital physiological functions. The CCO exists in two forms, namely, CCD (Carotenoid Cleavage Dioxygenase) and NCED (Nine-Cis Epoxycarotenoid Dioxygenase). This paper relates to a comparative study on CCD and NCED genes through phylogeny and codon usage analysis. The result of the phylogenetic analysis indicates a closer relationship between CCD and NCED subclass genes, while the RSCU values indicate a high preference for CUC codon in both CCD and NCED gene families. The mean ENc value of NCED genes was found to be 48.76, suggesting a higher codon bias compared to CCD genes. However, the ENc-GC3S plot suggests that both the gene families are under mutational pressure with variations according to their species-specific role. Similarly, the multivariate analysis also suggests that nucleotide mutation bias influences codon usage. Correlation analysis of Axis I and codon adaptation index values indicate a significant correlation between critical indices. Even though the prominence of the variations in codon usage between the two gene families, they are exerted towards the time-specific functional requirement for that plant species. This is evident from the cleaving roles of these enzymes against various carotenoids at different growth stages. The result of this investigation indicates that CCD and NCED genes are under mutational pressure. This codon bias study paves the way for increasing the production of apocarotenoids, which have a great significance in the industry.

Genetic Analysis and Evolutionary Changes of the Torque teno sus Virus.

The torque teno sus virus (TTSuV) is an emerging virus threating the Suidae species of unclear pathogenicity, although it was previously reported as a worsening factor of other porcine diseases, in particular, porcine circovirus associated disease (PCVAD). Here, a comprehensive codon usage analysis of the open reading frame 1 (ORF1), which encodes the viral capsid protein, was undertaken for the first time to reveal its evolutionary history. We revealed independent phylogenetic processes for the two genera during TTSuV evolution, which was confirmed by principal component analysis (PCA). A low codon usage bias was observed in different genera and different species, with Kappatorquevirus a (TTSuVk2a) displaying the highest, which was mainly driven by mutation pressure and natural selection, especially natural selection. Overall, ATs were more abundant than GCs, along with more A-ended synonymous codons in relative synonymous codon usage (RSCU) analysis. To further confirm the role of natural selection and TTSuV adaptation to the Suidae species, codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses were performed, which showed different adaptations for different TTSuVs. Importantly, we identified a more dominant role of Sus scrofa in the evolution of Iotatorquevirus (TTSuV1), with the highest CAI values and lowest RCDI values compared to Sus scrofa domestica. However, in TTSuVk2, the roles of Sus scrofa and Sus scrofa domestica were the same, regarding codon usage, with similar CAI and RCDI values. Our study provides a new perspective of the evolution of TTSuV and valuable information to develop control measures against TTSuV.

Codon Usage Bias Analysis of Citrus tristeza Virus: Higher Codon Adaptation to Citrus reticulata Host.

Citrus tristeza virus (CTV), a member of the aphid-transmitted closterovirus group, is the causal agent of the notorious tristeza disease in several citrus species worldwide. The codon usage patterns of viruses reflect the evolutionary changes for optimization of their survival and adaptation in their fitness to the external environment and the hosts. The codon usage adaptation of CTV to specific citrus hosts remains to be studied; thus, its role in CTV evolution is not clearly comprehended. Therefore, to better explain the host–virus interaction and evolutionary history of CTV, the codon usage patterns of the coat protein (CP) genes of 122 CTV isolates originating from three economically important citrus hosts (55 isolate from Citrus sinensis, 38 from C. reticulata, and 29 from C. aurantifolia) were studied using several codon usage indices and multivariate statistical methods. The present study shows that CTV displays low codon usage bias (CUB) and higher genomic stability. Neutrality plot and relative synonymous codon usage analyses revealed that the overall influence of natural selection was more profound than that of mutation pressure in shaping the CUB of CTV. The contribution of high-frequency codon analysis and codon adaptation index value show that CTV has host-specific codon usage patterns, resulting in higheradaptability of CTV isolates originating from C. reticulata (Cr-CTV), and low adaptability in the isolates originating from C. aurantifolia (Ca-CTV) and C. sinensis (Cs-CTV). The combination of codon analysis of CTV with citrus genealogy suggests that CTV evolved in C. reticulata or other Citrus progenitors. The outcome of the study enhances the understanding of the factors involved in viral adaptation, evolution, and fitness toward their hosts. This information will definitely help devise better management strategies of CTV.

Comparative genome‐wide analysis of codon usage of different bacterial species infecting Oryza sativa

Oryza sativa is vastly affected by microbial pathogen, causing blight-related diseases, which in turn deplete the growth and productivity of rice. In this study, we analyzed four bacterial rice pathogen genomes and reported on their codon usage that might have greater implication in mutation-related research. Differential codon usage indices, such as codon adaptation index (CAI), codon bias index (CBI), effective number of codons (ENc), relative synonymous codon usage (RSCU), correspondence analysis (COA), and parity plots, were applied on coding sequences of Pseudomonas fuscovaginae, Pseudomonas syringae, Xanthomonas oryzae, and Pseudomonas avenae speices. The RSCU results proposed a high-frequency usage of CUG and CGC that codes for leucine and arginine in all of the species. The CBI and CAI values between the genomes range from 0.17 to 0.3 and from 0.26 to 0.35, respectively, indicating a direct proportionality between these indexes. The mean ENc value of P. avenae coding sequence showed high codon bias compared with other genomes. The axis I variation from COA analysis shows a mean value of 42.28% codon variations in these bacterial species. Correlation studies between axis I and ENc-GC3, along with CAI and CBI, suggested the presence of nucleotide bias and mutational pressure as major forces for codon bias within these species. Hence, certain genes with high CAI-CBI have been correlated for better gene expression. Our study highlights the importance of nucleotide biasness, mutation pressure, and natural selection in shaping protein-coding genes in these four rice-affecting bacteria. This would further help in investigating the evolution of pathogenic gene families, which may direct research toward synthetic genes that could be suppressed or overrepresented based on their codon usage pattern toward pathogenicity.

Rice black-streaked dwarf virus Genome in China: Diversification, Phylogeny, and Selection.

Rice black-streaked dwarf virus (RBSDV), a Fijivirus, causes maize rough dwarf disease and rice black-streaked dwarf disease in the summer maize-growing regions of the Yellow and Huai rivers, respectively, in China. Nevertheless, the diversification and selection of the entire genome from S1 to S10 have not been illuminated. Molecular variation, evolution, conserved regions, and other genomic properties were analyzed in 21 RBSDV isolates from maize (Zea mays L.) and rice (Oryza sativa) hosts sampled from nine geographic locations in China. Low codon adaptation index values ranging from 0.1878 to 0.2918 indicated a low degree of codon-usage bias and low potential expression for all 13 RBSDV open reading frames (ORFs). ORF9-2 showed a stronger effect of codon usage bias than did other ORFs, as the majority of points for this ORF lay close to the standard curve in the Nc plot (the effective number of codons [Nc] versus the frequency of G+C at synonymous third-base positions [GC3]). A 9-bp deletion mutation was detected in the RBSDV genome in the 3' UTR of S8. Nucleotide diversity analysis indicated that the structural proteins of RBSDV, such as S2 and S4, were all more conserved than nonstructural proteins such as S9. Nucleotide diversity (π) was highest among S9 sequences (0.0656), and was significantly higher than among S4 sequences (0.0225, P < 0.01). The number of conserved regions among the 10 segments varied substantially. The highest number of conserved regions (5) was found in S5, whereas no conserved regions were identified in S9. Nucleotide diversity and the number of conserved regions were independent of the lengths of segments. Nucleotide diversity was also not correlated with the number of conserved regions in segments. Ten recombination events in 21 isolates were found in seven segments with breakpoint positions in UTRs, intergenic spacer regions, and gene coding regions. The number of recombination events was also independent of the lengths of segments. RBSDV isolates from China could be phylogenetically classified into two groups using either 10 segment sequences or the concatenated sequence of S1 through S10, regardless of host or geographical location. The phylogenetic tree generated from pairwise nucleotide identities of individual RBSDV segments such as S9 and S3, with nucleotide identity values of 93.74% and 95.86%, respectively, is similar to the tree constructed from the concatenated sequences of the entire RBSDV genome. The 13 RBSDV ORFs were under negative and purifying selection (Ka/Ks < 1). ORF5-2 was under the greatest selection pressure; however, ORF2, which encodes the core protein of RBSDV, was under the lowest selection pressure.

Novel methods to optimize gene and statistic test for evaluation \u2013 an application for Escherichia coli

BackgroundSince the recombinant protein was discovered, it has become more popular in many aspects of life science. The value of global pharmaceutical market was $87 billion in 2008 and the sales for industrial enzyme exceeded $4 billion in 2012. This is strong evidence showing the great potential of recombinant protein. However, native genes introduced into a host can cause incompatibility of codon usage bias, GC content, repeat region, Shine-Dalgarno sequence with host’s expression system, so the yields can fall down significantly. Hence, we propose novel methods for gene optimization based on neural network, Bayesian theory, and Euclidian distance.ResultThe correlation coefficients of our neural network are 0.86, 0.73, and 0.90 in training, validation, and testing process. In addition, genes optimized by our methods seem to associate with highly expressed genes and give reasonable codon adaptation index values. Furthermore, genes optimized by the proposed methods are highly matched with the previous experimental data.ConclusionThe proposed methods have high potential for gene optimization and further researches in gene expression. We built a demonstrative program using Matlab R2014a under Mac OS X. The program was published in both standalone executable program and Matlab function files. The developed program can be accessed from http://www.math.hcmus.edu.vn/~ptbao/paper_soft/GeneOptProg/.