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

Genomic Diversity of a Globally Used, Live Attenuated Mycoplasma Vaccine

The Mycoplasma synoviae live attenuated vaccine strain MS-H (Vaxsafe MS; Bioproperties Pty., Ltd., Australia) is commonly used around the world to prevent chronic infections caused by M. synoviae in birds and to minimize economic losses in the poultry industry. MS-H is a temperature-sensitive strain that is generated via the chemical mutagenesis of a virulent M. synoviae isolate, 86079/7NS. 32 single nucleotide polymorphisms have been found in the genome of MS-H compared to that of 86079/7NS, including 25 in predicted coding sequences (CDSs). There is limited information on the stability of these mutations in MS-H in vitro during the propagation of the vaccine manufacturing process or in vivo after the vaccination of chickens. Here, we performed a comparative analysis of MS-H genomes after in vitro and in vivo passages under different circumstances. Studying the dynamics of the MS-H population can provide insights into the factors that potentially affect the health of vaccinated birds. The genomes of 11 in vitro laboratory passages and 138 MS-H bird reisolates contained a total of 254 sequence variations. Of these, 39 variations associated with CDSs were detected in more than one genome (range = 2 to 62, median = 2.5), suggesting that these sequences are particularly prone to mutations. From the 25 CDSs containing previously characterized variations between MS-H and 86079/7NS, 7 were identified in the MS-H reisolates and progenies examined here. In conclusion, the MS-H genome contains individual regions that are prone to mutations that enable the restoration of the genotype or the phenotype of wild-type 86079/7NS in those regions. However, accumulated mutations in these regions are rare. IMPORTANCE Preventative measures, such as vaccination, are commonly used for the control of mycoplasmal infections in poultry. A live attenuated vaccine strain (Vaxsafe MS; MS-H; Bioproperties Pty. Ltd., Australia) is used for the prevention of disease caused by M. synoviae in many countries. However, information on the stability of previously characterized mutations in the MS-H genome is limited. In this study, we performed a comparative analysis of the whole-genome sequences of MS-H seeds used for vaccine manufacturing, commercial batches of the vaccine, cultures minimally passaged under small-scale laboratory and large-scale manufacturing conditions, MS-H reisolated from specific-pathogen-free (SPF) chickens that were vaccinated under controlled conditions, and MS-H reisolated from vaccinated commercial poultry flocks around the world. This study provides a comprehensive assessment of genome stability in MS-H after in vitro and in vivo passages under different circumstances and suggests that most of the mutations in the attenuated MS-H vaccine strain are stable.

Evolutionary dynamics of codon usages for peste des petits ruminants virus

Peste des petits ruminants virus (PPRV) is an important agent of contagious, acute and febrile viral diseases in small ruminants, while its evolutionary dynamics related to codon usage are still lacking. Herein, we adopted information entropy, the relative synonymous codon usage values and similarity indexes and codon adaptation index to analyze the viral genetic features for 45 available whole genomes of PPRV. Some universal, lineage-specific, and gene-specific genetic features presented by synonymous codon usages of the six genes of PPRV that encode N, P, M, F, H and L proteins reflected evolutionary plasticity and independence. The high adaptation of PPRV to hosts at codon usages reflected high viral gene expression, but some synonymous codons that are rare in the hosts were selected in high frequencies in the viral genes. Another obvious genetic feature was that the synonymous codons containing CpG dinucleotides had weak tendencies to be selected in viral genes. The synonymous codon usage patterns of PPRV isolated during 2007–2008 and 2013–2014 in China displayed independent evolutionary pathway, although the overall codon usage patterns of these PPRV strains matched the universal codon usage patterns of lineage IV. According to the interplay between nucleotide and synonymous codon usages of the six genes of PPRV, the evolutionary dynamics including mutation pressure and natural selection determined the viral survival and fitness to its host.

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.

A Crosstalk on Codon Usage in Genes Associated with Leukemia

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

Viral adaption of staphylococcal phage: A genome-based analysis of the selective preference based on codon usage Bias

Given the high therapeutic value of the staphylococcal phage, the genome co-evolution of the phage and the host has gained great attention. Though the genome-wide AT richness in staphylococcal phages has been well-studied with nucleotide usage bias, here we proved that host factor, lifestyle and taxonomy are also important factors in understanding the phage nucleotide usages bias using information entropy formula. Such correlation is especially prominent when it comes to the synonymous codon usages of staphylococcal phages, despite the overall scattered codon usage pattern represented by principal component analysis. This strong relationship is explained by nucleotide skew which testified that the usage biases of nucleotide at different codon positions are acting on synonymous codons. Therefore, our study reveals a hidden relationship of genome evolution with host limitation and phagic phenotype, providing new insight into phage genome evolution at genetic level.

Analyses of genetics and pathogenesis of Salmonella enterica QH with narrow spectrum of antibiotic resistance isolated from yak

Salmonella is an important pathogen for public health due to food poisoning and acute infectious intestinal disease by zoonotic trait. We isolated Salmonella enterica QH which represents the normal growth condition in Luria-Bertani culture and displays a wide range of susceptibility for multiple antibiotics. To further investigate genetic and pathogenic traits of S. enterica QH, the sequencing genome of S. enterica QH and oral Salmonella infection in mice were performed in this study. Compared with other Salmonella strains, several large sequences containing prophages and genomic islands were inserted into S. enterica QH genome. Furthermore, nucleotide and synonymous codon usage patterns display mutation pressure and natural selection serving as drivers for the evolutionary trend of S. enterica QH at gene level. The unique codon usage pattern of S. enterica QH probably contributes to adaptation to environmental/host niches and to pathogenicity. In an early oral S. enterica QH infection, the levels of CD4⁺ and CD8⁺ lymphocytes significantly reduce in peripheral blood of mice, but the increasing transcription levels of some cytokines (IFN-β1, IFN-γ and CXCL10) might have pleiotypic immune effects against S. enterica QH infection. Of note, IL10 displays significant enhancement at levels of transcription and translation, suggesting that immunosuppressive effects mediated by IL10 may function as an early oral S. enterica QH infection. The systemic investigations, including genomic and genetic characterizations and biological traits of S. enterica QH in vivo and in vitro may reflect the basic lifestyle of S. enterica QH, requiring intestine colonization, undergoing environmental stresses and performing dissemination.

Significance and roles of synonymous codon usage in the evolutionary process of Proteus

Proteus spp. bacteria frequently serve as opportunistic pathogens that can infect many animals and show positive survival and existence in various natural environments. The evolutionary pattern of Proteus spp. is an unknown topic, which benefits understanding the different evolutionary dynamics for excellent bacterial adaptation to various environments. Here, the eight whole genomes of different Proteus species were analyzed for the interplay between nucleotide usage and synonymous codon usage. Although the orthologous average nucleotide identity and average nucleotide identity display the genetic diversity of these Proteus species at the genome level, the principal component analysis further shows that these species sustain the specific genetic niche at the aspect of synonymous codon usage patterns. Interestingly, although these Proteus species have A/T rich genes with underrepresented G (guanine) or C (cytosine) at the third codon positions and overrepresented A or T at these positions, some synonymous codons with A or T end are obviously suppressed in usage. The overall codon usage pattern reflected by the effective number of codons (ENC) has a significantly positive correlation with GC3 content (GC content at the third codon position), and ENC has a significantly negative correlation with the adaptation index for these species. These results suggest that the mutation pressure caused by nucleotide composition constraint serves as a dominant evolutionary dynamic driving evolutionary trend of Proteus spp., along with other selections related to natural selection, replication and fine-tune translation, and so on. Taken together, the analyses help to understand the evolutionary interplay between nucleotide and codon usage at the gene level of Proteus.

Synonymous Codon Usages as an Evolutionary Dynamic for Chlamydiaceae

The family of Chlamydiaceae contains a group of obligate intracellular bacteria that can infect a wide range of hosts. The evolutionary trend of members in this family is a hot topic, which benefits our understanding of the cross-infection of these pathogens. In this study, 14 whole genomes of 12 Chlamydia species were used to investigate the nucleotide, codon, and amino acid usage bias by synonymous codon usage value and information entropy method. The results showed that all the studied Chlamydia spp. had A/T rich genes with over-represented A or T at the third positions and G or C under-represented at these positions, suggesting that nucleotide usages influenced synonymous codon usages. The overall codon usage trend from synonymous codon usage variations divides the Chlamydia spp. into four separate clusters, while amino acid usage divides the Chlamydia spp. into two clusters with some exceptions, which reflected the genetic diversity of the Chlamydiaceae family members. The overall codon usage pattern represented by the effective number of codons (ENC) was significantly positively correlated to gene GC3 content. A negative correlation exists between ENC and the codon adaptation index for some Chlamydia species. These results suggested that mutation pressure caused by nucleotide composition constraint played an important role in shaping synonymous codon usage patterns. Furthermore, codon usage of T3ss and Pmps gene families adapted to that of the corresponding genome. Taken together, analyses help our understanding of evolutionary interactions between nucleotide, synonymous codon, and amino acid usages in genes of Chlamydiaceae family members.

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.