Evolutionary and sequence-based relationships in bacterial AdoMet-dependent non-coding RNA methyltransferases

Complete list of canonical post-transcriptional modifications in the Bacillus subtilis ribosome and their link to RbgA driven large subunit assembly

Ribosomal RNA modifications in prokaryotes have been sporadically studied, but there is a lack of a comprehensive picture of modification sites across bacterial phylogeny. B. subtilis is a preeminent model organism for gram-positive bacteria, with a well-annotated and editable genome, convenient for fundamental studies and industrial use. Yet remarkably, there has been no complete characterization of its rRNA modification inventory. By expanding modern MS tools for the discovery of RNA modifications, we found a total of 25 modification sites in 16S and 23S rRNA of B. subtilis, including the chemical identity of the modified nucleosides and their precise sequence location. Furthermore, by perturbing large subunit biogenesis using depletion of an essential factor RbgA and measuring the completion of 23S modifications in the accumulated intermediate, we provide a first look at the order of modification steps during the late stages of assembly in B. subtilis. While our work expands the knowledge of bacterial rRNA modification patterns, adding B. subtilis to the list of fully annotated species after E. coli and T. thermophilus, in a broader context, it provides the experimental framework for discovery and functional profiling of rRNA modifications to ultimately elucidate their role in ribosome biogenesis and translation.

The Bacillus subtilis ywbD gene encodes RlmQ, the 23S rRNA methyltransferase forming m7G2574 in the A-site of the peptidyl transferase center

Ribosomal RNA contains many posttranscriptionally modified nucleosides, particularly in the functional parts of the ribosome. The distribution of these modifications varies from one organism to another. In Bacillus subtilis, the model organism for Gram-positive bacteria, mass spectrometry experiments revealed the presence of 7-methylguanosine (m7G) at position 2574 of the 23S rRNA, which lies in the A-site of the peptidyl transferase center of the large ribosomal subunit. Testing several m7G methyltransferase candidates allowed us to identify the RlmQ enzyme, encoded by the ywbD open reading frame, as the MTase responsible for this modification. The enzyme methylates free RNA and not ribosomal 50S or 70S particles, suggesting that modification occurs in the early steps of ribosome biogenesis.

DNA Methyltransferases: From Evolution to Clinical Applications

DNA methylation is an epigenetic mark that living beings have used in different environments. The MTases family catalyzes DNA methylation. This process is conserved from archaea to eukaryotes, from fertilization to every stage of development, and from the early stages of cancer to metastasis. The family of DNMTs has been classified into DNMT1, DNMT2, and DNMT3. Each DNMT has been duplicated or deleted, having consequences on DNMT structure and cellular function, resulting in a conserved evolutionary reaction of DNA methylation. DNMTs are conserved in the five kingdoms of life: bacteria, protists, fungi, plants, and animals. The importance of DNMTs in whether methylate or not has a historical adaptation that in mammals has been discovered in complex regulatory mechanisms to develop another padlock to genomic insurance stability. The regulatory mechanisms that control DNMTs expression are involved in a diversity of cell phenotypes and are associated with pathologies transcription deregulation. This work focused on DNA methyltransferases, their biology, functions, and new inhibitory mechanisms reported. We also discuss different approaches to inhibit DNMTs, the use of non-coding RNAs and nucleoside chemical compounds in recent studies, and their importance in biological, clinical, and industry research.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

Evolutionary History of RNA Modifications at N6-Adenosine Originating from the R-M System in Eukaryotes and Prokaryotes

Simple Summary

The m⁶A is the most abundant and well-studied modification of mRNA, and plays an important role in transcription and translation. It is known to be evolutionarily conserved machinery present in the last eukaryotic common ancestor (LECA). The writers and erasers responsible for adding or removing m⁶A belong to specific protein families, respectively, suggesting that these members are evolutionarily related. However, only some of these mRNA m⁶A modification-associated proteins have been studied from an evolutionary perspective, while there has been no comprehensive and systematic analysis of the distributions and evolutionary history of N6mA-associated proteins in the three kingdoms of life. In this study, we identified orthologues of all the reported N6mA-associated proteins in 88 organisms from three kingdoms of life and comprehensively reconstructed the evolutionary history of the RNA N6mA modification machinery. The results demonstrate that RNA N6mA-MTases are derived from at least two different types of prokaryotic DNA MTases (class α and β MTases). As the m⁶A reader, YTH proteins may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from the N-terminals of an R-M system endonuclease. In addition, the origin of eukaryotic ALKBH family proteins is inferred to be driven by at least two occasions of independent HTG from the bacterial ALKB family.

Abstract

Methylation at the N6-position of adenosine (N6mA) on mRNA (m⁶A) is one of the most widespread, highly selective and dynamically regulated RNA modifications and plays an important role in transcription and translation. In the present study, a comprehensive analysis of phylogenetic relationships, conserved domain sequence characteristics and protein structure comparisons were employed to explore the distribution of RNA N6mA modification (m⁶A, m^6,6A, m⁶Am, m^{6, 6}Am and m⁶t⁶A)-associated proteins (writers, readers and erasers) in three kingdoms of life and reveal the evolutionary history of these modifications. These findings further confirmed that the restriction-modification (R-M) system is the origin of DNA and RNA N6mA modifications. Among them, the existing mRNA m⁶A modification system derived from the last eukaryotic common ancestor (LECA) is the evolutionary product of elements from the last universal common ancestor (LUCA) or driven by horizontal gene transfer (HGT) from bacterial elements. The subsequent massive gene gains and losses contribute to the development of unique and diverse functions in distinct species. Particularly, RNA methyltransferases (MTases) as the writer responsible for adding N6mA marks on mRNA and ncRNAs may have evolved from class α and β prokaryotic “orphan” MTases originating from the R-M system. The reader, YTH proteins that specifically recognize the m⁶A deposit, may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from N-terminals of an R-M system endonuclease. The eraser, which emerged from the ALKB family (ALKBH5 and FTO) in eukaryotes, may be driven by independent HTG from bacterial ALKB proteins. The evolutionary history of RNA N6mA modifications was inferred in the present study, which will deepen our understanding of these modifications in different species.

The critical function of the plastid rRNA methyltransferase, CMAL, in ribosome biogenesis and plant development

Methylation of nucleotides in ribosomal RNAs (rRNAs) is a ubiquitous feature that occurs in all living organisms. The formation of methylated nucleotides is performed by a variety of RNA-methyltransferases. Chloroplasts of plant cells result from an endosymbiotic event and possess their own genome and ribosomes. However, enzymes responsible for rRNA methylation and the function of modified nucleotides in chloroplasts remain to be determined. Here, we identified an rRNA methyltransferase, CMAL (Chloroplast MraW-Like), in the Arabidopsis chloroplast and investigated its function. CMAL is the Arabidopsis ortholog of bacterial MraW/ RsmH proteins and accounts to the N4-methylation of C1352 in chloroplast 16S rRNA, indicating that CMAL orthologs and this methyl-modification nucleotide is conserved between bacteria and the endosymbiont-derived eukaryotic organelle. The knockout of CMAL in Arabidopsis impairs the chloroplast ribosome accumulation and accordingly reduced the efficiency of mRNA translation. Interestingly, the loss of CMAL leads not only to defects in chloroplast function, but also to abnormal leaf and root development and overall plant morphology. Further investigation showed that CMAL is involved in the plant development probably by modulating auxin derived signaling pathways. This study uncovered the important role of 16S rRNA methylation mediated by CMAL in chloroplast ribosome biogenesis and plant development.

Genome-Wide Identification and Transcriptional Expression of the METTL21C Gene Family in Chicken

The chicken is a common type of poultry that is economically important both for its medicinal and nutritional values. Previous studies have found that free-range chickens have more skeletal muscle mass. The methyltransferase-like 21C gene (METTL21C) plays an important role in muscle development; however, there have been few reports on the role of METTL21C in chickens. In this study, we performed a genome-wide identification of chicken METTL21C genes and analyzed their phylogeny, transcriptional expression profile, and real-time quantitative polymerase chain reaction (qPCR). We identified 10 GgMETTL21C genes from chickens, 11 from mice, and 32 from humans, and these genes were divided into six groups, which showed a large amount of variation among these three species. A total of 15 motifs were detected in METTL21C genes, and the intron phase of the gene structure showed that the METTL21C gene family was conservative in evolution. Further, both the transcript data and qPCR showed that a single gene’s (GgMETTL21C3) expression level increased with the muscle development of chickens, indicating that the METTL21C genes are involved in the development of chicken muscles. Our results provide some reference value for the subsequent study of the function of METTL21C.

RNA ribose methylation (2'-O-methylation): Occurrence, biosynthesis and biological functions

Methylation of riboses at 2'-OH group is one of the most common RNA modifications found in number of cellular RNAs from almost any species which belong to all three life domains. This modification was extensively studied for decades in rRNAs and tRNAs, but recent data revealed the presence of 2'-O-methyl groups also in low abundant RNAs, like mRNAs. Ribose methylation is formed in RNA by two alternative enzymatic mechanisms: either by stand-alone protein enzymes or by complex assembly of proteins associated with snoRNA guides (sno(s)RNPs). In that case one catalytic subunit acts at various RNA sites, the specificity is provided by base pairing of the sno(s)RNA guide with the target RNA. In this review we compile available information on 2'-OH ribose methylation in different RNAs, enzymatic machineries involved in their biosynthesis and dynamics, as well as on the physiological functions of these modified residues.

Bacillus subtilis exhibits MnmC-like tRNA modification activities

The MnmE-MnmG complex of Escherichia coli uses either ammonium or glycine as a substrate to incorporate the 5-aminomethyl or 5-carboxymethylaminomethyl group into the wobble uridine of certain tRNAs. Both modifications can be converted into a 5-methylaminomethyl group by the independent oxidoreductase and methyltransferase activities of MnmC, which respectively reside in the MnmC(o) and MnmC(m) domains of this bifunctional enzyme. MnmE and MnmG, but not MnmC, are evolutionarily conserved. Bacillus subtilis lacks genes encoding MnmC(o) and/or MnmC(m) homologs. The glycine pathway has been considered predominant in this typical gram-positive species because only the 5-carboxymethylaminomethyl group has been detected in tRNA^Lys_UUU and bulk tRNA to date. Here, we show that the 5-methylaminomethyl modification is prevalent in B. subtilis tRNA^Gln_UUG and tRNA^Glu_UUC. Our data indicate that B. subtilis has evolved MnmC(o)- and MnmC(m)-like activities that reside in non MnmC homologous protein(s), which suggests that both activities provide some sort of biological advantage.

Pangenome-wide and molecular evolution analyses of the Pseudomonas aeruginosa species

BACKGROUND

Drug treatments and vaccine designs against the opportunistic human pathogen Pseudomonas aeruginosa have multiple issues, all associated with the diverse genetic traits present in this pathogen, ranging from multi-drug resistant genes to the molecular machinery for the biosynthesis of biofilms. Several candidate vaccines against P. aeruginosa have been developed, which target the outer membrane proteins; however, major issues arise when attempting to establish complete protection against this pathogen due to its presumably genotypic variation at the strain level. To shed light on this concern, we proposed this study to assess the P. aeruginosa pangenome and its molecular evolution across multiple strains.

RESULTS

The P. aeruginosa pangenome was estimated to contain more than 16,000 non-redundant genes, and approximately 15 % of these constituted the core genome. Functional analyses of the accessory genome indicated a wide presence of genetic elements directly associated with pathogenicity. An in-depth molecular evolution analysis revealed the full landscape of selection forces acting on the P. aeruginosa pangenome, in which purifying selection drives evolution in the genome of this human pathogen. We also detected distinctive positive selection in a wide variety of outer membrane proteins, with the data supporting the concept of substantial genetic variation in proteins probably recognized as antigens. Approaching the evolutionary information of genes under extremely positive selection, we designed a new Multi-Locus Sequencing Typing assay for an informative, rapid, and cost-effective genotyping of P. aeruginosa clinical isolates.

CONCLUSIONS

We report the unprecedented pangenome characterization of P. aeruginosa on a large scale, which included almost 200 bacterial genomes from one single species and a molecular evolutionary analysis at the pangenome scale. Evolutionary information presented here provides a clear explanation of the issues associated with the use of protein conjugates from pili, flagella, or secretion systems as antigens for vaccine design, which exhibit high genetic variation in terms of non-synonymous substitutions in P. aeruginosa strains.