Integration of early disease-resistance phenotyping, histological characterization, and transcriptome sequencing reveals insights into downy mildew resistance in impatiens

Abstract

Downy mildew (DM), caused by obligate parasitic oomycetes, is a destructive disease for a wide range of crops worldwide. Recent outbreaks of impatiens downy mildew (IDM) in many countries have caused huge economic losses. A system to reveal plant–pathogen interactions in the early stage of infection and quickly assess resistance/ susceptibility of plants to DM is desired. In this study, we established an early and rapid system to achieve these goals using impatiens as a model. Thirty-two cultivars of Impatiens walleriana and I. hawkeri were evaluated for their responses to IDM at cotyledon, first/second pair of true leaf, and mature plant stages. All I. walleriana cultivars were highly susceptible to IDM. While all I. hawkeri cultivars were resistant to IDM starting at the first true leaf stage, many (14/16) were susceptible to IDM at the cotyledon stage. Two cultivars showed resistance even at the cotyledon stage. Histological characterization showed that the resistance mechanism of the I. hawkeri cultivars resembles that in grapevine and type II resistance in sunflower. By integrating full-length transcriptome sequencing (Iso-Seq) and RNA- Seq, we constructed the first reference transcriptome for Impatiens comprised of 48,758 sequences with an N50 length of 2060 bp. Comparative transcriptome and qRT-PCR analyses revealed strong candidate genes for IDM resistance, including three resistance genes orthologous to the sunflower gene RGC203, a potential candidate associated with DM resistance. Our approach of integrating early disease-resistance phenotyping, histological characterization, and transcriptome analysis lay a solid foundation to improve DM resistance in impatiens and may provide a model for other crops.

References

[1] Thines, M., Voglmayr, H. & Göker, M. in Oomycete Genetics and Genomics: Diversity, Interactions, and Research Tools. (eds Lamour, K. & Kamoun, S.) 47-75 (John Wiley & Sons, 2009).

[2] Delmas, C. E. et al. Adaptation of a plant pathogen to partial host resistance: selection for greater aggressiveness in grapevine downy mildew. Evol. Appl. 9, 709-725 (2016).

[3] Zhang, N., Lindhout, P., Niks, R. & Jeuken, M. Genetic dissection of Lactuca saligna nonhost resistance to downy mildew at various lettuce developmental stages. Plant Pathol. 58, 923-932 (2009).

[4] Meena, P., Thomas, L. & Singh, D. Assessment of yield losses in Brassica juncea due to downy mildew (Hyaloperonospora brassicae). J. Oilseed Brassica 1, 73-77 (2016).

[5] Chelpuri, D. et al. Mapping quantitative trait loci (QTLs) associated with resistance to major pathotype-isolates of pearl millet downy mildew pathogen. Eur. J. Plant Pathol. 154, 983-994 (2019).

[6] Lebeda, A., Křístková, E., Sedláková, B. & Widrlechner, M. P. Recent advances in cucurbit downy mildew research and their contribution to the development of integrated protection of cucurbits. Ecol. Evol. Org. Biol. Proc. 3. https://lib.dr.iastate.edu/eeob_ag_conf/3 (2019).

[7] Delmas, C. E. et al. Soft selective sweeps in fungicide resistance evolution: recurrent mutations without fitness costs in grapevine downy mildew. Mol. Ecol. 26, 1936-1951 (2017).

[8] Callan, B. & Carris, L. in Biodiversity of Fungi: Inventory and Monitoring Methods. (eds Mueller, G. M., Bills, G. F. & Foster, M. S.) 105-126 (Elsevier Academic Press, 2004).

[9] Pecrix, Y., Penouilh-Suzette, C., Muños, S., Vear, F. & Godiard, L. Ten broad spectrum resistances to downy mildew physically mapped on the sunflower genome. Front. Plant Sci. 9, 1780 (2018).

[10] Parra, L. et al. Rationalization of genes for resistance to Bremia lactucae in lettuce. Euphytica 210, 309-326 (2016).

[11] Vezzulli, S. et al. The Rpv3-3 haplotype and stilbenoid induction mediate downy mildew resistance in a grapevine interspecific population. Front. Plant Sci. 10, 234 (2019).

[12] Radwan, O., Bouzidi, M. F., Nicolas, P. & Mouzeyar, S. Development of PCR markers for the Pl5/Pl8 locus for resistance to Plasmopara halstedii in sunflower, Helianthus annuus L. from complete CC-NBS-LRR sequences. Theor. Appl. Genet. 109, 176-185 (2004).

[13] Michelmore, R. & Wong, J. Classical and molecular genetics of Bremia lactucae, cause of lettuce downy mildew. Eur. J. Plant Pathol. 122, 19 (2008).

[14] Hermanns, M., Slusarenko, A. J. & Schlaich, N. L. Organ-specificity in a plant disease is determined independently of R gene signaling. Mol. Plant-Microbe Interact 16, 752-759 (2003).

[15] She, H. et al. Fine mapping and candidate gene screening of the downy mildew resistance gene RPF1 in Spinach. Theor. Appl. Genet. 131, 2529-2541 (2018).

[16] Morgan, R. J. (ed.) Impatiens: the Vibrant World of Busy Lizzies, Balsams, and Touch-me-nots (Timber Press, 2007).

[17] Lim, T. (ed.) Edible Medicinal and Non-medicinal Plants (Springer, 2014).

[18] Uchneat, M. S. (ed.) Flower Breeding and Genetics (Springer, 2007).

[19] US Department of Agriculture, National Agricultural Statistics Service . https://www.nass.usda.gov (2019).

[20] Wang, W., He, Y., Cao, Z. & Deng, Z. Induction of tetraploids in impatiens (Impatiens walleriana) and characterization of their changes in morphology and resistance to downy mildew. HortScience 53, 925-931 (2018).

[21] Vajna, L. First report of Plasmopara obducens on impatiens (Impatiens walleriana) in Hungary. New Dis. Rep 24, 2044-0588 (2011).

[22] Cunnington, J., Aldaoud, R., Loh, M., Washington, W. & Irvine, G. First record of Plasmopara obducens (downy mildew) on impatiens in Australia. Plant Pathol. 57, 371 (2008).

[23] Crouch, J., Ko, M. & McKemy, J. First report of impatiens downy mildew outbreaks caused by Plasmopara obducens throughout the Hawai’ian islands. Plant Dis. 98, 696- 696 (2014).

[24] Palmateer, A., Lopez, P., Seijo, T. & Peres, N. Severe outbreak of downy mildew caused by Plasmopara obducens on Impatiens walleriana in Florida. Plant Dis. 97, 687- 687 (2013).

[25] US Department of Agriculture, National Agricultural Statistics Service . https://www.nass.usda.gov (2016).

[26] Conner, K., Olive, J., Hagan, A., Zhang, L. & Bloodworth, M. First report of impatiens downy mildew caused by Plasmopara obducens in Alabama. Plant Dis. 98, 1006- 1006 (2014).

[27] Warfield, C. Y. Downy mildew of impatiens. GrowerTalks 75, 78-86 (2012).

[28] Hansen, M. A., Bush, E. A., Latimer, J. G. & Hong, C. Impatiens Downy Mildew (Virginia Coop Extension Virgia State Univ., 2013).

[29] Satou, M., Sugawara, K., Nagashima, S., Tsukamoto, T. & Matsushita, Y. Downy mildew of busy lizzie caused by Plasmopara obducens in Japan. J. Gen. Plant Pathol. 79, 205-208 (2013).

[30] Jones, D. & O’Neill, T. Impatiens downy mildew. East Malling, UK: Horticultural Development Council FactSheet 5, 8 (2004).

[31] Choi, Y. J., Han, J. G., Park, M. J. & Shin, H. D. Downy mildew of Impatiens balsamina and I. walleriana in Korea. Plant Pathol. J. 25, 433- 433 (2009).

[32] Keach, J. E. & Bridgen, M. Towards improvement of Impatiens. Acta Hortic. 1140, 317-325 (2016).

[33] Eskandari, F. & Shishkoff, N. Systemic infection of Impatiens balsamina through inoculation of roots with vegetative sporangia of the Impatiens downy mildew (Plasmopara obducens). Phytopathol. 107, S4.3 (2017).

[34] Suarez, S. & Palmateer, A. Overview of impatiens downy mildew in Florida. Phytopathol. 104, 115 (2014).

[35] Suarez, S., Lopez, P., Chase, A. & Palmateer, A. Preventative fungicide applications in production and their impact on residual efficacy against impatiens downy mildew in the landscape. Phytopathol. 106, 68 (2016).

[36] Gessler, C., Pertot, I. & Perazzolli, M. Plasmopara viticola: a review of knowledge on downy mildew of grapevine and effective disease management. Phytopathol. Mediterr. 50, 3-44 (2011).

[37] Stassen, J. H. et al. Effector identification in the lettuce downy mildew Bremia lactucae by massively parallel transcriptome sequencing. Mol. Plant Pathol. 13, 719-731 (2012).

[38] Li, X. et al. Comparative transcriptome analysis reveals defense-related genes and pathways against downy mildew in Vitis amurensis grapevine. Plant Physiol. Biochem. 95, 1-14 (2015).

[39] Mestre, P. et al. Comparative analysis of expressed CRN and RXLR effectors from two Plasmopara species causing grapevine and sunflower downy mildew. Plant Pathol. 65, 767-781 (2016).

[40] Kandel, S. L. et al. Transcriptional analyses of differential cultivars during resistant and susceptible interactions with Peronospora effusa, the causal agent of spinach downy mildew. Sci. Rep. 10, 1-13 (2020).

[41] Bhattarai, K., Wang, W., Cao, Z. & Deng, Z. Comparative analysis of impatiens leaf transcriptomes reveal candidate genes for resistance to downy mildew caused by Plasmopara obducens. Int. J. Mol. Sci. 19, 2057 (2018).

[42] Suarez, S., Naveed, Z. A. & Ali, G. Transcriptional profiling of Impatiens walleriana genes through different stages of downy mildew infection reveals novel genes involved in disease susceptibility. Preprint at https://www.biorxiv.org/content/10.1101/622480v2 (2019).

[43] Rhoads, A. & Au, K. F. PacBio sequencing and its applications. Genomics Proteomics Bioinformatics 13, 278-289 (2015).

[44] Florida Automated Weather Network . https://fawn.ifas.ufl.edu (2015).

[45] Radwan, O., Bouzidi, M. F. & Mouzeyar, S. Molecular characterization of two types of resistance in sunflower to Plasmopara halstedii, the causal agent of downy mildew. Phytopathol. 101, 970-979 (2011).

[46] Bachlava, E. et al. Downy mildew (Pl₈ and Pl₁₄) and rust (R_Adv) resistance genes reside in close proximity to tandemly duplicated clusters of non-TIR-like NBS-LRR-encoding genes on sunflower chromosomes 1 and 13 . Theor. Appl. Genet. 122, 1211-1221 (2011).

[47] Song, J. et al. Gene RB cloned from Solanum bulbocastanum confers broad spectrum resistance to potato late blight. Proc. Natl Acad. Sci. USA 100, 9128-9133 (2003).

[48] Van Der Vossen, E. et al. An ancient R gene from the wild potato species Solanum bulbocastanum confers broad-spectrum resistance to Phytophthora infestans in cultivated potato and tomato. Plant J. 36, 867-882 (2003).

[49] Coelho, P. S., Valério, L. & Monteiro, A. A. Leaf position, leaf age and plant age affect the expression of downy mildew resistance in Brassica oleracea. Eur. J. Plant Pathol. 125, 179-188 (2009).

[50] Pyne, R. M., Koroch, A. R., Wyenandt, C. A. & Simon, J. E. A rapid screening approach to identify resistance to basil downy mildew (Peronospora belbahrii). HortScience 49, 1041-1045 (2014).

[51] Coates, M. E. & Beynon, J. L. Hyaloperonospora arabidopsidis as a pathogen model. Ann. Rev. Phytopathol. 48, 329-345 (2010).

[52] Bellin, D. et al. Resistance to Plasmopara viticola in grapevine ‘Bianca’is controlled by a major dominant gene causing localised necrosis at the infection site. Theor. Appl. Genet. 120, 163-176 (2009).

[53] Trigiano, R. N. in Plant Pathology Concepts and Laboratory Exercises (eds Ownley, B. H. & Trigiano, R. N. ) 349-356 (CRC Press, 2007).

[54] Lapin, D., Meyer, R. C., Takahashi, H., Bechtold, U. & Van den Ackerveken, G. Broad-spectrum resistance of Arabidopsis C24 to downy mildew is mediated by different combinations of isolate-specific loci. New Phytol. 196, 1171-1181 (2012).

[55] Govindarajulu, M., Epstein, L., Wroblewski, T. & Michelmore, R. W. Host‐induced gene silencing inhibits the biotrophic pathogen causing downy mildew of lettuce. Plant Biotechnol. J. 13, 875-883 (2015).

[56] Diez-Navajas, A., Wiedemann-Merdinoglu, S., Greif, C. & Merdinoglu, D. Non-host versus host resistance to the grapevine downy mildew, Plasmopara viticola, studied at the tissue level. Phytopathol. 98, 776-780 (2008).

[57] Yin, X. et al. Pathogen development and host responses to Plasmopara viticola in resistant and susceptible grapevines: an ultrastructural study. Hortic. Res. 4, 17033 (2017).

[58] Perazzolli, M. et al. Downy mildew resistance induced by Trichoderma harzianum T39 in susceptible grapevines partially mimics transcriptional changes of resistant genotypes. BMC Genomics 13, 660 (2012).

[59] Kunjeti, S. G. et al. RNA‐Seq reveals infection‐related global gene changes in Phytophthora phaseoli, the causal agent of lima bean downy mildew. Mol. Plant Pathol. 13, 454-466 (2012).

[60] Kulkarni, K. S. et al. De novo transcriptome sequencing to dissect candidate genes associated with pearl millet-downy mildew (Sclerospora graminicola Sacc.) interaction. Front. Plant Sci. 7, 847 (2016).

[61] Wang, M. et al. A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation. New Phytol. 217, 163-178 (2018).

[62] Zhang, G. et al. PacBio full‐length cDNA sequencing integrated with RNA‐seq reads drastically improves the discovery of splicing transcripts in rice. Plant J. 97, 296-305 (2019).

[63] Hu, Z. et al. Full-length transcriptome assembly of Italian ryegrass root integrated with RNA-Seq to identify genes in response to plant cadmium stress. Int. J. Mol. Sci. 21, 1067 (2020).

[64] Feng, S., Xu, M., Liu, F., Cui, C. & Zhou, B. Reconstruction of the full-length transcriptome atlas using PacBio Iso-Seq provides insight into the alternative splicing in Gossypium australe. BMC Plant Biol. 19, 1-16 (2019).

[65] Kuang, X., Sun, S., Wei, J., Li, Y. & Sun, C. Iso-Seq analysis of the Taxus cuspidata transcriptome reveals the complexity of Taxol biosynthesis. BMC Plant Biol. 19, 210 (2019).

[66] Tan, X. et al. Global expression analysis of nucleotide binding site-leucine rich repeat-encoding and related genes in Arabidopsis. BMC Plant Biol. 7, 1-20 (2007).

[67] Li, M. Y. et al. CRISPR/Cas9-mediated VvPR4b editing decreases downy mildew resistance in grapevine (Vitis vinifera L.). Hortic. Res. 7, 149 (2020).

[68] Qi, L., Ma, G., Li, X. & Seiler, G. Diversification of the downy mildew resistance gene pool by introgression of a new gene, Pl₃₅, from wild Helianthus argophyllus into oilseed and confection sunflowers (Helianthus annuus L.) . Theor. Appl. Genet. 132, 2553-2565 (2019).

[69] Ma, G. et al. Molecular dissection of resistance gene cluster and candidate gene identification of Pl₁₇ and Pl₁₉ in sunflower by whole-genome resequencing . Sci. Rep. 9, 1-10 (2019).

[70] Salmela, L. & Rivals, E. LoRDEC: accurate and efficient long read error correction. Bioinformatics 30, 3506-3514 (2014).

[71] Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658-1659 (2006).

[72] Alamancos, G. P., Pagès, A., Trincado, J. L., Bellora, N. & Eyras, E. Leveraging transcript quantification for fast computation of alternative splicing profiles. RNA 21, 1521-1531 (2015).

[73] Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[74] Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120 (2014).

[75] Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644-652 (2011).

[76] Pertea, G. et al. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19, 651-652 (2003).

[77] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[78] Conesa, A. et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674-3676 (2005).

[79] Finn, R. D., Clements, J. & Eddy, S. R. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 39, W29-W37 (2011).

[80] Delorenzi, M., Speed, T. & An, H. M. M. model for coiled-coil domains and a comparison with PSSM-based predictions. Bioinformatics 18, 617-625 (2002).

[81] Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313 (2014).

[82] Magalhães, D. M. et al. LRR-RLK family from two Citrus species: genome-wide identification and evolutionary aspects. BMC Genomics 17, 623 (2016).

[83] Fischer, S. et al. Using OrthoMCL to assign proteins to OrthoMCL‐DB groups or to cluster proteomes into new ortholog groups. Curr. Protoc. Bioinformatics 35, 6.12.1-6.12.19 (2011).

[84] Anders, S., Pyl, P. T. & Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166-169 (2015).

[85] Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).

[86] Vera Alvarez, R., Pongor, L. S., Mariño-Ramírez, L. & Landsman, D. TPMCalculator: one-step software to quantify mRNA abundance of genomic features. Bioinformatics 35, 1960-1962 (2019).