International Journal of Genomics

International Journal of Genomics / 2008 / Article

Research Article | Open Access

Volume 2008 |Article ID 571023 |

Zuoming Nie, Ruobing Zhou, Jian Chen, Dan Wang, Zhengbing Lv, Pingan He, Xuedong Wang, Hongdan Shen, Xiangfu Wu, Yaozhou Zhang, "Subcellular Localization and RNA Interference of an RNA Methyltransferase Gene from Silkworm, Bombyx Mori", International Journal of Genomics, vol. 2008, Article ID 571023, 7 pages, 2008.

Subcellular Localization and RNA Interference of an RNA Methyltransferase Gene from Silkworm, Bombyx Mori

Academic Editor: Jurg Bahler
Received01 Nov 2007
Accepted17 Mar 2008
Published21 May 2008


RNA methylation, which is a form of posttranscriptional modification, is catalyzed by S-adenosyl-L-methionone-dependent RNA methyltransterases (RNA MTases). We have identified a novel silkworm gene, BmRNAMTase, containing a 369-bp open reading frame that encodes a putative protein containing 122 amino acid residues and having a molecular weight of 13.88 kd. We expressed a recombinant His-tagged BmRNAMTase in E. coli BL21 (DE3), purified the fusion protein by metal-chelation affinity chromatography, and injected a New Zealand rabbit with the purified protein to generate anti-BmRNAMTase polyclonal antibodies. Immunohistochemistry revealed that BmRNAMTase is abundant in the cytoplasm of Bm5 cells. In addition, using RNA interference to reduce the intracellular activity and content of BmRNAMTase, we determined that this cytoplasmic RNA methyltransferase may be involved in preventing cell death in the silkworm.

1. Introduction

S-adenosyl-L-methionone-dependent RNA methyltrans-ferases (RNA MTases) are enzymes that use S-adenosyl-L-methionone (AdoMet) as the methyl donor group when catalyzing the methylation of RNA substrates [13]. Methylation of RNAs may modulate RNA maturation, alter RNA structure stability, and affect translation rates [4]. Several studies have indicated that almost all MTases have similar basic structures, however, the substrates targeted by various Mtases differ because of small differences in regions of the enzymes that contain the substrate binding site [5, 6]. The catalytic domain of MTases is highly conserved, including the common structural core and the AdoMet binding site, which accommodates the catalytic reaction and releases the final catalyzed product. The catalytic domain also plays an important role in stabilizing the conformation of the enzyme [7].

Due to variations in the binding sites for AdoMet, the binding methyl groups are allocated to different substrates; accordingly, AdoMet-dependent MTases are divided into five families (Class I–Class V) [8]. BmRNAMTase , which was derived from a Bombyx mori pupae stage cDNA library constructed in our laboratory, belongs to Class IV, otherwise known as the SPOUT (SpoU-TrmD) family of RNA MTases [9, 10]. The SpoU and TrmD families were previously considered to be independent families, however, amino acid sequence analyses have shown that members of these families share three common sequence motifs and also contain a “deep trefoil knot" motif in the C-terminus. Sequence analyses also indicate that mature SpoU proteins can potentially be targeted for several types of posttranslational modifications including N-myristoylation, cAMP- and cGMP-dependent protein kinase phosphorylation, and glycosylation.

In this study, we performed bioinformatics analysis of BmRNAMTase and expressed the BmRNAMTase open reading frame in E. coli BL21 (DE3). Rabbits inoculated with the recombinant protein generated high titer polyclonal antibodies, which we used to determine the subcellular localization of BmRNAMTase in a Bombyx mori cell line by confocal immunofluorescence microscopy. Finally, we also investigated the effect of BmRNAMTase dsRNA on Bm5 cell growth and determined that silencing BmRNAMTase activity resulted in increased cell death but did not induce apoptosis.

2. Materials and Methods

2.1. The Bioinformatic Analysis of BmRNAMTase

The similarity analysis of nucleotide and protein sequences was carried out using the BLASTN and BLASTP algorithms (NCBI). The deduced amino acid sequence and characteristics of BmRNAMTase were analyzed using the Expert Protein Analysis System ( The phylogenetic tree for BmRNAMTase and additional RNA methylases was generated using the MEGA 3.1 software.

2.2. Construction of Recombinant Plasmid

To amplify the open reading frame of BmRNAMTase from the cDNA fragment in a silkworm pupa cDNA library constructed by our laboratory [11], we performed the polymerase chain reaction (PCR) using upstream oligonucleotide primer P1 and downstream oligonucleotide primer P2 (Table 1). We digested the PCR product with SacI and SalI (Promega , Madison, USA) and ligated it into SacI- SalI- digested pET-28a expression vector. The recombinant plasmid was transformed into E. coli BL21 (DE3).

NameSequencesThe instead of underlineFunction

P1 SacI Amplifying the fragment of BmRNAMTase
P2 SalI
P3 T7 promoter Synthesizing the dsRNA of BmRNAMTase
P4 T7 promoter

2.3. Expression and Purification of Recombinant Protein

The transformed strain was grown at 37ºC in Luria-Bertani medium containing kanamycin (100  g/ l) until the optical density at 600 nm of the culture reached 0.6. Isopropyl-1-thio-β-D-galactopyranoside (IPTG) (Promega) was added to the culture to a final concentration of 1 mM, and cells were grown for an additional 5 hours before harvesting. The cells were pelleted by centrifugation; the cell pellet was suspended in phosphate-buffered saline (PBS), and the cells were lysed by sonication. The lysed cells were centrifuged and the supernatant was loaded onto an SP-Sepharose (G200) column; fractions collected from the Sepharose column that contained the recombinant protein were loaded onto an Ni2+-SephadexTM G-25 Superfine (HiTrapTM) desalting column (Amersham, Piscataway, USA) [12].Purified BmRNAMTase was analyzed by SDS-PAGE.

2.4. Preparation of Polyclonal Antibodies and Western Blot Analysis

Rabbits inoculated with purified recombinant protein generated high titer anti-BmRNAMTase polyclonal antibodies, as determined by ELISA [12]. We determined the specificity of the polyclonal antibodies for recombinant BmRNAMTase and for BmRNAMTase in silkworm pupae by Western blot analysis [12].

2.5. Subcellular Localization of BmRNAMTase

Bm5 cells were seeded into a dish for confocal microscopy. Cells were rinsed twice with 1mL PBS, fixed in 3.7% formaldehyde at 25ºC for 25 minutes, and rinsed three additional times with PBS. The fixed cells were blocked in a 3% BSA solution at 25ºC for 2 hours. The cells were then incubated with anti-BmRNAMTase serum (dilution, 1 : 50) at 4ºC for 12 hours. After three 10-minute washes with PBST (PBS % Tween-20), the cells were incubated with Cy3-labelled goat antirabbit antibody (dilution, 1 : 2000) (Promega) at 37ºC for 2 hours, and finally cells were incubated with 4-6-Diamidino-2-phenylindole (DAPI) (dilution, 1 : 2000) (Promega) at room temperature for 10 minutes. Following two 10-minute washes with PBST and another two 10-minute washes with PBS, cells were observed under a Nikon ECLIPSE TE2000-E confocal microscope and images were analyzed using EZ-C1 software.

2.6. RNAi

To determine BmRNAMTase function, we transfected Bm5 cells with a BmRNAMTase dsRNA. The BmRNAMTase ORF was amplified by PCR using primers containing T7 promoter sequences (P3 and P4) (Table 1). The primer used for the sense strand was P3 and for the antisense strand was P4. pET-28a -BmRNAMTase was used as the template for the PCR, and we amplified the two single strands using two primer pairs: P3 and P2, and P4 and P1 (Table 1). Using the T7 RiboMAXExpress RNAi System kit (Promega), we synthesized the dsRNA and purified PCR products according to the manufacturer’s instructions. The products were quantified spectrophotometrically at 260 nm and confirmed by agarose gel electrophoresis. Silkworm Bm5 cells were transfected with the BmRNAMTase dsRNA using liposomes (Invitrogen, Carlsbad, USA), and cytoactivity was determined after treatment with different concentrations of dsRNA (1.2, 2.4, 3.6, 4.8, 6.0, 7.2, and 8.4 g/ l) and the effectiveness of different durations of the RNAi treatment (12, 24, 36, 48, 60, 72, 84, and 96 hours) was determined using MTT [13]. Using a cell counter, we determined when the cell density reached about  units/mL after the RNAi treatment. At that time, we extracted the total protein from the cells and determined the change in BmRNAMTase content in the cells caused by the RNAi treatement using ELISAs. The effect of the treatment on cell growth and death was determined by examining the chromosomal DNA; the appearance of a characteristic DNA ladder signifies cell apoptosis [13].

3. Results

3.1. Bioinformatic Analysis of BmRNAMTase Sequence

We obtained a unique cDNA containing a 369 bp open reading frame that encodes an AdoMet-dependent RNA methyltransferase protein containing 122 amino acid and having a predicted molecular weight of 13.88 kd from a silkworm pupae cDNA library constructed in our laboratory (GenBank accession number DQ813500). This gene, BmRNAMTase, contains two exons and one intron (Figure 1(a)) and encodes an RNA methyltransferase that contains a domain that is highly conserved (SpoU) in similar genes in other species (Figure 1(b)).We have determined that BmRNAMTase is the first reported Bombyx mori gene to encode an RNA methyltransferase geneusing the BLAST software available on the NCBI web site.

The evolutionary relationships of a protein can provide important information about its functions. Therefore, we performed a phylogenetic tree analysis of RNA MTases to identify homologs of BmRNAMTase in other species (Figure 2). The phylogenetic tree shows the evolutionary relationships of the structural domains of RNA MTases from various species. It contains two primary branches: I and II. Branch I is primarily composed of vertebrate homologs, and branch II includes only invertebrate homologs. The vertebrate homologs are more highly evolutionary conserved than the invertebrate homologs. Therefore, we presumed that the functional domains in RNA MTases changed significantly during the evolutionary split between invertebrates and vertebrates. The species that are most closely related showed a high degree of homology in the structural domain within the RNA MTases.

3.2. Cloning, Expression, and Purification

We cloned the BmRNAMTase ORF into the expression vector pET-28a . We then transformed E. colito the recombinant plasmid, grew transformed cells in culture to the appropriate density, and used IPTG to induce expression of an 18.1 kd recombinant protein containing a His-tag (data not shown). We purified the recombinant protein using successive steps of SP-Sepharose (G200) column chromatography followed by metal-chelating affinity chromatography on a column containing Ni2+-SephadexTM G-25 Superfine. We obtained a yield of approximately 4.0 mg of purified protein from a one liter E. coli culture. The purity of the recombinant protein, as estimated by SDS-PAGE, was greater than 90% (Figure 3).

3.3. Immunohistochemistry Analysis

ELISAs showed that the titer of the polyclonal antibody is about , and Western blot analyses revealed that the antibody is highly specific for BmRNAMTase. We also determined that BmRNAMTase is highly expressed in silkworm pupae. However, Western blot analyses indicated an apparent molecular weight of about 20 kd, which is 6 kd greater than the predicted molecular weight (Figure 4). To address this inconsistency, we used to identify potential sites in the BmRNAMTase amino acid sequence that might be posttranslationally modified. We found that BmRNAMTase contains two potential sites for Casein kinase II phosphorylation: 5–8 (TRPE) and 21–24 (TVED); three potential sites for N-myristoylation: 40–45 (GVDYST), 73–78 (GLRLNI), and 90–95 (GMATAV); and two potential sites for Protein kinase C phosphorylation: 64–66 (SYK) and 69–71 (SSR). The amino acid sequence of BmRNAMTase also contains several potential glycosylation sites, which is consistent with predictions made using It is possible that posttranslational modifications are responsible for the increased apparent molecular weight of BmRNAMTase.

3.4. BmRNAMTase is Abundant in the Cytoplasm of Bm5 Cell

To determine the subcellular localization of BmRNAMTase, we performed immonohistochemistry, using the anti-BmRNAMTase polyclonal antibodies as the primary antibody and Cy3-conjugated goat antirabbit antibodies as the secondary antibody, to locate BmRNAMTase in Bm5 cells. We observed immunostained cells under a Nikon ECLIPSE TE2000-E confocal microscope and determined that BmRNAMTase is found mostly in the cytoplasm; it is particularly enriched in the inboard of karyotheca, and little BmRNAMTase is found in the karyon (Figure 5). The localization of BmRNAMTase suggests that it may play a role in posttranscriptional modification of RNAs.

3.5. Effect of RNAi on Cell Growth

MTT experiments demonstrated that the optimal concentration of the BmRNAMTase dsRNA for use in RNAi experiments is 6.0  g/ l and that the optimal duration of the interference is 72 hours (data not shown). We determined cytoactivity by MTT and determined BmRNAMTase content in the cell by ELISA after cells had been treated for 72 hours (Figure 6). We also determined whether the treated Bm5 cells had become apoptotic by examining chromosomal DNA for the appearance of a characteristic DNA ladder (Figure 7).

We concluded that the activity and content of BmRNAMTase in the Bm5 cells decreased significantly following the RNAi treatment and that silencing the BmRNAMTase gene resulted in increased cell death (Figure 6). However, we did not observe the formation of a DNA ladder in the chromosomal DNA, which is characteristic indicator of apoptotic cells [13]. The appearance of the chromosomal DNA in the RNAi-treated cells indicated that cell death did increase during the treatment, however, it also demonstrated that there was no significant increase in apoptotic cell death in the treated cells. The changes in the activity and content of BmRNAMTase in the cells were possibly due to necrotic cell death.

4. Discussion

In the present study, we discovered the first BmRNAMTase gene in a silkworm pupae cDNA library. We expressed a recombinant fusion BmRNAMTase protein in the E. coli BL21 (DE3) strain, purified the recombinant protein, and used it to generate anti-BmRNAMTase polyconal antibodies for use in determining the subcellular localization of BmRNAMTase. Immunohistochemistry indicated that the majority of BmRNAMTase is found in the cytoplasm of Bm5 cells. In addition, RNAi experiments indicated that silencing BmRNAMTase results in increased cell death.

RNA MTases collectively form an enzyme superfamily [14, 15] that has important biological functions. Methylation is a common post-transcriptional modification of RNAs; this modification can mediate control of translational processing. Prerequisites for RNA methylation are the availability of AdoMet-dependent MTases, which contain a common catalytic core and binding sites for the substrate, and both AdoMet and the substrate that corresponds to the binding site in the enzyme; binding sites vary slightly, depending on the substrate. Several studies on the mechanism of AdoMet-dependent MTase recognition of the substrate have been completed [16]. Recently, studies have focused on identifying catalytic and functional domains as well as determining the specific RNA sequences recognized in targeted substrates [17, 18]. The catalytic characteristics of several RNA AdoMet-dependent MTases have been investigated extensively, however, the mechanism responsible for recognition of substrates remains obscure.

In this study, we determined that inhibiting BmRNAMTase activity by depleting its content through RNA interference results in increased cell death (Figures 6 and 7). Analysis of the amino acid sequence revealed that the functions of BmRNAMTase might be similar to those of other members of the SpoU protein family. The key amino acid residues in SpoU family proteins are: R41, S150, P143, E124, N35, and S37 [7]. By comparing the BmRNAMTase sequence to those of other RNA MTases using CLUSTALW we determined that the key amino acid residues in BmRNAMTase may be R75, S61, P79, E91, N88, and S86 and that a possible function for BmRNAMTase is catalyzing the methylation of the 2-OH at guanosine 18 in tRNAs.

Subcellular localization of BmRNAMTase in Bombyxmori by confocal immunofluorescence microscopy indicated that the protein is concentrated in the cytoplasm (Figure 5), which differs significantly from the observed subcellular distribution of FTSJ2, a novel human RNA methyltransferase [19]. FTSJ2 shares significant sequence homology with FtsJ/RrmJ, an E. coli 23S rRNA uridine-2-O-methyltransferase [20] and is primarily found in the nucleolus. BmRNAMTase shares significant sequence homology with SpoU family proteins, which primarily catalyze the transfer of a methyl group from AdoMet to the 2-OH of guanosine 18 in tRNAs [21]. However, BmRNAMTase and other classes of RNA methyltransferases show different conserved sequences (Figure 8). The existence of differing conserved sequences in various classes of methylases is probably necessary to form different functional domains and binding sites for targeting and modifying different substrates. Alternative conserved amino acid sequences might also contribute to the differences in subcellular localization observed for various methylases.

We have studied the characteristics of BmRNAMTase and established a foundation for further functional studies on this enzyme. Western blot analyses have revealed an apparent molecular weight for BmRNAMTase in a silkworm pupa that is 6 kd greater than the predicted molecular weight; determining the underlying reason for this discrepancy requires further investigation. In addition, more study is necessary to determine the catalytic domain and the mechanism of substrate recognition of BmRNAMTase.


This work was supported by financial grants from the National Basic Research Program of China (973 Program, no. 2005CB121006), the National Natural Science Foundation of China (no. 30740015), and the Zhejiang Natural Science Foundation (no. Y205449, Y304122).


  1. X. Cheng and R. J. Roberts, “AdoMet-dependent methylation, DNA methyltransferases and base flipping,” Nucleic Acids Research, vol. 29, no. 18, pp. 3784–3795, 2001. View at: Publisher Site | Google Scholar
  2. T. G. Mosbacher, A. Bechthold, and G. E. Schulz, “Structure and function of the antibiotic resistance-mediating methyltransferase AviRb from Streptomyces viridochromogenes,” Journal of Molecular Biology, vol. 345, no. 3, pp. 535–545, 2005. View at: Publisher Site | Google Scholar
  3. M. Aittaleb, T. Visone, M. O. Fenley, and H. Li, “Structural and thermodynamic evidence for a stabilizing role of Nop5p in S-adenosyl-L-methionine binding to fibrillarin,” Journal of Biological Chemistry, vol. 279, no. 40, pp. 41822–41829, 2004. View at: Publisher Site | Google Scholar
  4. K. Datta, J. L. Fuentes, and J. R. Maddock, “The yeast GTPase Mtg2p is required for mitochondrial translation and partially suppresses an rRNA methyltransferase mutant, mrm2,” Molecular Biology of the Cell, vol. 16, no. 2, pp. 954–963, 2005. View at: Publisher Site | Google Scholar
  5. P. G. Foster, C. R. Nunes, P. Greene, D. Moustakas, and R. M. Stroud, “The first structure of an RNA m5C methyltransferase, Fmu, provides insight into catalytic mechanism and specific binding of RNA substrate,” Structure, vol. 11, no. 12, pp. 1609–1620, 2003. View at: Publisher Site | Google Scholar
  6. K. Ghoshal, S. Majumder, J. Datta et al., “Role of human ribosomal RNA (rRNA) promoter methylation and of methyl-CpG-binding protein MBD2 in the suppression of rRNA gene expression,” Journal of Biological Chemistry, vol. 279, no. 8, pp. 6783–6793, 2004. View at: Publisher Site | Google Scholar
  7. K. Watanabe, O. Nureki, S. Fukai et al., “Roles of conserved amino acid sequence motifs in the SpoU (TrmH) RNA methyltransferase family,” Journal of Biological Chemistry, vol. 280, no. 11, pp. 10368–10377, 2005. View at: Publisher Site | Google Scholar
  8. H. L. Schubert, R. M. Blumenthal, and X. Cheng, “Many paths to methyltransfer: a chronicle of convergence,” Trends in Biochemical Sciences, vol. 28, no. 6, pp. 329–335, 2003. View at: Publisher Site | Google Scholar
  9. K. Watanabe, O. Nureki, S. Fukai, Y. Endo, and H. Hori, “Functional categorization of the conserved basic amino acid residues in TrmH (tRNA (Gm18) methyltansferase) enzymes,” Journal of Biological Chemistry, vol. 281, no. 45, pp. 34630–34639, 2006. View at: Publisher Site | Google Scholar
  10. G. Michel, V. Sauvé, R. Larocque, Y. Li, A. Matte, and M. Cygler, “The structure of the RlmB 23S rRNA methyltransferase reveals a new methyltransferase fold with a unique knot,” Structure, vol. 10, no. 10, pp. 1303–1315, 2002. View at: Publisher Site | Google Scholar
  11. Y. Zhang, J. Chen, Z. Nie et al., “Expression of open reading frames in silkworm pupal cDNA library,” Applied Biochemistry and Biotechnology, vol. 136, no. 3, pp. 327–343, 2007. View at: Publisher Site | Google Scholar
  12. S. Joseph and W. R. David, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, USA, 3rd edition, 2002.
  13. L. S. David, D. G. Robert, and A. L. Leslie, Cells A Laboratory Manual, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, USA, 1998.
  14. D. Chen, P. A. Frey, B. W. Lepore, D. Ringe, and F. J. Ruzicka, “Identification of structural and catalytic classes of highly conserved amino acid residues in lysine 2,3-aminomutase,” Biochemistry, vol. 45, no. 42, pp. 12647–12653, 2006. View at: Publisher Site | Google Scholar
  15. S. Roje, “S-Adenosyl-L-methionine: beyond the universal methyl group donor,” Phytochemistry, vol. 67, no. 15, pp. 1686–1698, 2006. View at: Publisher Site | Google Scholar
  16. O. Nureki, K. Watanabe, S. Fukai et al., “Deep knot structure for construction of active site and cofactor binding site of tRNA modification enzyme,” Structure, vol. 12, no. 4, pp. 593–602, 2004. View at: Publisher Site | Google Scholar
  17. G. Maravić, M. Feder, S. Pongor, M. Flögel, and J. M. Bujnicki, “Mutational analysis defines the roles of conserved amino acid residues in the predicted catalytic pocket of the rRNA:m6A methyltransferase ErmC',” Journal of Molecular Biology, vol. 332, no. 1, pp. 99–109, 2003. View at: Publisher Site | Google Scholar
  18. E. Purta, F. van Vliet, K. L. Tkaczuk et al., “The yfhQ gene of Escherichia coli encodes a tRNA:Cm32/Um32 methyltransferase,” BMC Molecular Biology, vol. 7, article 23, 2006. View at: Publisher Site | Google Scholar
  19. Y.-P. Ching, H.-J. Zhou, J.-G. Yuan, B.-Q. Qiang, H.-F. Kung, and D.-Y. Jin, “Identification and characterization of FTSJ2, a novel human nucleolar protein homologous to bacterial ribosomal RNA methyltransferase,” Genomics, vol. 79, no. 1, pp. 2–6, 2001. View at: Publisher Site | Google Scholar
  20. T. Caldas, E. Binet, P. Bouloc, A. Costa, J. Desgres, and G. Richarme, “The FtsJ/RrmJ heat shock protein of Escherichia coli is a 23 S ribosomal RNA methyltransferase,” Journal of Biological Chemistry, vol. 275, no. 22, pp. 16414–16419, 2000. View at: Publisher Site | Google Scholar
  21. H. Hori, S. Kubota, K. Watanabe et al., “Aquifex aeolicus tRNA (Gm18) methyltransferase has unique substrate specificity: tRNA recognition mechanism of the enzyme,” Journal of Biological Chemistry, vol. 278, no. 27, pp. 25081–25090, 2003. View at: Publisher Site | Google Scholar

Copyright © 2008 Zuoming Nie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.