- Research article
- Open Access
Expression of geminiviral AC2 RNA silencing suppressor changes sugar and jasmonate responsive gene expression in transgenic tobacco plants
BMC Plant Biology volume 12, Article number: 204 (2012)
RNA-silencing is a conserved gene regulation and surveillance machinery, which in plants, is also used as major defence mechanism against viruses. Various virus-specific dsRNA structures are recognized by the silencing machinery leading to degradation of the viral RNAs or, as in case of begomoviruses, to methylation of their DNA genomes. Viruses produce specific RNA silencing suppressor (RSS) proteins to prevent these host defence mechanisms, and as these interfere with the silencing machinery they also disturb the endogenous silencing reactions. In this paper, we describe how expression of AC2 RSS, derived from African cassava mosaic geminivirus changes transcription profile in tobacco (Nicotiana tabacum) leaves and in flowers.
Expression of AC2 RSS in transgenic tobacco plants induced clear phenotypic changes both in leaves and in flowers. Transcriptomes of these plants were strongly altered, with total of 1118 and 251 differentially expressed genes in leaves and flowers, respectively. The three most up-regulated transcript groups were related to stress, cell wall modifications and signalling, whereas the three most down-regulated groups were related to translation, photosynthesis and transcription. It appears that many of the gene expression alterations appeared to be related to enhanced biosynthesis of jasmonate and ethylene, and consequent enhancement of the genes and pathways that are regulated by these hormones, or to the retrograde signalling caused by the reduced photosynthetic activity and sugar metabolism. Comparison of these results to a previous transcriptional profiling of HC-Pro RSS-expressing plants revealed that some of same genes were induced by both RSSs, but their expression levels were typically higher in AC2 than in HC-Pro RSS expressing plants. All in all, a large number of transcript alterations were found to be specific to each of the RSS expressing transgenic plants.
AC2 RSS in transgenic tobacco plants interferes with the silencing machinery. It causes stress and defence reactions for instance via induction of the jasmonate and ethylene biosynthesis, and by consequent gene expression alteration regulated by these hormones. The changed sugar metabolism may cause significant down-regulation of genes encoding ribosomal proteins, thus reducing the general translation level.
Gene silencing in plants
The conserved molecular machinery of RNA-silencing constitutes a very complex genetic regulatory network in all eukaryotes. The common features for these regulatory pathways are their induction by double stranded (ds) RNA sequences. In plants, these are cleaved by the RNAse III type DICER-LIKE (DCL) enzymes, assisted by HYPONASTI LEAVES1 (HYL1), SERRATE (SE) and DAWLE (DLL) proteins  into 21–25 nucleotide long RNAs, called either (micro-) miRNAs or (short interfering) siRNA. Upon cleavage, the short RNA fragments are methylated at their 3’ ends by the S-adenosyl methionine dependent methyltransferase HEN1 , and loaded into effector complexes, called RITC (RNA induced transcription silencing complex) or transported to cytoplasm apparently with the assistance of HASTY protein  and loaded RISC (RNA induced silencing complex) reviewed in . The RITS complexes contain siRNA, AGO4 and AGO6 to mediate transcriptional silencing (TGS) to repetitive or over-expressed DNA sequences via histone and DNA methylation, leading to heterochromatin maintenance, control of transposon mobility and transgene silencing. This RNA-dependent DNA methylation (RdDM) machinery utilizes at least the DNA-directed RNA polymerases PolIV and PolV, RNA-dependent RNA polymerase 2 (RDR2), Dicer-like 3 (DCL3) enzyme, AGO4, AGO5, AGO6 and DNA methyltransferase (DRM2) , chromatin-remodelling protein (DRD1), structural-maintenance-of-chromosomes-protein (DMS3) and RdDM effector molecule KTF1, which binds scaffold transcripts generated by PolV and recruits AGO4 bound siRNAs to form an RdDM effector complexes [6, 7].
The RISC complexes, targeted to RNA sequences, cause the post-transcriptional silencing (PTGS) of the target mRNAs either via their cleavage or translational arrest . MiRNA-mediated PTGS is used, for instance, to regulate the temporal and spatial expression of multiple different transcription factors, many of which are needed in the developmental differentiation of plant organs, and also regulate various components of the silencing pathways, thus back-regulating these pathways themselves. In addition, numerous miRNAs, as well as a plethora of various siRNAs are known to be involved in genetic responses and signalling cascades induced by various hormones and biotic and abiotic stresses (reviewed in [9–13]).
The small RNAs target these silencing complexes to partially or fully complementary DNA or RNA sequences, respectively. On these target molecules the effector complexes mediate multiple essential regulatory functions like transcriptional gene silencing (TGS) or post-transcriptional gene silencing (PTGS) or decreasing the translational rate (reviewed by [8, 14]). In plants RNA silencing is used as the major defence mechanism against viruses.
Suppression of RNA Silencing in plants
To counteract the silencing mediated host defence, plant viruses encode various viral RSSs . Functional suppression of RNA silencing appears to be required for any virus to successfully invade and accumulate in its host plant, and several of the RSSs have been identified as essential viral pathogenicity factors and symptom determinants already before their function in RNA silencing suppression was discovered . It is of interest that the RSS factors encoded by different viral families have no similarities with each other, and many viruses produce their RSS activity as a secondary function of some gene product that has also some other function in viral life cycle, e.g. in replication, cell-to cell or systemic movement or encapsulation (reviewed in ). The high complexity of the silencing machinery provides multiple steps where it can be disturbed, and the mechanisms of the different RSS are very diverse. The hosts’ silencing machineries and their molecular interactions with viral RSSs have been intensively studied over the last decade [8, 10, 17–20] and they appear to be excellent tools for interrupting and analysing of the RNA silencing pathways.
Begomoviruses and silencing suppression in plants
Begomoviruses, with circular ssDNA genomes, produce transcripts in opposite orientation. They form partly dsRNA structures and activate plant’s RNA silencing machinery, but can resist this suppression by producing active RSSs [21–23]. The silencing is not (at least totally) targeted against viral RNA transcripts, but rather against the viral replicative intermediate, i.e. the minichromosome structures formed of replicative viral dsDNA, combined with histone proteins . The silencing is mediated by the TGS, i.e. by methylation of these viral dsDNA sequences, which leads to repression of their replication and transcription . Total of three different types of RSS are known to be encoded by African cassava mosaic geminivirus (ACMV), which nowadays is classified as a begomovirus [23, 25]. The main RSS has been identified as the transcriptional activator protein AC2, encoded in opposite sense of the begomoviral minichromosome.
Here we have analysed the transcriptome and proteome of the transgenic tobacco plants expressing the AC2 RSS, derived from ACMV begomovirus. Alterations detected in these gene expression profiles were compared to the those that we have earlier detected in similar transgenic tobacco plants that express the HC-Pro RSS derived from Potato virus Y potyvirus, revealing the fundamental similarities and differences between the functional mechanisms of these two RSS in transgenic tobacco plants.
Phenotype of tobacco plants expressing AC2 RSS
A previously characterized transgenic tobacco line expressing Geminiviral AC2 RSS under a constitutive CaMV 35S promoter was used in this study . Transgenic tobacco plants expressing AC2 RSS (Additional file 1) were grown concurrently under the same growth conditions that were earlier used for HC-Pro expressing transgenic tobacco plants . Expression of nuclear AC2 protein in transgenic tobacco plants induced clear phenotypic changes both in leaves and flowers, including different malformations of veins and leaves of the young transgenic plants (Figure 1A and 1B). These malformations gradually disappeared as plants become older, but there still were some phenotypic changes in six week-old leaves; the leaves were narrower and paler green than leaves in wild type plants (Figure 1C). The growth habit was stubby. Lower leaves were often yellow, the root neck and root system was less developed than in wild type plants, and mature plants were slightly smaller in size compared to wild type tobacco plants. Flower buds were sometimes twisted round other flower buds (Figure 1E). The structure of flowers were often, but not always drastically malformed; petals were grown together having either round or more often a rectangular shape. Stamen filaments were also often changed to extra petals (Figure 1G-1I).
Experimental design and differential gene expression in AC2 expressing plants
The gene expression profiling of the AC2 RSS expressing plants was performed at the same time as that of the HC-Pro expressing transgenic tobacco plants by using the same wild type and empty vector (pBIN61) expressing control plants as were used for the transgenic HC-Pro expressing plants, using three biological replicates for each. The intensity values indicating differential gene expression between control and the AC2 expressing plants were normalized using the Chipster program (CSC, Espoo, Finland). The intensity values originating from the leaf and flower samples were normalized separately as was performed for the HC-Pro samples , (Additional file 2). Two-fold change in the expression between AC2 transgenic and control samples was regarded as significant level of differential expression. Statistical significance of differentially expressed genes was tested using Students t-test, p values less than 0.05 with False Discovery Rate (FDR) (See Additional file 3 Additional file 4 Additional file 5 Additional file 6). Microarray results were verified by using RT-qPCR. with a good correlation between RT-qPCR and array results both in up- and –down-regulated transcripts, and both in leaves and flowers (Table 1). Transcripts, whose expression did not change, were used as reference transcripts .
The numbers of differentially expressed transcripts, assigned to various functional groups, in leaves and flowers of transgenic AC2 expressing plants is presented in Table 2. The total number of altered transcripts was 1118 and 251 in leaves and flowers, respectively. In leaves total of 726 transcripts were up-regulated and 392 down-regulated, with the three biggest groups of up-regulated transcripts being related to signalling, cell wall modifications, and stress, whereas the three most down-regulated groups of transcripts being related to protein synthesis, photosynthesis and transcriptional regulation. The results demonstrated in Table 2 were based on counting of genes in functional categories similarly as it was performed in our previous publication  but not the overpresentation analysis.
About half of the transcripts up-regulated in leaves of the HC-Pro expressing plants were found up-regulated also in the leaves of the transgenic AC2 expressing plants (Figure 2), but more than five hundred additional transcripts were up-regulated only in the AC2 expressing plants. The same pattern was also found in the down-regulated transcripts of the leaves; more than three hundred transcripts were down-regulated only in AC2 expressing plants. The number of differentially expressed transcripts was much lower in flowers but these also included common and specific transcripts between AC2 and HC-Pro expressing plants (Figure 2). See also Additional file 7 for over and under presentation of functional categories between HC-Pro and AC2 RSS expressing plants.
Expression of AC2 in tobacco plants induces clear defence responses
The expression of AC2 RSS in transgenic plants induced many defence- and oxidative stress- related transcripts, including transcripts for systemic acquired resistance (SAR), ROS scavenging and respiratory burst oxidase (Table 3, Additional file 8) which were not detected in HC-Pro RSS expressing transgenic plants . The amount of oxidative stress was further studied using NBT staining for superoxide radical and DAB for H2O2 to indicate the amount different ROSes in leaves of AC2 expressing plants. The amount of ROSes was found increased only in one-month old AC2 expressing tobacco plants (Figure 3), while the amount of ROS accumulation in older and more developed leaves was similar as in wild-type leaves (data not shown).
Interestingly, hypersensitive reaction (HR, or necrotic lesions) was frequently detected in lower leaves of two month old transgenic AC2 expressing plants (Figure 3F-3H). These kinds of necrotic lesions were never detected in leaves of wild type tobaccos of similar age. This suggested that particularly the young developing leaves suffered from oxidative stress where also the enhanced levels of ROS-related transcripts were detected.
Jasmonic acid biosynthesis and jasmonate-responsive genes
Multiple genes related to jasmonic acid (JA) biosynthesis were up-regulated in leaves of transgenic AC2 expressing plants. Some of these genes were also up-regulated in HC-Pro expressing plants, but in these the level of up-regulation was much lower . Also genes encoding jasmonic acid-mediated regulatory proteins like NtJAZ1 and NtJAZ3 repressors were up-regulated (Table 3, Additional file 7). The presence of jasmonates (e.g.JA-Ile) induces proteasomal degradation of these JAZ repressors and thereby activates transcription e.g. from promoters containing MYC binding sites (reviewed in ). Many of the JA-responsive genes e.g. defence and stress-related transcripts were also up-regulated. Jasmonic acid signalling cascade also positively regulates the genes involved in biosynthesis of jasmonic acid itself (from chloroplast 18:3 fatty acids to JA) [29, 30], and in the biosynthesis of ethylene [31, 32]. Table 3 shows almost full induction of jasmonic acid biosynthesis related transcripts, jasmonic acid regulatory genes as well as some well-known jasmonic acid responsive genes, including the 1-aminocyclopropane-1-carboxylic acid oxidase (ACC oxidase) transcript, required for the biosynthesis of ethylene from its ACC precursor [31, 32] in leaves of transgenic AC2 expressing plants. In addition, a transcript of ethylene response factor 1 (ERF1), a transcription factor responsible for mediating ethylene responsive gene expression  was up-regulated about five-fold.
Chlorophyll biosynthesis and photosynthesis related genes were down–regulated in AC2 expressing transgenic plants
Most of the phenotypic changes occurred in young, less than two-month old plants, but not many differences were detected in older leaves of AC2 expressing plants. Leaves of young transgenic plants (up to two-month old) were more yellow than those of the wild type tobaccos, indicating changes in their pigment content. The pigment analysis revealed that the amount of both chlorophyll a and b were decreased in them (Additional file 9), and there was also a clear decrease in the accumulation of anthocyanin in the same plants, as compared to wild type plants. The decreased amount of chlorophyll in leaves was not surprising since many of the genes related to chlorophyll biosynthesis were clearly down-regulated (Table 4). Table 4 also depicts down-regulation of transcripts of photosystem II and I related genes, as well transcripts related to primary carbon metabolism.
Since there were clear changes in the genes encoding photosynthetic machinery proteins, the light-responsive oxygen evolution was measured from the thylakoid samples isolated from AC2 expressing and wild type leaves. Both photosynthetic oxygen evolution in freshly isolated thylakoids and the amount of starch were significantly decreased in transgenic AC2 expressing leaves, as compared to wild type tobacco leaves (Figures 4 and 5). The results were similar, but clearer than those published earlier for HC-Pro RSS expressing plants .
Down-regulated accumulation of translation factor and ribosomal protein transcripts
One of the most dramatic changes in the transcriptome of the transgenic AC2 expressing plants was the down regulation of genes encoding for translation factors and ribosomal protein subunits (Table 1 and Table 5, Additional file 7). Total of fifty two transcripts related to translation machinery, including genes encoding for the cytosolic large (60S) and small (40S) ribosomal protein subunits as well as chloroplast ribosomal protein subunits (50S and 30S), were down-regulated. There were also changes in the distribution of ribosomal RNAs (Figure 6), as the amount and ratio of 23S rRNA was clearly decreased in AC2 expressing tobacco plants. Consequently, the amount of total protein, as measured against fresh weight was decreased about 30% in AC2 expressing leaves compared to wild type tobacco leaves (Additional file 10). The proteomes of wild type and AC2 expressing tobacco leaves were compared by loading equal amounts of protein extracts (250μg) into the 2D-SDS-PAGE analysis (Figure 7), meaning that the analysis actually enhanced the protein spot intensities of the AC2-plant samples, as their total protein content was really lower than in the wild type plants. The AC2-plant proteome contained multiple changes, as compared to the wild type plant sample. To our surprise, the changes occurring in the proteomes of AC2 and HC-Pro expressing tobacco plants were quite similar to those observed in the HC-Pro expressing plants , although in general more proteins were down-regulated in AC2 expressing than in and HC-Pro expressing tobacco leaf.
Expression of silencing related transcripts
Interestingly, the microarray data indicated that some silencing related transcripts involved in the DNA methylation were altered (Table 6). Two of the up-regulated transcripts (KTF1 and AGO5) have been described to function in RNA-directed DNA methylation [33, 34]. Based on its ability to bind 24nt small RNAs, AGO5 could influence small viral RNA-directed DNA methylation (RdDM). The AGO5 has an effect on CG motifs rather than C residues in a CNG or CNN context. Counteracting the increased DNA methylation by these gene products, also Repressor of Silencing1 (NtROS1), functioning as a negative regulator of silencing was up-regulated. In addition, two chromatin methyl transferases and a structural maintenance of chromosomes, also suggested to function in de novo DNA methylation, were down-regulated.
To investigate possible alterations in the level of DNA methylation level, we also performed PCR-amplification experiments over some genomic DNA regions that contained at least three methylation sensitive restriction enzyme sites. These analyses indicated increased methylation at least in the coding regions of ERF1 and AP24 genes but not in 18S RNA (Figure 8). However, other gene fragments showed no change in the methylation status (data not shown).
Neither genes encoding S-adenosyl-L-methionine (SAM) biosynthesis nor genes encoding SAM-recycling were down regulated as much as in HC-Pro expressing plants . Thus, it appears that recycling of SAM-provided methyl groups for trans-methylation reactions was not affected in AC2-expressing plants. Interestingly, a large group of histone transcripts were down-regulated in AC2 expressing tobacco plants (Table 6, Additional file 12), while a similar group of histone transcripts were up-regulated in HC-Pro expressing tobacco plants. Also some histone acetyl-transferases were up-regulated, possibly leading to opening of the chromatin structures and thus increasing transcriptional activity. Another interesting phenomenon was the down-regulation of cell cycle-related transcripts in AC2 expressing plants, as similar transcripts were up-regulated also in the HC-Pro expressing plants (Additional file 12).
In this work we have analysed the transcriptome profile of transgenic tobacco plants expressing the AC2 RSS derived from ACMV, nowadays classified as a member of the genus Begomovirus, and shown that this single RSS protein causes massive changes in the gene expression. These changes may have been caused trough at least three identified regulatory nodes, which all interfere with large metabolic networks. The first node is the induction of plant hormones like jasmonates and ethylene, leading to jasmonate and ethylene responsive gene expression. Both jasmonates and ethylene are known to be highly involved in defence and stress regulated gene expression, and jasmonate hormones are also responsible for inducing oxidative stress [30, 35]. The second regulatory node appears to be related to photosynthetic end products leading to changes in the transcriptional and translational regulation. The role of chloroplasts in this retrograde signalling is crucial. Results demonstrated that photosynthetic light and dark reactions were clearly down-regulated. The third regulatory node appears to consist of silencing directed modifications in cell cycle including changes in transcripts involved in DNA replication and chromatin structures.
Defence and stress reactions in transgenic AC2 expressing tobacco plants
The expression of the AC2 RSS in transgenic tobacco plants induced more pronounced defence and stress related gene expression than was found in plants expressing HC-Pro RSS in our previous study . Although many similar defence and stress related transcripts were induced in these two types of transgenic plants, a high number of transcripts were enhanced only in the AC2 expressing plants, the most prominent difference being the induction of transcripts related to scavenging of ROS in these plants (Table 3, Additional file 8). These included transcripts coding for scavengers of hydrogen peroxide (H2O2) and of superoxide related oxygen radicals. The presence of ROS also was indicated by spontaneous HR lesions in the young leaves (Figure 3). Also, transcripts related to systemic acquired resistance (SAR) were induced, maybe via induction of avirulance (Avr) determinant related genes (Table 3) (reviewed in ). Balance in the redox homeostasis seems also be affected in AC2 RSS expressing transgenic tobacco plants, as several genes that function in balancing the redox state in the cells, like glutathione-S-transferases, glutaredoxins and thioredoxins were up-regulated in the leaves (Table 3). All in all, this may have affected the redox-responsive gene regulation in the AC2 RSS expressing tobacco cells .
Photosynthesis is also regulated by jasmonates
Photosynthesis is the key regulator of almost all metabolic and differentiation responses in the plant cells but it has also a major role in regulating gene expression by retrograde signalling. Our results clearly pointed out that many in photosynthesis related transcript (Table 4), and the photosynthetic oxygen evolution (Figure 4) were significantly decreased in AC2 expressing plants. Interestingly, jasmonate hormones have been shown to regulate photosynthesis related transcripts in Arabidopsis, to down-regulate the oxygen evolution (i.e. the Hill reaction activity) in isolated barley thylakoids , to up-regulate defence, stress, senescence and cell wall related transcripts and to causes oxidative stress in plants and thus also to change redox state of proteins . Thus all these observed effects in AC2 expressing plants could have been related to up-regulated jasmonate biosynthesis . Further on, the reduced photosynthesis may have changed the sugar metabolism and thus caused major secondary effects on retrograde signalling between the nucleus and the chloroplast. For instance, it may explain the changes in protein synthesis (e.g. ribosomal proteins, translation factors etc.), cell wall synthesis, and chromosome/histone modifications related transcripts (Tables 5 and 6, Additional file 7) in the AC2 expressing plants.
Several studies have shown that energy deprivation and different stresses change energy-dependent transcriptional regulation [40, 41] with drastic down-regulation of transcripts related to different biosynthetic processes and especially protein syntheses (e.g. translation factor and ribosomal protein encoding transcripts) and up-regulation of catabolic reactions to restore energy deprivation on the cells.
Protein synthesis affected in AC2 expressing plants
Results clearly indicated that expression of AC2 protein in transgenic plants had a major effect on transcripts related both protein synthesis and protein degradation (Figure 7). The down-regulation of protein synthesis caused to some extend by reduction of levels of multiple transcripts, but also the down-regulation of genes encoding both nuclear and plastid ribosomal proteins, ribosomal RNA or translation factors (Tables 45 and Additional file 7 and Additional file 12). Not only synthesis of proteins but also degradation of proteins was affected. Over and under presentation analysis (Additional file 7) indicated that both cysteine and aspartate proteases were up-regulated in leaves of AC2 RSS expressing plants. There were also indications of proteasomal protein degradation (Additional file 12). Interestingly, it has suggested that the geminiviral C2 protein might be the key player in proteasomal regulation of protein degradation [42, 43]. It was shown that C2 protein is capable of binding to proteins involved in proteasomal complexes e.g. SCF E3 ligase complexes and thus inhibiting jasmonate related gene expression in Arabidopsis (See also ). However, our results clearly indicate that the geminiviral protein AC2 has an opposite effect; enhancement of jasmonate signalling and jasmonate responsive gene regulation  and no effects to transcripts related to neither SA synthesis nor SA responsive gene expression .
Cell cycle, genome methylation and histone expression is altered in AC2 expressing transgenic plants
Begomoviruses are known to replicate their single-stranded DNA genomes through double-stranded DNA intermediates that are associated with cellular histone proteins to form minichromosomes . Thus it is important to Begomoviruses to induce the accumulation of the DNA replication machinery in mature plant cells. This is achieved most likely by modifying cell cycle and transcriptional regulation, for instance, the infection of cabbage leaf curl geminivirus (CaLCUV) has been shown to induce changes in the cell cycle in virus infected cells . Interestingly, our results also suggested that expression of AC2 RSS may have an effect on cell cycle (Additional file 12), providing replication-competent environment and preventing expression of genes needed for mitosis [48, 49].Our results demonstrated that the expression AC2 RRS in transgenic plants up-regulated transcripts involved in DNA methylation (KTF1 and AGO5) and down-regulated transcripts involved in histone synthesis (Table 6). This suggested that, although the transcriptional activator AC2 suppresses silencing (DNA methylation) in the native host(s) of ACMV, and apparently enhances the virus replication in these hosts, in tobacco (non-host to ACMV) it induces defence mechanisms, leading to opposite direction, i.e. to increase of methylation and decrease of the histone synthesis.
The expression of AC2 RSS in transgenic tobacco plants induced clear phenotypic changes in tobacco plants, even though tobacco is not considered as native host for ACMV. This study indicated several possible action mechanisms of the AC2 RSS. First both jasmonate and ethylene hormone responsive regulatory cascades were induced, secondly, the rate of photosynthesis and also the carbon metabolism were reduced. The stress and defence responsive gene expression was enhanced including both induction and response of reactive oxygen species (ROS). All these regulatory nodes may have been interconnected to one another. Also, changes occurred in regulation of cell cycle and transcription as well as changes in the genomic architecture like DNA methylation and abundance of histones. Thus it seems that the silencing suppression activity of AC2 in tobacco is indicated in the enhancement of the jasmonate and ethylene biosynthesis, known to be regulated and repressed via RNA silencing . Up-regulation of these plant hormones was related to induction of a large variety of stress-related genes, and severe developmental defects in the plants. Interestingly, the natural silencing suppression mechanism of the geminiviral RSS, i.e. reduction of the DNA methylation, did not occur in this non-host plant, but rather the opposite effect was observed, apparently related to the induced plant defences.
The wild type tobacco (Nicotiana tabacum), empty vector tobacco control (pBIN61), and transgenic tobacco plants expressing AC2 transgene were grown in greenhouse conditions at 60% relative humidity and 22°C, with a day/night regime of 16h light (150 μmol photons m-2s-1) and 8 h dark. Leaf samples (third leaf from the top) were taken from one-month-old plants, the plants were at that time about 20 centimetres of height. Leaf and flower samples were taken from the same plant. Flower samples were taken one day prior to opening. Both leaf and flower samples were directly frozen in liquid nitrogen and stored at −80°C.
RNA extraction, cDNA labelling, microarray hybridization and scanning of the microarray chips
Total RNA was isolated from leaves and flowers of wild type and transgenic plants using TRIsure-reagent (Bioline, UK) according to manufacturer’s recommendations. Total RNA was further purified using RNeasy clean up column (QIAGEN inc. USA). The cDNA labelling, the quality checking of total RNA and labelled cDNAs, hybridization on a Agilent’s 4x44K tobacco chip (Design ID 21113), washing and the scanning the chips was performed in a similar fashion, concurrently with already published HC-Pro microarray samples .
Microarray data analysis
All data handling was performed using Chipster, a visual program based on R Project for Statistical Computing program (Center of Scientific Calculating (CSC), Finland) (, Agi4x44k preprocess, Lopez-Romero, 2010). In order to compare intensity values of different samples, both wild type and empty vector control tobacco samples (6) vs. transgenic AC2 expressing leaf samples, (3), were normalized together. The flower samples were normalized together in a similar fashion. Normalization of three biological replicates was performed using median signal values and median background values. A background offset value (50) was added to prevent negative values during normalization. Normalization of the arrays was performed using a “quantile” parameter. The array results have been deposited into ArrayExpress with accession number E-MEXP-3724. Differentially expressed transcripts were re-annotated as was performed in earlier publication by Soitamo et al. 2011 . Microarray results were visualized by using MapMan program and enrichment analysis of functional categories by using PageMan program [52, 53].
Verification of differentially expressed genes
The array results were verified by using RT-qPCR according to MIQE guidelines . The RT-qPCR was performed from the same RNA samples as were previously used in microarray experiments. The cDNA was synthesized from 1 μg of purified leaf or flower total RNA using RevertAid H-Minus M-MuLV reverse transcriptase according to manufacturer’s recommendations (Product # EPO451, Fermentas). Produced cDNA was diluted and 3μl was used in RT-qPCR (Maxima SYBR Green /Fluorescein qPCR MasterMix (2X) (Product # KO242, Fermentas). The gene specific reference and sample primers used in RT-qPCR are listed in Additional file 11. For each three biological replicates, three-four technical replicates were run to minimize pipetting errors. RT-qPCR reactions were run in a 96-well plate containing both wild type (reference) and AC2 transgenic samples. The RT-qPCR was performed using Bio-RAD’s iQ5 machine. The results were calculated using the quantification cycle (Cq) method (delta delta Cq) according to Bio-RAD’s iQ5 default settings. All primer pairs produced only one peak in DNA melting curves indicating high specificity of primers. Standard error of mean (s.e) was also calculated of three to five biological replicates.
Detection of oxygen radicals in leaves
Hydrogen peroxide and superoxide anions were detected by using stains that react with these radicals. DAB (0.1mg/ml DAB pH to 3.8 using NAOH) was used to detect presence of hydrogen peroxide and NBT (0.1mg/ml NBT in 25mM HEPES/KOH pH 7.4) was used to detect presence of superoxide radicals . Staining whole leaves was performed in 20ml DAB or NBT solution in a small Petri dish. Leaf was cut from the stalk under staining solution to get a new surface area for stain infiltration. Leaves were then kept in the dark overnight. Next morning, chlorophyll was removed from the leaf using 96% ethanol. Ethanol was changed a few times, until all chlorophyll was removed (took at least 24h). Pictures of leaves were taken from which all chlorophyll was removed.
Photosynthetic and chlorophyll quantification measurements
Equal amount of intact wild type and HC-Pro transgenic tobacco leaves (1.0 g) were ground in an ice cold mortel in 4 ml of thylakoid isolation buffer (0.3 M sorbitol, 50 mM Hepes/KOH pH 7.4, 5 mM MgCl2, 1 mM EDTA and 1% BSA). Suspension was filtered through a Miracloth and 2 ml thylakoid suspension was pelleted in Eppendorf-centrifuge 12 000xg for 2 minutes. The pellet was resuspended into 100 μl of O2-electrode measuring buffer (0.3 M sorbitol, 50 mM Hepes/KOH pH 7.4, 5 mM MgCl2, 1 mM KH2PO4). Oxygen evolution was measured directly in a Clark type O2-electrode using 0.5 mM DCBQ as electron donor. The chlorophyll concentration was calculated according to Porra et al. . Samples in the cuvette were quantified based on equal amount of total chlorophyll. The anthocyanin concentration was measured from three 8 mm leaf discs of wild type and of AC2 expressing tobacco plants according to Neff and Chory .
Isolation of proteins and 2D-PAGE
Protein samples of leaves from wild type and HC-Pro expressing transgenic plants were isolated concurrently with the RNA isolation using TRIsure-reagent (Bioline). The protocol was adapted from TRIzol (Invitrogen inc. USA) and performed according to manufacturer’s recommendations. The protein concentration was measured using Lowry method. Proteins were first separated by Bio-Rad laboratories 7cm IPG strips pH 3–10 according to manufacturer’s recommendations. 250 μg of protein was loaded per a strip. Strips containing wild type and transgenic HC-Pro focused protein samples were then run simultaneously in a large gel in Protean II apparatus (Bio-Rad) to produce a similar mobility of focused proteins of both strips. Protein gels were then fixed and stained in colloidal Coomassie blue stain (PageBlue staining kit, Fermentas) according to manufacturer’s recommendations, destained and photographed.
Detection of DNA methylation using PCR amplification after DNA methylation sensitive restriction enzyme cutting
DNA was isolated from the leaves of wild type and AC2 expressing tobacco leaves according to protocol of Genomic DNA from Plant (Nucleospin Plant II, Macherye-Nagel Germany). DNA was cut with either BamHI or BamHI and HpaII restriction enzymes. The cut fragments were amplified using Taq polymerase with 2 mM MgCl2. Primers for amplification are given in Additional file 11. Primers were designed so that the fragments contained three CCGG sites inside amplified region. Methylation in the HpaII sites in the DNA fragments prevented the cleavage of the DNA and an amplified PCR product was detected.
Fang Y, Spector DL: Identification of nuclear dicing bodies containing proteins for microRNA biogenesis in living Arabidopsis plants. Curr Biol. 2007, 17 (9): 818-823. 10.1016/j.cub.2007.04.005.
Yang Z, Ebright YW, Yu B, Chen X: HEN1 recognizes 21–24 nt small RNA duplexes and deposits a methyl group onto the 2’ OH of the 3’ terminal nucleotide. Nucleic Acids Res. 2006, 34 (2): 667-675. 10.1093/nar/gkj474.
Park MY, Wu G, Gonzalez-Sulser A, Vaucheret H, Poethig RS: Nuclear processing and export of microRNAs in Arabidopsis. Proc Natl Acad Sci USA. 2005, 102 (10): 3691-3696. 10.1073/pnas.0405570102.
Vaucheret H: Plant ARGONAUTES. Trends Plant Sci. 2008, 13 (7): 350-358. 10.1016/j.tplants.2008.04.007.
Matzke M, Kanno T, Daxinger L, Huettel B, Matzke AJ: RNA-mediated chromatin-based silencing in plants. Curr Opin Cell Biol. 2009, 21 (3): 367-376. 10.1016/j.ceb.2009.01.025.
Wadsworth S, Dunoyer P: Plant RNA-silencing immunity and viral counter-defence strategies. pp. 1-35. In: Molecular Plant-Microbe Interactions. Eds. Bouarab et al. 2009, 1-35.
Navarro L, Jay F, Nomura K, He SY, Voinnet O: Suppression of the microRNA pathway by bacterial effector proteins. Science. 2008, 321 (5891): 964-967. 10.1126/science.1159505.
Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.
Jay F, Renou JP, Voinnet O, Navarro L: Biotic stress-associated microRNAs: identification, detection, regulation, and functional analysis. Methods Mol Biol. 2010, 592: 183-202. 10.1007/978-1-60327-005-2_13.
Ruiz-Ferrer V, Voinnet O: Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol. 2009, 60: 485-510. 10.1146/annurev.arplant.043008.092111.
Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8: 25-10.1186/1471-2229-8-25.
Voinnet O: Post-transcriptional RNA silencing in plant-microbe interactions: a touch of robustness and versatility. Curr Opin Plant Biol. 2008, 11 (4): 464-470. 10.1016/j.pbi.2008.04.006.
Chellappan P, Xia J, Zhou X, Gao S, Zhang X, Coutino G, Vazquez F, Zhang W, Jin H: siRNAs from miRNA sites mediate DNA methylation of target genes. Nucleic Acids Res. 2010, 38 (20): 6883-6894. 10.1093/nar/gkq590.
Chapman EJ, Carrington JC: Specialization and evolution of endogenous small RNA pathways. Nat Rev Genet. 2007, 8 (11): 884-896. 10.1038/nrg2179.
Moissiard G, Voinnet O: Viral suppression of RNA silencing in plants. Mol Plant Pathol. 2004, 5 (1): 71-82. 10.1111/j.1364-3703.2004.00207.x.
Brigneti G, Voinnet O, Li WX, Ji LH, Ding SW, Baulcombe DC: Viral pathogenicity determinants are suppressors of transgene silencing in Nicotiana benthamiana. EMBO J. 1998, 17 (22): 6739-6746. 10.1093/emboj/17.22.6739.
Voinnet O: Induction and suppression of RNA silencing: insights from viral infections. Nat Rev Genet. 2005, 6 (3): 206-220. 10.1038/nrg1555.
Lakatos L, Szittya G, Silhavy D, Burgyan J: Molecular mechanism of RNA silencing suppression mediated by p19 protein of tombusviruses. EMBO J. 2004, 23 (4): 876-884. 10.1038/sj.emboj.7600096.
Lakatos L, Csorba T, Pantaleo V, Chapman EJ, Carrington JC, Liu YP, Dolja VV, Calvino LF, Lopez-Moya JJ, Burgyan J: Small RNA binding is a common strategy to suppress RNA silencing by several viral suppressors. EMBO J. 2006, 25 (12): 2768-2780. 10.1038/sj.emboj.7601164.
Ding SW: RNA-based antiviral immunity. Nat Rev Immunol. 2010, 10 (9): 632-644. 10.1038/nri2824.
Vanitharani R, Chellappan P, Fauquet CM: Geminiviruses and RNA silencing. Trends Plant Sci. 2005, 10 (3): 144-151. 10.1016/j.tplants.2005.01.005.
Chellappan P, Masona MV, Vanitharani R, Taylor NJ, Fauquet CM: Broad spectrum resistance to ssDNA viruses associated with transgene-induced gene silencing in cassava. Plant Mol Biol. 2004, 56 (4): 601-611. 10.1007/s11103-004-0147-9.
Patil BL, Fauquet CM: Cassava mosaic geminiviruses: actual knowledge and perspectives. Mol Plant Pathol. 2009, 10 (5): 685-701. 10.1111/j.1364-3703.2009.00559.x.
Raja P, Sanville BC, Buchmann RC, Bisaro DM: Viral genome methylation as an epigenetic defense against geminiviruses. J Virol. 2008, 82 (18): 8997-9007. 10.1128/JVI.00719-08.
Fauquet CM, Bisaro DM, Briddon RW, Brown JK, Harrison BD, Rybicki EP, Stenger DC, Stanley J: Revision of taxonomic criteria for species demarcation in the family Geminiviridae, and an updated list of begomovirus species. Arch Virol. 2003, 148 (2): 405-421. 10.1007/s00705-002-0957-5.
Soitamo AJ, Jada B, Lehto K: HC-Pro silencing suppressor significantly alters the gene expression profile in tobacco leaves and flowers. BMC Plant Biol. 2011, 11: 68-10.1186/1471-2229-11-68.
Siddiqui SA, Sarmiento C, Truve E, Lehto H, Lehto K: Phenotypes and functional effects caused by various viral RNA silencing suppressors in transgenic Nicotiana benthamiana and N. tabacum. Mol Plant Microbe Interact. 2008, 21 (2): 178-187. 10.1094/MPMI-21-2-0178.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55 (4): 611-622. 10.1373/clinchem.2008.112797.
Memelink J: Regulation of gene expression by jasmonate hormones. Phytochemistry. 2009, 70 (13–14): 1560-1570.
Avanci NC, Luche DD, Goldman GH, Goldman MH: Jasmonates are phytohormones with multiple functions, including plant defense and reproduction. Genet Mol Res. 2010, 9 (1): 484-505. 10.4238/vol9-1gmr754.
Lin Z, Zhong S, Grierson D: Recent advances in ethylene research. J Exp Bot. 2009, 60 (12): 3311-3336. 10.1093/jxb/erp204.
Lorenzo O, Piqueras R, Sanchez-Serrano JJ, Solano R: ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell. 2003, 15 (1): 165-178. 10.1105/tpc.007468.
He XJ, Hsu YF, Zhu S, Wierzbicki AT, Pontes O, Pikaard CS, Liu HL, Wang CS, Jin H, Zhu JK: An effector of RNA-directed DNA methylation in arabidopsis is an ARGONAUTE 4- and RNA-binding protein. Cell. 2009, 137 (3): 498-508. 10.1016/j.cell.2009.04.028.
Takeda A, Iwasaki S, Watanabe T, Utsumi M, Watanabe Y: The mechanism selecting the guide strand from small RNA duplexes is different among argonaute proteins. Plant Cell Physiol. 2008, 49 (4): 493-500. 10.1093/pcp/pcn043.
Sun JQ, Jiang HL, Li CY: Systemin/Jasmonate-mediated systemic defense signaling in tomato. Mol Plant. 2011, 4 (4): 607-615. 10.1093/mp/ssr008.
Soosaar JL, Burch-Smith TM, Dinesh-Kumar SP: Mechanisms of plant resistance to viruses. Nat Rev Microbiol. 2005, 3 (10): 789-798. 10.1038/nrmicro1239.
Alvarez S, Zhu M, Chen S: Proteomics of Arabidopsis redox proteins in response to methyl jasmonate. J Proteomics. 2009, 73 (1): 30-40. 10.1016/j.jprot.2009.07.005.
Jung C, Lyou SH, Yeu S, Kim MA, Rhee S, Kim M, Lee JS, Choi YD, Cheong JJ: Microarray-based screening of jasmonate-responsive genes in Arabidopsis thaliana. Plant Cell Rep. 2007, 26 (7): 1053-1063. 10.1007/s00299-007-0311-1.
Maslenkova LT, Zanev Y, Popova LP: Oxygen-evolving activity of thylakoids from barley plants cultivated on different concentrations of jasmonic Acid. Plant Physiol. 1990, 93 (4): 1316-1320. 10.1104/pp.93.4.1316.
Baena-Gonzalez E, Rolland F, Thevelein JM, Sheen J: A central integrator of transcription networks in plant stress and energy signalling. Nature. 2007, 448 (7156): 938-942. 10.1038/nature06069.
Baena-Gonzalez E, Sheen J: Convergent energy and stress signaling. Trends Plant Sci. 2008, 13 (9): 474-482. 10.1016/j.tplants.2008.06.006.
Lozano-Duran R, Bejarano ER: Geminivirus C2 protein might be the key player for geminiviral co- option of SCF-mediated ubiquitination. Plant Signal Behav. 2011, 6 (7): 999-1001. 10.4161/psb.6.7.15499.
Lozano-Duran R, Rosas-Diaz T, Gusmaroli G, Luna AP, Taconnat L, Deng XW, Bejarano ER: Geminiviruses subvert ubiquitination by altering CSN-mediated derubylation of SCF E3 ligase complexes and inhibit jasmonate signaling in Arabidopsis thaliana. Plant Cell. 2011, 23 (3): 1014-1032. 10.1105/tpc.110.080267.
Yang JY, Iwasaki M, Machida C, Machida Y, Zhou X, Chua NH: Î2C1, the pathogenicity factor of TYLCCNV, interacts with AS1 to alter leaf development and suppress selective jasmonic acid responses. Genes Dev. 2008, 22 (18): 2564-2577. 10.1101/gad.1682208.
Lewsey MG, Murphy AM, Maclean D, Dalchau N, Westwood JH, Macaulay K, Bennett MH, Moulin M, Hanke DE, Powell G, Smith AG, Carr JP: Disruption of two defensive signaling pathways by a viral RNA silencing suppressor. Mol Plant Microbe Interact. 2010, 23 (7): 835-845. 10.1094/MPMI-23-7-0835.
Miao Y, Zentgraf U: The antagonist function of Arabidopsis WRKY53 and ESR/ESP in leaf senescence is modulated by the jasmonic and salicylic acid equilibrium. Plant Cell. 2007, 19 (3): 819-830. 10.1105/tpc.106.042705.
Ascencio-Ibanez JT, Sozzani R, Lee TJ, Chu TM, Wolfinger RD, Cella R, Hanley-Bowdoin L: Global analysis of Arabidopsis gene expression uncovers a complex array of changes impacting pathogen response and cell cycle during geminivirus infection. Plant Physiol. 2008, 148 (1): 436-454. 10.1104/pp.108.121038.
Nagar S, Hanley-Bowdoin L, Robertson D: Host DNA replication is induced by geminivirus infection of differentiated plant cells. Plant Cell. 2002, 14 (12): 2995-3007. 10.1105/tpc.005777.
Hanley-Bowdoin L, Settlage SB, Orozco BM, Nagar S, Robertson D: Geminiviruses: models for plant DNA replication, transcription, and cell cycle regulation. Crit Rev Biochem Mol Biol. 2000, 35 (2): 105-140.
Schommer C, Palatnik JF, Aggarwal P, ChÃÂ©telat A, Cubas P, Farmer EE, Nath U, Weigel D: Control of Jasmonate Biosynthesis and Senescence by miR319 Targets. PLoS Biol. 2008, 6 (9):
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Usadel B, Nagel A, Thimm O, Redestig H, Blaesing OE, Palacios-Rojas N, Selbig J, Hannemann J, Piques MC, Steinhauser D, Scheible WR, Gibon Y, Morcuende R, Weicht D, Meyer S, Stitt M: Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physiol. 2005, 138 (3): 1195-1204. 10.1104/pp.105.060459.
Usadel B, Nagel A, Steinhauser D, Gibon Y, Blasing OE, Redestig H, Sreenivasulu N, Krall L, Hannah MA, Poree F, Fernie AR, Stitt M: PageMan: an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments. BMC Bioinforma. 2006, 7: 535-10.1186/1471-2105-7-535.
Fryer MJ, Oxborough K, Mullineaux PM, Baker NR: Imaging of photo-oxidative stress responses in leaves. J Exp Bot. 2002, 53 (372): 1249-1254. 10.1093/jexbot/53.372.1249.
Porra RJ, Thompson WA, Kriedemann PE: Determination of accurate extinction coefficients and simultaneous equations for assaying chlorophyll a and b with four different solvents: verification of the concentration of chlorophyll by atomic absorption spectroscopy. Biochim Biophys Acta. 1989, 975: 384-394. 10.1016/S0005-2728(89)80347-0.
Neff MM, Chory J: Genetic interactions between phytochrome A, phytochrome B, and cryptochrome 1 during Arabidopsis development. Plant Physiol. 1998, 118 (1): 27-35. 10.1104/pp.118.1.27.
The Finnish Microarray and Sequencing Centre (FMSC) at Turku Centre for Biotecnology are acknowledged for labeling the cDNAs, hybridizations, scanning the chips and producing the raw microarray data. Dr. Mika Keränen is acknowledged for helping data analysis and Center of Scientific Calculating (CSC, Espoo, Finland) for Chipster program. Research was supported by The Academy of Finland, grant numbers 127203 and 128943.
The authors have no non-financial, financial or patent related competing interests.
AJS grew and collected the plant material for the experiments. The work was planned by KL. The experimental work and the reannotation of significantly altered transcripts were carried out by AJS and BJ together. AJS and KL wrote the manuscript. All authors have read and approved the manuscript.
Electronic supplementary material
Additional file 2: Quality control of intensity values in leaf and flower samples. Graphical Box Plot presentation of microarray data of intensity values originating from leaf and flower samples of wild type and transgenic AC2 expressing plants. (PDF 334 KB)
Over and under presentation analysis of functional categorization of leaves expressing AC2 or HC-Pro RSS.
Additional file 7: Analysis was performed using PageMan program (MapMan 3.5.1R2). Numbers after functional categorization indicate log2 values of differentially expressed genes and intensity of coloured boxes the z-scores of p values with FDR<0.05 . (PDF 796 KB)
Additional file 8: Visual presentation of transcripts involved in biotic stress. Data consists of up or down regulated transcripts with p-values less than 0.05 (FDR) in leaf and flower samples expressing AC2 or HC-Pro RSS. (PDF 575 KB)
Additional file 9: The amount of chlorophyll and anthocyanin in leaves of wild type and in AC2 expressing plants. (XLSX 16 KB)
Additional file 10: The amount of total protein against fresh weight (FW) and dry weight (DW) in wild type and AC2 expressing tobacco leaves. (XLSX 12 KB)
Additional file 12: Data consists of up or down regulated transcripts in leaf samples expressing AC2 or HC-Pro RSS. (PDF 427 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Soitamo, A.J., Jada, B. & Lehto, K. Expression of geminiviral AC2 RNA silencing suppressor changes sugar and jasmonate responsive gene expression in transgenic tobacco plants. BMC Plant Biol 12, 204 (2012). https://doi.org/10.1186/1471-2229-12-204
- Transgenic Tobacco Plant
- Transcriptional Gene Silence
- Express Plant
- Wild Type Tobacco
- Silence Suppression