Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans

Precise CRISPR/Cpf1 genome editing system in the Deinococcus radiodurans with superior DNA repair mechanisms

Deinococcus radiodurans, with its high homologous recombination (HR) efficiency of double-stranded DNA breaks (DSBs), is a model organism for studying genome stability maintenance and an attractive microbe for industrial applications. Here, we developed an efficient CRISPR/Cpf1 genome editing system in D. radiodurans by evaluating and optimizing double-plasmid strategies and four Cas effector proteins from various organisms, which can precisely introduce different types of template-dependent mutagenesis without off-target toxicity. Furthermore, the role of DNA repair genes in determining editing efficiency in D. radiodurans was evaluated by introducing the CRISPR/Cpf1 system into 13 mutant strains lacking various DNA damage response and repair factors. In addition to the crucial role of RecA-dependent HR required for CRISPR/Cpf1 editing, D. radiodurans showed higher editing efficiency when lacking DdrB, the single-stranded DNA annealing (SSA) protein involved in the RecA-independent DSB repair pathway. This suggests a possible competition between HR and SSA pathways in the CRISPR editing of D. radiodurans. Moreover, off-target effects were observed during the genome editing of the pprI knockout strain, a master DNA damage response gene in Deinococcus species, which suggested that precise regulation of DNA damage response is critical for a high-fidelity genome editing system.

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

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

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Analysis of codon usage patterns in Haloxylon ammodendron based on genomic and transcriptomic data

Haloxylon ammodendron, a xero-halophytic shrub of Chenopodiaceae, is a dominant species in deserts, which has a strong drought and salt tolerance and plays an important role in sand fixation. However, the codon usage bias (CUB) in H. ammodendron is still unclear at present. In this study, the codon usage patterns of 38,657 coding sequences (CDSs) in the newly released whole-genome sequence data of H. ammodendron and 3,948 CDSs in the previously obtained transcriptome sequencing data were compared and analyzed. The results showed that the CDSs with the total guanineandcytosine(GC)content in the range of 40% ∼ 45% was the most in the genome and transcriptome. Among which, the GC1, GC2, and GC3 contents of genomic CDSs were 50.83%, 40.56%, and 40.23%, respectively, and those of CDSs in the transcriptome were 47.16%, 39.02%, and 39.59%, respectively. Therefore, the bases in H. ammodendron were rich in adenine and thymine, and the overallcodonusage was biasedtoward A- and U-ending codons. The analysis of neutrality plot, effective number of codon (ENC) plot, and parity rule 2 (PR2) bias plot showed that both natural selection and mutation pressure had great influences on the CUB of H. ammodendron, but natural selection was the most important determinant. Besides, gene expression level and the function and protein length of some specific genes also had influences on the codon usage pattern. Finally, a total of 25 common optimal codons were found in the genomic and transcriptomic data, and AU/GC-ending codons ratio was 24:1. It should be noted that the salt-tolerant unigenes had similar codon usage, and the highly expressed genes had higher usage frequency of optimal codons and lower GC content than the lowly expressed genes. In addition, there was no difference in the ENC values of salt-tolerant unigenes in H. ammodendron, and the expression level of the genes had no correlation with CAI. This study will help to elucidate the formation mechanism of H. ammodendron codon usage bias, and make contributions to the identification of new genes and the genetic engineering study on H. ammodendron.

Whole genome sequencing and comparative genomic analyses of Pseudomonas aeruginosa strain isolated from arable soil reveal novel insights into heavy metal resistance and codon biology

Elevated concentration of non-essential persistent heavy metals and metalloids in the soil is detrimental to essential soil microbes and plants, resulting in diminished diversity and biomass. Thus, isolation, screening, and whole genomic analysis of potent strains of bacteria from arable lands with inherent capabilities of heavy metal resistance and plant growth promotion hold the key for bio remedial applications. This study is an attempt to do the same. In this study, a potent strain of Pseudomonas aeruginosa was isolated from paddy fields, followed by metabolic profiling using FTIR, metal uptake analysis employing ICP-MS, whole genome sequencing and comparative codon usage analysis. ICP-MS study provided insights into a high degree of Cd uptake during the exponential phase of growth under cumulative metal stress to Cd, Zn and Co, which was further corroborated by the detection of cadA gene along with czcCBA operon in the genome upon performing whole-genome sequencing. This potent strain of Pseudomonas aeruginosa also harboured genes, such as copA, chrA, znuA, mgtE, corA, and others conferring resistance against different heavy metals, such as Cd, Zn, Co, Cu, Cr, etc. A comparative codon usage bias analysis at the genomic and genic level, whereby several heavy metal resistant genes were considered in the backdrop of two housekeeping genes among 40 Pseudomonas spp. indicated the presence of a relatively strong codon usage bias in the studied strain. With this work, an effort was made to explore heavy metal-resistant bacteria (isolated from arable soil) and whole genome sequence analysis to get insight into metal resistance for future bio remedial applications.

Influencing elements of codon usage bias in Birnaviridae and its evolutionary analysis

Birnaviridae is a family of double stranded (ds) RNA virus with non-enveloped virions and 2-segmented genomes. These viruses are known to cause diseases in many hosts. Virus of this family has affected the fish and poultry economy in a wide sector. Unevenness in the use of synonymous codons for a particular amino acid in the coding strand of DNA is known as codon usage bias (CUB). Codons that code the same amino acid are used with variable frequency in a variety of life forms. To understand the pattern of CUB in Birnaviridae, we carried out bioinformatics study to understand the properties of coding sequences of proteins. ENC value of Birnaviridae suggested low CUB. Nucleotide analysis revealed high GC content. Parameters such as RSCU values, nucleotide skewness, translational selection, parity plot and neutrality plot were studied to investigate the pattern of codon use and it was clear that both mutational pressure and natural selection contributed to the designing of CUB in Birnaviridae family. The neutrality plot revealed natural selection to dominate the structuring of CUB and hence remained the major CUB determinant in Birnaviridae. Outcome of our study exemplified the pattern of codon use in the Birnaviridae genomes and contributed the basic primary data for fundamental evolutionary research on them.

Comparative analysis of transcriptomic data shows the effects of multiple evolutionary selection processes on codon usage in Marsupenaeus japonicus and Marsupenaeus pulchricaudatus

BACKGROUND

Kuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data.

RESULTS

Using CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω > 1) had a higher correlation coefficient than that of experienced purifying selection (ω < 1). Parity Rule 2 (PR2) and effective number of codons (ENc) plot analysis showed that the codon usage patterns of both species were influenced by both mutations and selection. Moreover, the average observed ENc value was lower than the expected value for both species, suggesting that factors other than GC may play roles in these phenomena. The results of multispecies clustering based on codon preference were consistent with traditional classification.

CONCLUSIONS

This study provides a relatively comprehensive understanding of the correlations among codon usage bias, gene expression, and selection pressures of CDSs for M. japonicus and M. pulchricaudatus. The genetic evolution was driven by mutations and selection pressure. Moreover, the results point out new insights into the specificities and evolutionary characteristics of the two Marsupenaeus species.

Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB

Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite-induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duodenalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G. duodenalis CUB formation. This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB. These results will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

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.

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.

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.

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.

Codon usage of human hepatitis C virus clearance genes in relation to its expression

Hepatitis C virus (HCV) infection is among the leading causes of hepatocellular carcinoma and liver cirrhosis globally, with a high economic burden. The disease progression is well established, but less is known about the spontaneous HCV infection clearance. This study tries to establish the relationship between codon biasness and expression of HCV clearance candidate genes in normal and HCV infected liver tissues. A total of 112 coding sequences comprising 151 679 codons were subjected to the computation of codon indices, namely relative synonymous codon usage, an effective number of codon (Nc), frequency of optimal codon, codon adaptation index, codon bias index, and base compositions. Codon indices report of GC3s, GC12, hydropathicity, and aromaticity implicates both mutational and translational selection in the candidate gene set. This was further correlated with the differentially expressed genes among the selected genes using BioGPS. A significant correlation is observed between the gene expression of normal liver and cancerous liver tissues with codon bias (Nc). Gene expression is also correlated with relative codon bias values, indicating that CCL5, APOA2, CD28, IFITM1, and TNFSF4 genes have higher expression. These results are quite encouraging in selecting the high responsive genes in HCV clearance. However, there could be additional genes which could also orchestrate the clearance role with the above mentioned first line of defensive genes.

Differences in Codon Usage Bias between Photosynthesis-Related Genes and Genetic System-Related Genes of Chloroplast Genomes in Cultivated and Wild Solanum Species

Solanum is one of the largest genera, including two important crops-potato (Solanum tuberosum) and tomato (Solanum lycopersicum). In this study we compared the chloroplast codon usage bias (CUB) among 12 Solanum species, between photosynthesis-related genes (Photo-genes) and genetic system-related genes (Genet-genes), and between cultivated species and wild relatives. The Photo-genes encode proteins for photosystems, the photosynthetic electron transport chain, and RuBisCO, while the Genet-genes encode proteins for ribosomal subunits, RNA polymerases, and maturases. The following findings about the Solanum chloroplast genome CUB were obtained: (1) the nucleotide composition, gene expression, and selective pressure are identified as the main factors affecting chloroplast CUB; (2) all these 12 chloroplast genomes prefer A/U over G/C and pyrimidines over purines at the third-base of codons; (3) Photo-genes have higher codon adaptation indexes than Genet-genes, indicative of a higher gene expression level and a stronger adaptation of Photo-genes; (4) gene function is the primary factor affecting CUB of Photo-genes but not Genet-genes; (5) Photo-genes prefer pyrimidine over purine, whereas Genet-genes favor purine over pyrimidine, at the third position of codons; (6) Photo-genes are mainly affected by the selective pressure, whereas Genet-genes are under the underlying mutational bias; (7) S. tuberosum is more similar with Solanum commersonii than with Solanum bulbocastanum; (8) S. lycopersicum is greatly different from the analyzed seven wild relatives; (9) the CUB in codons for valine, aspartic acid, and threonine are the same between the two crop species, S. tuberosum and S. lycopersicum. These findings suggest that the chloroplast CUB contributed to the differential requirement of gene expression activity and function between Photo-genes and Genet-genes and to the performance of cultivated potato and tomato.

Codon Usage Pattern of Genes Involved in Central Nervous System

Codon usage bias (CUB) is the non-uniform usage of synonymous codons in which some codons are more preferred to others in the transcript. Analysis of codon usage bias has applications in understanding the basics of molecular biology, genetics, gene expression, and molecular evolution. To understand the patterns of codon usage in genes involved in the central nervous system (CNS), we used bioinformatic approaches to analyze the protein-coding sequences of genes involved in the CNS. The improved effective number of codons (ENC) suggested that the overall codon usage bias was low. The relative synonymous codon usage (RSCU) revealed that the most frequently occurring codons had a G or C at the third codon position. The codons namely TCC, AGC, CTG, CAG, CGC, ATC, ACC, GTG, GCC, GGC, and CGG (average RSCU > 1.6) were over-represented. Both mutation pressure and natural selection might affect the codon usage pattern as evident from correspondence and parity plot analyses. The overall GC content (59.93) was higher than AT content, i.e., genes were GC-rich. The correlation of GC12 with GC3 suggested that mutation pressure might affect the codon usage pattern.

Comprehensive Analysis and Comparison on the Codon Usage Pattern of Whole Mycobacterium tuberculosis Coding Genome from Different Area

Phenomenon of unequal use of synonymous codons in Mycobacterium tuberculosis is common. Codon usage bias not only plays an important regulatory role at the level of gene expression, but also helps in improving the accuracy and efficiency of translation. Meanwhile, codon usage pattern of Mycobacterium tuberculosis genome is important for interpreting evolutionary characteristics in species. In order to investigate the codon usage pattern of the Mycobacterium tuberculosis genome, 12 Mycobacterium tuberculosis genomes from different area are downloaded from the GeneBank. The correlations between G₃, GC₁₂, whole GC content, codon adaptation index, codon bias index, and so on of Mycobacterium tuberculosis genomes are calculated. The ENC-plot, relationship between A₃/(A₃ + T₃) and G₃/(G₃ + C₃), GC₁₂ versus GC₃ plot, and the RSCU of overall/separated genomes all show that the codon usage bias exists in all 12 Mycobacterium tuberculosis genomes. Lastly, relationship between CBI and the equalization of ENC shows a strong negative correlation between them. The relationship between protein length and GC content (GC₃ and GC₁₂) shows that more obvious differences in the GC content may be in shorter protein. These results show that codon usage bias existing in the Mycobacterium tuberculosis genomes could be used for further study on their evolutionary phenomenon.

A modified GC-specific MAKER gene annotation method reveals improved and novel gene predictions of high and low GC content in Oryza sativa

BACKGROUND

Accurate structural annotation depends on well-trained gene prediction programs. Training data for gene prediction programs are often chosen randomly from a subset of high-quality genes that ideally represent the variation found within a genome. One aspect of gene variation is GC content, which differs across species and is bimodal in grass genomes. When gene prediction programs are trained on a subset of grass genes with random GC content, they are effectively being trained on two classes of genes at once, and this can be expected to result in poor results when genes are predicted in new genome sequences.

RESULTS

We find that gene prediction programs trained on grass genes with random GC content do not completely predict all grass genes with extreme GC content. We show that gene prediction programs that are trained with grass genes with high or low GC content can make both better and unique gene predictions compared to gene prediction programs that are trained on genes with random GC content. By separately training gene prediction programs with genes from multiple GC ranges and using the programs within the MAKER genome annotation pipeline, we were able to improve the annotation of the Oryza sativa genome compared to using the standard MAKER annotation protocol. Gene structure was improved in over 13% of genes, and 651 novel genes were predicted by the GC-specific MAKER protocol.

CONCLUSIONS

We present a new GC-specific MAKER annotation protocol to predict new and improved gene models and assess the biological significance of this method in Oryza sativa. We expect that this protocol will also be beneficial for gene prediction in any organism with bimodal or other unusual gene GC content.

Analysis of transcriptome data reveals multifactor constraint on codon usage in Taenia multiceps

Background

Codon usage bias (CUB) is an important evolutionary feature in genomes that has been widely observed in many organisms. However, the synonymous codon usage pattern in the genome of T. multiceps remains to be clarified. In this study, we analyzed the codon usage of T. multiceps based on the transcriptome data to reveal the constraint factors and to gain an improved understanding of the mechanisms that shape synonymous CUB.

Results

Analysis of a total of 8,620 annotated mRNA sequences from T. multiceps indicated only a weak codon bias, with mean GC and GC3 content values of 49.29% and 51.43%, respectively. Our analysis indicated that nucleotide composition, mutational pressure, natural selection, gene expression level, amino acids with grand average of hydropathicity (GRAVY) and aromaticity (Aromo) and the effective selection of amino-acids all contributed to the codon usage in T. multiceps. Among these factors, natural selection was implicated as the major factor affecting the codon usage variation in T. multiceps. The codon usage of ribosome genes was affected mainly by mutations, while the essential genes were affected mainly by selection. In addition, 21codons were identified as “optimal codons”. Overall, the optimal codons were GC-rich (GC:AU, 41:22), and ended with G or C (except CGU). Furthermore, different degrees of variation in codon usage were found between T. multiceps and Escherichia coli, yeast, Homo sapiens. However, little difference was found between T. multiceps and Taenia pisiformis.

Conclusions

In this study, the codon usage pattern of T. multiceps was analyzed systematically and factors affected CUB were also identified. This is the first study of codon biology in T. multiceps. Understanding the codon usage pattern in T. multiceps can be helpful for the discovery of new genes, molecular genetic engineering and evolutionary studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3704-8) contains supplementary material, which is available to authorized users.

Analysis of codon usage patterns in Morus notabilis based on genome and transcriptome data

Codons play important roles in regulating gene expression levels and mRNA half-lives. However, codon usage and related studies in multicellular organisms still lag far behind those in unicellular organisms. In this study, we describe for the first time genome-wide patterns of codon bias in Morus notabilis (mulberry tree), and analyze genome-wide codon usage in 12 other species within the order Rosales. The codon usage of M. notabilis was affected by nucleotide composition, mutation pressure, nature selection, and gene expression level. Translational selection optimal codons were identified and highly expressed genes of M. notabilis tended to use the optimal codons. Genes with higher expression levels have shorter coding region and lower amino acid complexity. Housekeeping genes showed stronger translational selection, which, notably, was not caused by the large differences between the expression level of housekeeping genes and other genes.

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

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

Codon Usage Bias and Determining Forces in Taenia solium Genome

The tapeworm Taenia solium is an important human zoonotic parasite that causes great economic loss and also endangers public health. At present, an effective vaccine that will prevent infection and chemotherapy without any side effect remains to be developed. In this study, codon usage patterns in the T. solium genome were examined through 8,484 protein-coding genes. Neutrality analysis showed that T. solium had a narrow GC distribution, and a significant correlation was observed between GC12 and GC3. Examination of an NC (ENC vs GC3s)-plot showed a few genes on or close to the expected curve, but the majority of points with low-ENC (the effective number of codons) values were detected below the expected curve, suggesting that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally. We also identified 26 optimal codons in the T. solium genome, all of which ended with either a G or C residue. These optimal codons in the T. solium genome are likely consistent with tRNAs that are highly expressed in the cell, suggesting that mutational and translational selection forces are probably driving factors of codon usage bias in the T. solium genome.

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

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

Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes

The sub-3 Mbp genomes from microsporidian species of the Encephalitozoon genus are the smallest known among eukaryotes and paragons of genomic reduction and compaction in parasites. However, their diminutive stature is not characteristic of all Microsporidia, whose genome sizes vary by an order of magnitude. This large variability suggests that different evolutionary forces are applied on the group as a whole. In this study, we have compared the codon usage bias (CUB) between eight taxonomically distinct microsporidian genomes: Encephalitozoon intestinalis, Encephalitozoon cuniculi, Spraguea lophii, Trachipleistophora hominis, Enterocytozoon bieneusi, Nematocida parisii, Nosema bombycis and Nosema ceranae. While the CUB was found to be weak in all eight Microsporidia, nearly all (98%) of the optimal codons in S. lophii, T. hominis, E. bieneusi, N. parisii, N. bombycis and N. ceranae are fond of A/U in third position whereas most (64.6%) optimal codons in the Encephalitozoon species E. intestinalis and E. cuniculi are biased towards G/C. Although nucleotide composition biases are likely the main factor driving the CUB in Microsporidia according to correlation analyses, directed mutational pressure also likely affects the CUB as suggested by ENc-plots, correspondence and neutrality analyses. Overall, the Encephalitozoon genomes were found to be markedly different from the other microsporidians and, despite being the first sequenced representatives of this lineage, are uncharacteristic of the group as a whole. The disparities observed cannot be attributed solely to differences in host specificity and we hypothesize that other forces are at play in the lineage leading to Encephalitozoon species.

Analysis of codon usage pattern in Taenia saginata based on a transcriptome dataset

BACKGROUND

Codon usage bias is an important evolutionary feature in a genome and has been widely documented in many genomes. Analysis of codon usage bias has significance for mRNA translation, design of transgenes, new gene discovery, and studies of molecular biology and evolution, etc. However, the information about synonymous codon usage pattern of T. saginata genome remains unclear. T. saginata is a food-borne zoonotic cestode which infects approximataely 50 million humans worldwide, and causes significant health problems to the host and considerable socio-economic losses as a consequence. In this study, synonymous codon usage in T. saginata were examined.

METHODS

Total RNA was isolated from T. saginata cysticerci and 91,487 unigenes were generated using Illumina sequencing technology. After filtering, the final sequence collection containing 11,399 CDSs was used for our analysis.

RESULTS

Neutrality analysis showed that the T. saginata had a wide GC3 distribution and a significant correlation was observed between GC12 and GC3. NC-plot showed most of genes on or close to the expected curve, but only a few points with low-ENC values were below it, suggesting that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally. We also identified twenty-three optimal codons in the T. saginata genome, all of which were ended with a G or C residue. These results suggest that mutational and selection forces are probably driving factors of codon usage bias in T. saginata genome. Meanwhile, other factors such as protein length, gene expression, GC content of genes, the hydropathicity of each protein also influence codon usage.

CONCLUSIONS

Here, we systematically analyzed the codon usage pattern and identified factors shaping in codon usage bias in T. saginata. Currently, no complete nuclear genome is available for codon usage analysis at the genome level in T. saginata. This is the first report to investigate codon biology in T. sagninata. Such information does not only bring about a new perspective for understanding the mechanisms of biased usage of synonymous codons but also provide useful clues for molecular genetic engineering and evolutionary studies.

A comparative analysis of synonymous codon usage bias pattern in human albumin superfamily

Synonymous codon usage bias is an inevitable phenomenon in organismic taxa across the three domains of life. Though the frequency of codon usage is not equal across species and within genome in the same species, the phenomenon is non random and is tissue-specific. Several factors such as GC content, nucleotide distribution, protein hydropathy, protein secondary structure, and translational selection are reported to contribute to codon usage preference. The synonymous codon usage patterns can be helpful in revealing the expression pattern of genes as well as the evolutionary relationship between the sequences. In this study, synonymous codon usage bias patterns were determined for the evolutionarily close proteins of albumin superfamily, namely, albumin, α-fetoprotein, afamin, and vitamin D-binding protein. Our study demonstrated that the genes of the four albumin superfamily members have low GC content and high values of effective number of codons (ENC) suggesting high expressivity of these genes and less bias in codon usage preferences. This study also provided evidence that the albumin superfamily members are not subjected to mutational selection pressure.

Analysis of Codon Usage Patterns in Toxic Dinoflagellate Alexandrium tamarense through Expressed Sequence Tag Data

We have analyzed synonymous codon usage in the genome of A. tamarense CCMP 1598 for protein-coding sequences from 10865 expressed sequence tags (ESTs). We reconstructed a total of 4284 unigenes, including 74 ribosomal protein and 40 plastid-related genes, from ESTs using FrameDP, an open reading frame (ORF) prediction program. Correspondence analysis of A. tamarense genes based on codon usage showed that the GC content at the third base of synonymous codons (GC3s) was strongly correlated with the first axis (r = 0.93 with P < .001). On the other hand, the second axis discriminated between presumed highly and low expressed genes, with expression levels being confirmed by the analysis of EST frequencies (r = −0.89 with P < .001). Our results suggest that mutational bias is the major factor in shaping codon usage in A. tamarense genome, but other factors, namely, translational selection, hydropathy, and aromaticity, also appear to influence the selection of codon usage in this species.

Synonymous codon usage bias dependent on local nucleotide context in the class Deinococci

To study the evolution of mutation biased synonymous codon usage, we examined nucleotide co-occurrence patterns in the Deinococcus radiodurans, D. geothermalis, and Thermus thermophilus genomes for nucleotide replacement dependent on the surrounding nucleotide context. Nucleotides on the third codon site were found to be strongly correlated with nucleotide sites at most six nucleotides away in all three species, where abundance patterns were dependent on whether two nucleotides share the same purine(R)/pyrimidine(Y) status. In the class Deinococci adjacent third site nucleotides were strongly correlated, where NNR|NNR and NNY|NNY codon pairs were overabundant while NNR|NNY and NNY|NNR codon pairs were underabundant. By far the largest deviations in all three species occur for NN(YR)|(YR)NN codon pairs. In the Thermus species, the NNY|YNN and NNR|RNN codon pairs were overabundant versus the underabundant NNY|RNN and NNR|YNN codon pairs, whereas in the Deinococcus species the opposite over-/underabundance relationship held for adjacent (GC) bases. We also observed a weaker overabundance of NNR|NRN and NNY|NYN codon pairs versus the underabundant NNR|NYN and NNY|NRN codon pairs. The perfect purine/pyrimidine symmetry of each of these cases, plus the lack of significant deviations for nucleotide pairs on other length scales up to 20 codons apart demonstrates that a pervasive pattern of nucleotide replacement dependent on local nucleotide context, and not codon bias, has occurred in these species. This nucleotide replacement has led to modified synonymous codon usage within the class Deinococci that affects which codons are positioned at particular codon sites dependent on the local nucleotide context.

The effects of differential gene expression on coding sequence features: Analysis by one-way ANOVA

It is well-established that non-random patterns in coding DNA sequence (CDS) features can be partially explained by translational selection. Recent extensions of microarray and proteomic expression data have stimulated many genome-wide investigations of the relationships between gene expression and various CDS features. However, only modest correlations have been found. Here we introduced the one-way ANOVA, a more powerful extension of previous grouping methods, to re-examine these relationships at the whole genome scale for Saccharomyces cerevisiae, where genome-wide protein abundance has been recently quantified. Our results clarify that coding sequence features are inappropriate for use as genome-wide estimators for protein expression levels. This analysis also demonstrates that one-way ANOVA is a powerful and simple method to explore the influence of gene expression on CDS features.