Novel Genomic and Evolutionary Perspective of Cyanobacterial tRNAs

Anticodon table of the chloroplast genome and identification of putative quadruplet anticodons in chloroplast tRNAs

The chloroplast genome of 5959 species was analyzed to construct the anticodon table of the chloroplast genome. Analysis of the chloroplast transfer ribonucleic acid (tRNA) revealed the presence of a putative quadruplet anticodon containing tRNAs in the chloroplast genome. The tRNAs with putative quadruplet anticodons were UAUG, UGGG, AUAA, GCUA, and GUUA, where the GUUA anticodon putatively encoded tRNA^Asn. The study also revealed the complete absence of tRNA genes containing ACU, CUG, GCG, CUC, CCC, and CGG anticodons in the chloroplast genome from the species studied so far. The chloroplast genome was also found to encode tRNAs encoding N-formylmethionine (fMet), Ile2, selenocysteine, and pyrrolysine. The chloroplast genomes of mycoparasitic and heterotrophic plants have had heavy losses of tRNA genes. Furthermore, the chloroplast genome was also found to encode putative spacer tRNA, tRNA fragments (tRFs), tRNA-derived, stress-induced RNA (tiRNAs), and the group I introns. An evolutionary analysis revealed that chloroplast tRNAs had evolved via multiple common ancestors and the GC% had more influence toward encoding the tRNA number in the chloroplast genome than the genome size.

Virtual 2D map of cyanobacterial proteomes

Cyanobacteria are prokaryotic Gram-negative organisms prevalent in nearly all habitats. A detailed proteomics study of Cyanobacteria has not been conducted despite extensive study of their genome sequences. Therefore, we conducted a proteome-wide analysis of the Cyanobacteria proteome and found Calothrix desertica as the largest (680331.825 kDa) and Candidatus synechococcus spongiarum as the smallest (42726.77 kDa) proteome of the cyanobacterial kingdom. A Cyanobacterial proteome encodes 312.018 amino acids per protein, with a molecular weight of 182173.1324 kDa per proteome. The isoelectric point (pI) of the Cyanobacterial proteome ranges from 2.13 to 13.32. It was found that the Cyanobacterial proteome encodes a greater number of acidic-pI proteins, and their average pI is 6.437. The proteins with higher pI are likely to contain repetitive amino acids. A virtual 2D map of Cyanobacterial proteome showed a bimodal distribution of molecular weight and pI. Several proteins within the Cyanobacterial proteome were found to encode Selenocysteine (Sec) amino acid, while Pyrrolysine amino acids were not detected. The study can enable us to generate a high-resolution cell map to monitor proteomic dynamics. Through this computational analysis, we can gain a better understanding of the bias in codon usage by analyzing the amino acid composition of the Cyanobacterial proteome.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Plastome evolution of Aeonium and Monanthes (Crassulaceae): insights into the variation of plastomic tRNAs, and the patterns of codon usage and aversion

This study reported 13 new plastomes from Aeonium and Monanthes, and observed new markers for phylogeny and DNA barcoding, such as novel tRNA structures and codon usage bias and aversion. The Macaronesian clade of Crassulaceae consists of three genera: Aichryson, with about 15 species; Monanthes, with about 10 species; Aeonium, with about 40 species. Within this clade, Aeonium, known as "the botanical equivalent of Darwin's finches", is regarded as an excellent model plant for researching adaptive evolution. Differing from the well-resolved relationships among three genera of the Macaronesian clade, the internal branching patterns within the genus Aeonium are largely unclear. In this study, we first reported 13 new plastomes from genus Aeonium and the closely related genus Monanthes. We further performed comprehensive analyses of the plastomes, with focuses on the secondary structures of pttRNAs and the patterns of codon usage and aversion. With a typical circular and quadripartite structure, the 13 plastomes ranged from 149,900 to 151,030 bp in size, and the unique pattern in IR junctions might become a family-specific marker for Crassulaceae species. Surprisingly, the π values of plastomes from Monanthes were almost twice those from Aeonium. Most importantly, we strongly recommend that highly polymorphic regions, novel putative pttRNA structures, patterns of codon usage bias and aversion derived from plastomes might have phylogenetic implications, and could act as new markers for DNA barcoding of plants. The results of phylogenetic analyses strongly supported a clear internal branching pattern in Macaronesian clade (represented by Aeonium and Monanthes), with higher nodal support values. The findings reported here will provide new insights into the variation of pttRNAs, and the patterns of codon usage and aversion of the family Crassulaceae.

PlantMWpIDB: a database for the molecular weight and isoelectric points of the plant proteomes

The molecular weight and isoelectric point of the proteins are very important parameters that control their subcellular localization and subsequent function. Although the genome sequence data of the plant kingdom improved enormously, the proteomic details have been poorly elaborated. Therefore, we have calculated the molecular weight and isoelectric point of the plant proteins and reported them in this database. A database, PlantMWpIDB, containing protein data from 342 plant proteomes was created to provide information on plant proteomes for hypothesis formulation in basic research and for biotechnological applications. The Molecular weight and isoelectric point (pI) are important molecular parameters of proteins that are useful when conducting protein studies involving 2D gel electrophoresis, liquid chromatography-mass spectrometry, and X-ray protein crystallography. PlantMWpIDB provides an easy-to-use and efficient interface for search options and generates a summary of basic protein parameters. The database represents a virtual 2D proteome map of plants, and the molecular weight and pI of a protein can be obtained by searching on the name of a protein, a keyword, or by a list of accession numbers. The PlantMWpIDB database also allows one to query protein sequences. The database can be found in the following link https://plantmwpidb.com/ . The individual 2D virtual proteome map of the plant kingdom will enable us to understand the proteome diversity between different species. Further, the molecular weight and isoelectric point of individual proteins can enable us to understand their functional significance in different species.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

Structural variation and evolution of chloroplast tRNAs in green algae

As one of the important groups of the core Chlorophyta (Green algae), Chlorophyceae plays an important role in the evolution of plants. As a carrier of amino acids, tRNA plays an indispensable role in life activities. However, the structural variation of chloroplast tRNA and its evolutionary characteristics in Chlorophyta species have not been well studied. In this study, we analyzed the chloroplast genome tRNAs of 14 species in five categories in the green algae. We found that the number of chloroplasts tRNAs of Chlorophyceae is maintained between 28-32, and the length of the gene sequence ranges from 71 nt to 91 nt. There are 23-27 anticodon types of tRNAs, and some tRNAs have missing anticodons that are compensated for by other types of anticodons of that tRNA. In addition, three tRNAs were found to contain introns in the anti-codon loop of the tRNA, but the analysis scored poorly and it is presumed that these introns are not functional. After multiple sequence alignment, the Ψ-loop is the most conserved structural unit in the tRNA secondary structure, containing mostly U-U-C-x-A-x-U conserved sequences. The number of transitions in tRNA is higher than the number of transversions. In the replication loss analysis, it was found that green algal chloroplast tRNAs may have undergone substantial gene loss during the course of evolution. Based on the constructed phylogenetic tree, mutations were found to accompany the evolution of the Green algae chloroplast tRNA. Moreover, chloroplast tRNAs of Chlorophyceae are consistent with those of monocotyledons and gymnosperms in terms of evolutionary patterns, sharing a common multi-phylogenetic pattern and rooted in a rich common ancestor. Sequence alignment and systematic analysis of tRNA in chloroplast genome of Chlorophyceae, clarified the characteristics and rules of tRNA changes, which will promote the evolutionary relationship of tRNA and the origin and evolution of chloroplast.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNA^Ile contained the anti-codon (C-A-U), which was generally the anti-codon of tRNA^Met. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ-stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNATyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences "U-U-x2-C" and "U-x-G-x2-T," thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.

Evolutionary analysis of chloroplast tRNA of Gymnosperm revealed the novel structural variation and evolutionary aspect

Gymnosperms such as ginkgo, conifers, cycads, and gnetophytes are vital components of land ecosystems, and they have significant economic and ecologic value, as well as important roles as forest vegetation. In this study, we investigated the structural variation and evolution of chloroplast transfer RNAs (tRNAs) in gymnosperms. Chloroplasts are important organelles in photosynthetic plants. tRNAs are key participants in translation where they act as adapter molecules between the information level of nucleic acids and functional level of proteins. The basic structures of gymnosperm chloroplast tRNAs were found to have family-specific conserved sequences. The tRNAΨ -loop was observed to contain a conforming sequence, i.e., U-U-C-N-A-N2. In gymnosperms, tRNAIle was found to encode a "CAU" anticodon, which is usually encoded by tRNAMet. Phylogenetic analysis suggested that plastid tRNAs have a common polyphyletic evolutionary pattern, i.e., rooted in abundant common ancestors. Analyses of duplication and loss events in chloroplast tRNAs showed that gymnosperm tRNAs have experienced little more gene loss than gene duplication. Transition and transversion analysis showed that the tRNAs are iso-acceptor specific and they have experienced unequal evolutionary rates. These results provide new insights into the structural variation and evolution of gymnosperm chloroplast tRNAs, which may improve our comprehensive understanding of the biological characteristics of the tRNA family.

Construction of anti-codon table of the plant kingdom and evolution of tRNA selenocysteine (tRNASec)

Background

The tRNAs act as a bridge between the coding mRNA and incoming amino acids during protein translation. The anti-codon of tRNA recognizes the codon of the mRNA and deliver the amino acid into the protein translation chain. However, we did not know about the exact abundance of anti-codons in the genome and whether the frequency of abundance remains same across the plant lineage or not.

Results

Therefore, we analysed the tRNAnome of 128 plant species and reported an anti-codon table of the plant kingdom. We found that CAU anti-codon of tRNA^Met has highest (5.039%) whereas GCG anti-codon of tRNA^Arg has lowest (0.004%) abundance. However, when we compared the anti-codon frequencies according to the tRNA isotypes, we found tRNA^Leu (7.808%) has highest abundance followed by tRNA^Ser (7.668%) and tRNA^Gly (7.523%). Similarly, suppressor tRNA (0.036%) has lowest abundance followed by tRNA^Sec (0.066%) and tRNA^His (2.109). The genome of Ipomoea nil, Papaver somniferum, and Zea mays encoded the highest number of anti-codons (isoacceptor) at 59 each whereas the genome of Ostreococcus tauri was found to encode only 18 isoacceptors. The tRNA^Sec genes undergone losses more frequently than duplication and we found that tRNA^Sec showed anti-codon switch during the course of evolution.

Conclusion

The anti-codon table of the plant tRNA will enable us to understand the synonymous codon usage of the plant kingdom and can be very helpful to understand which codon is preferred over other during the translation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07216-3.

Genome-wide analysis revealed novel molecular features and evolution of Anti-codons in cyanobacterial tRNAs

Transfer RNAs (tRNAs) play important roles to decode the genetic information contained in mRNA in the process of translation. The tRNA molecules possess conserved nucleotides at specific position to regulate the unique function. However, several nucleotides at different position of the tRNA undergo modification to maintain proper stability and function. The major modifications include the presence of pseudouridine (Ψ) residue instead of uridine and the presence of m⁵-methylation sites. We found that, Ψ13 is conserved in D-stem, whereas Ψ38 & Ψ39 were conserved in the anti-codon loop (AL) and anti-codon arm (ACA), respectively. Furthermore, Ψ55 found to be conserved in the Ψ loop. Although, fourteen possible methylation sites can be found in the tRNA, cyanobacterial tRNAs were found to possess conserved G9, m³C₃₂, C₃₆, A₃₇, m⁵C₃₈ and U₅₄ methylation sites. The presence of multiple conserved methylation sites might be responsible for providing necessary stability to the tRNA. The evolutionary study revealed, tRNA^Met and tRNA^Ile were evolved earlier than other tRNA isotypes and their evolution is date back to at least 4000 million years ago. The presence of novel pseudouridination and m⁵-methylation sites in the cyanobacterial tRNAs are of particular interest for basic biology. Further experimental study can delineate their functional significance in protein translation.

Analysis of genomic tRNA revealed presence of novel genomic features in cyanobacterial tRNA

Transfer RNAs (tRNA) are important molecules that involved in protein translation machinery and acts as a bridge between the ribosome and codon of the mRNA. The study of tRNA is evolving considerably in the fields of bacteria, plants, and animals. However, detailed genomic study of the cyanobacterial tRNA is lacking. Therefore, we conducted a study of cyanobacterial tRNA from 61 species. Analysis revealed that; cyanobacteria contain thirty-six to seventy-eight tRNA gens per genome that encodes for 20 tRNA isotypes. The number of iso-acceptors (anti-codons) ranged from thirty-two to forty-three per genome. tRNA^Ile with anti-codon AAU, GAU, and UAU was reported to be absent from the genome of Gleocapsa PCC 73,106 and Xenococcus sp. PCC 7305. Instead, they were contained anti-codon CAU that is common to tRNA^Met and tRNA^Ile as well. The iso-acceptors ACA (tRNA^Cys), ACC (tRNA^Gly), AGA, ACU (tRNA^Ser), AAA (tRNA^Phe), AGG (tRNA^Pro), AAC (tRNA^Val), GCG (tRNA^Arg), AUG (tRNA^His), and AUC (tRNA^Asp) were absent from the genome of cyanobacterial lineages studied so far. A few of the cyanobacterial species encode suppressor tRNAs, whereas none of the species were found to encode a selenocysteine iso-acceptor. Cyanobacterial species encode a few putative novel tRNAs whose functions are yet to be elucidated.

Comparative Analysis of the Complete Mitochondrial Genomes for Development Application

This present research work reports the comparative analysis of the entire nucleotide sequence of mitochondrial genomes of Serranochromis robustus and Buccochromis nototaenia and phylogenetic analyses of their protein-coding genes in order to establish their phylogenetic relationship within Cichlids. The mitochondrial genomes of S. robustus and B. nototaenia are 16,583 and 16,580 base pairs long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and one control region (D-loop) which is 888 and 887 base pairs long, respectively, showing the same gene order and identical number of gene or regions with other well-elucidated mitogenomes of Cichlids. However, with exception of cytochrome-c oxidase subunit-1 (COX-1) gene, all the identified PCGs were initiated by ATG-codons. Structurally, 11 tRNA genes in B. nototaenia species and 9 tRNA genes in S. robustus species, folded into typical clover-leaf secondary structure created by the regions of self-complementarity within tRNA. All the 22 tRNA genes in both species lack variable loop. Moreover, 28 genes which include 12-protein-coding genes are encoded on the H-strand and the remaining 9 genes including one protein-coding gene are encoded on the L-strand. Thirteen sequences of concatenated mitochondrial protein-coding genes were aligned using MUSCLE, and the phylogenetic analyses performed using maximum likelihood and Bayesian inference showed that S. robustus and B. nototaenia had a broad phylogenetic relationship. These results may be a useful tool in resolving higher-level relationships in organisms and a useful dataset for studying the evolution of the Cichlidae mitochondrial genome, since Cichlids are well-known model species in the study of evolutionary biology, because of their extreme morphological, biogeographical, parental care behavior for eggs and larvae and phylogenetic diversities.

Genomic and evolutionary aspects of chloroplast tRNA in monocot plants

BACKGROUND

Chloroplasts are one of the most indispensable organelles that make life forms on the earth possible by their capacity to photosynthesize. These organelles possess a circular genome with a number of coding genes responsible for self-regulation. tRNAs are an important evolutionary-conserved gene family that are responsible for protein translation. However, within the chloroplast genome, tRNA machinery are poorly understood.

RESULTS

In the present study, the chloroplast genome of six monocot plants, Oryza nivara (NC_005973), Oryza sativa (NC_001320), Sachharum officinarum (NC_006084), Sorghum bicolor (NC_008602), Triticum aestivum (NC_002762), and Zea mays (NC_001666) were downloaded and analyzed to identify tRNA sequences. Further analysis of the tRNA sequences in the chloroplast genomes of the monocot plants resulted in the identification of several novel features. The length of tRNAs in the chloroplast genome of the monocot plants ranged from 59 to 155 nucleotides. Pair-wise sequence alignment revealed the presence of a conserved A-C-x-U-A-x-U-A-x-U-x5-U-A-A nucleotide consensus sequence. In addition, the tRNAs in chloroplast genomes of the monocot plants also contain 21-28 anti-codons against 61 sense codons in the genome. They also contain a group I intron and a C-A-U anti-codon for tRNAIle, which is a common anti-codon of tRNAMet. Evolutionary analysis indicates that tRNAs in the chloroplast genome have evolved from multiple common ancestors, and tRNAMet appears to be the ancestral tRNA that underwent duplication and diversification to give rise to other tRNAs.

CONCLUSION

The results obtained from the study of chloroplast tRNA will greatly help to increase our understanding of tRNA biology at a new level. Functional studies of the reported novel aspects of the chloroplast tRNA of the monocot plants will greatly help to decipher their roles in diverse cellular processes.