Host-induced silencing of essential genes in Puccinia triticina through transgenic expression of RNAi sequences reduces severity of leaf rust infection in wheat Vinay Panwar1,2, Mark Jordan1, Brent McCallum1 and Guus Bakkeren2,* 1Agriculture and Agri-Food Canada, Morden Research and Development Centre, Morden, MB, Canada 2Agriculture and Agri-Food Canada, Summerland Research and Development Centre, Summerland, BC, Canada Received 15 February 2017; revised 18 August 2017; accepted 20 September 2017. *Correspondence (Tel 250-494-6368; fax 250-494-0755; email Guus.Bakkeren@agr.gc.ca) Keywords: Puccinia triticina, wheat leaf rust, host-induced gene silencing, RNAi, transgenic wheat, rust resistance. Summary Leaf rust, caused by the pathogenic fungus Puccinia triticina (Pt), is one of the most serious biotic threats to sustainable wheat production worldwide. This obligate biotrophic pathogen is prevalent worldwide and is known for rapid adaptive evolution to overcome resistant wheat varieties. Novel disease control approaches are therefore required to minimize the yield losses caused by Pt. Having shown previously the potential of host-delivered RNA interference (HD- RNAi) in functional screening of Pt genes involved in pathogenesis, we here evaluated the use of this technology in transgenic wheat plants as a method to achieve protection against wheat leaf rust (WLR) infection. Stable expression of hairpin RNAi constructs with sequence homology to Pt MAP-kinase (PtMAPK1) or a cyclophilin (PtCYC1) encoding gene in susceptible wheat plants showed efficient silencing of the corresponding genes in the interacting fungus resulting in disease resistance throughout the T2 generation. Inhibition of Pt proliferation in transgenic lines by in planta-induced RNAi was associated with significant reduction in target fungal transcript abundance and reduced fungal biomass accumulation in highly resistant plants. Disease protection was correlated with the presence of siRNA molecules specific to targeted fungal genes in the transgenic lines harbouring the complementary HD-RNAi construct. This work demonstrates that generating transgenic wheat plants expressing RNAi-inducing transgenes to silence essential genes in rust fungi can provide effective disease resistance, thus opening an alternative way for developing rust-resistant crops. Introduction Rust fungi, belonging to the genus Puccinia, are among the most economically destructive biotrophic phytopathogenic fungi infect- ing cereal crops. Members of this genus pose a serious threat to global wheat production which is a major staple food for mankind in many parts of the world. Wheat is a host to three important rust pathogen species: Puccinia graminis f. sp. tritici (Pgt), Puccinia striiformis f. sp. tritici (Pst) and Puccinia triticina (Pt), which are the causal agents of stem rust, stripe rust and leaf rust, respectively (Bolton et al., 2008; Chen et al., 2014; Singh et al., 2011). Of the three rust diseases, wheat leaf rust (WLR, caused by Pt) is comparatively more prevalent occurring regularly wherever wheat is grown (Kolmer, 2005). Because of its adaptation to diverse climatic conditions and widespread occur- rence, leaf rust results in greater total annual losses worldwide than stem and stripe rusts (Huerta-Espino et al., 2011). Under severe epidemic conditions, Pt can inflict yield losses ranging from 10% to 70% (Herrera-Foessel et al., 2006). Fungicides are usually applied for management of rust diseases, but financial costs and adverse environmental impacts are associated with chemical inputs. In addition, a major concern with excessive application of antimicrobial compounds is a high risk of reduced sensitivity or development of resistant pathogen populations. Genetic resistance remains the most economical, effective and ecologically sustainable approach to minimize the extent of crop damage from leaf rust disease (Draz et al., 2015). Current disease control measures against Pt rely predominantly on breeding and development of resistant cultivars. So far, more than 70 leaf rust resistance (R) genes have been reported, the majority of which confers protection to a particular genetic variant of pathogen species mediated by gene-for-gene-type interactions (Ellis et al., 2014; Flor, 1971; Singla et al., 2017). However, resistance alleles deployed into elite cultivars followed by monoculture farming instigate enormous selection pressures on rust pathogens. This often leads to new virulent races of the fungus, evolved to avoid recognition by these alleles and thereby rendering the resistant variety ineffective in a relatively short time (McCallum et al., 2016; Webb and Fellers, 2006). Further, finding and introgressing durable resistance genes from wild donor species without any undesirable trait drag remain a challenge to breeders. It is expected that rapid evolutionary dynamics of rust pathogens combined with anthro- pogenic climate changes will accelerate disease outbreaks in cultivatedwheat worldwide (Helfer, 2014). Therefore, a continuous search for developing alternative biotechnological solutions is required toalleviate thenegative impact of rustpathogens inwheat. Pt is a highly specialized macrocyclic, heteroecious fungus having a complex life cycle involving up to five different spore stages on two taxonomically different host plants with species in Triticeae being the primary hosts (Bolton et al., 2008). The fungus reproduces on wheat by the formation of asexual dikaryotic urediniospores, and under favourable conditions, this phase can lead to multiple infection cycles thereby causing crop failure within weeks. Pt infection in wheat occurs via urediniospore germ ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. 1013 Plant Biotechnology Journal (2018) 16, pp. 1013–1023 doi: 10.1111/pbi.12845 http://creativecommons.org/licenses/by/4.0/ tubes that differentiate to form appressoria over leaf stomata to penetrate into substomatal cavities. Subsequently, the plant tissue is colonized with an intercellular mycelium, breaching plant cell walls and invaginating the plasma membrane to form haustoria. This intimate connection surrounded by a zone of extrahaustorial matrix marks the basis of molecular crosstalk between the host and fungal cells (Garnica et al., 2014). Through haustoria, the fungus acquires nutrients from the host and there is evidence to suggest that these structures also play a crucial role in the delivery of virulence effectors to suppress host defence responses (Petre et al., 2015; Struck, 2015). The asexual cycle is completed within 7–10 days when sporogenous uredinia are formed containing thousands of infective propagules that are easily wind-dispersed over long distances to establish new infections. The availability of many cereal rust fungus genome sequences, along with the integration of comparative analyses, has con- tributed valuable insight into genome organization and annota- tion of genes involved in pathogenesis (Cantu et al., 2011; Cuomo et al., 2017; Duplessis et al., 2011; Wu et al., 2017; Zheng et al., 2013). However, attempts to translate this knowl- edge to decipher pathogen biology are hampered because these biotrophic fungi are not amenable to many of the traditional molecular genetic approaches, mainly due to difficulties in culturing these species outside of their host. In the postgenomics era, RNA interference (RNAi) has emerged as a powerful technique for de novo discovery of gene function. This technique, termed post-transcriptional gene silencing (PTGS) in plants, involves the endonucleolytic cleavage of double-stranded RNA (dsRNA) into small, ~20–23 nucleotide (nt)-long, small interfering RNA (siRNA) duplexes by the Dicer-like (DCL) enzymes which along with their associated proteins control the expression of genetic information (Baulcombe, 2004; Wilson and Doudna, 2013). An RNAi-based approach referred to as host-induced gene silencing (HIGS) or host-delivered RNAi (HD-RNAi) in which small RNAs are produced by the host plant to target parasite transcripts has provided a promising strategy for improving plant resistance against pathogens, including fungi and oomycetes (Andrade et al., 2015; Govindarajulu et al., 2015; Huang et al., 2006; Jahan et al., 2015; Koch and Kogel, 2014; Koch et al., 2013; Nowara et al., 2010; Panwar et al., 2013a, 2016; Song and Thomma, 2016; Yin et al., 2010). When targeted against crucial pathogenicity genes, HD-RNAi can potentially be developed as a genetic method to curb pathogen virulence for pesticide-free disease control in crop plants. Recently, functional genomics assays based on HIGS method- ology were developed for rust fungi, facilitating the character- ization of some rust fungi genes required for virulence (Panwar et al., 2013a,b; Yin et al., 2014; Zhang et al., 2012). These studies harnessed transient gene silencing approaches, either delivering the HD-RNAi constructs through Barley stripe mosaic virus (BSMV) or Agrobacterium tumefaciens-mediated transfer. Applying these methods, we reported that Pt Map-kinase (PtMAPK1), cyclophilin (PtCYC1) and calcineurin B (PtCNB1) encoding genes were required for full pathogenicity on the wheat host (Panwar et al., 2013a,b). Here, we test the efficacy of stably integrated HD-RNAi constructs in transgenic wheat plants to achieve resistance against WLR disease by silencing the PtMAPK1 or PtCYC1 genes in Pt. We show that transgenic expression of HIGS cassettes targeting Pt genes provides a viable means to reduce WLR infection in otherwise highly susceptible wheat plants, providing an attractive and precise genetic engineering strategy for combating aggressive rust pathogens. Results Generation and analysis of transgenic wheat lines for Pt resistance The HD-RNAi constructs hp-PtMAPK1RNAi and hp-PtCYC1RNAi designed to express the 520-nt and 501-nt long 30-end coding segments of the fungal PtMAPK1 and PtCYC1 genes, respec- tively, as hairpin loop dsRNA (Experimental procedures) were used for wheat transformation. These constructs were intro- duced into scutellar tissue of immature embryos of Pt susceptible wheat cultivar Fielder by microprojectile bombardment. The HD- RNAi cassettes were highly specific to the target fungal genes with no discernible homology to off-target sequences in avail- able wheat genomic resources (Panwar et al., 2013a). Fourteen to twenty independent transgenic T0 wheat lines harbouring the hp-PtMAPK1RNAi and hp-PtCYC1RNAi construct, respectively, were generated and selfed to obtain the T1 seeds. Integration of the PtMAPK1 or PtCYC1 gene segment in T1 lines was verified by polymerase chain reaction (PCR), and expression of the dsRNA was confirmed by reverse-transcription PCR (RT-PCR), using Pt gene-specific primers (Figure S1, Table S1). No visible physiological or growth alterations were observed in the transgenic wheat lines expressing the respective HD-RNAi construct when compared with the nontransformed control plants of the same genotype, confirming no unintended RNAi effect in the host. To examine whether the presence of hp-PtMAPK1RNAi or hp- PtCYC1RNAi constructs in wheat lines affected the susceptibility to infection by Pt, 2-week-old transgenic T1 plants along with nontransformed control plants were challenged with fungal urediniospores and visually observed for WLR development. Analysis of disease phenotype identified six hp-PtMAPK1RNAi and seven hp-PtCYC1RNAi lines that were variably effective in restricting the magnitude of WLR symptoms as compared to controls. Segregation analysis revealed that the fungal transgenes segregated in a Mendelian ratio in these lines, indicating that integration of functional inserts occurred at one (3:1 ratio) or two (15:1 ratio) genetic loci (Table S2). The disease response of Pt- inoculated wheat plants in terms of infection type (IT) was scored at 12 days postinfection (dpi) using the 0–4 scale (Dakouri et al., 2013) to reflect uredinia size and abundance. Resistance to Pt in T1 plants expressing the respective HD-RNAi constructs was characterized by the appearance of fewer and smaller uredinia with repressed sporulation, corresponding to low infection types varying from 1 to 3. In contrast, Pt-challenged control wheat plants consistently exhibited vigorously sporulating uredinia resembling high ITs, indicative of highly susceptible reactions (Figure 1). Interestingly, the suppression of infection by Pt in transgenic wheat plants with low ITs (IT < 2) was also notable for delayed appearance of symptoms as compared to controls. WLR symptoms in control plants became apparent at 5–6 dpi which eventually progressed to large and more distinct uredinia, whereas in transgenic T1 plants displaying a low IT, appearance of disease symptoms was delayed by 1–3 days, compared to controls, with markedly reduced uredinia size and density (Figure S2). Taken together, a quantitative reduction of WLR disease symptoms was observed in transgenic wheat lines carrying a HIGS construct targeting the fungal PtCYC1 or PtMAPK1 gene. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Vinay Panwar et al.1014 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense Transgenic wheat lines expressing the HD-RNAi constructs impede Pt development We next explored whether suppression of Pt development in transgenic wheat T1 lines was associated with impairment of fungal growth inside leaf tissue by measuring the fungal biomass accumulation in leaves challengedwith Pt. Determination of fungal biomass showed that the ratio of fungal-to-plant DNA was clearly, but varyingly, reduced in transgenic lines, compared to nontrans- formed controls (Table 1). The hp-PtMAPK1RNAi expressing plants of lines MAPK1-2163, MAPK1-2169 and MAPK1-2166 were the most effective in restricting fungal development, showing biomass reductions by up to 77% in comparison with controls (Figures 1a and 2a). Similarly, fungal biomass was reduced by as much as 79% in hp-PtCYC1RNAi expressing plants of lines CYC1-2236, CYC1- 2248 and CYC1-2224 as compared to controls (Figures 1b and 2a). This marked reduction in fungal load in T1 lines suggested that the presence of the HD-RNAi constructs in transgenic wheat plants interfered with Pt virulence and thereby attenuating the capability of the fungus to colonize the host tissue. Generally, resistance to rust pathogens is not an apparent case of immunity versus susceptibility but, more often, is perceived as a quantitative difference in fungal growth between more or less resistant and susceptible plants which makes visual disease scoring systems sporadically prone to variability and error. To circumvent this possibility, we grouped transgenic T1 HD-RNAi plants into four different resistance types based on the quantifi- cation of fungal biomass in Pt-infected tissue by qPCR. The rationale behind such classification was to ensure accurate and unbiased disease resistance pattern detection in transgenic lines. Transgenic plants were considered as highly resistant (HR) to Pt if the HIGS-derived reduction in fungal biomass was more than 75% relative to control plants, reduction in fungal biomass by 50%–75% was designated as a resistant (R) phenotype, 25%– 50% reduction in fungal biomass was classified as moderately resistant (MR), whereas plants in which fungal biomass accumu- lation was suppressed by <25% as compared to controls were categorized as susceptible (S; Table 1). Overall, we observed a close relationship between WLR disease severity and the amount of fungal biomass detected in Pt-challenged wheat plants. For example, wheat plants with visibly high disease symptoms (high ITs) also showed higher fungal biomass in infected leaves than those suppressing fungal development (low ITs) that had less fungal biomass (Figure 2a). HD-RNAi results in reduction of targeted gene transcript levels in Pt infecting transgenic wheat lines In order to substantiate the hypothesis that protection to WLR infection in hp-PtMAPK1RNAi or hp-PtCYC1RNAi lines was Figure 1 Leaf rust disease resistance in transgenic wheat T1 lines. Pt- inoculated leaves of representative transgenic wheat T1 lines expressing the hp-PtMAPK1RNAi (a) or hp-PtCYC1RNAi (b) construct show suppressed disease development with low infection types (IT), compared to nontransformed control (NTC) plants. %, represents per cent fungal biomass reduction by HD-RNAi in a specific transgenic plant selected from the lines shown compared to control. Photographs were taken at 12 dpi. Table 1 Response of transgenic wheat T1 lines to Pt infection HD-RNAi construct and line number Disease intensity based on reduction in fungal biomass n S MR R HR hp-PtMAPK1RNAi MAPK1-2159 15 8 5 2 – MAPK1-2162 12 7 4 1 – MAPK1-2163 12 8 2 1 1 MAPK1-2166 23 9 9 5 – MAPK1-2169 18 9 4 4 1 MAPK1-2170 16 7 6 3 – Fielder-NTC 15 15 – – – hp-PtCYC1RNAi CYC1-2216 14 6 5 3 – CYC1-2224 15 10 2 2 1 CYC1-2226 6 4 2 – – CYC1-2227 11 3 5 3 – CYC1-2228 14 9 2 3 – CYC1-2236 15 9 3 2 1 CYC1-2248 13 6 4 3 – Fielder-NTC 15 15 – – – n, number of plants tested (PCR positive for transgenic lines); NTC, nontrans- genic control. HR, highly resistant, fungal biomass reduced by more than 75%; R, resistant, fungal biomass reduction by 50%–75%; MR, moderately resistant, reduction in fungal biomass by 25%–50%; S, susceptible, reduction in fungal biomass by <25% as compared to control plants. Ratio of fungal-to-wheat DNA analysed by qPCR using total genomic DNA extracted from Pt-inoculated leaves at 12 dpi. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Leaf rust resistance in transgenic wheat by HD-RNAi 1015 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense caused by in planta-derived silencing of the target PtMAPK1 or PtCYC1 genes, we quantified the transcript levels of the respec- tive genes in Pt infecting transgenic T1 and nontransformed control wheat plants at 5 dpi. Transcript levels of the PtMAPK1 and PtCYC1 genes were reduced by approximately 41%–65% in Pt invading transgenic wheat lines expressing the complementary HD-RNAi constructs compared to Pt-infected controls (Figure 3). This indicated that integration of HD-RNAi expression cassettes in transgenic wheat was sufficient to interfere with transcription of the targeted fungal genes during infection. The endogenous silencing of the target fungal genes in transgenic wheat plants was specific to the presence of the corresponding HD-RNAi construct because no effect was observed on the expression of the PtCYC1 gene in the hp-PtMAPK1RNAi lines and vice versa (Figure 3). These findings support the conclusion that PtMAPK1 and PtCYC1 contribute to Pt virulence and knockdown of their expression by HIGS compromised the fungus ability to proliferate in transgenic wheat. Next, we asked whether there was an association between HIGS-derived reduction of PtCYC1 or PtMAPK1 transcript levels in Pt and the growth penalty imposed on the fungus in transgenic wheat HD-RNAi lines. Two lines, MAPK1-2169 and CYC1-2236, one from each hp-PtMAPK1RNAi and hp-PtCYC1RNAi transgenic events, were selected, and the relative abundance of fungal PtCYC1 or PtMAPK1 transcripts was measured in transgenic plants representing the range of four different resistance types and compared to levels present in Pt infecting control plants. The analysis showed a close relationship between the endogenous PtMAPK1 and PtCYC1 transcript levels in the fungus invading transgenic wheat plants expressing the corresponding HD-RNAi construct and the amount of fungal biomass accumulating within the infected host tissues (Figure 2b). The highest reduction in fungal PtMAPK1 or PtCYC1 transcript abundance by HD-RNAi was observed in transgenic plants that also displayed the lowest fungal biomass accumulation (rated as HR) followed by resistant and moderately resistant plants. These results indicate that the Figure 2 HD-RNAi of PtMAPK1 and PtCYC1 genes in Pt infecting transgenic wheat plants. (a) Disease symptoms on transgenic wheat T1 plants expressing the hp-PtMAPK1RNAi (line 2169) or hp-PtCYC1RNAi (line 2236) construct and nontransformed control (NTC) plants following challenge with Pt. Representative leaves are shown from transgenic plants of selected lines classified into highly resistant (HR), resistant (R), moderately resistant (MR) and susceptible (S) categories based on per cent fungal biomass reduction (%) in infected tissue compared to controls, along with the scored infection type (IT). Photographs were taken at 12 dpi, the same time point when fungal biomass ratios were analysed by qPCR. (b) Quantitation of fungal PtMAPKI or PtCYC1 transcript abundance in Pt-inoculated leaves of transgenic T1 plants of lines MAPK1-2169 and CYC1-2236 representing different resistance types (as presented in a) by RT-qPCR at 5 dpi. The expression levels of PtMAPK1 or PtCYC1 were normalized to the levels of the Pt succinate dehydrogenase gene and are reported as a percentage, relative to that of the NTC plants (set at 100%). A correlation was observed between HIGS-derived relative reduction in PtMAPK1 and PtCYC1 transcript levels and per cent fungal biomass reduction (% in a). Bars represent mean values � SDs of two independent experiments. Values with different letters indicate statistically significant differences (Student’s t-test; P < 0.05). ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Vinay Panwar et al.1016 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense effectiveness of HD-RNAi to thwart WLR infection in transgenic wheat plants was associated with the silencing efficiency of the targeted genes in the colonizing fungus. Resistance to Pt in transgenic wheat persists in the advanced generation To investigate the impact of HD-RNAi on suppression of Pt infection in advanced generations, three transgenic wheat lines for each of the hp-PtMAPK1RNAi (MAPK1-2163, MAPK1-2169 and MAPK1-2166) and hp-PtCYC1RNAi (CYC1-2236, CYC1- 2248 and CYC1-2224) construct that showed enhanced disease protection were selected and selfed to obtain T2 progeny. PCR and RT-PCR screening revealed the stable inheritance and expression of respective HD-RNAi constructs in transgenic wheat T2 plants (Figure S3). Resistance to Pt was assessed by inoculating 2-week-old transgenic and control wheat plants with fungal urediniospores. Disease phenotype analyses showed that WLR development was radically altered in transgenic wheat T2 lines, as indicated by fewer pustules of much smaller size often sur- rounded by chlorotic or necrotic halos. In contrast, leaves of Pt- infected control plants were completely overgrown by large sporulating uredinia (Figure 4). The disease lesions evaluated at 12 dpi were mainly scored as low IT (IT 1–3) in the T2 plants, reflecting restricted Pt proliferation, which was substantially less severe than the highly susceptible phenotype (IT4) that was consistently observed in control plants. Also, a delay in both the onset and progression of WLR symptoms was seen in transgenic T2 plants with low ITs (IT < 2) as compared to control plants (Figure S4). These results suggested that HIGS of PtMAPK1 or PtCYC1 genes was effective in limiting the growth of pathogenic Pt in transgenic wheat T2 lines. Figure 3 HIGS-derived silencing of targeted Pt genes in transgenic wheat T1 lines. RT-qPCR analyses show significant reduction in transcript abundance of fungal PtMAPK1 or PtCYC1 genes in selective transgenic wheat T1 lines expressing the respective hp-PtMAPK1RNAi (a) and hp- PtCYC1RNAi (b) construct as compared to nontransformed control (NTC) plants. cDNA was generated from total RNA isolated from Pt-challenged leaves at 5 dpi. The expression levels were normalized to the level of the endogenous Pt reference succinate dehydrogenase gene with control set at 100%. No HD-RNAi effect is seen on the other, sequence-unrelated target fungal gene. Bars represent mean values � SDs of three independent sample collections, each using leaves pooled from an average 20–25 infection sites of 4–7 plants. *Statistically significant difference with controls, at P < 0.05 (Student’s t-test). Figure 4 Leaf rust disease resistance in transgenic wheat T2 lines. Pt- inoculated leaves of representative transgenic wheat T2 lines expressing the hp-PtMAPK1RNAi (a) or hp-PtCYC1RNAi (b) construct show suppressed disease symptoms with low infection types (IT), compared to nontransformed control (NTC) plants. Photographs were taken at 12 dpi. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Leaf rust resistance in transgenic wheat by HD-RNAi 1017 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense To quantify the inhibitory effect of HD-RNAi on fungal development, the amount of fungal DNA in Pt-inoculated leaf samples of transgenic T2 and control wheat plants was measured by qPCR. As expected, transgenic T2 lines accumulated substan- tially reduced fungal biomass averaging between 44% and 64%, compared to control plants (Figure 5a). RT-qPCR analysis con- firmed that in planta inhibition of fungal growth in transgenic T2 lines was due to silencing of targeted genes in Pt. Silencing induced by HD-RNAi significantly reduced PtMAPK1 and PtCYC1 relative transcript levels by approximately 44%–58% in Pt colonizing the corresponding transgenic wheat T2 line, compared to the normal expression levels detected in fungi infecting control plants (Figure 5b). Microscopic examination of Pt development in infected tissue revealed normal mycelial morphology with exten- sive colonization of host cells in control plants, whereas hyphal growth was significantly compromised in transgenic lines (Fig- ure 5c). This suggests that depletion of target fungal gene transcripts by HD-RNAi imposed a growth arrest on invading Pt in its transgenic host. Taken together, these results showed that HD- RNAi constructs targeting the PtMAPK1 or PtCYC1 gene in Pt were able to confer genetically stable resistance to WLR infection in transgenic wheat. Production of siRNAs in transgenic wheat We examined whether the silencing of the targeted PtMAPK1 and PtCYC1 genes in Pt by HD-RNAi was associated with the presence of homologous siRNAs in transgenic wheat plants. Northern blot analysis of total RNA extracted from transgenic T2 wheat lines expressing the respective hp-PtMAPK1RNAi or hp-PtCYC1RNAi construct prior to Pt inoculation revealed the presence of corre- sponding siRNA molecules (Figure 6). No detectable hybridization signals were observed in nontransformed control plants. The findings showed that hp-PtMAPK1RNAi or hp-PtCYC1RNAi sequences expressed in transgenic wheat plants were processed by the host silencing machinery into siRNA molecules, and these probably are translocated to fungal cells during infection to act upon the target endogenous mRNAs in a nucleotide-specific manner thereby inhibiting their function in the fungus. Figure 5 Quantitative protection from Pt infection in transgenic wheat T2 lines carrying HD-RNAi constructs. (a) Fungal biomass measurements determined by qPCR in Pt-challenged leaves of hp-PtMAPK1RNAi or hp-PtCYC1RNAi expressing transgenic wheat lines and nontransformed control (NTC) plants at 12 dpi. (b) Relative transcript abundance of target PtMAPK1 or PtCYC1 genes analysed by RT-qPCR in Pt-challenged transgenic wheat T2 lines expressing the corresponding hp-PtMAPK1RNAi or hp-PtCYC1RNAi constructs and NTC plants. cDNA was generated from total RNA isolated at 5 dpi. The expression levels were normalized to the level of the endogenous Pt reference succinate dehydrogenase gene, with control set at 100%. Bars represent mean values � SDs of three independent sample collections, each using leaves pooled from an average 20–25 infection sites of 4–7 plants. *Statistically significant difference with controls, at P < 0.05 (Student’s t-test). (c) Confocal microscopy of Pt infection in wheat leaves at 5 dpi, showing restricted mycelial development in transgenic samples expressing the respective HD-RNAi construct as indicated, compared to the control plant. Bars: 70 lm. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Vinay Panwar et al.1018 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense Discussion In this study, we have shown that developing transgenic wheat plants stably expressing HD-RNAi constructs designed to target Pt genes that are crucial for fungal pathogenicity is an effective strategy to obtain enhanced resistance against the destructive WLR disease. The constitutive expression of HD-RNAi constructs in transgenic wheat lines induced efficient silencing of the targeted PtMAPK1 or PtCYC1 genes in the interacting Pt resulting in subsequent impaired fungal development and a reduction in disease severity and progression in otherwise highly susceptible wheat plants. The engineered resistance trait was heritable and stable in the T2 generation. Suppression of WLR development was correlated with the presence of siRNA molecules specific to the fungal PtMAPK1 or PtCYC1 genes in the transgenic wheat lines expressing the complementary HD-RNAi construct. The in vivo reduction in PtMAPK1 and PtCYC1 transcript levels in Pt feeding on transgenic wheat leaves expressing the homologous siRNA molecules indicated the manifestation of PTGS by HIGS. The MAP-kinase and cyclophilin-encoding genes were selected for the transgenic HD-RNAi approach because they are implicated in playing fundamental roles in the regulation of various physiological and pathogenesis processes in eukaryotes (Wang and Heitman, 2005). In addition, we demonstrated previously by transient HIGS assays that PtMAPK1 and PtCYC1 genes are important for pathogenicity in Pt (Panwar et al., 2013a,b). Cyclophilins have been reported to act as virulence determinants in diverse plant pathogenic fungi such as Magnaporthe grisea (Viaud et al., 2002), Botrytis cinerea (Viaud et al., 2003) and Cryphonectria parasitica (Chen et al., 2011). MAP-kinases from phytopathogenic fungi have also been extensively studied and shown to regulate cell wall biogenesis, morphogenesis, patho- genic development and disease (Mart�ınez-Soto and Ruiz-Herrera, 2015; Zhao et al., 2007). Earlier, it was reported that the PtMAPK1 gene could complement an Ustilago maydis MAPK (kpp2 kpp6) double deletion mutant to restore mating and full virulence (Hu et al., 2007a). Similarly, Fus3/Kss1-type MAPK from Pst, when expressed in ascomycete Fusarium graminearum and M. oryzae, complemented the corresponding kss1 and pmk1 mutants to partially restore vegetative growth and ability to infect plants (Guo et al., 2011). Consistent with the well-defined functions of MAP-kinases and cyclophilins in ensuring pathogen ability to cause disease, we have shown that silencing their homologs in Pt by transgenic HD-RNAi significantly protected host wheat plants from the disease. An obligatory requirement for HIGS is the mobility of silencing signals from host plant species into the invading pathogen. The reduction in endogenous transcript abundance of PtMAPK1 and PtCYC1 genes in Pt by the respective HD-RNAi construct was associated with the presence of complementary siRNAs in the transgenic wheat plants. This indicated that these could be the silencing molecules that initiated in vivo RNAi of the correspond- ing PtMAPK1 and PtCYC1 genes in the colonizing fungus. Recently, siRNAs have been reported to move across interacting organisms, encompassing cross-species silencing and communi- cation between host and pathogens or parasites (Knip et al., 2014; Weiberg et al., 2015). A correlation between HIGS in phytopathogenic fungi and the presence of pathogen gene- specific siRNAs in host plants was also observed in other diverse interactions, suggesting that there is a successful translocation of siRNAs from host plants to parasitic fungi (Andrade et al., 2015; Ghag et al., 2014; Govindarajulu et al., 2015; Nowara et al., 2010; Tinoco et al., 2010). In the wheat–rust pathosystem, the HD-RNAi signals may traverse from host to fungus via the haustorial interphase (Panwar et al., 2013a) utilizing the mech- anism of vesicle-mediated (exosomes) transport (Knip et al., 2014). Despite the progress being made, several aspects of the trans-kingdom RNAi are still poorly understood and more extensive research is required to validate the nature of silencing signals (whether siRNA or precursor dsRNA; e.g., Koch et al., 2016) and their route of delivery during plant–pathogen interac- tions. Inhibition of WLR infection in transgenic wheat lines demon- strated that the degree of silencing achieved by HIGS of the target PtMAPK1 and PtCYC1 genes was sufficient to attenuate fungal pathogenicity. Although the severity of disease was significantly compromised in transgenic wheat leaves, we did not observe complete protection against Pt as few pustules continued to develop from initial infection sites. This can be attributed to partial silencing of the PtMAPK1 or PtCYC1 gene as seen in incomplete knockdown of the target fungal transcripts (Figure 3 and 5b); it is possible that the presence of residual transcripts was sufficient to confer some level of virulence. While RNAi can inhibit transcript expression, it normally does not completely eliminate gene function (Eamens et al., 2008). The fact that we observed variability from R to S in the resistance response within transgenic lines in the T1 generation (Fig. 2 and Table 1) is most likely due to segregation of the transgene locus or loci in this generation. Gene expression levels would be expected to vary between plants where the transgene locus is homozygous vs heterozygous which could explain why some plants are scored as R and some as MR, and segregants lacking the transgene will be scored as S. Other factors possibly affecting HIGS efficiency could be the abundance and turnover rate of target mRNAs. Nevertheless, we observed that the efficacy of HD- RNAi constructs in reducing fungal growth in transgenic wheat lines correlated well with the reduced levels of the corresponding PtMAPK1 and PtCYC1 gene transcripts in the invading fungus. A similar cause–effect relationship between HIGS-derived depletion of pathogen mRNA levels and impairment of host infection competency was also reported in different fungal–plant pathosys- tems (Andrade et al., 2015; Koch et al., 2013). These results indicate that, when exploring HIGS to generate effective plant disease resistance, the extent of target gene silencing in the Figure 6 siRNA detection in transgenic wheat lines. RNA blots showing the presence of PtMAPK1 (a) and PtCYC1 (b) sequence-specific siRNA molecules in transgenic wheat T2 lines expressing the respective hp- PtMAPK1RNAi or hp-PtCYC1RNAi construct. Total RNA was extracted from wheat leaf tissue prior to fungal inoculations. RNA blots were hybridized with target fungal gene-specific DIG-labelled probes. Ethidium bromide-stained rRNA serve as loading controls (LC). The arrow indicates oligonucleotide marker size. NTC, nontransformed control. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Leaf rust resistance in transgenic wheat by HD-RNAi 1019 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense pathogen is likely to be important. However, the choice of effective candidate genes is arguably the most crucial parameter for the success of this technology (Panwar et al., 2016). The WLR disease protection in the transgenic wheat lines was marked by reduced uredinia density and size and by a general delay in symptom appearance in infected host tissue. This supports the notion that the HD-RNAi effect against Pt was quantitative in nature. From a breeding point of view, partial resistance is often regarded as more durable, as it can potentially limit the trade-offs between virulence and reproductive fitness in pathogen population, thereby preventing the emergence of a virulent race (Burdon et al., 2014). Indeed, in our transgenic plants, reduced sporulation would decrease fungal population size in the field. In contrast, rust disease control measures that rely on exploiting major resistance genes that recognize pathogen effectors or their actions have often proved to be ephemeral owing to rapid adaptations by these fungi to escape recognition in resistant plants (Dangl et al., 2013; Stokstad, 2007). In this regard, targeting pathogen genes that are indispensable for morphogenesis and/or pathogenesis by HD-RNAi can presumably provide long-lasting protection because modifying or deleting such essential genes would likely result in major consequences to the pathogen. To increase the degree of functional constraints facing the pathogen upon infection, an effective HIGS design can be silencing multiple genes that are required for pathogenicity rather than relying on single gene targets. In a recent study, simultaneous silencing of three cytochrome P450 lanosterol C- 14a-demethylase (CYP51) genes in F. graminearum by HD-RNAi was shown to completely restrict fungal infection in transgenic Arabidopsis and barley (Koch et al., 2013), although multiple targets may not always be more effective (Chen et al., 2016). Therefore, one approach would be crossing the hp-PtMAPK1R- NAi and hp-PtCYC1RNAi lines to investigate whether or not a cumulative effect would result in more Pt-resistant lines. Earlier, transient HD-RNAi assays were utilized for analysing gene function in the wheat rust pathogens (Panwar et al., 2013a, b; Yin et al., 2015; Zhang et al., 2012). A major limitation of transient RNAi techniques, despite being useful in fast forward functional genomic platforms, is that they elicit uneven and short- lived silencing responses (Bhaskar et al., 2009; Holzberg et al., 2002) which are not suitable for sustainable applications of HIGS in agriculture. Stable HD-RNAi on the other hand should offer many advantages such as the establishment of more uniform gene silencing and the transmission of disease resistance traits to next generations. Production of transgenic crop plants will therefore provide a genetic resource to harvest the full potential of HIGS for engineering long-lasting resistance against rust diseases. However, a drawback of transforming a monocotyle- donous species, such as wheat, is that it is an arduously lengthy and random process, making it practically unsuitable for screen- ing large number of genes to identify efficient HIGS targets. Therefore, a rational approach would be to first validate the function of potential candidate genes identified in the pathogen by transient HD-RNAi approaches and those showing the ability to deter parasitism in the host can then be tested for the development of transgenic disease-resistant varieties. The ability to express HIGS constructs in transgenic plants for control of plant pathogens offers an attractive strategy for managing pests of significant importance (Knip et al., 2014). However, the effectiveness of HD-RNAi to impart genetically stable resistance against rust pathogens has not been explored. An important issue to investigate in the future is whether in subsequent wheat generations RNA-mediated DNA methylation of the silencing constructs may occur, thereby reducing transgene expression levels. Our study provides proof of concept that HIGS against rust fungi can be made operational in stable transgenic plants and results are in line with, and extend beyond, technical advancements made in recent years over the potential application of this technology for fighting rust infections in cereal crop plants (Panwar et al., 2013a,b; Yin et al., 2014). The lack of homology between the chosen Pt sequences and published wheat genomic sequences ensured that the likelihood of any off-target silencing effect against a host gene was low; indeed, we did not observe any sort of deviation from a normal growth phenotype in the transgenic wheat plants. By the same analogy though not directly tested here, the transgenic wheat hp-PtMAPK1RNAi and hp- PtCYC1RNAi lines will likely also prove effective against many other leaf (Pt), stem (Pgt) and stripe (Pst) rust isolates as the sequences used in the HD-RNAi constructs have high homology to the corresponding gene sequences in other sequenced isolates and species. Earlier, we showed through extensive molecular studies using the same constructs in the transient system that the suppression of disease symptoms was highly correlated with the silencing of the target genes and reduction in fungal biomass for various different isolates from all three rust species (Panwar et al., 2013a). Due to complexities associated with conducting reverse genetics studies in rust pathogens, to date, the role of only a few fungal genes involved in the infection processes has been identified. Given that HIGS has emerged as a promising method to analyse gene function in otherwise genetically recalcitrant rust fungi, more studies are needed to address the roles of the large repertoires of predicted pathogenicity-related genes available from genome sequencing projects. The use of such virulence- determining or essential genes that might be functionally difficult to overcome by pathogens promises a potential avenue towards sustainable rust disease protection in wheat cultivars. Experimental procedures Plant materials, Pt inoculations and disease assessment Leaf inoculations were performed using Pt isolate CCDS. Uredin- iospores were increased on susceptible wheat cultivar Fielder plants. The inoculum was prepared by suspending urediniospores in light mineral oil (Baytol) at a rate of 12 mg/mL. Inoculum was applied to the leaves of 2-week-old wheat plants using an air brush. The inoculated plants were transferred to a dark dew chamber preconditioned at 22 °C and 100% relative humidity. After 18 h of incubation, plants were transferred to the green- house and grown under 16:8 h day: night regime at 22 °C. Disease development in terms of infection types (IT) was assessed when rust symptoms were fully expressed at 12 dpi using the seedling 0–4 scale (Dakouri et al., 2013). In short, infection types 0: no flecks or uredinia; IT 1: small uredinia with necrosis; IT 2: small-to-medium uredinia with necrosis: IT 3: moderate-to-large size uredinia with/without chlorosis; IT 4: very large uredinia. The variation within each class is indicated by the use of – (less than average) and + (more than average). Construction of HD-RNAi vectors The Pt genes encoding a MAP-kinase (PtMAPK1) and a cyclophilin (PtCYC1) that we earlier characterized to have an essential role in pathogenicity were used in this study (Panwar et al., 2013a). Binary vectors for wheat transformation were constructed as described earlier (Panwar et al., 2013a). Briefly, a 520-bp ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Vinay Panwar et al.1020 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense fragment representing the 30-end coding region of the PtMAPK1 gene (Hu et al., 2007a) GenBank accession #68303937) was PCR-amplified from cDNA of Pt Race 1 isolate BBBD using primers PtMAPK1-F and PtMAPK1-R (Table S1). The amplified segment was cloned into pENTERTM vector using the D-TOPO� Cloning kit (Invitrogen, Carlsbad, CA) and the recombinant entry module recombined in the cereal-specific binary vector pIPK007 (Him- melbach et al., 2007) to generate hp-PtMAPK1RNAi. The result- ing RNAi vector allowed expression of hpRNA molecules under the control of the Zea mays ubiquitin promoter and terminator sequences. To generate hp-PtCYC1RNAi, a 501-bp 30-end coding region of PtCYC1 (Hu et al., 2007b) GenBank accession #BU672663) was amplified from the EST sequence (PtCon- tig6674) using primers PtCYC1-F and PtCYC1-R (Table S1). The amplified fragment was sequentially cloned into plasmid pUB- leX1-RNAi to generate an intron-spliced hpRNA segment that was subsequently cloned into the binary vector pMCG161 (http:// www.chromdb.org/rnai/vector_info.html) under the control of the Cauliflower mosaic virus 35S promoter and octopine synthase terminator. All final constructs were verified by sequencing. Genomic DNA extraction Fungal genomic DNA was extracted from germinating uredin- iospores as described earlier (Panwar et al., 2013a). Uredin- iospores were germinated by dusting them over sterile water in dishes for over 8 h using a volatile nonanol solution [1.56 lL nonanol (Sigma-Aldrich, ON, Canada), 1 mL acetone, 19 mL of ddH20] spotted on filter paper suspended in the lids to stimulate germination. Plant genomic DNA was extracted from leaf tissue using the DNeasy plant Mini Kit (Qiagen, Mississauga, ON, Canada) as described by the manufacturer. Total RNA extraction and cDNA synthesis Total RNA was extracted using TRIzol reagent (Invitrogen) as described by the manufacturer. The quality and concentration of total RNA were assessed via denaturing 1.2% agarose gel and 260/280ABS measurements on a NanoDrop spectrophotometer (ThermoFisher Scientific, Mississauga, ON, Canada). Total RNA was treated with DNase I using the TURBO DNA-freeTM kit (ThermoFisher Scientific, Mississauga, ON, Canada) following manufacturer instructions to remove DNA contamination prior to cDNA synthesis. The first stand of cDNA was synthesized using 1.0 lg of RNA sample with a poly(A)18 oligonucleotide primer using the Reverse Transcriptase III (ThermoFisher Scientific, Mississauga, ON, Canada) according to the manufacturer’s protocol. Wheat transformation Plasmid suspensions at a concentration of 1 lg/lL in 10 mM Tris (pH 8.0) were used for wheat transformation by particle bombardment as described earlier (Jordan, 2000) with minor modifications. Immature embryos of wheat cv. Fielder were cultured on Murashige and Skoog medium supplemented with 0.5 mg/L thiamine, 100 mg/L glutamine, 1 mM niacinamide, 2 mg/L 2,4-D and solidified with 2.5 g/L phytagel (Sigma- Aldrich). A 2-day preculture was used prior to bombardment using 0.6-lm gold, 650 psi rupture discs and a 3 cm target distance. After callus growth on RG5 medium for 14 days, calli were transferred to 1/2 strength MS with 25 mg/L of hygromycin B (Sigma-Aldrich) for 21 days under a 16:8 h light: dark regime prior to subculture every 14 days on the same selectable medium, followed by planting in soil in growth cabinets to grow to maturity. Integration and expression of HD-RNAi transgene in transgenic wheat plants were confirmed by PCR and RT-PCR analysis using candidate Pt gene-specific primers (Table S1). PCRs were subjected to 30 cycles of denaturation (94 °C, 1 min), annealing (58 °C, 30 s) and extension (72 °C, 30 s). The good- ness of fit of the observed segregation ratio for the transgene was tested against the Mendelian segregation ratio of 3:1 and 15:1 using the chi-square (v2) test. Quantitative real-time PCR analysis Quantitative real-time PCR analysis was carried out using a CFX96TM Real-Time PCR machine (Bio-Rad, Mississauga, ON, Canada). Transcript levels of the target Pt genes, PtMAPK1 and PtCYC1, were measured using cDNA prepared from total RNA isolated from wheat leaf tissue inoculated with Pt by RT-qPCR as described previously (Panwar et al., 2013a). Specific primers for each gene were designed using the Primer 3.0 program and are shown in Table S1. Quantification results were analysed using the comparative 2�DDCT method. To assess fungal transcript levels, the Pt succinate dehydrogenase (SDH) gene was used as normalizing reference gene. To quantify fungal biomass, the ratio of single copy PtRTP1 and the wheat TaEF1 genes was assessed in genomic DNA isolated from infected leaves using gene-specific primers (Panwar et al., 2013a). The relative amounts of PCR product of PtRTP1 and EF1 in Pt-infected samples were calculated using generated gene-specific standard curves to quantify the Pt and wheat gDNA, respectively. Confocal microscopy Microscopic analysis of mycelium growth in Pt-inoculated wheat leaves was carried out using a Leica SP2-AOBS laser scanning confocal microscope (Leica, Mannheim, Germany) at 5 dpi. Fungal structures were stained with Uvitex-2B (Sigma-Aldrich), and Acridine Orange (Sigma-Aldrich) was used to stain plant cell walls following procedures as described previously (Panwar et al., 2013b). The fluorescence of Uvitex-2B and Acridine Orange was detected by excitation at 405 and 514 nm, and scanning was performed with filter setting at 411–485 and 550–560 nm, respectively. siRNA detection Small RNA detection was performed using total RNA extracted from wheat leaves prior to Pt inoculation. Equal amounts of RNA (15 lg) from different samples were resolved on denaturing 15% urea–polyacrylamide gels and transferred to Hybond-N+ mem- branes (Amersham Biosciences, Mississauga, ON, Canada) using a semidry blotter (Bio-Rad). Transferred RNA was then cross-linked to the membrane by UV with an energy output of 1200 lJ/cm2 for 5 min. DIG-labelled DNA probes were synthesized using target fungal gene sequences as template. Prehybridization was performed using DIG Easy Hybridization buffer (Sigma-Aldrich) followed by hybridization at 42 °C overnight. The membranes were washed twice with 19 SSC, 0.1% SDS for a total of 15 min at 42 °C. Detection of the hybridized probe was performed using the DIG Nucleic Acid Detection Kit (Sigma-Aldrich) following the manufacturer instructions. Photoemissions were captured using a ChemiDoc-IT Imaging System (Bio-Rad). Acknowledgements This work was supported by grants from the Western Grains Research Foundation and the Alberta Wheat Commission (Project ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Leaf rust resistance in transgenic wheat by HD-RNAi 1021 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense http://www.ncbi.nlm.nih.gov/nuccore/BU672663 http://www.chromdb.org/rnai/vector_info.html http://www.chromdb.org/rnai/vector_info.html #A08853) and from the Canadian Wheat Improvement Flagship as part of the Canadian Wheat Alliance. The authors declare no conflict of interest. References Andrade, C.M., Tinoco, M.L.P., Rieth, A.F., Maia, F.C.O. and Arag~ao, F.J.L. (2015) Host-induced gene silencing in the necrotrophic fungal pathogen Sclerotinia sclerotiorum. Plant. Pathol. 65, 626–632. Baulcombe, D. (2004) RNA silencing in plants. Nature, 431, 356–363. Bhaskar, P.B., Venkateshwaran, M., Wu, L., An�e, J.-M. and Jiang, J. (2009) Agrobacterium-mediated transient gene expression and silencing: a rapid tool for functional gene assay in potato. PLoS ONE, 4, e5812. Bolton, M.D., Kolmer, J.A. and Garvin, D.F. (2008) Wheat leaf rust caused by Puccinia triticina. Mol. Plant Pathol. 9, 563–575. Burdon, J.J., Barrett, L.G., Rebetzke, G. and Thrall, P.H. (2014) Guiding deployment of resistance in cereals using evolutionary principles. Evol. Appl. 7, 609–624. Cantu, D., Govindarajulu, M., Kozik, A., Wang, M., Chen, X., Kojima, K.K., Jurka, J. et al. (2011) Next generation sequencing provides rapid access to the genome of Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe rust. PLoS ONE, 6, e24230. Chen, M.-M., Jiang, M., Shang, J., Lan, X., Yang, F., Huang, J., Nuss, D.L. et al. (2011) CYP1, a hypovirus-regulated cyclophilin, is required for virulence in the chestnut blight fungus. Mol. Plant Pathol. 12, 239–246. Chen, W., Wellings, C., Chen, X., Kang, Z. and Liu, T. (2014) Wheat stripe (yellow) rust caused by Puccinia striiformis f. sp. tritici. Mol. Plant Pathol. 15, 433–446. Chen, W., Kastner, C., Nowara, D., Oliveira-Garcia, E., Rutten, T., Zhao, Y., Deising, H.B. et al. (2016) Host-induced silencing of Fusarium culmorum genes protects wheat from infection. J. Exp. Bot. 67, 4979–4991. Cuomo, C.A., Bakkeren, G., Khalil, H.B., Panwar, V., Joly, D., Linning, R., Sakthikumar, S. et al. (2017) Comparative analysis highlights variable genome content of wheat rusts and divergence of the mating loci. G3 (Bethesda), 7, 361–376. Dakouri, A., McCallum, B.D., Radovanovic, N. and Cloutier, S. (2013) Molecular and phenotypic characterization of seedling and adult plant leaf rust resistance in a world wheat collection. Mol. Breed. 32, 663–677. Dangl, J.L., Horvath, D.M. and Staskawicz, B.J. (2013) Pivoting the plant immune system from dissection to deployment. Science, 341, 746–751. Draz, I.S., Abou-Elseoud, M.S., Kamara, A.-E.M., Alaa-Eldein, O.A.-E. and El- Bebany, A.F. (2015) Screening of wheat genotypes for leaf rust resistance along with grain yield. Ann. Agric. Sci. 60, 29–39. Duplessis, S., Cuomo, C., Lin, Y.-C., Aerts, A., Tisserant, E., Veneault- Fourrey, C., Joly, D.L. et al. (2011) Obligate biotrophy features unraveled by the genomic analysis of rust fungi. Proc. Nat. Acad. Sci. USA, 108, 9166–9171. Eamens, A., Wang, M.-B., Smith, N.A. and Waterhouse, P.M. (2008) RNA silencing in plants: yesterday, today, and tomorrow. Plant Physiol. 147, 456– 468. Ellis, J.G., Lagudah, E.S., Spielmeyer, W. and Dodds, P.N. (2014) The past, present and future of breeding rust resistant wheat. Front. Plant Sci. 5, 641. Flor, H. (1971) Current status of the gene-for-gene concept. Annu. Rev. Phytopathol. 9, 275–296. Garnica, D.P., Nemri, A., Upadhyaya, N.M., Rathjen, J.P. and Dodds, P.N. (2014) The ins and outs of rust haustoria. PLoS Pathog. 10, e1004329. Ghag, S.B., Shekhawat, U.K.S. and Ganapathi, T.R. (2014) Host-induced post- transcriptional hairpin RNA-mediated gene silencing of vital fungal genes confers efficient resistance against Fusariumwilt in banana. Plant Biotechnol. J. 12, 541–553. Govindarajulu, M., Epstein, L., Wroblewski, T. and Michelmore, R.W. (2015) Host-induced gene silencing inhibits the biotrophic pathogen causing downy mildew of lettuce. Plant Biotechnol. J. 13, 875–883. Guo, J., Dai, X., Xu, J.-R., Wang, Y., Bai, P., Liu, F., Duan, Y. et al. (2011) Molecular characterization of a Fus3/Kss1 Type MAPK from Puccinia striiformis f. sp. tritici, PsMAPK1. PLoS ONE, 6, e21895. Helfer, S. (2014) Rust fungi and global change. New Phytol. 201, 770–780. Herrera-Foessel, S.A., Singh, R.P., Huerta-Espino, J., Crossa, J., Yuen, J. and Djurle, A. (2006) Effect of leaf rust on grain yield and yield traits of durum wheats with race-specific and slow-rusting resistance to leaf rust. Plant Dis. 90, 1065–1072. Himmelbach, A., Zierold, U., Hensel, G., Riechen, J., Douchkov, D., Schweizer, P. and Kumlehn, J. (2007) A set of modular binary vectors for transformation of cereals. Plant Physiol. 145, 1192–1200. Holzberg, S., Brosio, P., Gross, C. and Pogue, G.P. (2002) Barley stripe mosaic virus-induced gene silencing in a monocot plant. Plant J. 30, 315–327. Hu, G., Kamp, A., Linning, R., Naik, S. and Bakkeren, G. (2007a) Complementation of Ustilago maydis MAPK mutants by a wheat leaf rust, Puccinia triticina homolog: potential for functional analyses of rust genes. Mol. Plant Microbe Interact. 20, 637–647. Hu, G.G., Linning, R., Mccallum, B., Banks, T., Cloutier, S. and Butterfield, Y. (2007b) Generation of a wheat leaf rust, Puccinia triticina, EST database from stage-specific cDNA libraries. Mol. Plant Pathol. 8, 451–467. Huang, G., Allen, R., Davis, E.L., Baum, T.J. and Hussey, R.S. (2006) Engineering broad root-knot resistance in transgenic plants by RNAi silencing of a conserved and essential root-knot nematode parasitism gene. Proc. Natl Acad. Sci. USA, 103, 14302–14306. Huerta-Espino, J., Singh, R.P., Germ�an, S., McCallum, B.D., Park, R.F., Chen, W.Q., Bhardwaj, S.C. et al. (2011) Global status of wheat leaf rust caused by Puccinia triticina. Euphytica, 179, 143–160. Jahan, S.N., �Asman, A.K.M., Corcoran, P., Fogelqvist, J., Vetukuri, R.R. and Dixelius, C. (2015) Plant-mediated gene silencing restricts growth of the potato late blight pathogen Phytophthora infestans. J. Exp. Bot. 66, 2785– 2794. Jordan, M.C. (2000) Green fluorescent protein as a visual marker for wheat transformation. Plant Cell Rep. 19, 1069–1075. Knip, M., Constantin, M.E. and Thordal-Christensen, H. (2014) Trans-kingdom cross-talk: small RNAs on the move. PLoS Genet. 10, e1004602. Koch, A. and Kogel, K.-H. (2014) New wind in the sails: improving the agronomic value of crop plants through RNAi-mediated gene silencing. Plant Biotechnol. J. 12, 821–831. Koch, A., Kumar, N., Weber, L., Keller, H., Imani, J. and Kogel, K.-H. (2013) Host-induced gene silencing of cytochrome P450 lanosterol C14a- demethylase–encoding genes confers strong resistance to Fusarium species. Proc. Natl Acad. Sci. USA, 110, 19324–19329. Koch, A., Biedenkopf, D., Furch, A., Weber, L., Rossbach, O., Abdellatef, E., Linicus, L. et al. (2016) An RNAi-based control of Fusarium graminearum infections through spraying of long dsRNAs involves a plant passage and is controlled by the fungal silencing machinery. PLoS Pathog. 12, e1005901. Kolmer, J.A. (2005) Tracking wheat rust on a continental scale. Curr. Opin. Plant Biol. 8, 441–449. Mart�ınez-Soto, D. and Ruiz-Herrera, J. (2015) Regulation of the expression of the whole genome of Ustilago maydis by a MAPK pathway. Arch. Microbiol. 197, 575–588. McCallum, B.D., Hiebert, C.W., Cloutier, S., Bakkeren, G., Rosa, S.B., Humphreys, D.G., Marais, G.F. et al. (2016) A review of wheat leaf rust research and the development of resistant cultivars in Canada. Can. J. Plant Pathol. 38, 1–18. Nowara, D., Gay, A., Lacomme, C., Shaw, J., Ridout, C., Douchkov, D., Hensel, G. et al. (2010) HIGS: host-induced gene silencing in the obligate biotrophic fungal pathogen Blumeria graminis. Plant Cell, 22, 3130–3141. Panwar, V., McCallum, B. and Bakkeren, G. (2013a) Endogenous silencing of Puccinia triticina pathogenicity genes through in planta-expressed sequences leads to the suppression of rust diseases on wheat. Plant J. 73, 521–532. Panwar, V., McCallum, B. and Bakkeren, G. (2013b) Host-induced gene silencing of wheat leaf rust fungus Puccinia triticina pathogenicity genes mediated by the Barley stripe mosaic virus. Plant Mol. Biol. 81, 595–608. Panwar, V., McCallum, B., Jordan, M., Loewen, M., Fobert, P., McCartney, C.A. and Bakkeren, G. (2016) RNA silencing approaches for identifying pathogenicity and virulence elements towards engineering crop resistance to plant pathogenic fungi. CAB Rev. 11, 1–13. Petre, B., Saunders, D.G.O., Sklenar, J., Lorrain, C., Win, J., Duplessis, S. and Kamoun, S. (2015) Candidate effector proteins of the rust pathogen ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Vinay Panwar et al.1022 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense Melampsora larici-populina target diverse plant cell compartments.Mol. Plant Microbe Interact. 28, 689–700. Singh, R.P., Hodson, D.P., Huerta-Espino, J., Jin, Y., Bhavani, S., Njau, P., Herrera-Foessel, S. et al. (2011) The emergence of Ug99 races of the stem rust fungus is a threat to world wheat production. Annu. Rev. Phytopathol. 49, 465–481. Singla, J., L€uthi, L., Wicker, T., Bansal, U., Krattinger, S.G. and Keller, B. (2017) Characterization of Lr75: a partial, broad-spectrum leaf rust resistance gene in wheat. Theor. Appl. Genet. 130, 1–12. Song, Y. and Thomma, B.P.H.J. (2016) Host-induced gene silencing compromises Verticillium wilt in tomato and Arabidopsis. Mol. Plant Pathol. https://doi.org/10.1111/mpp.12500. Stokstad, E. (2007) Deadly wheat fungus threatens world’s breadbaskets. Science, 315, 1786–1787. Struck, C. (2015) Amino acid uptake in rust fungi. Front. Plant Sci. 6, 40. Tinoco, M., Dias, B., Dall’Astta, R., Pamphile, J. and Arag~ao, F. (2010) In vivo trans-specific gene silencing in fungal cells by in planta expression of a double-stranded RNA. BMC Biol. 8, 1–11. Viaud, M.C., Balhad�ere, P.V. and Talbot, N.J. (2002) A Magnaporthe grisea cyclophilin acts as a virulence determinant during plant infection. Plant Cell, 14, 917–930. Viaud, M., Brunet-Simon, A., Brygoo, Y., Pradier, J.-M. and Levis, C. (2003) Cyclophilin A and calcineurin functions investigated by gene inactivation, cyclosporin A inhibition and cDNA arrays approaches in the phytopathogenic fungus Botrytis cinerea. Mol. Microbiol. 50, 1451–1465. Wang, P. and Heitman, J. (2005) The cyclophilins. Genome Biol. 6, 226. Webb, C.A. and Fellers, J.P. (2006) Cereal rust fungi genomics and the pursuit of virulence and avirulence factors. FEMS Microbiol. Lett. 264, 1–7. Weiberg, A., Bellinger, M. and Jin, H. (2015) Conversations between kingdoms: small RNAs. Curr. Opin. Biotechnol. 32, 207–215. Wilson, R.C. and Doudna, J.A. (2013) Molecular mechanisms of RNA interference. Annu. Rev. Biophys. 42, 217–239. Wu, J.Q., Sakthikumar, S., Dong, C., Zhang, P., Cuomo, C.A. and Park, R.F. (2017) Comparative genomics integrated with association analysis identifies candidate effector genes corresponding to Lr20 in phenotype-paired Puccinia triticina isolates from Australia. Front. Plant Sci. 8, 148. Yin, C., Jurgenson, J.E. and Hulbert, S.H. (2010) Development of a host-induced RNAi system in the wheat stripe rust fungus Puccinia striiformis f. sp. tritici. Mol. Plant Microbe Interact. 24, 554–561. Yin, C., Park, J.J., Gang, D.R. and Hulbert, S.H. (2014) Characterization of a tryptophan 2-monooxygenase gene from Puccinia graminis f. sp. tritici involved in auxin biosynthesis and rust pathogenicity. Mol. Plant Microbe Interact. 27, 227–235. Yin, C., Downey, S.I., Klages-Mundt, N.L., Ramachandran, S., Chen, X., Szabo, L.J., Pumphrey, M. et al. (2015) Identification of promising host-induced silencing targets among genes preferentially transcribed in haustoria of Puccinia. BMC Genom. 16, 579. Zhang, H., Guo, J., Voegele, R.T., Zhang, J., Duan, Y., Luo, H. and Kang, Z. (2012) Functional characterization of calcineurin homologs PsCNA1/PsCNB1 in Puccinia striiformis f. sp. tritici using a host-induced RNAi system. PLoS ONE, 7, e49262. Zhao, X., Mehrabi, R. and Xu, J.-R. (2007) Mitogen-activated protein kinase pathways and fungal pathogenesis. Eukaryot. Cell, 6, 1701–1714. Zheng, W., Huang, L., Huang, J., Wang, X., Chen, X. and Zhao, J. (2013) High genome heterozygosity and endemic genetic recombination in the wheat stripe rust fungus. Nat. Commun. 4, 2673. Supporting information Additional Supporting Information may be found online in the supporting information tab for this article: Figure S1 Molecular analysis of transgenic wheat T1 lines. Integration of hp-PtMAPK1RNAi or hp-PtCYC1RNAi construct in transgenic plants analyzed by PCR-amplification of the corre- sponding PtMAPK1 (a) or PtCYC1 (b) transgenes. Genomic DNA extracted from transgenic and non-transformed control (NTC) plants prior to Pt inoculation was used as template in PCR reactions. Expression of PtMAPK1RNAi or hp-PtCYC1RNAi in the same plants analysed by RT-PCR using primers specific to the corresponding PtMAPK1 (c) or PtCYC1 (d) transgene. cDNA prepared from total RNA extracted from transgenic and control plants prior to Pt inoculations was used as template in RT-PCR reactions. M indicates 1 kb DNA ladder (Invitrogen); plasmid DNA served as positive control (PC). The PCR products were fraction- ated on a 1% agarose gel. Results from representative plants of each selected line are shown (lanes 1–20). Figure S2 Delay in progression of symptoms in leaves of hp- PtMAPK1RNAi or hp-PtCYC1RNAi expressing wheat T1 plants. Pt inoculated transgenic leaves are characterized by a slower and more restrictive uredinia development as compared to non- transformed control Fielder plants. Representative leaves from transgenic T1 lines MAPK1-2166 and CYC1-2236 are shown. Photographs were taken at 6 dpi. Figure S3 Molecular analysis of transgenic wheat T2 lines. Integration of hp-PtMAPK1RNAi or hp-PtCYC1RNAi constructs in transgenic plants analyzed by PCR-amplification of the corre- sponding PtMAPK1 (a) and PtCYC1 (b) transgenes. Genomic DNA extracted from transgenic and non-transformed control (NTC) plants prior to Pt inoculation was used as template in PCR reactions. Expression of PtMAPK1RNAi or hp-PtCYC1RNAi in the same plants was analysed by RT-PCR using primers specific to the corresponding PtMAPK1 (c) or PtCYC1 (d) transgene. cDNA prepared from total RNA extracted from transgenic and control plants prior to Pt inoculations was used as template in RT-PCR reactions. M indicates 1 kb DNA ladder (Invitrogen); BC, buffer control; plasmid DNA served as positive control (PC). The PCR products were fractionated on a 1% agarose gel. Results from representative plants of each selected line are shown (lanes 1–3). Figure S4 Delay in progression of symptoms in leaves of hp- PtMAPK1RNAi or hp-PtCYC1RNAi expressing wheat T2 plants. Pt inoculated transgenic leaves are characterized by a slower and more restrictive uredinia development as compared to non- transformed control Fielder plants (NTC). Representative leaves from transgenic T2 lines MAPK1-2166 and CYC1-2236 are shown. Photographs were taken at 6 dpi. Table S1 Primers used in the research. Table S2 Segregation analysis of transgenic wheat T1 lines. ª 2017 Her Majesty the Queen in Right of Canada. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 16, 1013–1023 Leaf rust resistance in transgenic wheat by HD-RNAi 1023 14677652, 2018, 5, D ow nloaded from https://onlinelibrary.w iley.com /doi/10.1111/pbi.12845 by C ochrane C anada Provision, W iley O nline L ibrary on [06/04/2023]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense https://doi.org/10.1111/mpp.12500