The genome of Mycoplasma capricolum

Comparative Analysis of Codon Usage Bias in Six Eimeria Genomes

The codon usage bias (CUB) of genes encoded by different species' genomes varies greatly. The analysis of codon usage patterns enriches our comprehension of genetic and evolutionary characteristics across diverse species. In this study, we performed a genome-wide analysis of CUB and its influencing factors in six sequenced Eimeria species that cause coccidiosis in poultry: Eimeria acervulina, Eimeria necatrix, Eimeria brunetti, Eimeria tenella, Eimeria praecox, and Eimeria maxima. The GC content of protein-coding genes varies between 52.67% and 58.24% among the six Eimeria species. The distribution trend of GC content at different codon positions follows GC1 > GC3 > GC2. Most high-frequency codons tend to end with C/G, except in E. maxima. Additionally, there is a positive correlation between GC3 content and GC3s/C3s, but a significantly negative correlation with A3s. Analysis of the ENC-Plot, neutrality plot, and PR2-bias plot suggests that selection pressure has a stronger influence than mutational pressure on CUB in the six Eimeria genomes. Finally, we identified from 11 to 15 optimal codons, with GCA, CAG, and AGC being the most commonly used optimal codons across these species. This study offers a thorough exploration of the relationships between CUB and selection pressures within the protein-coding genes of Eimeria species. Genetic evolution in these species appears to be influenced by mutations and selection pressures. Additionally, the findings shed light on unique characteristics and evolutionary traits specific to the six Eimeria 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.

In Silico Analyses of Burial Codon Bias Among the Species of Dipterocarpaceae Through Molecular and Phylogenetic Data

Introduction:

DNA barcode, a molecular marker, is used to distinguish among the closely related species, and it can be applied across a broad range of taxa to understand ecology and evolution. MaturaseK gene (matK) and rubisco bisphosphate carboxylase/oxygenase form I gene (rbcL) of the chloroplast are highly conserved in a plant system, which are used as core barcode. This present endeavor entails the comprehensive examination of the under threat plant species based on success of discrimination on DNA barcode under selection pressure.

Result:

The family Dipterocarpaceae comprising of 15 genera is under threat due to some factors, namely, deforestation, habitat alteration, poor seed, pollen dispersal, etc. Species of this family was grouped into 6 clusters for matK and 5 clusters and 2 sub-clusters for rbcL in the phylogenetic tree by using neighbor-joining method. Cluster I to cluster VI of matK and cluster I to cluster V of rbcL genes were analyzed by various codon and substitution bias tools. Mutational pressure guided the codon bias which was favored by the avoidance of higher GC content and significant negative correlation between GC12 and GC3 (in sub-cluster I of cluster I [0.03 < P], cluster I [0.00001 < P], and cluster II [0.01 < P] of rbcL, and cluster IV [0.013 < P] of matK). After refining the results, it could be speculated that the lower null expectation values (R = 0.5 or <0.5) were less divergent from the evolutionary perspective. Apart from that, the higher null expectation values (R = >0.85) also showed the same result, which possibly could be due to the negative impact of very high and low transition rate than transversion.

Conclusion:

Through the analysis of inter-generic, inter/intra-specific variation and phylogenetic data, it was found that both selection and mutation played an important role in synonymous codon choice in these genes, but they acted inconsistently on the genes, both matK and rbcL. In vitro stable proteins of both matK and rbcL were selected through natural selection rather than mutational selection. matK gene had higher individual discrimination and barcode success compared with rbcL. These discriminatory approaches may describe the problem related to the extinction of plant species. Hence, it becomes very imperative to identify and detect the under threat plant species in advance.

Compositional bias coupled with selection and mutation pressure drives codon usage in Brassica campestris genes

The plant Brassica campestris includes the vegetables turnip and Chinese cabbage, important plants of economic importance. Here, we have analysed the codon usage bias of B. campestris for 116 protein coding genes. Neutrality analysis showed that B. campestris had a wide range of GC3s, and a significant correlation was observed between GC12 and GC3. Nc versus GC3s plot showed a few genes on or proximate to the expected curve, but the majority of points were found to be scattered distantly from the expected curve. Correspondence analysis on codon usage revealed that the position preference of codons on multidimensional space totally depends on the presence of A and T at synonymous third codon position. These results altogether suggest that composition bias along with selection (major) and mutation pressure (minor) affects the codon usage pattern of the protein coding genes in Brassica campestris.

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

Because of the low GC content of the gene population, amino acids of the two mycoplasmas tend to be encoded by synonymous codons with an A or T end. Compared with the codon usage of ovine, Mycoplasma capricolum and M. agalactiae tend to select optimal codons, which are rare codons in ovine. Due to codon usage pattern caused by genes with key biological functions, the overall codon usage trends represent a certain evolutionary direction in the life cycle of the two mycoplasmas. The overall codon usage trends of a gene population of M. capricolum subsp. capricolum can be obviously separated from other mycoplasmas, and the overall codon usage trends of M. agalactiae are highly similar to those of M. bovis. These results partly indicate the independent evolution of the two mycoplasmas without the limits of the host cell's environment. The GC and AT skews estimate nucleotide composition bias at different positions of nucleotide triplets and the protein consideration caused by the nucleotide composition bias at codon positions 1 and 2 largely take part in synonymous codon usage patterns of the two mycoplasmas. The correlation between the codon adaptation index and codon usage variation indicates that the effect of codon usage on gene expression in M. capricolum subsp. capricolum is opposite to that of M. agalactiae, further suggesting independence of the evolutionary process influencing the overall codon usage trends of gene populations of mycoplasmas.

Analysis of Codon Usage Patterns in Herbaceous Peony (Paeonia lactiflora Pall.) Based on Transcriptome Data

Codon usage bias, which exists in many genomes, is mainly determined by mutation and selection. To elucidate the genetic features and evolutionary history of herbaceous peony (Paeonia lactiflora), a well-known symbol of prosperity in China, we examined synonymous codon usage in 24,216 reconstructed genes from the P. lactiflora transcriptome. The mean GC content was 44.4%, indicating that the nucleotide content of P. lactiflora genes is slightly AT rich and GC poor. The P. lactiflora genome has a wide range of GC3 (GC content at the third synonymous codon position) distribution, with a significant correlation between GC12 and GC3. ENC (effective number of codons) analysis suggested that mutational bias played a major role in shaping codon usage. Parity Rule 2 (PR2) analysis revealed that GC and AU were not used proportionally. We identified 22 “optimal codons”, most ending with an A or U. Our results suggested that nucleotide composition mutation bias and translational selection were the main driving factors of codon usage bias in P. lactiflora. These results lay the foundation for exploring the evolutionary mechanisms and heterologous expression of functionally-important proteins in P. lactiflora.

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.

Molecular biology of mycoplasmas

Although mycoplasmas lack cell walls, they are in many respects similar to the gram-positive bacteria with which they share a common ancestor. The molecular biology of mycoplasmas is intriguing because the chromosome is uniquely small (< 600 kb in some species) and extremely A-T rich (as high as 75 mol% in some species). Perhaps to accommodate DNA with a lower G + C content, most mycoplasmas do not have the "universal" genetic code. In these species, TGA is not a stop codon; instead it encodes tryptophan at a frequency 10 times greater than TGG, the usual codon for this amino acid. Because of the presence of TGA codons, the translation of mycoplasmal proteins terminates prematurely when cloned genes are expressed in other eubacteria, such as Escherichia coli. Many mycoplasmas possess strikingly dynamic chromosomes in which high-frequency changes result from errors in DNA repair or replication and from highly active recombination systems. Often, high-frequency changes in the mycoplasmal chromosome are associated with antigenic and phase variation, which regulate the production of factors critical to disease pathogenesis.

RNase P RNA of Mycoplasma capricolum

There are at least six small stable RNAs in Mycoplasma capricolum cells besides tRNAs and rRNAs. One of them, MCS5 RNA, is a homolog of RNase P RNA. The predicted secondary structure of this RNA is essentially the same as that of other eubacterial RNase P RNAs. MCS5 RNA is more similar to the RNase P RNA of B. Subtilis than to that of E. coli. This is consistent with previous conclusions that mycoplasmas are phylogenetically related to the low G + C Gram-positive bacterial group. The major substrates for MCS5 RNA must be the precursors of tRNAs. The precursor of MCS6 RNA, which is a homolog of the E. coli 10Sa RNA, may also be a substrate for the MCS5 RNA because this RNA has a tRNA-like structure at its 5' and 3' ends.

A small RNA of Mycoplasma capricolum that resembles eukaryotic U6 small nuclear RNA

Mycoplasma capricolum, a parasitic prokaryote, contains several small stable RNAs, besides rRNAs and tRNAs. One of them, designated MCS4 RNA (125 nucleotides in length), has been isolated and sequenced. This RNA is abundant in the cell, and is encoded by two genes. Unexpectedly, MCS4 RNA has been found to reveal extensive sequence similarity to eukaryotic U6 snRNAs. This finding suggests that MCS4 and U6 snRNAs are derived from a common ancestral RNA that has existed before the divergence of prokaryotes and eukaryotes.

Hydrolysis of urea by Ureaplasma urealyticum generates a transmembrane potential with resultant ATP synthesis

When urea is added to Ureaplasma urealyticum, it is hydrolysed internally by a cytosolic urease. Under our measuring conditions, and at an external pH of 6.0, urea hydrolysis caused an ammonia chemical potential equivalent to almost 80 mV and, simultaneously, an increase in proton electrochemical potential (delta p) of about 24 mV with resultant de novo ATP synthesis. Inhibition of the urease with the potent inhibitor flurofamide abolished both the chemical potential and the increase of delta p such that ATP synthesis was reduced to approximately 5% of normally obtained levels. Uncouplers of electrochemical gradients had little or no effect on these systems. The electrochemical parameters and ATP synthesis were measured similarly at three other external pH values. Any change in delta p was primarily via membrane potential (delta psi), and the level of de novo ATP synthesis was related to the increase in delta p generated upon addition of urea and more closely to the ammonia chemical potential. Although the organisms lack an effective mechanism for internal pH homeostasis, they maintained a constant delta pH. The data reported are consistent with, and give evidence for, the direct involvement of a chemiosmotic mechanism in the generation of around 95% of the ATP by this organism. Furthermore, the data suggest that the ATP-generating system is coupled to urea hydrolysis by the cytosolic urease via an ammonia chemical potential.

'boxA'-like sequence between the 16 S/23 S spacer in rRNA operon of mycoplasmas

We have found that a boxA-like sequence is conserved in the 16 S and 23 S rRNA intergenic spacer regions of mycoplasmas, and that it always locates on loop regions of the hypothetical secondary stem-loop structures. A nucleotide sequence similar to the '-10' box of prokaryotic promoters was identified at upstream sites of the boxA-like sequence in the 16 S/23 S spacer regions. These structures may represent an internal promoter between the 16 S and 23 S rRNA genes in mycoplasmas.

Levels of tRNAs in bacterial cells as affected by amino acid usage in proteins

Transfer RNAs of Mycoplasma capricolum were separated by two-dimensional polyacrylamide gel electrophoresis, and the relative abundance of each of the 28 known tRNA species was measured. There existed a correlation between the relative amount of isoacceptor tRNAs and the frequency in choosing synonymous codons that could be translated by the isoacceptors. Furthermore, it was observed that the total amount of tRNAs for a particular amino acid was paralleled by the composition of the amino acid in ribosomal proteins. A similar relationship was obtained from reexamination of the previous data on Escherichia coli tRNAs, suggesting that the amount of tRNAs for an amino acid is affected by the usage of the amino acid in proteins.

Isolation and detection of urease genes in Ureaplasma urealyticum

Urease from ureaplasmas was purified by immunoaffinity chromatography, and the N-terminal amino acid sequence was determined for two of the three subunits. These sequences were used to design primers for a polymerase chain reaction (PCR) that amplified most of the gene coding for one of the subunits. By using a novel "PCR walking" technique, we synthesized almost the complete locus on two overlapping PCR products. We present here a partial nucleotide sequence of the urease locus from Ureaplasma urealyticum (serotype 8), which agrees with our N-terminal amino acid data but differs slightly from the sequence previously reported (A. Blanchard, Mol. Microbiol. 4:669-676, 1990). Also described are PCR primers, intended for diagnostic use, that amplify a sequence from all Ureaplasma strains tested but not from any other mycoplasmas or urease-positive bacteria.

CGG: an unassigned or nonsense codon in Mycoplasma capricolum

CGG is an arginine codon in the universal genetic code. We previously reported that in Mycoplasma capricolum, a relative of Gram-positive eubacteria, codon CGG did not appear in coding frames, including termination sites, and tRNA(ArgCCG) pairing with codon CGG, was not detected. These facts suggest that CGG is a nonsense (unassigned and untranslatable) codon--i.e., not assigned to arginine or to any other amino acid. We have investigated whether CGG is really an unassigned codon by using a cell-free translation system prepared from M. capricolum. Translation of synthetic mRNA containing in-frame CGG codons does not result in "read-through" to codons beyond the CGG codons--i.e., translation ceases just before CGG. Sucrose-gradient centrifugation profiles of the reaction mixture have shown that the bulk of peptide that has been synthesized is attached to 70S ribosomes and is released upon further incubation with puromycin. The result suggests that the peptide is in the P site of ribosome in the form of peptidyl-tRNA, leaving the A site empty. When in-frame CGG codons are replaced by UAA codons in mRNA, no read-through occurs beyond UAA, just as in the case of CGG. However, the synthesized peptide is released from 70S ribosomes, presumably by release factor 1. These data suggest strongly that CGG is an unassigned codon and differs from UAA in that CGG is not used for termination.

The organization and evolution of transfer RNA genes in Mycoplasma capricolum

The genes for presumably all the tRNA species in Mycoplasma capricolum, a derivative of Gram-positive eubacteria, have been cloned and sequenced. There are 30 genes encoding 29 tRNA species. This number is the smallest in all the known genetic systems except for mitochondria. The sequences of 9 tRNA genes of them have been previously reported (1-3). Twenty-two genes are organized in 5 clusters consisting of nine, five, four and two genes (2 sets), respectively. The other eight genes exist as a single transcription unit. All the tRNAs are encoded each by a single gene, except for the occurrence of two tRNA(Lys)(TTT) genes. The arrangement of tRNA genes in the 9-gene cluster, the 5-gene cluster, the 4-gene cluster and one of the 2-gene clusters reveals extensive similarity with a part of the 21-tRNA gene cluster and/or the 16-tRNA gene cluster in Bacillus subtilis, respectively. The results suggest that the present M. capricolum tRNA genes have evolved from large tRNA gene clusters in the ancestral Gram-positive bacterial genome common to M. capricolum and B. subtilis, by discarding genes for redundant as well as non-obligate tRNAs, so that all the codons may be translated by as small a number of tRNAs as possible.

Location of sites of transposon Tn916 insertion in the Mycoplasma mycoides genome

The genome of Mycoplasma mycoides subsp. mycoides GC1176-2 was cleaved into six large fragments by the endonuclease KpnI which also cleaved the transposon Tn916 once. This has allowed genomic mapping of insertion sites for 50 transformants of GC1176-2 containing Tn916. Almost all of the mapped sites were clearly separate. The transformants provide a bank of genomes each with a KpnI site at a different position to facilitate mapping of gene loci.

Spectinomycin operon of Micrococcus luteus: evolutionary implications of organization and novel codon usage

The complete DNA sequence of the Micrococcus luteus spectinomycin (spc) operon and its adjacent regions has been determined. The sequence has revealed the presence of genes that are homologous to those of the Escherichia coli ribosomal and related proteins, L14, L24, L5, S8, L6, L18, S5, L30, L15, and secretion protein Y (sec Y), and the gene for adenylate kinase (adk). The gene arrangement in the spc operon is essentially the same as that of E. coli except for the absence in the M. luteus spc operon of the genes for S14 and X protein that exist in the E. coli spc operon. SecY and adk seem to be composed of another operon (adk operon) with at least an open reading frame. The deduced amino acid sequences for these ribosomal proteins are well conserved among the two species (40-65% identity). Reflecting the high genomic guanine and cytosine (GC) content of M. luteus (74%), the codon usage of the genes is extremely biased toward use of G and C, about 94% of the codon third positions being G or C. Seven codons, AUA, AAA, AGA, UUA, GUA, CUA, and CAA, all of which have A at the codon third positions, are completely absent in the M. luteus genes examined. Out of 11 genes in the M. luteus spc and adk operons, 5 (10) use GUG (UGA) and 6 (1) use AUG (UAA) as an initiation (termination) codon.

Codon recognition patterns as deduced from sequences of the complete set of transfer RNA species in Mycoplasma capricolum. Resemblance to mitochondria

The nucleotide sequences of the complete set of tRNA species in Mycoplasma capricolum, a derivative of Gram-positive eubacteria, have been determined. This bacterium represents the first genetic system in which the sequences of all the tRNA species have been determined at the RNA level. There are 29 tRNA species: three for Leu, two each for Arg, Ile, Lys, Met, Ser, Thr and Trp, and one each for the other 12 amino acids as judged from aminoacylation and the anticodon nucleotide sequences. The number of tRNA species is the smallest among all known genetic systems except for mitochondria. The tRNA anticodon sequences have revealed several features characteristic of M. capricolum. (1) There is only one tRNA species each for Ala, Gly, Leu, Pro, Ser and Val family boxes (4-codon boxes), and these tRNAs all have an unmodified U residue at the first position of the anticodon. (2) There are two tRNAThr species having anticodons UGU and AGU; the first positions of these anticodons are unmodified. (3) There is only one tRNA with anticodon ICG in the Arg family box (CGN); this tRNA can translate codons CGU, CGC and CGA. No tRNA capable of translating codon CGG has been detected, suggesting that CGG is an unassigned codon in this bacterium. (4) A tRNATrp with anticodon UCA is present, and reads codon UGA as Trp. On the basis of these and other observations, novel codon recognition patterns in M. capricolum are proposed. A comparatively small total, 13, of modified nucleosides is contained in all M. capricolum tRNAs. The 5' end nucleoside of the T psi C-loop (position 54) of all tRNAs is uridine, not modified to ribothymidine. The anticodon composition, and hence codon recognition patterns, of M. capricolum tRNAs resemble those of mitochondrial tRNAs.

The ribosomal protein gene cluster of Mycoplasma capricolum

The DNA sequence of the part of the Mycoplasma capricolum genome that contains the genes for 20 ribosomal proteins and two other proteins has been determined. The organization of the gene cluster is essentially the same as that in the S10 and spc operons of Escherichia coli. The deduced amino acid sequence of each protein is also well conserved in the two bacteria. The G + C content of the M. capricolum genes is 29%, which is much lower than that of E. coli (51%). The codon usage pattern of M. capricolum is different from that of E. coli and extremely biased to use of A and U(T): about 91% of codons have A or U in the third position. UGA, which is a stop codon in the "universal" code, is used more abundantly than UGG to dictate tryptophan.