Resistance To Wheat Streak Mosaic Virus Generated by Expression of An Artificial Polycistronic Microrna in Wheat
Resistance To Wheat Streak Mosaic Virus Generated by Expression of An Artificial Polycistronic Microrna in Wheat
Resistance To Wheat Streak Mosaic Virus Generated by Expression of An Artificial Polycistronic Microrna in Wheat
ª 2011 CSIRO
150 Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd
Virus resistance via multiple artificial miRNAs 151
utility of other approaches is worth exploring. One such Petrie, P. Shrestha, M. Fahim, Q. Liu, C. C. Wood and S. P.
approach, amiRNA-mediated gene silencing, has recently been Singh, unpublished).
developed specifically to address the risk of off-target effects Here, we expressed five pre-amiRNA, potentially generating
and transgene–virus recombination to form new biotypes ten amiRNA species, to different conserved regions of the
(Schwab et al., 2006). WSMV genome from a modified version of rice miRNA precur-
The amiRNA approach utilizes a naturally occurring miRNA sor miR395. The polycistronic amiRNA strategy is able to pro-
precursor as a backbone, with the mature miRNA sequence duce marker-free transgenic wheat plants immune to WSMV,
being replaced to gain new targeting ability (Ossowski et al., demonstrating that this is a viable strategy for a major crop spe-
2008; Vaucheret et al., 2004). In plants, amiRNAs have been cies. Moreover, it alleviates the concerns of recombinant novel
successfully used to down-regulate endogenous genes (Alvarez viral entities forming and also produces plants predicted to
et al., 2006; Khraiwesh et al., 2008; Molnar et al., 2009; avoid the loss of resistance caused by virus mutation. The reso-
Schwab et al., 2006; Warthmann et al., 2008) and also for lution of technological problems and concerns implies this strat-
developing transgenic virus resistance against Turnip yellow egy has a strong biotechnological potential for agriculture.
mosaic virus (TYMV), Turnip mosaic virus (TuMV) (Niu et al.,
2006), Cucumber mosaic virus (CMV) (Duan et al., 2008), Results
Potato virus X (PVX), and Potato virus Y (PVY) (Ai et al., 2011)
Design of polycistronic amiRNA construct
in Arabidopsis; against CMV (Qu et al., 2007) in tobacco;
against CMV (Zhang et al., 2011) in tomato; and against We chose pre-miR395 as the backbone for simultaneous
Cassava brown streak virus (CBSV) and Cassava brown streak expression of multiple amiRNAs targeting various conserved
Uganda virus (CBSUV) in cassava (Wagaba et al., 2010). regions in the WSMV genome. We combined published ami-
It has been argued that the use of short viral sequences in RNA selection criterion (http://wmd3.weigelworld.org/cgi-bin/
this amiRNA approach is less likely to enable the emergence of webapp.cgi) into a software application we call ‘miR Mate’.
novel viral entities through recombination and trans-encapsida- Applying this to the available full WSMV genome sequences
tion (Schnippenkoetter et al., 2001). However, when only a (five at the time) identified approximately 120 potential target
small viral sequence is used, the virus is more likely to evolve in sites for amiRNA. As a final selection filter, we searched Wheat
the amiRNA target sequence via transition mutation and enable TIGR mRNA databases for potential off-targets of these poten-
avoidance of amiRNA complementarity and defence (Simon- tial WSMV amiRNAs using the online miRU program (http://
Mateo and Garcia, 2006; Lin et al., 2009). Other examples plantgrn.noble.org/psRNATarget/). No expressed sequences
include HIV escape mutants to avoid RNAi (Boden et al., 2003; were found as potential targets when three mismatches were
Das et al., 2004; Westerhout et al., 2005). To substantially allowed. Through this process, five amiRNA were chosen and
reduce the risk of viruses evolving to avoid degradation, a strat- designated amiRNA-1, amiRNA-2, amiRNA-3, amiRNA-4, and
egy would be very useful where an amiRNA precursor gene amiRNA-5; their WSMV genome targets and their target coordi-
expressed multiple amiRNAs targeting different conserved struc- nates are given in Table 1 and Figure 1. The target of amiRNA-
tural and functional portions of the viral genome. Resistance to 1 lies in 5¢ UTR region, amiRNA-2 targets the newly described
this protection would require simultaneous mutations to avoid open reading frame (ORF) pipo region of P3 cistron (Chung
all the amiRNA sequences. A similar rationale was invoked in et al., 2008), amiRNA-3 targets P1 gene, amiRNA-4 targets P3
the work of Israsena et al. (2009) where they designed precur- cistron (upstream of pipo), and amiRNA-5 targets the HCpro
sor genes encoding three amiRNA against rabies virus and gene on WSMV genome (Figure 1). We deliberately chose a
tested them in cell culture. Likewise, multiple siRNA were devel- mix of targets on the genomic and replicative strands of the
oped from a multiplex miRNA directed against HIV in cultured virus in case one strand was more available for the amiRNA sur-
cells (Liu et al., 2008; ter Brake et al., 2006). veillance than the other. While further good targets could be
miRNA precursors that have been used for the delivery of identified in the 3¢ genes, the bias to the 5¢ genes is simply a
amiRNAs in plants include miR159a (Niu et al., 2006); miR171a result of beginning the screening for potential off-targets from
(Qu et al., 2007), miR172a (Schwab et al., 2006), miR30 (Zeng that end. Surprisingly, no 21-nt sequence in the extreme 3¢
et al., 2002), miR528 (Warthmann et al., 2008), and miR167b region of the virus could be identified, which was both con-
(Ai et al., 2011). In plants, Niu et al. (2006) demonstrated that served and met the design rules. The endogenous miRNAs and
a dimeric amiRNA precursor in Arabidopsis could be effective miRNA* that are derived from the miR395 precursor were
against two different viruses. Others have successfully targeted replaced with these amiRNA and amiRNA* sequences to con-
two endogenous transcripts with dimeric amiRNA precursors in serve the secondary structure of the transcript. The predicted
Arabidopsis (Park et al., 2009) and Chlamydomonas (Zhao secondary structure of the polycistronic miR395 and the modi-
et al., 2009) using different miRNA precursors. Here, in our fied artificial miR395 were almost identical (Figure 2), presum-
studies, we used the multiplex precursor of rice miR395 family ably enhancing the prospect of the predicted biogenesis of
of miRNAs that was identified in both Arabidopsis thaliana and mature amiRNAs. This artificial polycistronic precursor was
Oryza sativa computationally and was later experimentally veri- named FanGuard395 (FGmiR395).
fied (Guddeti et al., 2005; Kawashima et al., 2009). OsmiR395
Generation of wheat carrying the FanGuard395
targets ATP sulphurylases that are involved in sulphate assimila-
transgene
tion (Rotte and Leustek, 2000) and is induced in sulphur starva-
tion to regulate a low-affinity sulphate transporter and two ATP FGmiR395 was synthesized by Geneart and cloned into pWubi
sulphurylases (Allen et al., 2005; Jones-Rhoades and Bartel, vector, to generate FGWS-pWubi, where FGmiR395 was behind
2004). The rice miR395 is a single 1-kb transcript that gener- a constitutive maize polyubiquitin promoter separated from the
ates a convoluted RNA structure that generates seven fully pro- transgene by a spliceable ubiquitin intron. FGWS-pWubi was
cessed miRNA (Jones-Rhoades and Bartel, 2004; S. Belide, J. R. cobombarded into wheat immature embryos along with
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Table 1 Conservation of amiRNA targets in WSMV Genome. The alignment used to design amiRNA against conserved targets in WSMV genome. AlignX was used with default settings using Vec-
tor NTI 10. (a) Alignment of the five chosen target regions in the five published WSMV genomes at the time of the design of the amiRNA. (b) Alignment in the five target regions in the ten new
WSMV genomes that became available subsequent to amiRNA design. (c) Alignment of the amiRNA target regions in the WSMV-ACT isolate (unpublished and obtained subsequent to the design
152 Muhammad Fahim et al.
of FGmiR395 and production of the transgenic plants). The mismatched nucleotides with other isolates are highlighted. amiRNA-1 and amiRNA-2 target the replicating strand, signified by number-
ing 1–21 from left to right; amiRNA-3, amiRNA-4 and amiRNA-5 target the genomic strand of the WSMV, signified by numbering 1–21 from right to left. The chosen targets remain absolutely
conserved in all five target regions of the WSMV-ACT isolate, when virus population RNA was resequenced from inoculated plants in the study including )S, +R, and the transgenic breakdown
plants +S
amiRNA-1 Target sequence amiRNA-2 Target sequence amiRNA-3 Target sequence amiRNA-4 Target sequence amiRNA-5 Target sequence
AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
Accession Description Origin 1 21 1 21 21 1 21 1 21 1
(a)
AF057533 Sydney 81 Nebraska AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
1AF285169 Type Strain Kansas AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF285170 El Batan Mexico AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATACCACTATGTACCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF454454 Czech Czech AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTTCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF454455 Turkish Turkey AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
(b)
AF511614 H95S Kansas AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGGTTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF511615 H98 Kansas AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CTGCGAACGTCTTGCAAGTTA
AF511618 ID96 Idaho AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF511619 ID99 Idaho AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
AF511630 Mon96 Montana AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
EU914917 Naghadeh Iran AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTTCATTATGTACCAA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
EU914918 Sadat-Saher Iran AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
F511643 WA99 WA, USA AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTTCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
FJ348358 WA94 WA, USA AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTTCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCAAACGTCTTGCAAGTTA
FJ348359 ARG2 Argentina AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
(c)
Unpublished ACT Australia AGCTCTCGCATAGAGATAAGC TCGAGCAAGATCTTTCACACG GAAGATTCCATTATGTGCCGA CCAGGAAGCATTTTCTGGTCA CCGCGAACGTCTTGCAAGTTA
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 153
Figure 1 Structure of the Wheat streak mosaic virus (WSMV) genome (approximately 9400 nt), the target sites for amiRNAs and the FGmiR395
transgene. (a) Genome map of WSMV showing the five conserved regions (indicated by scissors) targeted by amiRNAs, amiRNA-1 to amiRNA-5. (b)
Design of FGmiR395 construct (1400 nt) used to transform wheat using biolistics; shown are the probe region for Southern blot and primer sequences
FgPf1 and M13RevP used in PCR. (c) Diagram of pCMneoSTLS2 containing the nptII gene for geneticin resistance, used in the cotransformation of
immature wheat embryos and showing the position of the PCR primers pNeo3 and pNeo5.
(a) (b)
1
1–5. These are the predicted fold structures of
transcripts using RNAfold. Bars showing regions
1 2
corresponding to amiRNA guide strand. Second-
0
0
ary structure probabilities are indicated by heat
map (Blue, weak; Red, strong).
plasmid pCMneoSTLS2 (Maas et al., 1997) that contained neo- within the promoter region, full ubiquitin intron and all of the
mycin phosphotransferase gene (nptII), conferring resistance to FGmiR395 transgene including the nos terminator (Figure 1b).
geneticin (G418). From a wheat transformation of 379 PCR analysis of T1 families confirmed that 14 of the 23 families
embryos, a total of 23 T0 transgenic wheat plants were gener- had FGmiR395 as well as nptII (Figure 3a, PCR data not shown)
ated from 16 different embryos; therefore, there were at least and 10 of these were from different embryos and were there-
16 independent events. fore independent transgenic events. Southern blot analysis of T2
The transgenic lines were designated FanGuard plasmid (FGP) plants from a subset of these lines confirmed the presence of
plus a number corresponding to the bombarded embryo. transgene (Figure 4). Nine families were negative for FGmiR395
Where multiple T0 plants were obtained from a single embryo, and positive for nptII. The latter lines were discarded after the
they were distinguished with lower-case letters, e.g. FGP1a and preliminary assessment for resistance (next section).
FGP1b. All T0, T1, and T2 transgenic plants were morphologi-
Preliminary assessment of FGmiR395 transgenic wheat
cally indistinguishable from the wild-type parental cultivar, Bob
in T1
White selection 26 (BW26), implying that the FGmiR395 does
not influence growth or development (Figure 3b). The T1 generation was subsequently challenged with WSMV
Plants recovered from the transformation selection cultures through mechanical inoculation at the three-leaf stage (4–17
were screened through PCR for the nptII selectable marker from plants per family) using the spray gun. Wheat streak mosaic dis-
the pCMneoSTLS2 cotransformation vector (Figure 1c). This ease is characterized by light-green-to-faint-yellow streaks in
confirmed that all 23 plants coming through the antibiotic wheat leaves parallel to the veins. The virus arrests growth, and
selection were true transgenics carrying the selectable marker. plants show moderate-to-severe stunting with prostrated tillers
Genomic PCR screening for the presence of FGmiR395 was car- often with empty spikes or spikes with shrivelled kernels.
ried out using primers, fgPf1 and M13RevP, that span from Serological characterization of the transgenic families involved
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
154 Muhammad Fahim et al.
Figure 4 Southern blot analysis of families in T2 generation. (a) BamHI digested wheat DNA for transgene copy number. (b) HindIII-digested wheat
DNA for transgene copy size. FC is full copy insert size 3 KB; TC is truncated copy inserts. Plant FGP18.6 is labelled )S by PCR and bioassay; it has FG
sequence present, but that sequence is truncated so that it is not detected by the PCR.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 155
phenotypes to be studied in greater numbers and detail in T1 inoculated transgenic individuals and ELISA readings very simi-
and T2 generations. lar to the ELISA readings of uninoculated controls, suggesting
the complete absence of virus from the inoculated transgenic
Analysis of transgenic family FGP6 for resistance
plants. Experiments were conducted to see whether infectious
to WSMV
virus or viral RNA could be recovered from the resistant inocu-
The preliminary analysis indicated that all FGmiR395-positive lated transgenic plants. Leaf sap from inoculated plants in
segregants in FGP1b, FGP4b, FGP5a, and FGP6 were resistant three transgenic families (FGP8c, FGP13b, and 15a) was
(Figure 3a). Subsequently, a bioassay on 26 T1 individuals of extracted and inoculated onto test plants of control BW26 at
FGP6 showed that all 13 FGmiR395-carrying segregants were 1 ⁄ 10 (w ⁄ v) dilution to investigate the presence of any infec-
symptom free at all four time points (7, 14, 21, and 28 d p.i.) tious WSMV particles. Results from these test inoculation
and had background ELISA ratio at 14 and 28 d p.i. One of the experiments revealed that all FGmiR395+R phenotypes were
major effects of virus infection on plant physiology is the severe immune to WSMV, as no infectious virus could be recovered
stunting and extreme reduction in plant height. When we plot- and carried over to control wheat through mechanical inocula-
ted the virus concentration (ELISA ratio of virus inoculated tion with the most concentrated leaf extract inoculum. Sap
plant) in segregating FGP6 progeny, we found that they clearly from inoculated segregants with FGmiR395 transgenes failed
grouped into two clusters (Figure 5b), showing that the plant to transmit infection to susceptible BW26 as judged by symp-
height is inversely related to the virus concentration. toms and ELISA, whereas sap from segregants with no trans-
PCR analysis of selectable marker nptII in this family revealed gene and non-transformed controls (BW26) did transmit
that 17 of 26 plants carried the selectable marker (Table 2). One infection in every case. Examples of these tests are shown in
immune FGmiR395 transgenic plant FGP6.22 was negative for Figure 6.
selectable marker nptII, while five FGmiR395-negative segregants
The formation of amiRNAs from the FGmiR395
were carrying nptII, and eight segregants were negative for both
transcript
transgenes. It is worth noting that the approach and technique
utilized in this work can yield marker-free immune transgenics. The expression of the amiRNAs was analysed in virus-free
Southern blot analysis carried out on T2 segregants confirmed transgenic wheat leaves and detected by splinted ligation using
stable integration of FGmiR395 transgene(s) into the wheat miRtect IT (Maroney et al., 2007). Potentially, from the five
genome (Table 2) and apparently multiple copies of the trans- duplex arms of the precursor FGmiR395, five guide strands
gene (BamHI digest: Figure 4a). (amiRNA) could be produced; moreover, if one also considers
the loading of passenger strand (amiRNA*) into RNA Induced
Efficacy of viral suppression in transgenic families
Silencing Complex (RISC), then a total of ten amiRNA could
The analyses reported thus far for the resistant transgenic potentially be produced against the virus. To detect both
families showed the complete absence of symptoms in amiRNA and amiRNA* sequences using splinted ligation, we
Figure 5 Segregation of resistance in FGP families expressing FGmiR395 transgene. Virus levels were detected by double-antibody sandwich enzyme-
linked immunosorbent assay at 28 d p.i. postinoculation and expressed as ratio of inoculated and healthy control. The plant height (in cm) was
measured at the heading stage and plotted against the corresponding ELISA ratio. )S indicated the FGmiR395 transgene–negative susceptible segre-
gant, +S indicated transgene-carrying susceptible segregant, +MR indicated segregants that carried the transgene but accumulated virus titre owing to
resistance breakdown, while +R indicated the transgene-carrying resistant segregants. Transgenic families are shown as (a) FGP4a, (b) FGP6, (c) FGP8c,
(d) FGP15a, and (e) FGP18. MR, moderate resistance.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
156 Muhammad Fahim et al.
Table 2 Segregation of transgene and resistance in T1 and T2. Analysis of six selected FGmiR395 families in T1 and T2. Recorded for T1 are
the results of Southern hybridisation, segregation for genomic PCR, and resistance. Segregation for resistance in T2 is also shown for a series
of selected derivative families. Selected T2 families were derived as shown from T1 individuals that were transgene carrying and resistant (+R);
transgene-carrying susceptible (+S); and transgene-negative susceptible ()S). Southern blots are shown in Figure 4
T0 (NT) T1 T2
FGP4a 2 1 21 : 11 24 : 8 15 : 9 4a.18 +S 20 : 15
4a.22 +R 29 : 0
4a.31 )S 0:8
FGP6 4 – 17 : 9 13 :173 13 : 13 6.3 )S 0 : 14
6.10 +R 22 : 15
FGP8c 3 1 49 : 21 51 : 19 30 : 21 8c.10 +S 14 : 21
8c.11 +R 24 : 8
8c.27 )S 0 : 19
FGP15a 1 1 39 : 13 39 : 13 32 : 7 15a.1 +R 34 : 1
15a.2 +S 24 : 4
15a.11 )S 0:6
FGP18 3 – 26 : 12 17 : 21 0 : 38 18.6 )S 0 : 15
18.10 +S 12 : 32
designed the ten bridging oligos accordingly (Table S1). All five
designed anti-WSMV miRNAs, amiRNA-1, amiRNA-2, amiRNA-
3, amiRNA-4, and amiRNA5, were expressed in the plant. Using
an immune plant in T2 generation, seven of the potential ten
amiRNA were generated from FGmiR395 and accumulated suf-
ficiently to be readily detected (Figure 7). At least these seven
amiRNAs would be loaded into RISC and would be expected to
target the invading virus at both the genomic RNA and replica-
tive RNA levels in four viral regions, the 5¢ UTR, P1, P3, and
HCpro regions.
Investigating the resistance breakdown phenotype
In addition to stable immunity FGmiR395+R, other phenotypes
were observed in some transgenic events described previously
as FGmiR395+MR and FGmiR395+S. These also were studied
in more detail in the five families chosen for closer examina-
tion. Negative segregants FGmiR395)S along with wild-type
BW26 were used as negative (susceptible) controls in all
assays.
Transgenic event FGP4a
In the T1 generation of FGP4a, 24 of 32 plants carried the
transgene; however, only 15 of these showed strong resistance
throughout the experiment. The resistance in nine individuals
Figure 6 Virus transmission from the inoculated transgenic plants onto
broke down with time, and symptoms began to appear as faint
Bobwhite26. Sap was extracted from inoculated transgenic plants at
streaks at 14 d p.i. Two of these nine individuals showed a
1 ⁄ 10 (w ⁄ v). Each sap extract was inoculated onto three wild-type BW26
plants. At 28 d p.i., samples were collected and processed for Wheat
moderate resistance (MR) as judged by no symptoms at 14 d
streak mosaic virus ELISA. (a) Plotted is the average ELISA ratio (inocu- p.i., but the gradual development of leaf symptoms, growth to
lated divided by healthy) for the three test-inoculated plants. (b) Also only intermediate height, and accumulated virus titre at 28 d
tabulated are the molecular analyses of the inoculated T1 individuals p.i.; these segregants were characterized as FGmiR395+MR
whose sap extracts were being tested. Shown for each sap donor is the (Figure 5a). PCR, ELISA ratio, and plant height plot suggested
FGmiR395 PCR result, the bioassay result for that donor plant (as 28 d four groups in this family. The third cluster (FGmiR395+S) accu-
p.i. ELISA ratio), and the symptom score at 28 d p.i. (0–4 scale). Control mulated less virus than the null segregants but was equally
BW(i) is inoculated BW26. BW(h) is the uninoculated BW26. affected in height.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 157
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
158 Muhammad Fahim et al.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 159
Resistance cosegregated with the transgene in most of the online (NCBI) had grown from five to 13. An alignment of all 13
transgenic events of this study. However, analysis of T1 families revealed that three of five chosen amiRNA targets are still
revealed a range of phenotypes. (i) Immune individuals were completely conserved in the WSMV genome (Table 1). This
obtained, which remained symptomless and with no viral coat highlights the importance of having multiple targets in polycis-
protein accumulation all the way through to maturity. Sap pre- tronic amiRNA and the importance of aligning as many virus ge-
pared from these inoculated plants failed to transmit infection nomes as possible to select highly conserved regions.
to susceptible controls. Such immunity was evident across two One of the targets in FGmiR395 was the functional region
generations, (ii) In other cases, the plants were resistant without called pipo (pretty interesting Potyviridae ORF) within the gene
symptoms or ELISA detected coat protein at 14 d p.i., but sub- encoding the P3 protein. Pipo was initially identified as a small
sequently, the resistance broke down, allowing virus to accumu- ORF embedded in the P3 cistron of Turnip mosaic virus (TuMV;
late by 28 d p.i., (iii) A third transgenic phenotype may be genus Potyvirus; family Potyviridae) (Chung et al., 2008), and its
called plant recovery, in which early susceptibility is followed by presence was confirmed in 48 viruses representing all genera in
full recovery and resistance, (iv) The forth phenotype is where the family Potyviridae, including WSMV. Mutation in pipo hin-
the presence of multiple copies of the transgene confers only ders various important functions; in the case of WSMV that
moderate resistance. We challenged the transgenic families with includes effects on replication and movement in the plant.
high-titre virus inoculum at two time points 0 and 10 d p.i. and When mutations are introduced into the pipo region of P3,
therefore suggest that the breakdown of resistance in some without affecting the amino acid sequence of the translated
events might result from an excessive virus pressure that over- protein, the virus loses the ability to replicate in protoplasts
whelms the amiRNA-mediated resistance. Such a phenomenon (Chung et al., 2008) or it is restricted to only a few cells in inoc-
has been observed in transgenic barley against Barley yellow ulated plants (Wen and Hajimorad, 2010).
dwarf virus (BYDV) where an increase in viruliferous aphid infes- As one would expect, the presence of one full copy of the
tation resulted in breakdown of RNAi resistance to BYDV (M.-B. transgene can be enough to confer resistance (immunity)
Wang, pers. commun.). Ai et al. (2011) showed that although against WSMV. However, we saw evidence in some transgenic
amiRNAs were detectable and resistance to PVX or PVY evident events of additional truncated FGmiR395 insertions and
in their transgenic plants, the resistance was overcome by rein- behaviour in T1 and T2 generations suggestive that the
oculation at 35 d p.i., resulting in increased viral pressure. The truncated copy may be interfering with the expression of the
pressure of inoculation in our experiments was very high and full-length copy; subsequent loss of the truncated copy restores
may be responsible for some of the more complex resistance effective resistance. This interpretation will require further
phenotypes. It would be interesting to investigate the efficiency experimentation to confirm. We sometimes observed resistance
of the FGmiR395-expressing wheat plants in the field under the segregation ratios inconsistent with Mendelian expectation
milder pressures expected from natural wheat curl mite infesta- based on the inferred insert locus number. Usually, this involved
tions. lower than expected numbers of resistant individuals and is
It could also be informative to quantify the level of amiRNA likely associated with some inserts of the transgene being inef-
expression in the various transgenic phenotypes to observe fective or conditionally effective (Matzke et al., 2009).
whether there is a correlation with the degree or stability of Previously, we have reported the use of long hairpin dsRNA-
resistance. Such an analysis would need to follow all ten of the mediated WSMV immunity in wheat (Fahim et al., 2010). The
potential amiRNA species and would ideally include the various comparison to the present study is a rare opportunity to con-
phases of resistance breakdown and various phases of recovery. trast the two approaches in the same genetic background
Previous studies have shown some degree of correlation against the same virus. Compared to amiRNA, long hairpin RNA
between amiRNA expression and virus resistance (PVX and PVY had the advantage of a very high frequency of insertion events
in Arabidopsis, Ai et al., 2011; CMV in tobacco, Qu et al., with stable heritable immunity. In this respect, long hairpin
2007). When there are at least five and potentially ten species RNAi is very attractive from a biotechnology application per-
of amiRNA attacking the virus, as with FGmi395, the analysis of spective. However, long hpRNA would have a higher probability
correlation with resistance will be complex. of unintended silencing of off-target genes in the host (Xu
One potential risk with amiRNA-mediated resistance is the et al., 2006; Duan et al., 2008; Khraiwesh et al., 2008). Fur-
generation of virus mutants that escape the amiRNA surveillance thermore, long hpRNA approaches are seen by some as posing
(Simon-Mateo and Garcia, 2006; Lin et al., 2009). We addressed a risk in the field of heterologous recombination with other
this issue in our studies by selecting amiRNA targets based upon virus genomes and resulting in new virus biotypes. Low temper-
conserved regions in five full WSMV genome sequences avail- atures can also compromise the efficacy of RNAi silencing strat-
able. We observed no evidence that the virus was evolving egies (Szittya et al., 2003). On the other hand, miRNAs appear
(mutating) during the course of the experiment. This was evi- to be completely temperature independent, and the transgenic
dent because virus populations collected from FGmiR395+S T1 lines expressing virus-derived amiRNA retain their resistance at
plants was used to inoculate T2 progeny. Not only did the resis- low temperatures (Niu et al., 2006; Szittya et al., 2003). We
tant transgenic continue to be resistant to these inoculations, were able to achieve immunity in wheat to WSMV from both
but progeny from some +S T1 parents included +R and +MR the long hairpin dsRNA and the amiRNA strategies. We con-
phenotypes. We take this as evidence that the original +S sus- clude that amiRNA-based viral resistance, especially polycistronic
ceptibility was not the result of virus mutation; otherwise, resis- amiRNA as advocated here, deserves and needs further in-depth
tance would not have emerged in T2 when challenged with the studies to improve the amiRNA efficiency. An even better com-
putatively mutant virus preparation. Furthermore, deep sequenc- parison would be achieved if the long hpRNA was designed to
ing of virus populations from infected transgenic plants con- cover the same regions as the amiRNA used in this study.
firmed no mutations in the target sequences. By the end of this The work described here exemplifies the utility of miR395
study, the number of full-length WSMV genomes available and similar miRNA clusters as a carrier of multiple amiRNAs.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
160 Muhammad Fahim et al.
They can be used to target multiple regions of the one virus (as The stemloop backbone
here), multiple viruses, or multiple endogenous mRNA species
For the delivery of the final five amiRNAs as a polycistronic
(S. Belide, J. R. Petrie, P. Shrestha, M. Fahim, Q. Liu, C. C.
transgene, we selected the precursor of rice miR395 that is
Wood and S. P. Singh, unpublished). Mixed viral infections are
expressed under sulphur stress conditions (Guddeti et al., 2005;
common in the field especially in fruits and vegetables. Using
Jones-Rhoades and Bartel, 2004; Kawashima et al., 2009).
the polycistronic amiRNA, it will be possible to target highly
A synthetic gene called FanGuard (FG) was designed by replac-
conserved regions of multiple viruses. Similarly, polycistronic
ing the five native miRNA sequences in the first five duplex
pre-amiRNA genes will be effective in targeting multiple plant
arms of native miR395 with five amiRNA designed to target
endogenous genes for functional genomics and in applications
WSMV. In a parallel study, a similar construct was used to
such as redirecting plant metabolism into novel products.
simultaneously silence five endogenous genes in Arabidopsis
We conclude that polycistronic amiRNAs can be utilized to
(Belide et al. submitted). The designed FGmiR395 was synthe-
induce virus resistance in commercially valuable plants, where
sized through GENEART GmbH (http://www.geneart.com)
there are limited options of natural resistance. We anticipate
flanked by restriction sites for BamHI and KpnI in the carrier
ongoing improvements in the understanding of miRNA biogene-
plasmid. The FG gene was excised from the carrier plasmid
sis and design of amiRNA to further enhance the utility for virus
using appropriate restriction enzymes and ligated between the
resistance and engineering other agronomically important traits.
Ubiquitin promoter and tm1 terminator of vector pWubi-tm1
Furthermore, the expression of multiple amiRNAs from a single
vector (Wang and Waterhouse, 2000) generating cereal trans-
precursor transgene will minimize the difficulties of repeated
formation plasmid FG-pWubi.
transformations, need for multiple selectable markers, and the
constraint of breeding with multiple independent loci. Wheat transformation
Transgenic wheat plants were generated following microparticle
Experimental procedures bombardment of 186 immature cv. Bob White 26 (BW26)
wheat embryos. The embryos were cobombarded with two
Designing WSMV-specific amiRNA
plasmids, FG-pWubi and a selectable marker plasmid
To select conserved regions in WSMV genome as targets for arti- pCMneoSTLS2, as described previously (Fahim et al., 2010;
ficial miRNAs, full-genome sequences of WSMV were retrieved Pellegrineschi et al., 2002).
from NCBI (Table 1). The sequences were aligned with Clustal
Analysis of T0 transgenic plants—PCR
W ⁄ AlignX (a component of Vector NTI Advance 10.3.0) to
screen for highly conserved regions in the viral genome. The DNA extraction was carried out using DNAeasy Plant Mini Kit
possible amiRNA sequences were generated from the highly following manufacturer’s instructions (Qiagen Inc, Valencia,
conserved regions (20 nt or more in length) using the basic crite- CA). For PCR-based genotyping of the FGmiR395 transgene,
ria defined at WMD3 (http://wmd3.weigelworld.org/, a web DNA was amplified (Figure 1b) using FgPf1 5¢-TGCAGCATC-
microRNA design tool) and incorporated into a software algo- TATTCATATGC-3¢ and M13RevP 5¢-CATGGTCATAGCTGTT-3¢,
rism called miR Mate developed specifically for this study. The that generated approximately 1.4 kb of FG-pWubi amplicon
algorithm was developed using Microsoft. NET Framework and covering promoter, transgene, and terminator regions, under
also incorporates the Vienna RNA Package 1.7 algorithm RNA- the following thermocycler conditions 94 C for 30 s, 60 C for
fold.exe (Hofacker et al., 1994; McCaskill, 1990; Zuker and Stie- 45 s, 72 C for 60 s for 35 cycles with a final extension at
gler, 1981). The miRNA design criteria used include A ⁄ U at 72 C for 10 min. For the selectable marker nptII (Figure 1c), a
position 1 (Mi et al., 2008; Eamens et al., 2009; Takeda et al., 700-bp nptII fragment was amplified using the forward primer
2008), A at position 10 (Reynolds et al., 2004; Mallory et al., Neo3 5¢-TACGGTATCGCCGCTCCCGAT-3¢ and reverse primer
2004), and G ⁄ C at position 21 (P. Waterhouse, pers. commun.). Neo5 5¢-GGCTATTCGGCTATGACTG-3¢, both sequences being
miR Mate utilizes the RNAfold algorithm to calculate minimal in the nptII coding region, using the following thermocycler
free energy (mfe) values for the formation of the candidate miR- conditions: 94 C for 30 s, 55 C for 30 s, 72 C for 60 s for
NA’s folded structure; values of £)30 kcal ⁄ mol represent opti- 40 cycles with a final extension at 72 C for 10 min.
mal stability. The negative values reflect the fact that stored
Analysis of T1 and T2 transgenic plants—Virus bioassay
energy is released during the formation of the structure;
the more negative the value, the more energy is released and Virus inoculum was prepared by grinding WSMV-infected tissue
the more favourable is the formation of the structure. Candidate in a mortar and pestle at a 1 : 10 w ⁄ v ratio in 0.02 M potas-
amiRNAs with the lowest mfe value (the highest stability) were sium phosphate buffer (pH 7). After filtering through four lay-
then assessed for potential off-targets in wheat and barley. ers of Miracloth (Calbiochem, La Jolla, CA), abrasive celite
The set of potential virus target sequences were used to (Johns-Manville, Denver, CO) was added at 2% w ⁄ v to serve
search for genes that may be potential off-targets in wheat or as an abrasive. For the analysis in T1 plants, the inoculum was
barley, using miRU: Plant microRNA Potential Target Finder prepared from virus-infected non-transgenic BW26; for the
http://bioinfo3.noble.org/miRNA/miRU.htm (a recent version analysis of T2 plants, sap from a mixture of FGmiR395 carrying
psRNATarget: A Plant Small RNA Regulator Target Analysis Ser- susceptible and FGmiR395-negative segregant susceptible
ver is available at http://bioinfo3.noble.org/psRNATarget/) (Bren- plants was used. The sap–celite mixture was first applied with
necke et al., 2005; Jones-Rhoades and Bartel, 2004; Lim et al., an air-powered spray gun, and then, leaves were gently rubbed
2005; Mallory et al., 2004). WSMV-derived amiRNAs were with fingers to ensure the infection of plant by the virus. At
selected having the least probability of targeting any sequence 10 d p.i., plants were reinoculated with the virus-infected sap
in the gene or EST databases of wheat or barley. to ensure high inoculum pressure.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 161
The plants were scored for symptoms at 7, 14, 21, and 28 d 15–50 bp was dissected and small RNA was extracted from
p.i. on a scale of 0–4 with 0 as healthy, 1 as mild with very excised gel overnight in 4 M NaCl. The RNA concentration was
few streaks, 2 as moderate with streaks that coalesce, 3 as measured in 1 lL of solution using Nanodrop (Thermo Scien-
severe with approximately 50% leaf area with streaks, and 4 as tific, Wilmington). The splinted ligation was performed on the
the most severe or lethal symptoms where the streaks develop purified fraction with miRtect-IT miRNA Labeling and Detec-
into chlorosis of more than 70% of leaf area based on visual tion Kit (USB, Cleveland, OH) (Maroney et al., 2007). Specific
observation. WSMV-specific ELISA was performed on leaf sam- bridge oligonucleotides (Table S1) were designed according to
ples collected at 14 d p.i. and 28 d p.i., using Agdia reagents the manufacturer’s directions. Using 50 ng of enriched smRNA
(Agdia, Elkhart, IN) following the manufacturer’s instructions. per reaction, amiRNAs were captured by a specific bridge oligo-
Plates were read at A405nm in ELISA Reader Spectra Max 340 nucleotide and ligated to the P32-labelled detection oligonucleo-
PC (Molecular Devices, Sunnyvale, CA) 60 min after the addi- tide with T4 DNA ligase. Ligated products were separated on
tion of substrates. Healthy controls were included on every 17% urea–polyacrylamide gel and visualized using Fujifilm FLA-
plate, every sample was duplicated, and duplicate value means 5000 phosphor imager.
were used in calculating the ELISA value ratio between inocu-
Small RNA library preparation and deep sequencing
lated and healthy controls.
Small RNAs were enriched using the mirVana miRNA Isolation
Segregation analysis of FanGuard transgene and
Kit (Invitrogen) following manufacturer instructions. Small RNA-
resistance in T1 and T2 generations
Seq libraries were prepared based on Illumina’s alternative v1.5
For detailed segregation analysis of selected events in T1 gen- protocol and a published method (Lu et al., 2007) and run on
eration, 25–35 seeds were germinated in pots. Leaf samples the Illumina’s GAIIx platform at the Genome Discovery Unit of
were collected, and DNA was extracted as described previ- Australian National University. WSMV sequences were assem-
ously, and genomic PCR was conducted to detect both bled with assistance from Dr Stephen Ohms, JCSMR, The
FGmiR395 and nptII amplicons. The cosegregation of resistance Australian National University.
with the transgene was assessed by challenging with WSMV
as described earlier. ELISA was performed 14 and 28 d p.i. in
Acknowledgements
T1 and only 28 d p.i. for T2 generations on inoculated plants.
Plant heights and symptoms were recorded at 7, 14, 21, and We are thankful to Dr Peter Waterhouse and Dr Ming-Bo Wang
28 d p.i. for their technical advice at various stages of these studies,
Mr. Zarman Mazhar Rizvi for his assistance with putting
Test inoculation to detect infectious virus in leaf sap
together various algorithms under one miR Mate, Drs. James
Sap was extracted at 28 d p.i. from each inoculated transgenic Petrie and Qing Liu for technical supervision in amiRNA detec-
plant to be tested, using 0.02 M potassium phosphate buffer at tion, and Dr Jun Fan, The John Curtin School of Medical
1 : 10 (1 g leaf per 10 mL buffer) concentration and mixed with Research, The Australian National University, for providing
celite abrasive and then inoculated onto three control BW26 excellent training and technical supervision in small RNA library
plants. This test inoculation (or back-inoculation) method was preparation. Dr Stephen Ohms is gratefully acknowledged for
used to evaluate the effectiveness of the FanGuard (FGmiR395) his assistance with bioinformatic analysis of small RNAs. The
transgene in eliminating viral replication and preventing the for- first author thankfully acknowledges AusAID for financial
mation of infectious particles. Symptoms were scored on the assistance as PhD studentship.
test-inoculated plants and leaf samples collected 28 d p.i. for
ELISA as described previously. References
Analysis of transgenic plants—Southern hybridization
Abbott, D., Wang, M.-B. and Waterhouse, P.M. (2002) A single copy of
Southern hybridization was carried out as described previously virus-derived transgene-encoding hairpin RNA confers BYDV immunity. In
(Fahim et al., 2010; Lagudah et al., 1991). Instead of T0 plants, Barley Yellow Dwarf Disease: Recent Advances and Future Strategies
a pool of 8–44 T1 individuals per family were used for the anal- (Henry, M. and McNab, A., eds), pp. 22–26. Mexico, DF: CIMMYT.
Abel, P.P., Nelson, R.S., De, B., Hoffmann, N., Rogers, S.G., Fraley, R.T.
ysis of transgene copy number and copy size. This method was
and Beachy, R.N. (1986) Delay of disease development in transgenic plants
used so as not to compromise the initial transgenic T0 plant yet
that express the tobacco mosaic virus coat protein gene. Science, 232, 738–
to capture all the insertion events likely to have been present in 743.
the T0 plant. DNA was digested with BamHI to determine the Ai, T., Zhang, L., Gao, Z., Zhu, C.X. and Guo, X. (2011) Highly efficient virus
number of independent insertions; there is only one site for resistance mediated by artificial microRNAs that target the suppressor of
BamH1 in FGmi395 (Figures 1b and 4). DNA was digested with PVX and PVY in plants. Plant Biol. (Stuttg), 13, 304–316.
HindIII to assess whether the inserted copies were full length or Allen, E., Xie, Z.X., Gustafson, A.M. and Carrington, J.C. (2005) microRNA-
truncated. directed phasing during trans-acting siRNA biogenesis in plants. Cell, 121,
207–221.
amiRNA analysis in immune transgenic plants Alvarez, J.P., Pekker, I., Goldshmidt, A., Blum, E., Amsellem, Z. and Eshed, Y.
(2006) Endogenous and synthetic microRNAs stimulate simultaneous,
The small RNA fraction was enriched from a fraction of total
efficient, and localized regulation of multiple targets in diverse species.
RNA extracted from 100 mg of transgene-carrying immune
Plant Cell, 18, 1134–1151.
FGP4a.22 T2 plants, using miRvana Kit (Ambion, Austin, TX) Baulcombe, D. (2004) RNA silencing in plants. Nature, 431, 356–363.
following manufacturer instructions. The extracted total RNA Boden, D., Pusch, O., Lee, F., Tucker, L. and Ramratnam, B. (2003) Human
was run on denaturing 17% PAGE and stained with EtBr. Using immunodeficiency virus type 1 escape from RNA interference. J. Virol. 77,
100-bp RNA ladder as reference, the region corresponding to 11531–11535.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
162 Muhammad Fahim et al.
ter Brake, O., Konstantinova, P., Ceylan, M. and Berkhout, B. (2006) Lu, C., Meyers, B.C. and Green, P.J. (2007) Construction of small RNA cDNA
Silencing of HIV-1 with RNA Interference: A Multiple shRNA Approach. libraries for deep sequencing. Methods, 43, 110–117.
Molecular Therapy, 14, 883–892. Maas, C., Simpson, C.G., Eckes, P., Schickler, H., Brown, J.W.S., Reiss, B.,
Brennecke, J., Stark, A., Russell, R.B. and Cohen, S.M. (2005) Principles of Salchert, K., Chet, I., Schell, J. and Reichel, C. (1997) Expression of intro
microRNA-target recognition. PLoS Biol. 3, 404–418. modified nptII genes in monocotylenonous and dicotyledonous plants cells.
Chung, B.Y.W., Miller, W.A., Atkins, J.F. and Firth, A.E. (2008) An Mol. Breed. 3, 15–28.
overlapping essential gene in the Potyviridae. Proc. Natl Acad. Sci. USA, Mallory, A.C., Reinhart, B.J., Jones-Rhoades, M.W., Tang, G.L., Zamore, P.D.,
105, 5897–5902. Barton, M.K. and Bartel, D.P. (2004) MicroRNA control of PHABULOSA in
Das, A.T., Brummelkamp, T.R., Westerhout, E.M., Vink, M., Madiredjo, M., leaf development: importance of pairing to the microRNA 5¢region. EMBO
Bernards, R. and Berkhout, B. (2004) Human immunodeficiency virus type 1 J. 23, 3356–3364.
escapes from RNA interference-mediated inhibition. J. Virol. 78, 2601–2605. Maroney, P.A., Chamnongpol, S., Souret, F. and Nilsen, T.W. (2007) A rapid,
Duan, C.G., Wang, C.H., Fang, R.X. and Guo, H.S. (2008) Artificial quantitative assay for direct detection of microRNAs and other small RNAs
microRNAs highly accessible to targets confer efficient virus resistance in using splinted ligation. RNA, 13, 930–936.
plants. J. Virol. 82, 11084–11095. Matzke, M., Kanno, T., Daxinger, L., Huettel, B. and Matzke, A.J. (2009)
Dwyer, G.I., Gibbs, M.J., Gibbs, A.J. and Jones, R.A.C. (2007) Wheat streak RNA-mediated chromatin-based silencing in plants. Curr. Opin. Cell Biol.
mosaic virus in Australia: relationship to isolates from the Pacific Northwest 21, 367–376.
of the USA and its dispersion via seed transmission. Plant Dis. 91, 164– McCaskill, J.S. (1990) The equilibrium partition-function and base pair
170. binding probabilities for RNA secondary structure. Biopolymers, 29, 1105–
Eamens, A.L., Smith, N.A., Curtin, S.J., Wang, M.B. and Waterhouse, P.M. 1119.
(2009) The Arabidopsis thaliana double-stranded RNA binding protein Mi, S., Cai, T., Hu, Y., Chen, Y., Hodges, E., Ni, F., Wu, L., Li, S., Zhou, H.,
DRB1 directs guide strand selection from microRNA duplexes. RNA, 15, Long, C., Chen, S., Hannon, G.J. and Qi, Y. (2008) Sorting of small RNAs
2219–2235. into Arabidopsis argonaute complexes is directed by the 5¢ terminal
Fahim, M., Ayala-Navarrete, L., Millar, A.A. and Larkin, P.J. (2010) Hairpin nucleotide. Cell, 133, 116–127.
RNA derived from viral NIa gene confers immunity to Wheat streak Molnar, A., Bassett, A., Thuenemann, E., Schwach, F., Karkare, S., Ossowski,
mosaic virus infection in transgenic wheat plants. Plant Biotechnol. J. 8, S., Weigel, D. and Baulcombe, D. (2009) Highly specific gene silencing by
821–834. artificial microRNAs in the unicellular alga Chlamydomonas reinhardtii.
Guddeti, S., Zhang, D.C., Li, A.L., Leseberg, C.H., Kang, H., Li, X.G., Zhai, Plant J. 58, 165–174.
W.X., Johns, M.A. and Mao, L. (2005) Molecular evolution of the rice Niu, Q.W., Lin, S.S., Reyes, J.L., Chen, K.C., Wu, H.W., Yeh, S.D. and Chua,
miR395 gene family. Cell Res. 15, 631–638. N.H. (2006) Expression of artificial microRNAs in transgenic Arabidopsis
Harvey, T.L. and Seifers, D.L. (1991) Transmission of wheat streak mosaic thaliana confers virus resistance. Nat. Biotechnol. 24, 1420–1428.
virus to sorghum by the wheat curl mite (Acari, Eriophyidae). J. Kans. Ossowski, S., Schwab, R. and Weigel, D. (2008) Gene silencing in plants
Entomol. Soc. 64, 18–22. using artificial microRNAs and other small RNAs. Plant J. 53, 674–690.
Hofacker, I.L., Fontana, W., Stadler, P.F., Bonhoeffer, L.S., Tacker, M. and Park, W., Zhai, J.X. and Lee, J.Y. (2009) Highly efficient gene silencing using
Schuster, P. (1994) Fast folding and comparison of RNA secondary perfect complementary artificial miRNA targeting AP1 or heteromeric
structures. Monatshefte Für Chemie, 125, 167–188. artificial miRNA targeting AP1 and CAL genes. Plant Cell Rep. 28, 469–
Israsena, N., Supavonwong, P., Ratanasetyuth, N., Khawplod, P. and 480.
Hemachudha, T. (2009) Inhibition of rabies virus replication by multiple Pellegrineschi, A., Noguera, L.M., Skovmand, B., Brito, R.M., Velazquez, L.,
artificial microRNAs. Antiviral Res. 84, 76–83. Salgado, M.M., Hernandez, R., Warburton, M. and Hoisington, D. (2002)
Jackson, A.L., Bartz, S.R., Schelter, J., Kobayashi, S.V., Burchard, J., Mao, M., Identification of highly transformable wheat genotypes for mass production
Li, B., Cavet, G. and Linsley, P.S. (2003) Expression profiling reveals of fertile transgenic plants. Genome, 45, 421–430.
off-target gene regulation by RNAi. Nat. Biotechnol. 21, 635–637. Qu, J., Ye, J. and Fang, R.X. (2007) Artificial microRNA-mediated virus
Jones, R.A.C., Coutts, B.A., Mackie, A.E. and Dwyer, G.I. (2005) Seed resistance in plants. J. Virol. 81, 6690–6699.
transmission of Wheat streak mosaic virus shown unequivocally in wheat. Reynolds, A., Leake, D., Boese, Q., Scaringe, S., Marshall, W.S. and
Plant Dis. 89, 1048–1050. Khvorova, A. (2004) Rational siRNA design for RNA interference. Nature
Jones-Rhoades, M.W. and Bartel, D.P. (2004) Computational identification of Biotechnology, 22, 326–330.
plant MicroRNAs and their targets, including a stress-induced miRNA. Mol. Rotte, C. and Leustek, T. (2000) Differential subcellular localization and
Cell, 14, 787–799. expression of ATP sulfurylase and 5¢-adenylylsulfate reductase during
Kawashima, C.G., Yoshimoto, N., Maruyama-Nakashita, A., Tsuchiya, Y.N., ontogenesis of arabidopsis leaves indicates that cytosolic and plastid forms
Saito, K., Takahashi, H. and Dalmay, T. (2009) Sulphur starvation induces of ATP sulfurylase may have specialized functions. Plant Physiol. 124, 715–
the expression of microRNA-395 and one of its target genes but in 724.
different cell types. Plant J. 57, 313–321. Schnippenkoetter, W.H., Martin, D.P., Willment, J.A. and Rybicki, E.P. (2001)
Khraiwesh, B., Ossowski, S., Weigel, D., Reski, R. and Frank, W. (2008) Forced recombination between distinct strains of Maize streak virus.
Specific gene silencing by artificial MicroRNAs in Physcomitrella patens: an J. Gen. Virol. 82, 3081–3090.
alternative to targeted gene knockouts. Plant Physiol. 148, 684–693. Schwab, R., Ossowski, S., Riester, M., Warthmann, N. and Weigel, D. (2006)
Lagudah, E.S., Appels, R. and Mcneil, D. (1991) The nor-D3 locus of Triticum Highly specific gene silencing by artificial microRNAs in Arabidopsis. Plant
tauschii – natural variation and genetic-linkage to markers in chromosome- Cell, 18, 1121–1133.
5. Genome, 34, 387–395. Seifers, D.L., Harvey, T.L., Martin, T.J. and Jensen, S.G. (1998) Partial host
Lim, L.P., Lau, N.C., Garrett-Engele, P., Grimson, A., Schelter, J.M., Castle, J., range of the high plains virus of corn and wheat. Plant Dis. 82, 875–879.
Bartel, D.P., Linsley, P.S. and Johnson, J.M. (2005) Microarray analysis Simon-Mateo, C. and Garcia, J.A. (2006) MicroRNA-Guided processing
shows that some microRNAs downregulate large numbers of target impairs Plum pox virus replication, but the virus readily evolves to escape
mRNAs. Nature, 433, 769–773. this silencing mechanism. J. Virol. 80, 2429–2436.
Lin, S.S., Wu, H.W., Elena, S.F., Chen, K.C., Niu, Q.W., Yeh, S.D., Chen, C.C. Sivamani, E., Brey, C.W., Dyer, W.E., Talbert, L.E. and Qu, R.D. (2000)
and Chua, N.H. (2009) Molecular evolution of a viral non-coding sequence Resistance to Wheat streak mosaic virus in transgenic wheat expressing the
under the selective pressure of amiRNA-mediated silencing. Plos Pathog. 5, viral replicase (NIb) gene. Mol. Breed. 6, 469–477.
e1000312. Sivamani, E., Brey, C.W., Talbert, L.E., Young, M.A., Dyer, W.E., Kaniewski,
Liu, Y.P., Haasnoot, J., ter Brake, O., Berkhout, B. and Konstantinova, P. W.K. and Qu, R.D. (2002) Resistance to Wheat streak mosaic virus in
(2008) Inhibition of HIV-1 by multiple siRNAs expressed from a single transgenic wheat engineered with the viral coat protein gene. Transgenic
microRNA polycistron. Nucleic Acids Res. 36, 2811–2824. Res. 11, 31–41.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163
Virus resistance via multiple artificial miRNAs 163
Slykhuis, J.T. (1955) Aceria tulipae Keifer (Acarina, Eriophyidae) in relation to Wilson, J.A. and Richardson, C.D. (2005) Hepatitis C Virus replicons escape
the spread of wheat streak mosaic. Phytopathology, 45, 116–128. RNA interference induced by a short interfering RNA directed against the
Smith, N.A., Singh, S.P., Wang, M.B., Stoutjesdijk, P.A., Green, A.G. and NS5b coding region. J. Virol. 79, 7050–7058.
Waterhouse, P.M. (2000) Gene expression – total silencing by intron- Xu, P., Zhang, Y., Kang, L., Roossinck, M.J. and Mysore, K.S. (2006)
spliced hairpin RNAs. Nature, 407, 319–320. Computational estimation and experimental verification of off-target
Szittya, G., Silhavy, D., Molnar, A., Havelda, Z., Lovas, A., Lakatos, L., silencing during posttranscriptional gene silencing in plants. Plant Physiol.
Banfalvi, Z. and Burgyan, J. (2003) Low temperature inhibits RNA silencing- 142, 429–440.
mediated defence by the control of siRNA generation. EMBO J. 22, 633– Zeng, Y., Wagner, E.J. and Cullen, B.R. (2002) Both natural and designed
640. micro RNAs technique can inhibit the expression of cognate mRNAs when
Takeda, A., Iwasaki, S., Watanabe, T., Utsumi, M. and Watanabe, Y. (2008) expressed in human cells. Mol. Cell, 9, 1327–1333.
The mechanism selecting the guide strand from small RNA duplexes is Zhang, X., Li, H., Zhang, J., Zhang, C., Gong, P., Ziaf, K., Xiao, F. and Ye, Z.
different among argonaute proteins. Plant Cell Physiol. 49, 493–500. (2011). Expression of artificial microRNAs in tomato confers efficient and
Vaucheret, H., Vazquez, F., Crete, P. and Bartel, D.P. (2004) The action of stable virus resistance in a cell-autonomous manner. Transgenic Res. 20,
ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA 569–581.
pathway are crucial for plant development. Genes Dev. 18, 1187–1197. Zhao, T., Wang, W., Bai, X. and Qi, Y. (2009) Gene silencing by artificial
Wagaba, H., Basavaprabhu, P.L., Jitender, Y.S., Nigel, T.J., Alicai, T., Baguma, microRNAs in Chlamydomonas. Plant J. 58, 157–164.
Y., Settumba Mukasa, B. and Fauquet, C.M. (2010). Testing the efficacy of Zuker, M. and Stiegler, P. (1981) Optimal computer folding of large RNA
artificial microRNAs to control cassava brown streak disease. Second sequences using thermodynamics and auxiliary information. Nucleic Acids
RUFORUM Biennial Meeting 20–24 September 2010, Entebbe, Uganda, Res. 9, 133–148.
287–291.
Wang, M.B. and Waterhouse, P.M. (2000) High-efficiency silencing of a beta-
glucuronidase gene in rice is correlated with repetitive transgene structure Supporting information
but is independent of DNA methylation. Plant Mol. Biol. 43, 67–82.
Additional Supporting information may be found in the online
Warthmann, N., Chen, H., Ossowski, S., Weigel, D. and Herve, P. (2008) Highly
version of this article:
specific gene silencing by artificial miRNAs in rice. PLoS ONE, 3, e1829.
Waterhouse, P.M., Graham, H.W. and Wang, M.B. (1998) Virus resistance Table S1 Oligo sequences used in miRTect analysis of amiRNA
and gene silencing in plants can be induced by simultaneous expression of accumulation in immune transgenics.
sense and antisense RNA. Proc. Natl Acad. Sci. USA, 95, 13959–13964.
Table S2 Families used in segregation analysis in T2.
Wen, R.H. and Hajimorad, M.R. (2010) Mutational analysis of the putative
pipo of soybean mosaic virus suggests disruption of PIPO protein impedes Please note: Wiley-Blackwell are not responsible for the content
movement. Virology, 400, 1–7. or functionality of any supporting materials supplied by the
Westerhout, E.M., ter Brake, O., Ooms, M., Vink, M., Das, A.T. and authors. Any queries (other than missing material) should be
Berkhout, B. (2005) HIV-1 can evade RNAi-mediated inhibition by altering
directed to the corresponding author for the article.
the secondary structure of its RNA genome. J. Biotechnol. 118, S69–S69.
ª 2011 CSIRO
Plant Biotechnology Journal ª 2011 Society for Experimental Biology, Association of Applied Biologists and Blackwell Publishing Ltd, Plant Biotechnology Journal, 10, 150–163