Skip to main content


N6-methyl-adenosine level in Nicotiana tabacum is associated with tobacco mosaic virus



N6-methyl-adenosine (m6A) is a prevalent RNA modification in many species. Abnormal m6A methylation levels can lead to RNA dysfunction and can cause diseases. Tobacco mosaic virus (TMV) is one of the most devastating viruses for agricultural plants. It has many hosts, particularly including tobacco and other members the family Solanaceae. However, it remains unclear whether the abnormal growth induced by TMV is associated with the m6A level.


A rapid and accurate analytical method using ultra-high-performance liquid chromatography coupled with high-resolution tandem mass spectrometry (UHPLC−HR − MS/MS) was developed to analyse the adenosine (A), cytidine (C), guanosine (G), uridine (U), and m6A contents in the tobacco leaf, and the m6A/G ratio was used to evaluate the m6A level. Subsequent protein sequence alignments were used to find the potential methylases and demethylases in Nicotiana tabacum (N. tabacum). Finally, reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) was used to analyse the gene expression levels of the potential methylases and demethylases in the N. tabacum leaf.


The results showed that TMV reduced the m6A level. Moreover, protein sequence alignments revealed partial homology among human ALKBH5, Arabidopsis (NP_001031793), and Nicotiana sylvestris (XP_009800010). The gene expression level of the potential demethylase XM_009801708 increased at 14 and 21 days in N. tabacum infected with TMV, whereas all of the potential methylases decreased.


The reversible m6A modification in N. tabacum mRNA might represent a novel epigenetic mechanism involved in TMV.


M6A is a prevalent RNA modification in species such as viruses, bacteria [1], yeasts [2], plants [3], and mammals [4, 5]. The modification of m6A is mediated by N6-adenosine methyltransferase complexes such as 70 kD SAM (S-adenosylmethionine)-binding subunit methyltransferase-like 3 (METTL3, also called MT-A70), methyltransferase-like 14 (METTL14), and Wilms tumour 1-associated protein (WTAP) [6, 7]. METTL3 and METTL14 serve as the catalytic subunits, and WTAP serves as the regulatory subunit. The discovery of m6A demethylases of FTO (fat mass and obesity-associated protein) [8] and ALKBH5 [9] reveals that RNA methylation is reversible and can dynamically the mRNA metabolism. m6A methylation plays a broad and important role in the functional interplay among m6A methylases, demethylases, and m6A binding proteins [3]. Abnormal m6A methylation levels, which are induced by defects in any factor in this network, can lead to RNA dysfunction and can cause diseases [10,11,12,13,14]. Zhong et al. [11] report that the inactivation of mRNA adenosine methylase (MTA) prevents the progression of the developing embryo from entering the globular stage, leading to an embryo-lethal phenotype with seed arrestment in plants. Bodi et al. [15] report that the reduced expression of MTA in Arabidopsis thaliana (A. thaliana) reduces the m6A level in mRNA and induces abnormal growth, including a reduction of apical dominance, abnormal organ definition, and an increase in trichome branching. Shen et al. [16] show that FIP37 is a core component of the m6A methyltransferase complex and determines Arabidopsis m6A mRNA modification pattern. FIP37 mediates m6A modification of key shoot meristem regulator mRNAs and, through the control of their stability, prevents overproliferation of the shoot apical meristem that produces all aerial organs. However, it remains unclear whether the abnormal growth induced by plant viruses is associated with the m6A levels.

Nicotiana tabacum (N. tabacum) is the most commonly grown of all plants in the Nicotiana genus. The three most commonly used tobacco types are Flue-Cured (or Virginia), Burley and Oriental. N. tabacum K326 belongs to type of Flue-Cured. Tobacco is used as a model for plant disease susceptibility, which it shares with other Solanaceae plants including potato, tomato and pepper. Recently research confirmed that Nicotiana sylvestris and Nicotiana tomentosiformis are the most likely progenitors of N. tabacum. [17]. Tobacco mosaic virus (TMV) is a positive-sense single stranded RNA virus, genus tobamovirus. It has a wide host range, including vegetables, ornamentals, legumes, and other important crops [18]. Although TMV has been studied for more than 110 years [19], the association between m6A and TMV has not been reported.

The traditional methods for studying RNA modifications are 32P-labelling, 2D thin-layer chromatography [11, 20, 21], dot-blot [8], high-performance liquid chromatography (HPLC) [10], and capillary electrophoresis coupled with laser-induced fluorescence detection [22]. However, these methods are tedious, semi-quantitative, and low-throughput, and they require a large amount of RNA. HPLC coupled with a triple-quadruple mass spectrometry using external standards was recently used to measure the m6A levels in mRNA [9, 23], other RNA modifications in tRNA and small RNA, and a femtomole level of sensitivity was achieved [24]. Furthermore, a method of stable isotope-labelled internal standards was used to measure these post-transcriptionally modified nucleosides in the cellular RNA species [25]. Although this method is more accurate than a triple-quadruple mass spectrometer using external standards, the production of stable isotope-labelled compounds is complicated and time-consuming. Therefore, a rapid and accurate analytical method is required to better understand the relationship between TMV and m6A methylation.

A rapid method to determine the levels of adenosine (A), cytidine (C), guanosine (G), uridine (U), and m6A in mRNA in Nicotiana tabacum (N. tabacum) using ultra-high-performance liquid chromatography coupled with high-resolution tandem mass spectrometry (UHPLC−HR − MS/MS) was developed in this study. This assay was used to determine the effect of TMV on m6A mRNA methylation in N. tabacum. Subsequent protein sequence alignments were applied to find the potential methylases and demethylases in N. tabacum. Finally, RT − qPCR was used to investigate the gene expression of potential methylases and demethylases in N. tabacum.


Plant materials and growth conditions

N. tabacum K326 seeds were purchased from the Guizhou Academy of Tobacco Science. TMV was obtained from the Wuhan Institute of Virology at the Chinese Academy of Sciences, and its source was preserved in a N. tabacum K326 plant over an extended period. The extraction and purification of TMV were conducted as previously described by Gooding and Iwata [26]. N. tabacum K326 seeds were grown in a disease-free greenhouse, and 12-week-old plants were used for the study. The temperature and relative humidity levels in the greenhouse were maintained close to 25 °C and 60%, respectively.

The leaf samples were divided into two groups: a control group of healthy N. tabacum leaves and a treatment group of N. tabacum leaves infected with TMV. N. tabacum were mechanically inoculated by dusting carborundum on the leaves, wetting a brush with the inoculum and rubbing the leaves of selected plants. The inoculated leaves were rinsed with sterile distilled water 30 min after inoculation. The controls were prepared in a similar fashion using sterile PBS [27]. PCR was used to detect whether the N. tabacum infected TMV (see Additional file 1). The leaf samples were collected at 3, 5, 7, 10, 14, and 21 days after inoculation with TMV, removed from the main leaf vein, frozen in liquid nitrogen, and stored at − 80 °C before use. The collected leaf samples in the treatment group were not inoculated by TMV directly. Each experiment was repeated three times.

Sample preparation

To prepare mRNA from N. tabacum leaves, 0.50 g of N. tabacum leaves was prepared, and the total RNA was extracted using TRIzol reagent (Takara Biotechnology Co. Ltd., Otsu, Japan) according to the manufacturer’s protocol. RNA pellets were dissolved in 200 μL of 1% diethylpyrocarbonate H2O (DEPC-water), and the quality was tested using an ASP-3700 ultraviolet/visible (UV-vis) spectrophotometer (Xiamen biotech co., LTD, Xiamen, China) and gel electrophoresis. The total mRNA was isolated using poly(dT) oligo magnetic beads, followed by rRNA depletion. The mRNA concentrations and qualities were measured using UV spectrophotometry.

To digest mRNAs to nucleosides, 0.5 unit of nuclease P1 (Sigma-Aldrich, St. Louis, MO) and 2.5 μL of 1 M NH4Ac (pH = 5.3) were added to 50 ng of mRNA. DEPC-water was added to obtain a 20 μL volume. The reaction mixture was incubated at 42 °C for 8 h. Before adding 1 unit of alkaline phosphatase (Sigma-Aldrich, St. Louis, MO, USA), 2.5 μL of NH4HCO3 was added. After digestion at 37 °C for 8 h, the resulting digestion mixture was filtered through a 0.22-μm nylon syringe filter for UHPLC−HR − MS/MS analysis.

Condition of UHPLC−HR − MS/MS

UHPLC−HR − MS/MS analyses of A, C, G, U, and m6A were conducted on an UltiMate 3000 UHPLC coupled with a Q Exactive high resolution mass spectrometer (Thermo Fisher Scientific, Foster City, CA, USA).

Chromatographic separation of A, C, G, U, and m6A was performed on an UltiMate 3000 UHPLC system, a binary solvent manager, a column oven, a solvent degasser, and an autosampler equipped with an Agilent ZORBOX SB-Aq column (4.6 mm × 50 mm, 1.8 μm particle size; Santa Clara, CA, USA). Solutions of 0.1% (v/v) formic acid in water (solution A) and 0.1% (v/v) formic acid in methanol (solution B) were used as the mobile phases. The nucleosides were separated using a gradient of 20% B for 1 min, 20–40% B for 2 min, 40% B for 3 min, 40–20% B for 0.01 min, and 20% B for 2 min. The flow rate was 0.3 mL/min. The column oven temperature was maintained at 40 °C to reduce the viscosity, and the temperature of the sample vial holder was set at 4 °C. The injected sample volume was maintained at 2 μL. The target compounds were all eluted within 8 min.

A Q Exactive high-resolution mass spectrometer equipped with an electrospray ionization source was used to quantify the tested analytes. Nitrogen was used as the sheath gas and aux gas. The electrospray ionization source-dependent parameter settings were as follows: sheath gas flow rate, 35 psi (nitrogen); aux gas flow rate, 8 psi (nitrogen); spray voltage, 3.5 kV; capillary temperature, 320 °C; aux gas heater temperature, 300 °C; and S-lens RF level, 50. The target selected ion monitoring (SIM) in positive mode was used for the quantitative determination of A, C, G, U, and m6A. The resolution of the mass spectrometer was 70,000. m/z values of 245.07681 (U), 282.11968 (m6A), 284.09895 (G), 244.09289 (C), and 268.10403 (A) were used for the qualitative detection and quantification.

Method validation

To ensure that the method was accurate and reliable, the method was validated to evaluate its performance according to the Food and Drug Administration Guidance for Industry Bioanalytical Method Validation [28]. The validation parameters included the linear range, limit of detection (LOD), limit of quantification (LOQ), precision, and accuracy.

The LOD and LOQ were 3 and 10, respectively, as determined using S/N. The intra- and inter-day accuracy and precision were assessed by analysing the standard solutions of ribonucleosides at three different concentrations. The inter-day results were obtained from analyses of five replicates conducted on 3 separate days. Accuracy and precision were expressed as the percent recovery and relative standard deviation (RSD), respectively. The stabilities of analytes present in mRNA were evaluated after three cycles of freeze (at − 20 °C for 24 h) and thaw (at room temperature). The percent recovery was evaluated by comparing the levels before and after the freeze-and-thaw cycles, and the RSD was calculated based on the present recovery obtained from three aliquots of RNA samples.

Protein sequence alignment

To discover the potential methylases and demethylases in N. tabacum and investigate conservation of a family of ALKB in diverse species, amino acid sequences of human ALKBH5, closely related proteins from A. thaliana, and Nicotiana sylvestris (N. sylvestris) were identified using database searches and aligned using AlignX (a component of Vector NTI Suite 11.0 software, InforMax Inc., Rockville, MD, USA), and the protein sequences were obtained from NCBI information (

Method of RT–qPCR

To investigate the reason for the m6A level decrease, RT − qPCR was used to study the gene expression level of the potential methylases and demethylases in N. tabacum.

The gene sequences analysed by RT-qPCR in N. tabacum were obtained from the BLAST searches in NCBI ( and compared to the methylase and potential demethylase of A. thaliana. The detailed RT-qPCR operation was performed as described below. First, 0.15 g of N. tabacum leaves was prepared, and the total RNA was extracted using TRIzol reagent according to the manufacturer’s protocol. RNA pellets were dissolved in 30 μL of RNA-free water with concentration and purity determined by an ASP-3700 UV spectrophotometer. Reverse transcription (RT) was performed to synthesize cDNA using the PrimeScript™ RT reagent Kit (Takara Biotechnology Co. Ltd., Otsu, Shiga, Japan) in C1000™ Thermal Cycler (Bio-Rad, Hercules, CA, USA). The qPCR primers were designed using the online NCBI sequence information. The target genes were XM_009766347, XM_009766348, XM_009775897, XM_009770927, XM_009801708, and XM_009801663, and the reference gene was β-actin. The primers of target genes were described in [see Additional file 1: Table S1]. Plasmid transformation and cloning was conducted using a pUCm-T DNA Cloning Vector Kit (Sangon Biotech, Shanghai, China). The sequences of target genes were validated by sequencing using a pUCm-T DNA Cloning Vector Kit, and the sequences agreed with the NCBI sequence information. qPCR analysis was conducted using SYBR Green technology with the SYBR® Premix Ex Taq™ II (TliRNaseH Plus) (Takara Biotechnology Co. Ltd., Otsu, Shiga, Japan) in an iCycleriQ™ instrument according to the manufacturer’s instructions (Bio-Rad, Hercules, CA, USA). The optimized RT − qPCR parameters included the annealing temperature, number of cycles, and cDNA concentration. The thermal cycling profile for qPCR amplifications was as follows: 95 °C for 30 s, followed by 35 cycles at 95 °C for 30 s, 58 °C for 30 s, and 72 °C for 30 s. The relative quantification of the expression levels of selected genes was performed using the 2-ΔΔCt method.

Statistical analysis

Data were expressed as the means ± standard deviation (SD) of at least triplicate assays. Statistical differences of the data were evaluated by using IBM SPSS Statistics 19. The differences between the measurement data were analysed using the t-test, and differences were considered significant at P ≤ 0.05 or P ≤ 0.01.


Optimization of MS/MS

The performances of three types of columns, including the Agilent ZORBOX SB-Aq column (4.6 mm × 50 mm, 1.8 μm particle size, Santa Clara, CA, USA), Thermo Scientific Hypersil GOLD column (2.1 mm × 150 mm, 1.9 μm particle size, Foster City, CA, USA), and Phenomenex Desc. Kinetex C18 (2.1 mm × 100 mm, 2.6 μm particle size, Los Angeles, CA, USA) were tested. The other two columns required a longer running time than did the Agilent ZORBOX SB-Aq column. Thus, A, C, G, U, and m6A could be well-separated on Agilent ZORBOX SB-Aq column in 8 min.

The linearity, LOD, and LOQ were obtained using the peak areas of the product ion with the target SIM mode. As shown in Table 1, the linearity was evaluated by preparing different calibration curves (A, C, G, U, and m6A). The regression equations and determination coefficients (R2) of all the standard solution curves indicated that all the target compounds showed excellent linearity (R2 ≥ 0.9989 in all cases) [see Additional file 1: Figures S1–S5]. The LODs for A, C, G, U, and m6A were 0.0010, 0.0045, 0.0030, 0.0476, and 0.0003 μM, respectively. The LOQs for A, C, G, U, and m6A were 0.0032, 0.0148, 0.0100, 0.1587, and 0.0009 μM, respectively.

Table 1 Calibration equations, determination coefficients (R2), LODs, and LOQs of A, C, G, U, and m6A in the solvent

The intra- and inter-day accuracy and precision were evaluated by analysing three different concentrations of standard solutions of nucleosides. As shown in Table 2, the method was satisfactorily accurate and precise in measuring the nucleosides. The intra-day recoveries for the five ribonucleosides ranged from 86.9 to 114.8%, and the RSDs for all the analytes ranged from 1.2 to 9.5%. We further investigated the stabilities of analytes present in mRNA after three cycles of freeze (at − 20 °C for 24 h) and thaw (to room temperature). The analytes were reasonably stable under freezing/thaw conditions, as reflected by the inter-day recovery of 89.3%–106.3% and the RSD range of 1.0 to 9.0%.

Table 2 Precision and accuracy for the measurements of A, C, G, U, and m6A from three validation runs (2 μL of aliquot was injected in each run)

An UHPLC−HR − MS/MS method was developed for the accurate assessment of A, C, G, U, and m6A in mRNA isolated from N. tabacum. Figure 1 and [see Additional file 1: Figures S6–S8] showed the representative UHPLC−HR − MS/MS results for the quantifications of A, C, G, U, and m6A in N. tabacum with selected-ion chromatograms. As shown in Fig. 1 and [see Additional file 1: Figures S6–S8], the characteristic ions at m/z 282.11969, m/z 284.09897, m/z 268.10406, m/z 244.09282, and m/z 245.07687 were assigned to m6A [C11H15O4N5 + H] +, G [C11H15O5N5 + H] +, A [C10H13O4N5 + H] +, C [C9H13O5N3 + H] +, U [C9H12O6N2 + H] +, respectively. All the mass errors of precursors were less than 0.25 ppm.

Fig. 1

Representative UHPLC-HR-MS/MS results for G (a) and m6A (b) in N. tabacum. TIC MS shows the separated components. FTMS + p ESI SIM MS shows the chromatogram of objective component and its precise mass, resolution, and mass error

TMV induced a decrease in m6A level in N. tabacum

In this study, the ratio of m6A/G in N. tabacum was recorded during TMV infection. TMV reduced the m6A level in N. tabacum. As shown in Fig. 2, the m6A methylation levels decreased by 12.30, 12.66, 23.82, 12.02, 37.34, and 40.00% for 3, 5, 7, 10, 14, and 21 days, respectively, compared to the control. The level of m6A decreased most significantly at 21 days after TMV infection. In addition, the tobacco leaves presented dysmorphosis as time went on, especially the occurring characteristic mosaic at 21 days. This finding suggested that m6A played an important role in the growth of N. tabacum and that TMV reduced the m6A level in N. tabacum.

Fig. 2

m6A levels after infection with TMV after 14 and 21 days compared with control. The tobacco leaves present clear dysmorphosis at day 14 and appear characteristically mosaic at day 21. Each experiment represents the mean ± SD of at least three duplicates (n = 3).*P ≤ 0.05, **P ≤ 0.01 vs. control

Conservation of a family of ALKB in diverse species

To analyse the conservation of ALKB gene family, we compared amino acid sequences from the ALKB family. The sequence alignment results showed that the potential methylases of Nicotiana sylvestris were methyltransferase-like protein 1 (LOC104232562) (XM_009785798), methyltransferase-like protein 1 (LOC104216329) (XM_009766348), and FKBP12-interacting protein of 37 kDa-like (LOC104224278) (XM_009775897). The potential demethylases of Nicotiana sylvestris were uncharacterized LOC104220112 (LOC104220112) (XM_009770927), uncharacterized LOC104245988 (LOC104245988) (XM_009801708), and uncharacterized LOC104245950 (LOC104245950) (XM_009801663).

As shown in Fig. 3, database analysis and protein sequence alignments revealed partial homology between human ALKBH5 (NP_060228), Arabidopsis (NP_001031793, putative demethylase), and N. sylvestris (XP_009800010, putative demethylase, which was encoded by XM_009801708).

Fig. 3

Conserved protein sequences of putative mRNA N6-methyl-adenosine demethylases. The red letters represent the consensus residues derived from the conserved amino acid at a given position. Identical residues are shown in purple, and blocks of similar residues are shown in blue. Weak similarities are depicted in green. The overall sequence identity was ~ 40% across all proteins, and the sequence similarity was almost 95% across this region

ALKB family expression analysis

In the RT − qPCR experiments, the concentration and 260/280 and 260/230 ratios of RNA were 1244.9–2397.8 ng/μL, 1.92–2.28, and 1.54–2.26, respectively.

As shown in Fig. 4, the gene expression level of the XM_009801708 increased at 14 and 21 days, and the expression levels of the tested potential methylases decreased. These results suggested that XP_009800010 was the most likely demethylase of N. tabacum.

Fig. 4

Gene expression levels of putative methylases and demethylases. The potential methylases are XM_009766347, XM_009766348, and XM_009775897. The potential demethylases are XM_009770927, XM_009801708, and XM_009801663 in N. tabacum.*P ≤ 0.05, **P ≤ 0.01 vs. control


LC − MS/MS was successfully used to identify and quantify post-transcriptionally modified nucleosides in RNA [9, 23,24,25]. In this study, an UHPLC−HR − MS/MS method was developed to detect A, C, G, U, and m6A in mRNA that was isolated from N. tabacum. This method showed several advantages compared with previously reported methods. First, only 0.50 g of N. tabacum leaves was needed to prepare RNA because of the high sensitivity. Second, the analytes could be simultaneously detected by UHPLC−HR − MS/MS under the target SIM mode within 8 min, which significantly increased its efficiency compared with previous methods [8, 10, 23, 25, 29]. Moreover, SIM mode could avoid any interference from the matrix. Third, more accurate detection results could be obtained by UHPLC−HR − MS/MS due to its high resolution.

In this experiment, we used the modified method of Desrosiers et al. [30, 31] for mRNA preparation. In this method, the poly(A) tract in eukaryotic mRNA was used to study mRNA methylation. The reduction of m6A in TMV-infected N. tabacum could occur in mRNA with poly(A) structure can be tested and not in rRNA and tRNA. In addition, no poly(A) was found in the 3′ ends of TMV mRNA [32]. Therefore, the prepared mRNA did not contain TMV mRNA.

m6A was the most prevalent methylated nucleoside in mRNA. The regulatory role of m6A in RNA metabolism included mRNA transcription, splicing, nuclear export, and translation ability and stability. The discovery of m6A demethylases of FTO and ALKBH5 revealed that RNA methylation was reversible [8, 9]. Abnormal m6A methylation levels might lead to RNA dysfunction and cause diseases [33]. For example, in Arabidopsis, the disruption of MTA led to an arrest at the globular stage of embryo development when the m6A level decreased [12]. In humans, a loss-of-function mutation in the FTO gene was responsible for a recessive lethal syndrome [34]. In mice, the loss of the FTO gene caused post-natal growth retardation, and FTO overexpression led to increased food intake and obesity [35, 36]. m6A mRNA methylation affected the mammalian circadian clock. The specific inhibition of m6A methylation elongates the circadian period [37], and considerable epidemiological evidence showed an association of circadian disruption and human diseases like obesity, diabetes, and cancer [38]. Corresponding to this result, we found that the m6A level decreased after TMV infection. It was interesting that m6A decreased significantly as time passed. For example, the m6A level exhibited no significant decrease at day 3 but had a large one at day 21. This may have occurred because there was a low TMV concentration day 3, and the TMV concentration subsequently increased with time.

Considering the close association of m6A with methylases and demethylases, we used RT-qPCR to study the gene expression levels and investigate two potential enzymes in N. tabacum leaves. This work revealed that TMV increased the expression level of XM_009801708 in N. tabacum. In addition, the resulting protein sequence alignments revealed that it was partially homologous to human ALKBH5 (NP_060228) and Arabidopsis (NP_001031793, putative demethylase). As shown in Fig. 3, the residues highlighted within the black frame were the candidate determinants for the m6A recognition, and the catalysis of ALKBH5 was relatively conserved in the ALKB family [39]. Our result agreed with this report, except for the second amino acid in these residues. This difference existed between human ALKBH5 and the other two. This finding might be due to the different classifications of the species: humans are mammals, but Arabidopsis and N. sylvestris are plants. These results suggested that TMV was associated with m6A. We thus hypothesise that TMV directly or indirectly affects the demethylases in N. tabacum. Demethylases regulate the m6A level, and m6A regulates mRNA metabolism. Consequently, the tobacco mRNA metabolism is disrupted and causes disease. Thus, we will attempt to prove this hypothesis as a next step.

In this study, the findings that mRNA expression of the putative demethylase is increased after TMV infection, and the putative methylatransferase is decreased, may support the finding that levels of m6A in N. tabacum appear to be reduced by TMV infection. However, many functional experiments are still required to show that changes in the m6A levels after infection contribute to disease. Furthermore, the expression of the putative m6A methylases and demethylases should be genetically modulated to show that these genes are indeed part of the m6A machinery and control the m6A levels in N. tabacum. In addition, actin composes microfilament which plays an important role in cytoskeleton. β-actin is one of the actin, serve as a recognized reference gene in PCR in general. Some researchers proposed that when the tobacco was infected with TMV, actin cytoskeleton may play a role in delivering and attaching the MP to the putative cell wall adhesion sites and the plasmodesmata [40]. This conclusion showed that actin could not be a good reference gene when tobacco infected with TMV. However, Hofmann, et al. demonstrated that TMV movement did not require an intact actomyosin system [41]. Moreover, Baeka, et al. reported that actin still was one of the most stable genes as inference gene post inoculation with TMV, CMV, PVX, and PVY in tobacco [42]. Furthermore, some researchers still use actin as reference gene in tobacco infected TMV [43]. In future study, we will compare the difference of reference gene in TMV infection.


We developed a rapid method, UHPLC−HR − MS/MS, to determine the A, C, G, U, and m6A levels in mRNA in N. tabacum. TMV promotes ALKBH5-dependent m6A demethylation. The reversible m6A modification in N. tabacum mRNA might represent a novel epigenetic mechanism involved in TMV.




A. thaliana :

Arabidopsis thaliana




Diethylpyrocarbonate H2O


Fat mass and obesity-associated protein




High-performance liquid chromatography


Limit of detection


Limit of quantification




Methyltransferase like 14


Methyltransferase like 3

N. tabacum :

Nicotiana tabacum

R 2 :

Regression equations and determination coefficients


Relative standard deviation


Reverse transcription quantitative real-time polymerase chain reaction


Tobacco mosaic virus




Ultra-high-performance liquid chromatography coupled with high resolution tandem mass spectrometry


Wilms tumour 1 associated protein


  1. 1.

    Deng X, Chen K, Luo GZ, Weng XC, Ji QJ, Zhou TH, He C. Widespread occurrence of N 6-methyladenosine in bacterial mRNA. Nucleic Acids Res. 2015;43:6557–67.

  2. 2.

    Schwartz S, Agarwala SD, Mumbach MR, Jovanovic M, Mertins P, Shishkin A, et al. High-resolution mapping reveals a conserved, widespread, dynamic mRNA methylation program in yeast meiosis. Cell. 2013;155:1409–21.

  3. 3.

    Luo GZ, MacQueen A, Zheng GQ, Duan HC, Dore LC, Lu ZK, et al. Unique features of the m6A methylome in Arabidopsis thaliana. Nat Commun. 2014;5:5630.

  4. 4.

    Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3' UTRs and near stop codons. Cell. 2012;149:1635–46.

  5. 5.

    Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485:201–6.

  6. 6.

    Liu JZ, Yue YN, Han DL, Wang X, Fu Y, Zhang L, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N 6-adenosine methylation. Nat Chem Biol. 2014;10:93–5.

  7. 7.

    Ping XL, Sun BF, Wang L, Xiao W, Yang X, Wang WJ, et al. Mammalian WTAP is a regulatory subunit of the RNA N 6-methyladenosine methyltransferase. Cell Res. 2014;24:177–89.

  8. 8.

    Jia GF, Fu Y, Zhao X, Dai Q, Zheng GQ, Yang Y, et al. N 6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO. Nat Chem Biol. 2011;7:885–7.

  9. 9.

    Zheng GQ, Dahl JA, Niu YM, Fedorcsak P, Huang CM, Li CJ, et al. ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol Cell. 2013;49:18–29.

  10. 10.

    Clancy MJ, Shambaugh ME, Timpte CS, Bokar JA. Induction of sporulation in Saccharomyces cerevisiae leads to the formation of N 6-methyladenosine in mRNA: a potential mechanism for the activity of the IME4 gene. Nucleic Acids Res. 2002;30:4509–18.

  11. 11.

    Zhong SL, Li HY, Bodi Z, Button J, Vespa L, Herzog M, et al. MTA is an Arabidopsis messenger RNA adenosine methylase and interacts with a homolog of a sex-specific splicing factor. Plant Cell. 2008;20:1278–88.

  12. 12.

    Hongay CF, Orr-Weaver TL. Drosophila inducer of MEiosis 4 (IME4) is required for notch signaling during oogenesis. Proc Natl Acad Sci U S A. 2011;108:14855–60.

  13. 13.

    Horiuchi K, Umetani M, Minami T, Okayama H, Takada S, Yamamoto M, et al. Wilms’ tumor 1-associating protein regulates G2/M transition through stabilization of cyclin A2 mRNA. Proc Natl Acad Sci U S A. 2006;103:17278–83.

  14. 14.

    Fukusumi Y, Naruse C, Asano M. Wtap is required for differentiation of endoderm and mesoderm in the mouse embryo. Dev Dyn. 2008;237:618–29.

  15. 15.

    Bodi Z, Zhong SL, Mehra S, Song J, Graham N, Li HY, et al. Adenosine methylation in Arabidopsis mRNA is associated with the 3′ end and reduced levels cause developmental defects. Front Plant Sci. 2012;3:48–56.

  16. 16.

    Shen LS, Liang Z, Gu XF, Chen Y, Teo ZWN, Hou XL, et al. N 6-methyladenosine RNA modification regulates shoot stem cell fate in Arabidopsis. Dev Cell. 2016;38:186–200.

  17. 17.

    Sierro N, Batte JND, Ouadi S, Bakaher N, Bovet L, Willig A, et al. The tobacco genome sequence and its comparison with those of tomato and potato. Nat Commun. 2014;5:3833.

  18. 18.

    Hu Q, Niu YB, Zhang K, Liu Y, Zhou XP. Virus-derived transgenes expressing hairpin RNA give immunity to Tobacco mosaic virus and Cucumber mosaic virus. Virol J. 2011;8:41.

  19. 19.

    Scholthof KG, Adkins S, Czosnek H, Palukaitis P, Jacquot E, Hohn T, et al. Top 10 plant viruses in molecular plant pathology. Mol Plant Pathol. 2011;12:938–54.

  20. 20.

    Keith G. Mobilities of modified ribonucleotides on two-dimensional cellulose thin-layer chromatography. Biochimie. 1995;77:142–4.

  21. 21.

    Bodi Z, Button JD, Grierson D, Fray RG. Yeast targets for mRNA methylation. Nucleic Acids Res. 2010;38:5327–35.

  22. 22.

    Cornelius MG, Schmeiser HH. RNA analysis by MEKC with LIF detection. Electrophoresis. 2007;28:3901–7.

  23. 23.

    Chen BE, Ye F, Yu L, Jia GF, Huang XT, Zhang XJ, et al. Development of cell-active N 6-methyladenosine RNA demethylase FTO inhibitor. J Am Chem Soc. 2012;134:17963–71.

  24. 24.

    Yan MH, Wang YG, Hu YN, Feng Y, Dai CG, Wu JX, et al. A high-throughput quantitative approach reveals more small RNA modifications in mouse liver and their correlation with diabetes. Anal Chem. 2013;85:12173–81.

  25. 25.

    Fu LJ, Amato NJ, Wang PC, McGowan SJ, Niedernhofer LJ, Wang YS. Simultaneous quantification of methylated cytidine and adenosine in cellular and tissue RNA by nano-flow liquid chromatography-tandem mass spectrometry coupled with the stable isotope-dilution method. Anal Chem. 2015;87:7653–9.

  26. 26.

    Midoh N, Iwata M. Expression of defense-related genes by probenazole or 1, 2-benziothiazole-3 (2H) -one 1, 1-dioxide. J Pestic Sci. 1997;22:45–7.

  27. 27.

    Zhang Y, Wen LY, Yang AG, Luo CG, Cheng LR, Jiang CH, et al. Development and application of a molecular marker for TMV resistance based on N gene in tobacco (Nicotiana tabacum). Euphytica. 2017;213:259–67.

  28. 28.

    Guidance for industry bioanalytical method validation. 2017. Accessed 20 Dec 2017.

  29. 29.

    Rana AP, Tuck MT. Analysis and in vitro localization of internal methylated adenine residues in dihydrofolate reductase mRNA. Nucleic Acids Res. 1990;18:4803–8.

  30. 30.

    Desrosies R, Frideric K, Rottman F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc Natl Acad Sci U S A. 1974;71:3971–5.

  31. 31.

    Jelinek W, Adesnik M, Salditt M, Sheiness D, Wall R, Molloy G, et al. Further evidence on the nuclear origin and transfer to the cytoplasm of polyadenylic acid sequences in mammalian cell RNA. J Mol Biol. 1973;75:515–32.

  32. 32.

    Li WM, Zhang YQ, Zhang C, Pei XW, Wang ZX, Jia SR. Presence of poly(A) and poly(A)-rich tails in a positive-strand RNA virus known to lack 3′ poly(A) tails. Virology. 2014;454-455:1–10.

  33. 33.

    Niu YM, Zhao X, Wu YS, Li MM, Wang XJ, Yang YG. N 6-methyl-adenosine (m6A) in RNA: an old modification with a novel epigenetic function. Genom Proteom Bioinform. 2013;11:8–17.

  34. 34.

    Boissel S, Reish O, Proulx K, Kawagoe-Takaki H, Sedgwick B, Yeo GS, et al. Loss-of-function mutation in the dioxygenase-encoding FTO gene causes severe growth retardation and multiple malformations. Am J Hum Genet. 2009;85:106–11.

  35. 35.

    Fischer J, Koch L, Emmerling C, Vierkotten J, Peters T, Bruning JC, et al. Inactivation of the FTO gene protects from obesity. Nature. 2009;458:894–8.

  36. 36.

    Church C, Moir L, McMurray F, Girard C, Banks GT, Teboul L, et al. Overexpression of FTO leads to increased food intake and results in obesity. Nat Genet. 2010;42:1086–92.

  37. 37.

    Fustin JM, Doi M, Yamaguchi Y, Hida H, Nishimura S, Yoshida M, et al. RNA-methylation-dependent RNA processing controls the speed of the circadian clock. Cell. 2013;155:793–806.

  38. 38.

    Zelinski EL, Deibel SH, McDonald RJ. The trouble with circadian clock dysfunction: multiple deleterious effects on the brain and body. Neurosci Biobehav Rev. 2014;40:80–101.

  39. 39.

    Feng C, Liu Y, Wang GQ, Deng ZQ, Zhang Q, Wu W, et al. Crystal structures of the human RNA demethylase ALKBH5 reveal basis for substrate recognition. J Biol Chem. 2014;289:11571–83.

  40. 40.

    Heinlein M, Padgett HS, Gens JS, Pickard BG, Casper SJ, Epel BL, et al. Changing patterns of localization of the tobacco mosaic virus movement protein and replicase to the endoplasmic reticulum and microtubules during infection. Plant Cell. 1998;10:1107–20.

  41. 41.

    Hofmann C, Niehl A, Sambade A, Steinmetz A, Heinlein M. Inhibition of tobacco mosaic virus movement by expression of an actin-binding protein. Plant Physiol. 2009;149:1810–23.

  42. 42.

    Baeka E, Yoon JY, Palukaitis P. Validation of reference genes for quantifying changes in gene expression in virus-infected tobacco. Virology. 2017;510:29–39.

  43. 43.

    Chen Z, Zeng MJ, Song BA, Hou CR, Hu DY, Fan HT, et al. Dufulin activates HrBP1 to produce antiviral responses in tobacco. PLoS One. 2012;7:e37944.

Download references


We are grateful to the group of Guifang Jia (Department of Chemical Biology and Synthetic and Functional Biomolecules Center, College of Chemistry and Molecular Engineering, Peking University) for their assistance with the analysis of m6A and Collaborative Innovation Center for Natural Products and Biological Drugs of Yunnan.


This research was financially supported by the National Natural Science Foundation of China (Nos. 21562013 and 21732002) and Subsidy Project for Outstanding Key Laboratory of Guizhou Province in China (No. 20154004).

Availability of data and materials

All data generated or analysed during this study are available from the corresponding author on reasonable request.

Author information

BAS and DYH designed this study. ZRL and XZZ collected the samples; ZRL, XZZ and LY performed the analyst tests; ZRL, JS, and LLR performed the RT-PCR tests; ZRL and LY analysed the data; ZRL wrote the paper. All authors read and approved the final manuscript.

Correspondence to Deyu Hu or Baoan Song.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Figure S1. Calibration curves for the quantification of A in mRNA. The amounts of A ranged from 0.1–2.5 μmol. Figure S2. Calibration curves for the quantification of C in mRNA. The amounts of C ranged from 0.5–25 μmol. Figure S3. Calibration curves for the quantification of G in RNA. The amounts of G ranged from 0.25–12.5 μmol. Figure S4. Calibration curves for the quantification of U in RNA. The amounts of U ranged from 5 to 250 μmol. Figure S5. Calibration curves for the quantification of m6A in RNA. The amounts of m6A ranged from 0.005–0.250 μmol. Figure S6. Representative UHPLC-HR-MS/MS results for the quantification of A in N. tabacum K326 with selected-ion chromatograms. Figure S7. Representative UHPLC-HR-MS/MS results for the quantification of C in N. tabacum K326 with selected-ion chromatograms. Figure S8. Representative UHPLC-HR-MS/MS results for the quantification of U in N. tabacum K326 with selected-ion chromatograms. Table S1. Primer sequences for RT-qPCR analysis. (PDF 1250 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Z., Shi, J., Yu, L. et al. N6-methyl-adenosine level in Nicotiana tabacum is associated with tobacco mosaic virus. Virol J 15, 87 (2018) doi:10.1186/s12985-018-0997-4

Download citation


  • N. tabacum
  • TMV
  • m6A
  • RT-qPCR