Nucleotide composition bias and codon usage trends of gene populations in Mycoplasma capricolum subsp. capricolum and M. Agalactiae

Codon usage bias and nucleotide bias are not influenced by the 5' flanking but by 3' and intronic region composition in SCID-associated genes

Compositional constraints, selectional and mutational forces, nucleotide skews, RNA folding free energy, gene expression, protein properties, and differential expression are a few attributes that define the molecular signatures of any gene. The absence of information regarding these attributes for genes associated with severe combined immunodeficiency disorder (SCID) prompted us to take up this study. The compositional bias influenced codon bias. Overall percent T composition is the lowest among all nucleotides. However, its distribution varies markedly across different codon positions and is not lowest at all codon positions. We, for the first time, determined the influence of intergenic elements and introns on nucleotide and codon bias on genes and found that in SCID-associated genes, the 5' flanking region neither influences the nucleotide nor codon bias, contrary to the intronic and 3' flanking region, which both influence nucleotide and codon bias in SCID associated genes. Codon usage in the SCID-associated gene set significantly differs from the codon usage present in overall human codon usage for 33 out of 59 codons (excluding start, stop, and trp encoding). Analysis of differentially expressed genes revealed that out of the 10 most differentially expressed genes, 07 genes are Zn finger proteins (ZNF728, ZNF726, ZNF676, ZNF667, ZNF439, ZNF257, and ZNF208). Applying the knowledge of codon bias, rare codons, minimum free energy, and codon adaptation index, codon deoptimization was carried out, and ZNF208 was the best suitable candidate. The study opened the area for the identification of peculiar molecular features and the development of more candidates for gene therapy purposes.

Evolution of codon usage in Taenia saginata genomes and its impact on the host

The beef tapeworm, also known as Taenia saginata, is a zoonotic tapeworm from the genus Taenia in the order Cyclophyllidea. Taenia saginata is a food-borne zoonotic parasite with a worldwide distribution. It poses serious health risks to the host and has a considerable negative socioeconomic impact. Previous studies have explained the population structure of T. saginata within the evolutionary time scale and adaptive evolution. However, it is still unknown how synonymous codons are used by T. saginata. In this study, we used 90 T. saginata strains, applying the codon usage bias (CUB). Both base content and relative synonymous codon usage (RSCU) analysis revealed that AT-ended codons were more frequently used in the genome of T. saginata. Further low CUB was observed from the effective number of codons (ENC) value. The neutrality plot analysis suggested that the dominant factor of natural selection was involved in the structuring of CUB in T. saginata. Further analysis showed that T. saginata has adapted host-specific codon usage patterns to sustain successful replication and transmission chains within hosts (Bos taurus and Homo sapiens). Generally, both natural selection and mutational pressure have an impact on the codon usage patterns of the protein-coding genes in T. saginata. This study is important because it characterized the codon usage pattern in the T. saginata genomes and provided the necessary data for a basic evolutionary study on them.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Analysis of codon usage bias of lumpy skin disease virus causing livestock infection

Lumpy skin disease virus (LSDV) causes lumpy skin disease (LSD) in livestock, which is a double-stranded DNA virus that belongs to the genus Capripoxvirus of the family Poxviridae. LSDV is an important poxvirus that has spread out far and wide to become distributed worldwide. It poses serious health risks to the host and causes considerable negative socioeconomic impact on farmers financially and on cattle by causing ruminant-related diseases. Previous studies explained the population structure of the LSDV within the evolutionary time scale and adaptive evolution. However, it is still unknown and remains enigmatic as to how synonymous codons are used by the LSDV. Here, we used 53 LSDV strains and applied the codon usage bias (CUB) analysis to them. Both the base content and the relative synonymous codon usage (RSCU) analysis revealed that the AT-ended codons were more frequently used in the genome of LSDV. Further low codon usage bias was calculated from the effective number of codons (ENC) value. The neutrality plot analysis suggested that the dominant factor of natural selection played a role in the structuring of CUB in LSDV. Additionally, the results from a comparative analysis suggested that the LSDV has adapted host-specific codon usage patterns to sustain successful replication and transmission chains within hosts (Bos taurus and Homo sapiens). Both natural selection and mutational pressure have an impact on the codon usage patterns of the protein-coding genes in LSDV. This study is important because it has characterized the codon usage pattern in the LSDV genomes and has provided the necessary data for a basic evolutionary study on them.

Analysis of codon usage bias in mitochondrial CO gene among platyhelminthes

The phenomenon of non-uniform usage of the synonymous codons, where some codons are given more preference to others, is known as codon usage bias (CUB). CUB is known to be determined by two major evolutionary forces i.e. mutation pressure and selection. We used various approaches to understand the codon usage pattern in mitochondrial CO (MT-CO) genes involved in complex IV of the respiratory chain (RC) as no work was reported yet. Our present study revealed that CUB was relatively high and the coding sequences were rich in A and T. Correspondence analysis further indicated that A/T compositional properties under mutational pressure might be affecting the codon usage pattern and was different in different classes for MT-CO gene. A highly significant correlation between A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% in all the classes indicated that compositional constraints under mutational pressure and natural selection might affect the CUB. Neutrality plot indicated that both natural selection and mutational bias affected the CUB, where, natural selection played the major role as compared to mutational pressure.

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.

Revisiting the Taxonomy of the Genus Arcobacter: Getting Order From the Chaos

Since the description of the genus Arcobacter in 1991, a total of 27 species have been described, although some species have shown 16S rRNA similarities below 95%, which is the cut-off that usually separates species that belong to different genera. The objective of the present study was to reassess the taxonomy of the genus Arcobacter using information derived from the core genome (286 genes), a Multilocus Sequence Analysis (MLSA) with 13 housekeeping genes, as well as different genomic indexes like Average Nucleotide Identity (ANI), in silico DNA–DNA hybridization (isDDH), Average Amino-acid Identity (AAI), Percentage of Conserved Proteins (POCPs), and Relative Synonymous Codon Usage (RSCU). The study included a total of 39 strains that represent all the 27 species included in the genus Arcobacter together with 13 strains that are potentially new species, and the analysis of 57 genomes. The different phylogenetic analyses showed that the Arcobacter species grouped into four clusters. In addition, A. lekithochrous and the candidatus species ‘A. aquaticus’ appeared, as did A. nitrofigilis, the type species of the genus, in separate branches. Furthermore, the genomic indices ANI and isDDH not only confirmed that all the species were well-defined, but also the coherence of the clusters. The AAI and POCP values showed intra-cluster ranges above the respective cut-off values of 60% and 50% described for species belonging to the same genus. Phenotypic analysis showed that certain test combinations could allow the differentiation of the four clusters and the three orphan species established by the phylogenetic and genomic analyses. The origin of the strains showed that each of the clusters embraced species recovered from a common or related environment. The results obtained enable the division of the current genus Arcobacter in at least seven different genera, for which the names Arcobacter, Aliiarcobacter gen. nov., Pseudoarcobacter gen. nov., Haloarcobacter gen. nov., Malacobacter gen. nov., Poseidonibacter gen. nov., and Candidate ‘Arcomarinus’ gen. nov. are proposed.

Comparative genomic analysis for nucleotide, codon, and amino acid usage patterns of mycoplasmas

The evolutionary factors in influencing the genetic characteristics of nucleotide, synonymous codon, and amino acid usage of 18 mycoplasma species were analyzed. The nucleotide usage at the 1st and 2nd codon position which determines amino acid composition of proteins has a significant correlation with the total nucleotide composition of gene population of these mycoplasma species, however, the nucleotide usage at the 3rd codon position which affects synonymous codon usage patterns has a slight correlation with either the total nucleotide composition or the nucleotide usage at the 1st and 2nd codon position. Other evolutionary factors join in the evolutionary process of mycoplasma apart from mutation pressure caused by nucleotide usage constraint based on the relationships between effective number of codons/codon adaptation index and nucleotide usage at the 3rd codon position. Although nucleotide usage of gene population in mycoplasma dominates in forming the overall codon usage trends, the relative abundance of codon with nucleotide context and amino acid usage pattern show that translation selection involved in translation accuracy and efficiency play an important role in synonymous codon usage patterns. In addition, synonymous codon usage patterns of gene population have a bigger power to represent genetic diversity among different species than amino acid usage. These results suggest that although the mycoplasmas reduce its genome size during the evolutionary process and shape the form, which is opposite to their hosts, of AT usages at high levels, this kind organism still depends on nucleotide usage at the 1st and 2nd codon positions to control syntheses of the requested proteins for surviving in their hosts and nucleotide usage at the 3rd codon position to develop genetic diversity of different mycoplasma species. This systemic analysis with 18 mycoplasma species may provide useful clues for further in vivo genetic studies on the related species.

Unique features of nucleotide and codon usage patterns in mycoplasmas revealed by information entropy

Currently, the comparison between GC usage pattern at the 3rd codon position and codon usage index is commonly used to estimate the roles of evolutionary forces in shaping synonymous codon usages, however, this kind of analysis often losses the information about the role of A/T usage bias in shaping synonymous codon usage bias. To overcome this limitation and better understand the interplay between nucleotide and codon usages for the evolution of bacteria at gene levels, in this study, we employed the information entropy method with some modification to estimate roles of nucleotide compositions in the overall codon usage bias for 18 mycoplasma species in combination with Davies-Bouldin index. At gene levels, the overall nucleotide usage bias represents A content as the highest, followed by T, G and C for mycoplasmas, resulting in a low GC content. This feature is universal across these species derived from different hosts, suggesting that the hosts have the limited impact on nucleotide usage bias of mycoplasmas. Information entropy and Davies-Bouldin index can better reveal that the nucleotide usage bias at the 3rd codon position is essential in shaping the overall nucleotide bias for all given mycoplasmas except M. pneumoniae M129. Davies-Bouldin index revealed that the 1st and 2nd codon position play more important role in synonymous codon usage bias than that of the 3rd position at gene levels. To our knowledge, this is the first comprehensive investigation into cooperation between nucleotide and codon usages for mycoplasma and extends our knowledge of the mechanisms that contribute to codon usage and evolution of this microorganism.

Analyses of nucleotide, codon and amino acids usages between peste des petits ruminants virus and rinderpest virus

Peste des petits ruminants virus (PPRV) and rinderpest virus (RPV) are two causative agents of an economically important disease for ruminants (i.e., sheep, cattle and goat). In this study, the nucleotide, codon and amino acid usages for PPRV and RPV have been analyzed by multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis represents that ACG for Thr and GCG for Ala are selected with under-representation in both PPRV and RPV, and AGA for Arg in PPRV and AGG for Arg in RPV are used with over-representation. The usage of nucleotide pair (CpG) tends to be removed from viral genes of the two viruses, suggesting that other evolutionary forces take part in evolutionary processes for viral genes in addition to mutation pressure from nucleotide usage at the third codon position. The overall nucleotide usage of viral gene is not major factor in shaping synonymous codon usage patterns, while the nucleotide usages at the third codon position and the nucleotide pairs play important roles in shaping synonymous codon usage patterns. Although PPRV and RPV are closely related antigenically, the codon and amino acid usage patterns for viral genes represent a significant genetic diversity between PPRV and RPV. Moreover, the overall codon usage trends for viral genes between PPRV and RPV are mainly influenced by mutation pressure from nucleotide usage at the third codon position and translation selection from hosts. Taken together, this is first comprehensive analyses for nucleotide, codon and amino acid usages of viral genes of PPRV and RPV and the findings are expected to increase our understanding of evolutionary forces influencing viral evolutionary pathway and adaptation toward hosts.

Genome-Wide Analysis of the Synonymous Codon Usage Patterns in Riemerella anatipestifer

Riemerella anatipestifer (RA) belongs to the Flavobacteriaceae family and can cause a septicemia disease in poultry. The synonymous codon usage patterns of bacteria reflect a series of evolutionary changes that enable bacteria to improve tolerance of the various environments. We detailed the codon usage patterns of RA isolates from the available 12 sequenced genomes by multiple codon and statistical analysis. Nucleotide compositions and relative synonymous codon usage (RSCU) analysis revealed that A or U ending codons are predominant in RA. Neutrality analysis found no significant correlation between GC12 and GC₃ (p > 0.05). Correspondence analysis and ENc-plot results showed that natural selection dominated over mutation in the codon usage bias. The tree of cluster analysis based on RSCU was concordant with dendrogram based on genomic BLAST by neighbor-joining method. By comparative analysis, about 50 highly expressed genes that were orthologs across all 12 strains were found in the top 5% of high CAI value. Based on these CAI values, we infer that RA contains a number of predicted highly expressed coding sequences, involved in transcriptional regulation and metabolism, reflecting their requirement for dealing with diverse environmental conditions. These results provide some useful information on the mechanisms that contribute to codon usage bias and evolution of RA.